题目链接:http://acm.nyist.net/JudgeOnline/problem.php?pid=479
解题思路:
欧拉函数打表即可。
n的欧拉函数计算方法如下:
先将n素分解,之后取全部素因子。则n的欧拉函数就是n*(1-/1p1)*(1-1/p2)*.....
这样,当n比较小的时候,我们可以通过打表的方法计算,先将f[i]初始化为i,然后从2开始向后遍历,如果该数是素数,就把这个数和这个数的全部倍数更新为f[j]*(1-1/该素因子),因为我们对所有的素数只遍历一次,所以我们需要遇到一个素数时,就把1-n中所有有这个素因子的数全部更新一次。
代码如下:
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int N = 10005;
long long euler[N];
void Fast_Euler()
{
for(int i = 1; i < N; ++i) //初始化
euler[i] = i;
for(int i = 2; i < N; ++i)
{
if(euler[i] == i)
{
for(int j = i; j < N; j += i)
euler[j] = euler[j] / i * (i - 1);
}
}
}
int main()
{
int n;
Fast_Euler();
while(~scanf("%d", &n))
printf("%lld\n", euler[n]);
}