欧拉函数 素数线性筛法模板

/**Eular素数筛*/
const int PN = 1000000;
int phi[PN+10], prime[PN+10], tot;
bool isNotPrime[PN+10];
void getPhiPrime()
{
    phi[1]=1;
    isNotPrime[1] = true;
    for(int i=2; i<=PN; i++) //相当于分解质因式的逆过程
    {
        if(!isNotPrime[i])
        {
            prime[++tot]=i;//筛素数的时候首先会判断i是否是素数。
            phi[i]=i-1;//当 i 是素数时 phi[i]=i-1
        }
        for(int j=1; j<=tot; j++)
        {
            if(i*prime[j]>PN)  break;
            isNotPrime[i*prime[j]]=1;//确定i*prime[j]不是素数
            if(i%prime[j]==0)//接着我们会看prime[j]是否是i的约数
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else  phi[i*prime[j]]=phi[i]*(prime[j]-1);//其实这里prime[j]-1就是phi[prime[j]]，利用了欧拉函数的积性
        }
    }
  //  printf("tot = %d\n", tot);
  //  for(int i = 1; i <= 100; i++) printf("%d ", prime[i]); puts("");
}