计算欧拉函数的原理

求欧拉函数的一般方法：

1.我们知道一个素数p的欧拉函数f(p)=p-1；那么p的k次幂，即n=p^k，则容易证明：f(n)=p^k-p^(k-1)；

 

证明：已知少于p^k的数有p^k-1，其中与p^k不互质的数有p^(k-1)-1个，分别为（p*1，p*2，p*3，····，p*（p^(k-1)-1））。所以f(n)=p^(k)-1-(p^(k-1)-1)=p^(k)-p^(k-1)。

 

2.假设p，q为两个互质的数，则：f(p*q)=f(p)*f(q)；

      

 原来这里一直写着一个略字，倒不是因为太简单，而是我没证出来，一开始自己写了一个证明，但后来发现不对，然后一直也没想出如何证，今天看初等数论第三章同余的时候找到了证明，所以赶紧过来整理一下。

 

证明：这里主要用到了简化剩余系，x1，x2分别通过模m1，m2的简化剩余系，则m2*x1+m1*x2通过模m1*m2的简化剩余系。由于m2*x1+m1*x2通过f（m1*m2）个整数，x1通过f(m1)个整数，x2通过f（m2）个整数，所以f（m1*m2）=f（m1）*f（m2）；

有了以上两个性质，欧拉函数就非常简单了，先对一个数进行整数分解，然后···

下面介绍一种快速的方法：

 

因为任意正整数都可以唯一表示成如下形式:
          k=p1^a1*p2^a2*……*pi^ai;(即分解质因数形式)                可以推出:E(k)=(p1-1)(p2-1)……(pi-1)*(p1^(a1-1))(p2^(a2-1))……(pi^(ai-1))
                          =k*(p1-1)(p2-1)……(pi-1)/(p1*p2*……pi);
                          =k*(1-1/p1)*(1-1/p2)....(1-1/pk)
      ps：在程序中利用欧拉函数如下性质，可以快速求出欧拉函数的值(a为N的质因素)                                                                                       若(N%a==0 && (N/a)%a==0) 则有:E(N)=E(N/a)*a;          若(N%a==0 && (N/a)%a!=0) 则有:E(N)=E(N/a)*(a-1);