ACM 第十九天

积性函数

积性函数线性筛，筛素数,u(n),欧拉函数：

vis[0]=vis[1]=1,mu[1]=1,phi[1]=1;
for (RG int i=2;i<=N;++i){
    if (!vis[i]) mu[i]=-1,phi[i]=i-1,prime[++cnt]=i;
    for (RG int j=1,k=i*prime[j];j<=cnt && k<=N;++j,k=i*prime[j]){
    vis[k]=1;
    if (!(i%prime[j])){ mu[k]=0,phi[k]=phi[i]*prime[j]; break; }
    else mu[k]=-mu[i],phi[k]=phi[i]*phi[prime[j]];
    }
}

可以发现，线性筛分为3部分：

1.n本身是素数，这个根据积性函数的定义可得，很容易求。

2.i%prime[j]!=0，这个也是根据积性函数的性质可得。f(a)f(b)=f(a

3.i%prime[j]==0，可能需要找规律。据ljh2000神犇的说法，可以用2,4,8或3,9,27这些数来找。

μ(n),ϕ