POJ 1284 Primitive Roots(欧拉函数)

Description
定义一个数p的原根x:{x^i=1,i=1,2,…,p-1}={1,2,…,p-1},先给出素数p,求p的原根个数
Input
多组输入,每组用例占一行为一素数p,以文件尾结束输入
Output
对于每组用例,输出p的原根个数
Sample Input
23
31
79
Sample Output
10
8
24
Solution
此处用到一个定理:如果p有原根,则它有euler[euler[p]]个原根,p是素数时,由于euler[p]=p-1,所以有euler[p-1]个原根
Code

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define maxn 66666
#define INF 0x3f3f3f3f
typedef long long ll;
int euler[maxn],prime[maxn],cnt;
void Get_Euler(int n)//求n以内所有数的欧拉函数值 
{
    memset(euler,0,sizeof(euler));
    euler[1]=1;
    cnt=0; 
    for(int i=2;i<=n;i++)
    {
        if(!euler[i])
        {
            euler[i]=i-1;
            prime[cnt++]=i;
        } 
        for(int j=0;j<cnt&&prime[j]*i<=n;j++)
        {
            if(i%prime[j]) 
                euler[prime[j]*i]=euler[i]*(prime[j]-1);
            else
             {
                euler[prime[j]*i]=euler[i]*prime[j];
                break;
            }
        }
    }
}
int main()
{
    Get_Euler(maxn);
    int n;
    while(~scanf("%d",&n))
        printf("%d\n",euler[n-1]);
    return 0;
}

你可能感兴趣的:(POJ 1284 Primitive Roots(欧拉函数))