poj 2407 Relatives

题目链接：http://poj.org/problem?id=2407

 

解题思路：

φ函数的值  Euler函数
通式：φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn)
其中p1, p2……pn为x的所有质因数，x是不为0的整数
φ(1)=1（唯一和1互质的数就是1本身）【注意：每种质因数只一个。比如12=2*2*3】

定理： 若n是素数p的k次幂，φ(n)=p^k-p^(k-1)=(p-1)p^(k-1)，因为除了p的倍数外，其他数都跟n互质
欧拉函数是积性函数——若m,n互质，φ(mn)=φ(m)φ(n)

代码：

#include<stdlib.h>
#include<cstdio>
#include<math.h>
int phi(int a)
{
	int temp=a;
	for(int i=2;i*i<=temp;i++)
		if(a%i==0)
		{
			while(a%i==0)                                        
				a/=i;
			temp=temp/i*(i-1);
		}
		if(a>1)
			temp=temp/a*(a-1);
	return temp;
}
int main()
{
	int N;
	while(scanf("%d",&N),N)
	{
		printf("%d\n",phi(N));
	}
	return 0;
}