莫比乌斯反演:
看不懂挺难得,看懂了就还好,不敢说什么都会,但是博客上的基础题gcd这种类型的会了;
莫比乌斯反演:就是已知和函数求分函数,表达可能不恰当,但就是这个意思;
莫比乌斯函数u(d):
性质一:d=1时,u(1)=1;
性质二:d=p1p2···pk,Pi均为素数时,u(d)=(-1)^k,注意素数的指数都是一次的;
性质三:其他情况u(d)=0;
这样用代码实现就好了;
线性筛选莫比乌斯反演函数的代码:
int prim[maxn],cnt,mu[maxn];
bool v[maxn];
void init()
{
cnt=0;
memset(v,0,sizeof(v));
mu[1]=1;//对应性质之一
for(int i=2;i<maxn;i++)
{
if(!v[i])
{
prim[cnt++]=i;
mu[i]=-1;//对应性质二
}
for(j=0;j<cnt&&i*prim[j]<maxn;j++)
{
v[i*prim[j]]=1;
if(i%prim[j])
mu[i*prim[j]]=-mu[i];//对应性质二
else
{
mu[i*prim[j]]=0;//对应性质三
break;
}
}
}
}
f(k):表示gcd(x,y)=k的对数的个数;
F(k):表示gcd(x,y)=k的倍数的个数(和);
F(k)=(m/k)*(n/k);
F(k)=signmaf(p*k);p=1,2,···min(m,n);这和莫比乌斯反演不大像,另d=p*k;
那么f(k)=sigamau(p)F(p*k)p=1,2,···,min(m,n);