JOJ1874:Relatives

传送门：http://acm.jlu.edu.cn/joj/showproblem.php?pid=1874

求欧拉函数。

分析：

在数论，对正整数n，欧拉函数是少于或等于n的数中与n互质的数的数目。此函数以其首名研究者欧拉命名，它又称为Euler's totient function、φ函数、欧拉商数等。

　　例如φ(8)=4，因为1,3,5,7均和8互质。

　　从欧拉函数引伸出来在环论方面的事实和拉格朗日定理构成了欧拉定理的证明。

　　φ函数的值

　　φ(1)=1（唯一和1互质的数就是1本身）。

　　若n是质数p的k次幂，φ(n)=p^k-p^(k-1)=(p-1)p^(k-1)，因为除了p的倍数外，其他数都跟n互质。

　　欧拉函数是积性函数——若m,n互质，φ(mn)=φ(m)φ(n)。

　　证明：设A, B, C是跟m, n, mn互质的数的集，据中国剩余定理，A*B和C可建立一一对应的关系。因此φ(n)的值使用算术基本定理便知，

　　若

　　n= ∏p^(α（下标p）)

　　p|n

　　则φ(n)=∏(p-1)p^(α（下标p）-1)=n∏(1-1/p)

　　p|n p|n

　　例如φ(72)=φ(2^3×3^2)=(2-1)2^(3-1)×(3-1)3^(2-1)=24

　　与欧拉定理、费马小定理的关系

　　对任何两个互质的正整数a, m, m>=2有

　　a^φ(m)≡1(mod m)

　　即欧拉定理

　　当m是质数p时，此式则为：

　　a^(p-1)≡1(mod m)

　　即费马小定理。

　　欧拉函数的编程实现：

　　利用欧拉函数和它本身不同质因数的关系，用筛法计算出某个范围内所有数的欧拉函数值。

　　欧拉函数和它本身不同质因数的关系：欧拉函数ψ（Ｎ）＝Ｎ｛∏ｐ｜Ｎ｝（１－１/ｐ）。（Ｐ是数Ｎ的质因数）

　　如：

　　ψ（１０）＝１０×（１－１/２）×（１－１/５）＝４；

　　ψ（３０）＝３０×（１－１/２）×（１－１/３）×（１－１/５）＝８；

　　ψ（４９）＝４９×（１－１/７）＝４２。

代码：

#include <cstdio>
using namespace std;

int main()
{
    int n, eular, i, j;
    while (scanf("%d", &n) != EOF)
    {
        if (n == 0)
            break;
        for (i = 2, eular = n; n > 1; ++i) //一个数，如果不能被a整除，绝对不可能被a的倍数整除。
        {
            if (n % i == 0)
            {
                eular -= eular / i;
                do
                {
                    n /= i;
                }
                while (n % i == 0);
            }
        }
        printf("%d\n", eular);
    }
    return 0;
}

n一直都大于等于i，所以循环结束的条件是n == 1；欧拉公式为： φ(n) = n * (1 - 1 / p1) * (1 - 1 / p2)....；p为n的质数因子。