关于质数

欧拉函数

求phi(n)

//phi(1)==1
int phi(int k){
    int ans=k,a=k;
    for(int i=2;i*i<=a;i++)
        if(a%i==0){
            ans-=ans/i;
            while(a%i==0)a/=i;//算数基本定理
        }
    if(a>1)ans-=ans/a;//存在大于sqrt(k)的质因子
    return ans;
}

欧拉线性筛

筛质数

优秀！O(n)时间复杂度

每个合数仅会被它的最小质因数筛去1次

为什么是最小质因数？

假设A=x*y（x是质数，y是合数），y=m*n（m是质数，显然m<y），那么A= x*m*n= m*Z。可得，A可以表示为一个更大的合数和一个更小的质数的乘积，所以A一定会被它最小的质因数筛掉。

比如24，它只会被2*12筛掉。

void getPrime(int n){//筛n以内的质数
    int cnt=0;
    for(int i=2;i<=n;i++){
        if(Mark[i]==false)Prime[++cnt]=i;//没被筛过，是质数
        for(int j=1;j<=cnt && i*Prime[j]<=n;j++){
            Mark[i*Prime[j]]=true;//标记为已经被筛
            if(i%Prime[j]==0)break;//保证后面的数不会重复被筛
        }
    }
}

求phi(1)+phi(2)+…+phi(n)

欧拉函数的性质（p为质数）：

1.phi(p)=p-1;

2.(k%p=0)phi(k*p)=p*phi(k);

3.(k%p!=0,即二者互质)phi(k*p)=(p-1)*phi(k)（欧拉函数是积性函数）.

优秀！同样O(n)的时间复杂度

void phi(int n){
    int cnt=0;
    for(int i=2;i<=n;i++){
        if(Mark[i]==false)Prime[++cnt]=i,Phi[i]=i-1;//性质1
        for(int j=1;j<=cnt && i*Prime[j]<=n;j++){
            Mark[i*Prime[j]]=true;
            if(i%Prime[j]==0){//性质2
                Phi[i*Prime[j]]=Phi[i]*Prime[j];break;
            }else Phi[i*Prime[j]]=Phi[i]*(Prime[j]-1);//性质3
        }
    }
}