51nodoj 1135 （原根）

1135 原根
基准时间限制：1 秒 空间限制：131072 KB 分值: 0 难度：基础题 收藏 关注
设m是正整数，a是整数，若a模m的阶等于φ(m)，则称a为模m的一个原根。（其中φ(m)表示m的欧拉函数）
给出1个质数P，找出P最小的原根。
Input
输入1个质数P(3 <= P <= 10^9)
Output
输出P最小的原根。
Input示例
3
Output示例
2

推荐大神博客好了， 自己不会写，也不算太懂这里写链接内容

#include <cstdio>#include <cstring>#define M 10010#define LL long longint prime[M], cnt;void getprime(int n){    cnt = 1;    for(int i=2; i<=n; i++)    {        if(n % i == 0)        {            prime[cnt++] = i;            while(n % i == 0)            {                n /= i;            }        }    }    if(n > 1)    {        prime[cnt++] = n;    }}LL quickpow(LL n, LL m, LL mod){    n %= mod;    LL res = 1;    while(m)    {        if(m & 1)        {            res = res * n % mod;        }        n = n * n % mod;        m >>= 1;    }    return res;}int main(){    LL n;    bool ok;    while(scanf("%lld", &n) != EOF)    {        getprime(n-1);//求出 n-1 的所有的质因数        for(LL i=2; i<=n; i++)        {            ok = 1;            for(LL j=1; j<cnt; j++)            {                if(quickpow(i, (n-1)/prime[j], n) == 1)                //快速幂，判断是否有不符合的                {                    ok = 0;                    break;                }            }            if(ok)//有成立的就直接输出            {                printf("%lld\n", i);                break;            }        }    }    return 0;}