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;
}