51Nod 1135 原根

http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1135

分析：

  

     元根的定义

设m是正整数，a是整数，若a模m的阶等于φ(m)，则称a为模m的一个原根。（其中φ(m)表示m的欧拉函数）

假设一个数g是P的原根，那么g^i mod P的结果两两不同,且有 1<g<P, 0<i<P,归根到底就是g^(P-1) = 1 (mod P)当且仅当指数为P-1的时候成立.(这里P是素数).

简单来说，g^i mod p ≠ g^j mod p （p为素数）

其中i≠j且i, j介于1至(p-1)之间

则g为p的原根。

求原根目前的做法只能是从2开始枚举，然后暴力判断g^(P-1) = 1 (mod P)是否当且仅当指数为P-1的时候成立

而由于原根一般都不大，所以可以暴力得到.

元根的性质

1）可以证明，如果正整数(a,m) = 1和正整数 d 满足a^d≡1(mod m)，则 d 整除 φ(m)。因此Ordm(a)整除φ(m)。在例子中，当a= 3时，我们仅需要验证 3 的 1 、2、3 和 6 次方模 7 的余数即可。

2）记δ = Ordm(a)，则a^1，……a^(δ-1)模 m 两两不同余。因此当a是模m的原根时，a^0,a^1，……a^(δ-1)构成模 m 的简化剩余系。

3）模m有原根的充要条件是m= 1,2,4,p,2p,p^n，其中p是奇质数，n是任意正整数。

4）对正整数(a,m) = 1，如果 a 是模 m 的原根，那么 a 是整数模n乘法群（即加法群 Z/mZ的可逆元，也就是所有与 m 互素的正整数构成的等价类构成的乘法群）Zn的一个生成元。由于Zn有 φ(m)个元素，而它的生成元的个数就是它的可逆元个数，即 φ(φ(m))个，因此当模m有原根时，它有φ(φ(m))个原根。（来自百度百科）

AC代码：

#include <iostream>
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <vector>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include<list>
#include <bitset>
#include <climits>
#include <algorithm>
#define gcd(a,b) __gcd(a,b)
#define mset(a,n) memset(a,n,sizeof(a))
#define FIN	freopen("input.txt","r",stdin)
#define FOUT 	freopen("output.txt","w",stdout)
typedef long long LL;
const LL mod=1e9+7;
const int INF=0x3f3f3f3f;
const double PI=acos(-1.0);
using namespace std;
LL prime[100005];
int len;
LL p;
void prime_div(LL n){//          质因数
    for (int i=2;i<=(LL)sqrt(n);i++){
        if (n%i==0){
            prime[len++]=i;
            while (n%i==0) n/=i;
        }
    }
    if (n>1) prime[len++]=n;
}
LL Power(LL a,LL b){
    LL temp=1;
    while (b){
        if(b%2)
            temp=(temp*a)%p;
        b/=2;
        a=a*a%p;
    }
    return temp;
}
int solve(int p){
    len=0;
    prime_div(p-1);
    for (int i=2;i<=(int)(p);i++){
        if (Power(i,p-1)%p==1){
            int flag=1;
            for (int j=0;j<len;j++){
                if (Power(i,(p-1)/prime[j])%p==1){
                    flag=0;break;
                }
            }
            if (flag){
                return i;
            }
        }
    }
}
int main (){
    scanf ("%lld",&p);
    printf ("%d\n",solve(p));
    return 0;
}