51nod 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 <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;
#define LL __int64
#define PB push_back
LL mul(LL a,LL b,LL m)
{
LL ret = 0;
a %= m;
while(b)
{
if(b & 1) ret = (ret + a) % m;
a = (a + a) % m;
b >>= 1;
}
return ret;
}
LL a_b_MOD_c(LL a,LL b,LL m)
{
LL ret = 1;
a %= m;
while(b)
{
if(b&1) ret = mul(ret,a,m);
a = mul(a,a,m);
b >>= 1;
}
return ret;
}

vector<LL> a;
bool g_test(LL g,LL p)
{
for(LL i=0; i<a.size(); ++i)
if(a_b_MOD_c(g,(p-1)/a[i],p)==1)
return 0;
return 1;
}
LL pri_root(LL p)
{
a.clear();
LL tmp=p-1;
for(LL i=2; i<=tmp/i; ++i)
if(tmp%i==0)
{
a.push_back(i);
while(tmp%i==0)
tmp/=i;
}
if(tmp!=1)
a.push_back(tmp);
LL g=1;
while(true)
{
if(g_test(g,p))
return g;
++g;
}
}

int main()
{
LL x;
while(~scanf("%lld",&x))
cout<<pri_root(x)<<endl;
}