莫比乌斯反演 学习笔记

预备知识

枚举除法

⌊ni⌋只有O(n√)种取值
并且对于i，⌊n⌊ni⌋⌋是i被n除并下取整取值相同的一段区间的右端点

一个非常有用性质：
⌊nab⌋=⌊⌊na⌋b⌋=⌊⌊nb⌋a⌋

积性函数

f(ab)=f(a)f(b),(a,b)=1
完全积性函数：不要求ab互质。

若函数f(n)为积性函数，那么
f(n)=∏if(pkii)
并且由于f(a∗1)=f(a)∗f(1)
可以得出f(1)=1

常见函数

id(n)=n
e(n)=[n=1]
1(n)=1
d(n)=n的约数数量
σ(n)=n所有约数的和
除了最后一个其余都是积性函数

狄利克雷卷积

对于算术函数f(n),g(n)，定义其狄利克雷卷积
(f×g)(n)=∑d|nf(d)g(nd)
比如说一个公式n=∑d|nφ(d)就可以表示为：id=φ×1
狄利克雷卷积的运算满足：
（1）交换律：f×g=g×f
（2）结合律：(f×g)×h=f×(g×h)
（3）分配律：f×(g+h)=f×g+f×h

几个性质：
1、存在单位函数e满足f=f×e=e×f
可以根据卷积的定义理解一下，枚举n的约数时，只有当d=n时e(nd)=1
2、如果f,g都为积性函数，那么f×g也为积性函数

莫比乌斯函数μ

定义

（1）若d=1，那么μ(d)=1
（2）若d=p1p2..pk，pi均为互异素数，那么μ(d)=(−1)k
（3）其他情况下，μ(d)=0

显然，μ为积性函数。

性质

[n=1]=∑d|nμ(d)
即1×μ=e
不会证明…

莫比乌斯反演

公式

F(n)和f(n)是定义在非负整数集合上的两个函数
若满足条件F(n)=∑d|nf(d)
那么可以得出结论f(n)=∑d|nμ(d)F(nd)
即已知F=f×1
则f=μ×F

证明

利用狄利克雷卷积证明莫比乌斯反演公式
1° 已知F=f×1，两边同乘μ得F×μ=f×1×μ
2° 1×μ=e，那么F×μ=f×e
3° 又因为f×e=f，得出f=μ×F
证毕

另一个公式

F(n)和f(n)是定义在非负整数集合上的两个函数
若满足条件F(n)=∑n|df(d)
那么可以得出结论f(n)=∑n|dμ(dn)F(d)

似乎这个公式在做题的时候更常用些
因为反演的题gcd用的比较多…

化式子常用技巧

∑i=1n∑d|i=∑d=1n⌊nd⌋

杜教筛

栗子：求∑i=1nφ(i)，∑i=1nμ(i)，1≤n≤231−1

关于杜教筛：
求F(n)=∑if(i),存在g=f×I, I表示恒等函数，即g=∑d|nf(d)
我们定义G(n)=∑ig(i) ,就可以得到F(n)=G(n)−∑iF(⌊ni⌋)
如果G(n)可以在一定时间内求解，那么我们最好情况下可以做到O(n23)的时间求解F(n)

假设计算出F(n)的复杂度为T(n)，则有T(n)=O(n√)+∑i=1n√T(i)+T(ni)，这里只展开一层就可以了，更深层的复杂度是高阶小量，所以有T(n)=∑∑i=1n√O(i√)+O(ni−−√)=O(n34)。
如果用筛法预处理前k个F(n)，且k≥n√，则复杂度变为T(n)=∑i=1nkni−−√=O(nk√)，当k=O(n23)时可以取到较好的复杂度T(n)=O(n23)。

实际上上面的复杂度证明我也不懂…只是搬运了神犇的讲解…

φ

利用一个反演公式n=∑d|nφ(d)=∑d|n,d<nφ(d)+φ(n)
得出φ(n)=n−∑d|n,d<nφ(d)
那么设P(n)=∑i=1nφ(n)
P(n)=∑i=1nφ(n)
=∑i=1n(i−∑d|i,d<iφ(d))
=n(n+1)2−∑i=1n∑d|i,d<iφ(d)
令i=id
那么原式=n(n+1)2−∑id=1n∑d|id,d<idφ(d)
=n(n+2)2−∑i=2n∑d=1⌊ni⌋φ(d)=n(n+1)2−∑i=2nP(⌊ni⌋)
也就是说P(n)=n(n+1)2−∑i=2nP(⌊ni⌋)

μ

设M(n)=∑i=1nμ(i)
同样利用一个反演公式1(n)=∑d|nμ(d)来化简
1(n)=∑i=1n[i=1]=∑i=1n∑d|iμ(d)=∑i=1n∑d=1n[d|i]μ(id)
令i=id
那么原式=∑id=1n∑d=1nμ(i)=∑d=1n∑i=1⌊nd⌋μ(i)=∑i=1n∑d=1⌊ni⌋μ(d)=∑i=1nM(⌊ni⌋)
也就是说∑i=1nM(⌊ni⌋)=1
又因为M(n)=∑i=1nM(⌊ni⌋)−∑i=2nM(⌊ni⌋)
所以M(n)=1−∑i=2nM(⌊ni⌋)
咦怎么感觉就是省了点常数…

加几个链接：
http://blog.csdn.net/skywalkert/article/details/50500009
http://jiruyi910387714.is-programmer.com/posts/195270.html
http://www.cnblogs.com/abclzr/p/6242020.html

例题