莫比乌斯反演入门--化简比较烦

最新推荐文章于 2022-04-15 10:14:32 发布

置顶 zjyang12345

最新推荐文章于 2022-04-15 10:14:32 发布

阅读量635

点赞数 1

分类专栏： —————数论————— 莫比乌斯&欧拉文章标签： ACM

本文链接：https://blog.csdn.net/zjyang12345/article/details/89143466

版权

—————数论————— 同时被 2 个专栏收录

31 篇文章 1 订阅

订阅专栏

莫比乌斯&欧拉

15 篇文章 1 订阅

订阅专栏

一.莫比乌斯函数

特点是，分解的素因子中有平方因子则为0，否则素因子个数奇偶性决定正负1。线性筛大全，含莫比乌斯

二.莫比乌斯反演公式

用的最多是第二条性质：

已知 $F(n)=\sum_{n|d}f(d)$ ，

可得 $f(n)=\sum_{n|d}\mu(\frac{d}{n})F(d)$

①注意反演的意义：个人拙见是类似反函数把Y=F(x)写成X=F(y)的形式，因为有时F(Y）比F(X)更容易求。

为啥引入莫比乌斯函数 $\mu$ ：实际上消去同类项，比如f(2)出现3次，必然要两项为﹢，一项为-；如果出现两次，有一项为0一项为+

具体见用容斥实例理解莫比乌斯反演过程

②注意公式符号的意义，n|d代表d是n的倍数（n是d的因子），实际上求的是n的所有倍数和

~~这里被符号绕晕了一天（还是太菜了）~~

三.反演实例

P2257 YY的GCD

①列出与GCD有关的莫比乌斯反演，都是套路！

设 $f(d)$ 为gcd(i,j)=d的个数，形式： $\sum_{?=1}^{max} (gcd=d)$

$F(n)$ 为gcd(i,j)=d的倍数的个数，

如F(2)=f(2*1)+f(2*2)+f(2*3)+... F(n)=f(d)+f(2d)+...

这样就满足反演前提： $F(n)=\sum_{n|d}f(d)$

②反演得到： $f(n)=\sum_{n|d}\mu(\frac{d}{n})F(d)$ （其实这步才刚开始，重点是后面的化简...）

③列出ANS=xxx，用f(a)代入，a为所求量（注意符号的意义、范围，别被绕晕）

这题得到 $Ans=\sum_{p\in prim}\sum_{p|d}\mu(\lfloor\frac{d}{p}\rfloor)F(d)$ （设p是d的质因子，即题目要求的量）

④化简，最难的一步，为了缩小范围，降低复杂度（推荐多看大佬过程，这里就引用大佬的化简过程，还有gcd化简全总结）

注意符号的意义范围

a.F(d)显然可以直接求得

b.令 $d=ip$ ， d是1~ min(n,m) ，则i范围为1~ $\dpi{100} {\left \lfloor \frac{min(m,n)}{p} \right \rfloor}$ (一般想求x|y之类的范围，我喜欢令y=kx之类，好理解)

$\sum_{isprime(p)}\sum_{i=1}^{\left \lfloor \frac{min(n,m)}{p} \right \rfloor}\mu(i){\left \lfloor \frac{n}{pi} \right \rfloor}{\left \lfloor \frac{m}{pi} \right \rfloor}$

c.如果这里停手的话，复杂度仍不满足需要，因为每次分块要枚举质数p（详见bzoj2301/hdu1695 求gcd=k代码），需要令 $T=di$ ，t=p,继续化简，这样分块的时候不用枚举p

$Ans=\sum_{T=1}^{min(n,m)}\sum_{t|T,t\in prime}\mu(\lfloor\frac{T}{t}\rfloor)\lfloor\frac{n}{T}\rfloor\lfloor\frac{m}{T}\rfloor$

由和式的性质，因为t是T的因数才会有贡献，所以可以提出来

$Ans=\sum_{T=1}^{min(n,m)}\lfloor\frac{n}{T}\rfloor\lfloor\frac{m}{T}\rfloor(\sum_{t|T,t\in prime}\mu(\lfloor\frac{T}{t}\rfloor))$

⑤ $\sum_{t|T,t\in prime}\mu(\lfloor\frac{T}{t}\rfloor)$ 可以预处理线性筛大全，含莫比乌斯

⑥后两项可以除法分块，记录前缀和即可，见链接

四.莫比乌斯反演性质

① $[n = 1] = \sum_{i|n} \mu(i)$

② $[gcd(i, j)=1] = \sum_{d|i, d|j}\mu(d)$

zjyang12345

关注

1
点赞
踩
2

收藏

觉得还不错? 一键收藏
0
评论
莫比乌斯反演入门--化简比较烦

一.莫比乌斯函数特点是，分解的素因子中有平方因子则为0，否则素因子个数奇偶性决定正负1。线性筛大全，含莫比乌斯二.莫比乌斯反演公式用的最多是第二条性质：已知，可得①注意反演的意义：个人拙见是类似反函数把Y=F(x)写成X=F(y)的形式，因为有时F(Y）比F(X)更容易求。为啥引入莫比乌斯函数：实际上消去同类项，比如f(2)出现3次，必然要两项为...
复制链接

扫一扫