51nod 1135 原根 (数论)

题目链接:
原根例题

求模素数 P 原根的方法:对这里写图片描述素因子分解,即这里写图片描述 P1 的标准分解式,若恒有这里写图片描述
成立,则这里写图片描述就是这里写图片描述的原根。(对于合数求原根,只需把这里写图片描述换成这里写图片描述即可)。

代码:

#include 
using namespace std;
int p[100010];
int t = 0;
//http://blog.csdn.net/acdreamers/article/details/8883285
void divide(long long n)
{
    for(int i = 2; i * i <= n; i++){
        if(n % i == 0){
            p[t ++] = i;
            while(n % i == 0) n /= i;
        }
    }
}

int solve(long long x,long long n,long long m)
{
    long long res = 1;
    while(n > 0){
        if(n & 1) res = (res * x) % m;
        x = (x * x) % m;
        n >>= 1;
    }
    return res;
}

int main()
{
    long long a;
    cin >> a;
    divide(a - 1);//求a-1所有的质因子 
    for(int i = 2; i < a; i++)
    {
        bool flag = false;
        for(int j = 0; j < t; j++)
        {
            if(solve(i,(a - 1) / p[j], a ) == 1)//判断是否是原根 
            {
                flag = true;
                break;
            }
        }
        if(!flag){
            cout << i << endl;
            break;
        }
    }
    return 0;
}

你可能感兴趣的:(ACM_快速幂,ACM_其他数论,51NOD)