Relatives 【欧拉函数】

1214: Relatives

NK

Given n, a positive integer, how many positive integers less than n are relatively prime to n? Two integers a and b are relatively prime if there are no integers x > 1, y > 0, z > 0 such that a = xy and b = xz.

Input

There are several test cases. For each test case, standard input contains a line with n <= 1,000,000,000. A line containing 0 follows the last case.

Output

For each test case there should be single line of output answering the question posed above.

Sample Input

7
12
0

Sample Output

6
4

 

 

分析：

在数论，对正整数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<stdio.h> #include<string.h> #define size 32700 int hk,y[size],p[size]; void prime() { int i,j; hk=0; for(i=2;i*i<=size;i++) { if(!y[i]){ for(j=i*i;j<=size;j+=i){y[j]=1;} } } for(i=2;i<=size;i++) { if(!y[i]){p[hk]=i;hk++;}} } void solve() { int n,i,flag,m,s1,s2,t; double sum; while(scanf("%d",&n),n) { sum=1; m=n; for(i=0;i<hk&&n>1;i++) { flag=0; t=1; while((0==n%p[i])&&n>1) { flag=1; n/=p[i]; t*=p[i]; } if(flag){sum*=(t/p[i])*(p[i]-1);} } if(n>1){sum*=n-1;} printf("%d/n",(int)sum); } } int main() { memset(y,0,sizeof(y)); prime(); solve(); return 0; }