欧拉函数的简单求法及欧拉打表法

欧拉函数的意义为，对于一个正整数n来说，euler（n）即为小于等于n，且与n互斥的数的个数

根据容斥原理可得：

euler（n）=n-n/p1-n/p2....-n/pn+n/p1p2+n/p1p3+n/p2p3...+n/pn-1pn-n/p1p2p3.......

=n(1-1/p1)(1-1/p2)....(1-1/pn)

同样根据积性原理可得:

若m，n互斥则

euler（m*n）=euler(m)*euler(n);

因为n=p1^k1+p2^k2....+pn^kn

所以euler(n)=euler(p1^k1)*euler(p2^k2)*....euler(pn^kn)

因为p为素数的话euler(pn^kn)=pn^kn-pn^(kn-1)=pn^(kn-1)(1-1/pn);

所以euler(n)=n(1-1/p1)(1-1/p2)....(1-1/pn)

计算欧拉函数:

 int euler_phi(int n)
 {
      int m=(int)sqrt(n+0.5);
      int ans=n;
      for(int i=2;i<=m;i++)
          if(n%i==0)//找素因子
          {
              ans=ans/i*(i-1);
              while(n%i==0)n/=i;//除尽
          }
     if(n>1)ans=ans/n*(n-1);
     return ans; 
 }

打表求欧拉函数:

int phi[maxn];
void phi_table(int n)
{
    for(int i=2;i<=n;i++)phi[i]=0;
    phi[1]=1;
    for(int i=2;i<=n;i++)
        if(!phi[i])
            for(int j=i;j<=n;j+=i)
            {
                if(!phi[j])phi[j]=j;
                phi[j]=phi[j]/i*(i-1);
            }
}