首页 > 其他分享 >最小公倍佩尔数 题解

最小公倍佩尔数 题解

时间:2023-02-04 08:11:06浏览次数:59  
标签:subset frac gcd 题解 sum 公倍 佩尔 ne empty

首先需要知道 \(f\) 是个啥

这里直接给出结论,过程可以看大佬的博客

\(f(n) = 2f(n - 1) + f(n - 2)\)

\(f(0)=0\)

\(f(1)=1\)

这种类似 斐波那契数列的递推式有结论 \(gcd(f(x), f(y)) = f(gcd(x, y))\) 通过辗转相减证明

那么一个集合的 \(gcd\) 是好求的,所以考虑进行 \(min-max\) 容斥

\(lcm(S)=\prod_{T \subset S,T \ne \not\empty}gcd(T)^{(-1)^{|T| - 1}}\)

那我们求的就是

\[ \prod_{T \subset S,T \ne \not\empty}f(gcd(T))^{(-1)^{|T| - 1}} \]

\[ = \prod_{d= 1}^{n}f(d)^{\sum_{T \subset S,T \ne \not\empty}[gcd(T) = d](-1)^{|T| - 1}} \]

现在只处理角标

\[ \sum_{T \subset S,T \ne \not\empty}[gcd(T) = d](-1)^{|T| - 1} \]

\[= \sum_{T \subset S,T \ne \not\empty}\sum_{x|\frac{gcd(T)}{d}}\mu(x)(-1)^{|T| - 1} \]

\[= \sum_{T \subset S,T \ne \not\empty}\sum_{xd | gcd(T)}\mu(x)(-1)^{|T| - 1} \]

\[= \sum_{x= 1}^{\frac{n}{d}}\mu(x)\sum_{T \subset S`}(-1)^{|T| - 1} \]

其中 \(s`\) 为 \(s\) 中 \(xd\) 的倍数除以 \(xd\),即\({1, 2, .. \frac{n}{xd}}\)

\[\sum_{T \subset S`,T \ne \not\empty}(-1)^{|T| - 1} \]

\[=\sum_{i = 1}^{size}(-1)^{i + 1} \]

\[=-\sum_{i = 1}^{size}(-1)^i \]

\[ = 1- \sum_{i = 0}^{size}(-1)^i \]

\[ =1-[T = \empty] \]

那么带回去

\[\sum_{x= 1}^{\frac{n}{d}}\mu(x)\sum_{T \subset S`}(-1)^{|T| - 1} \]

因为 \(dx <= n\) 所以\(S`\)非空

\[ =\sum_{x= 1}^{\frac{n}{d}}\mu(x) \]

于是我们要求的变成

\[ = \prod_{d= 1}^{n}f(d)^{\sum_{x= 1}^{\frac{n}{d}}\mu(x)} \]

考虑 \(dx <= n\) 的都会产生贡献,我们枚举 \(i = dx\),式子就是

\[= \prod_{i= 1}^{n}\prod_{d|i}f(\frac{i}{d})^{\mu(d)} \]

标签:subset,frac,gcd,题解,sum,公倍,佩尔,ne,empty
From: https://www.cnblogs.com/Chencgy/p/17090829.html

相关文章

  • P3119 [USACO15JAN]Grass Cownoisseur G 题解
    做过的原题,模拟赛时PDF里的题面实在有点难受。首先有显然结论:在一个环上反走一定是不值的,因为环上的点本来就相互可达。所以考虑缩点。缩点后的问题可以看成:求对于每一......
  • P8368 [LNOI2022] 串 题解
    题目链接题目分析题目要求我们构造一个最长的\(T\)序列,我们首先从每个\(T_i\)入手,思考如何安排才能合法。容易观察到对于每个\(T_i\),合法的\(T_{i-1}\)有两种方......
  • 【题解】P5219 无聊的水题 I
    思路prufer序列+卷积优化dp.首先考虑到令\(a\)为原树的prufer序列,则\(\sum\limits_{i=1}^{n-2}[a_i=k]=\operatorname{deg}(k)\),其中\(\operatorname......
  • IIS WordPress 单站点,多站点 中文URL乱码和重定向多次问题解决方法
    前提是需要安装IISURL重写组件(安装方法这里不说明,搜一下资料就有)1、站点网站根目录新增一个chineseurl.php文件用来处理中文url问题<?php//IISMod-Rewriteif(is......
  • Maven Use STAR or POSIX extensions to overcome this limit 报错问题解决
     问题:Mavenassembly-plugin在打包的时候如果出现groupid'707420648'istoobig(>2097151). UseSTARorPOSIXextensionstoovercomethislimit解决:只需在......
  • Hibernate下分页语句第一页和第二页sql语句不一致问题解决方法
    <propkey="hibernate.dialect">新增类OracleDialect</prop>OracleDialectextendsorg.hibernate.dialect.Oracle10gDialect重写getLimitString方法sql......
  • LeetCode 对称二叉树算法题解 All In One
    LeetCode对称二叉树算法题解AllInOne对称二叉树原理图解101.SymmetricTree对称二叉树https://leetcode.com/problems/symmetric-tree/https://leetcode.c......
  • 【题解】P3750 [六省联考 2017] 分手是祝愿
    出题老哥收收味吧,阿米奈!记一下几个常用的手段。昨晚CF的D是差不多的思路吧,不然不会来做。思路期望dp.先做一些准备工作,求一下逆元和每个数的因数,复杂度\(O(n\l......
  • 【题解】[CSP-S2019] Emiya 家今天的饭
    题目分析:(我竟然可以独立做出来正赛的题,表示震惊)这个题面显然就很神仙,不好分析,我们进行转化一下题意:给定一个\(n\timesm\)的矩阵,对于每一行我们可以选择一个数也可以......
  • CF1311E 题解
    题意:构造一个节点数为\(n\)的树,使得节点深度之和为\(d\).(根节点为\(1\)且深度为\(0\))显然,不断构造二叉树并检查是否为答案是行不通的,只能在\(O(n)\)的......