歐拉函數 素數線性篩法模板

/**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("");
}