POJ 2407 Relatives(欧拉函数)

Description
求一个数的欧拉函数值
Input
多组输入,每组一个正整数n,以0结束输入
Output
对于每组输入,输出n的欧拉函数值
Sample Input
7
12
0
Sample Output
6
4
Solution
直接套欧拉函数公式,即将n素分解后有n=p1^k1*p2^k2*…pm^km,则euler(n)=n(1-1/p1)(1-1/p2)…*(1-1/pm)
Code

#include<cstdio>
#include<iostream>
using namespace std;
#define maxn 1111111111
typedef long long ll;
ll Get_Euler(ll n)
{
    ll ans=n;
    for(ll i=2;i*i<=n;i++)
    {
        if(n%i==0)
        {
            ans=ans/i*(i-1);
            while(n%i==0)
                n/=i;
        }
    }
    if(n>1)
        ans=ans/n*(n-1);
    return ans;
}
int main()
{
    ll n;
    while(scanf("%lld",&n),n)
        printf("%lld\n",Get_Euler(n));
    return 0;
}

你可能感兴趣的:(POJ 2407 Relatives(欧拉函数))