莫比乌斯反演.pptVIP

莫比乌斯反演

什么是莫比乌斯反演？介绍这个之前，我们先来看一个函数根据F(n)的定义，我们显然有：F(1)=f(1)F(2)=f(1)+f(2)F(3)=f(1)+f(3)F(4)=f(1)+f(2)+f(4)F(5)=f(1)+f(5)F(6)=f(1)+f(2)+f(3)+f(6)F(7)=f(1)+f(7)F(8)=f(1)+f(2)+f(4)+f(8)

于是我们可以通过F(n)推导出f(n)f(1)=F(1)f(2)=F(2)-F(1)f(3)=F(3)-F(1)f(4)=F(4)-F(2)f(5)=F(5)-F(1)f(6)=F(6)-F(3)-F(2)+F(1)f(7)=F(7)-F(1)f(8)=F(8)-F(4)在推导的过程中我们是否发现了一些规律？

莫比乌斯函数的性质(1)对于任意正整数n有：证明：①当时显然②当时，将分解可以得到在的所有因子中，值不为零的只有所有质因子次数都为1的因子，其中质因数个数为个的因子有个那么显然有：

莫比乌斯函数的性质只需证明即可二项式定理：令，代入即可得证。

莫比乌斯函数的性质(2)对于任意正整数n有：其实没什么用，结论挺好玩的可以记一下只需要令，代入莫比乌斯反演的公式即可至于为什么就留给大家做思考题了

莫比乌斯函数的性质(3)积性函数数论上积性函数的定义：积性函数的性质：①②积性函数的前缀和也是积性函数

BZOJ2440完全平方数题目大意：求第k个无平方因子数无平方因子数(Square-FreeNumber)，即分解之后所有质因数的次数都为1的数首先二分答案问题转化为求[1,x]之间有多少个无平方因子数根据容斥原理可知对于sqrt(x)以内所有的质数有x以内的无平方因子数=0个质数乘积的平方的倍数的数的数量(1的倍数)-每个质数的平方的倍数的数的数量(9的倍数,25的倍数,...)+每2个质数乘积的平方的倍数的数的数量(36的倍数,100的倍数,...)-...

现在我们来证明莫比乌斯反演定理证明：这里利用到了这条性质形式二：证明同理一般要用到的都是这种形式

BZOJ2301Problembn次询问，每次询问有多少个数对(x,y)满足a=x=b,c=y=d且gcd(x,y)=kN=5W,1=a=b=5W,1=c=d=5W首先利用容斥原理将一个询问拆分成四个，每次询问有多少个数对(x,y)满足1=x=n,1=y=m且gcd(x,y)=k这个问题等价于询问有多少个数对(x,y)满足1=x=floor(n/k),1=y=floor(m/k)且x与y互质利用NOI2010能量采集中的方法，我们可以得到一个O(nlogn)的算法，但是这显然不能胜任此题的数据范围这时候我们就可以考虑莫比乌斯反演了

BZOJ2301Problemb由于之前的结论，我们可以令f(i)为1=x=n,1=y=m且gcd(x,y)=i的数对(x,y)的个数，F(i)为1=x=n,1=y=m且i|gcd(x,y)的数对(x,y)的个数那么显然有反演后可得枚举原题中k的每一个倍数，我们就可以O(n)时间处理每个询问了但是O(n)还是不能胜任本题的数据范围考虑进一步优化

BZOJ2301Problemb观察式子，发现最多有个取值那么就至多有个取值枚举这个取值，对莫比乌斯函数维护一个前缀和，就可以在时间内出解总时间复杂度枚举除法的取值这种方法在莫比乌斯反演的应用当中非常常用，且代码并不难写不难写？怎么写？

BZOJ2301Problembif(nm)s);for(i=1;i=n;i=last+1){ last=min(n/(n/i),m/(m/i)); re+=(n/i)*(m/i)*(sum[last]-sum[i-1]);}returnre;超级好写不是么？

BZOJ2820YGY的GCD题目大意：求有多少数对(x,y)(1=x=n,1=y=m)满足gcd(x,y)为质数n,m=1000W数据组数=1W首先我们枚举每一个质数那么答案就是直接做显然TLE考虑优化令，那么有

BZOJ3529数表题目大意：令F(i)为i的约数和，q次给定n,m,a