线性求欧拉函数

欧拉函数：在数论中，对正整数n，欧拉函数是少于或等于n的数中与n互质的数的数目。
先给出一个结论：
设P是素数，
若p是x的约数，则E(x*p)=E(x)*p.
若p不是x的约数，则E(x*p)=E(x)E(p)=E(x)(p-1).
证明如下：
E(x)表示比x小的且与x互质的正整数的个数。
*若p是素数，E(p)=p-1。
*E(p^k)=p^k-p^(k-1)=(p-1)*P^(k-1)
证：令n=p^k,小于n的正整数数共有n-1即(p^k-1)个,其中与p不质的数共[p^(k-1)-1]个(分别为1*p,2*p,3*p…p(p^(k-1)-1))。
所以E(p^k)=(p^k-1)-(p^(k-1)-1)=p^k-p^(k-1).得证。
*若ab互质，则E(a*b)=E(a)*E(b)，欧拉函数是积性函数.
对任意数n都可以唯一分解成n=p1^a1*p2^a2*p3^a3…*pn^an(pi为素数).
则E(n)=E(p1^a1)E(p2^a2)*E(p3^a3)…*E(pn^an)
=(p1-1)p1^(a1-1)(p2-1)p2^(a2-1)…*(pn-1)*pn^(an-1)
=(p1^a1*p2^a2*p3^a3*…pn^an)[(p1-1)(p2-1)(p3-1)…(pn-1)]/(p1*p2*p3*…*pn)
=n*(1-1/p1)(1-1/p2)…*(1-1/pn)
* E(p^k) =(p-1)*p^(k-1)=(p-1)*p^(k-2)*p
E(p^(k-1))=(p-1)*p^(k-2)
->当k>1时，E(p^k)=E(p*p^(k-1))=E(p^(k-1))*p.
(当k=1时，E(p)=p-1.)
由上式: 设P是素数，
若p是x的约数，则E(x*p)=E(x)*p.
若p不是x的约数，则E(x*p)=E(x)E(p)=E(x)(p-1). 证明结束。

代码

void Get_phi()///筛法求欧拉函数  
{  
    cnt = 0;  
    memset(flag, true, sizeof(flag));  
    phi[1] = 1;  
    for(int i=2; i<MAXN; i++)///线性筛法  
    {  
        if(flag[i])///素数  
        {  
            p[cnt++] = i;  
            phi[i] = i-1;///素数的欧拉函数值是素数 - 1  
        }  
        for(int j=0; j<cnt; j++)  
        {  
            if(i*p[j] > MAXN)  
                break;  
            flag[i*p[j]] = false;///素数的倍数，所以i*p[j]不是素数  
            if(i%p[j] == 0)///性质：i mod p == 0, 那么 phi(i * p) == p * phi(i)  
            {  
                phi[i*p[j]] = p[j] * phi[i];  
                break;  
            }  
            else  
                phi[i*p[j]] = (p[j]-1) * phi[i];///i mod p != 0, 那么 phi(i * p) == phi(i) * (p-1)   
        }  
    }  
}