莫比乌斯反演

莫比乌斯反演

 

略微学了一点皮毛，先记录下来以备复习。

反演分为两种形式：

①约数型

对于等式      

可以得到  

②倍数型

对于等式  

可以得到  

 

u函数的性质决定了莫比乌斯反演的本质是容斥原理。

一般情况下，我们要求的 f 函数很复杂，不便直接求出，但我们可以构造出一个与其有关的 F 函数，而 F 函数相比来说比较容易求出，这样我们可以通过F函数的反演求出f函数。

 

u函数可以通过线性筛求出

bool vis[max_];
int prime[max_];
int pl=0;
int mu[max_];
void getmu()
{
	mu[1]=1;
	for(int i=2;i<max_;i++)
	{
		if(vis[i]==false)
		{
			prime[++pl]=i;
			mu[i]=-1;
		}
		for(int j=1;j<=pl&&prime[j]*i<max_;j++)
		{
			vis[i*prime[j]]=true;
			if(i%prime[j]==0)
			{
				mu[i*prime[j]]=0;
				break;
			}
			else
			mu[i*prime[j]]=-mu[i];
		}
	}
}

莫比乌斯反演的难点，在于对当前要求的 f 函数,如何构造出适当的 F 函数，不但要满足上式的性质，还要容易求。

入门题基本都是计算有关gcd(x,y)=k的个数,这个构造非常典型，

当1<=x<=n，1<=y<=m时，求gcd(x,y)=1的对数。

设f(d)为满足gcd(x,y)=t的对数，

则可以构造F(d)为满足gcd(x,y)%d=0的对数，那么F(d)=(n/d)*(m/d)，可以直接求出，

那么由莫比乌斯反演公式，我们就可以反演出要求的f(d)

具体可以看题目：