莫比乌斯反演
略微学了一点皮毛,先记录下来以备复习。
反演分为两种形式:
①约数型
对于等式
可以得到
②倍数型
对于等式
可以得到
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)
具体可以看题目: