mdd的烦恼

mdd的烦恼

时间限制： 1000 ms  |  内存限制： 65535 KB

难度： 3

描述

今天mdd看到这么一段话：在数论，对正整数n，欧拉函数是少于或等于n的数中与n互质的数的数目。此函数以其首名研究者欧拉命名，它又称为Euler's totient function、φ函数、欧拉商数等。 例如φ(8)=4，因为1,3,5,7均和8互质。于是他想用计算机实现欧拉函数的功能，但是他又不想去写，你能帮帮他吗？

ps: 互质（relatively primeì）又叫互素。若N个整数的最大公因数是1，则称这N个整数互质。

输入

有多组测试数据组数小于1003，
每组测试数据有一个整数n(0<n<=65535^2+1).

输出

输出欧拉函数φ(n)的值。

样例输入

2
6
46

样例输出

1
2
22

欧拉函数问题：

通式：&(n)=n*(1-1/p1)*(1-1/p2)*(1-1/p3).....(1-1/pn);p是n的素因子

#include<stdio.h>
#include<math.h>
int b[1000];
int main()
{
    long long int n,i,j;
    while(scanf("%lld",&n)!=-1)
    {
        long long int g;
        g=n;
         int m=0;//求素因子
		for(i=2;i*i<=n;i++)
        {
			if(n&&n%i==0)
			{
				b[m++]=i;
				while(n&&n%i==0)
					n/=i;
			}
        }
        if(n!=1)
        {
            b[m++]=n;
        }
        long long int sum=1,ni=1;
        for(long long int i=0; i<m; i++)
        {
            sum=sum*(b[i]-1);//把上式画简了
            ni=ni*b[i];
        }
        printf("%lld\n",sum*g/ni);
    }
}