poj1284 Primitive Roots（原根，欧拉函数）

http://poj.org/problem?id=1284

题意：求一个奇素数的原根个数。

思路：原根定义：here。

感觉和费马小定理有点像，定义里都说“归根到底就是g^(P-1) = 1 (mod P)当且仅当指数为P-1的时候成立.(这里P是素数).”。

说白了就是用了一个定理：

一个数n有原根，那么他有phi(phi(n))个模n不同余的原根(n是否素数都可用)。

由于是素数，所以phi(n)=n-1，于是求原根就是求phi(n-1)。

貌似有证明过程，有空看看。。

#include <stdio.h>
#include <algorithm>
#include <string.h>

using namespace std;

typedef long long ll;
const int N = 66000;

int prime[N], num[N], euler[N], cnt = 0;

void table()
{
    memset(prime, 1, sizeof(prime));
    euler[1] = 1;
    for(int i = 2; i <= N; i++)
    {
        if(prime[i])//如果是素数
        {
            num[++cnt] = i;
            euler[i] = i-1;
        }
        for(int j = 1; j<=cnt && num[j]*i<N; j++)
        {
            prime[num[j]*i] = 0;
            if(i%num[j] == 0) euler[i*num[j]] = euler[i]*num[j];
            else euler[i*num[j]] = euler[i]*(num[j]-1);
        }
    }
}

int main()
{
  //  freopen("in.txt", "r", stdin);
    int n;
    table();
    while(~scanf("%d", &n))
    {
        if(n == 0) break;
        printf("%d\n", euler[n-1]);
    }
    return 0;
}