SGU 102 Coprimes（欧拉函数）

Description
给出一整数n，问1~n中有多少个数与n互素
Input
一整数n(1<=n<10000)
Output
输出1~n中与n互素的数的个数
Sample Input
9
Sample Output
6
Solution
求n的欧拉函数值，对n质因子分解得n=p1^a1*p2^a2*…pm^am，则euler(n)=n * (1-1/p1) * (1-1/p2) * … * (1-1/pm)
Code

#include<cstdio>
#include<iostream>
using namespace std;
int get_euler(int n)
{
    int ans=n;
    for(int 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()
{
    int n;
    while(~scanf("%d",&n))
        printf("%d\n",get_euler(n));
    return 0;
}