原根

首先，笔者要给大家介绍一下，什么叫原根。

设g是奇素数m的一个原根，则g¹ mod m , g² mod m , g³ mod m , … , g^m-1 mod m 两两不同。并且其结果恰好为1 ~ m-1。

---

下面是原根的求法以及其性质：

• 模数m有原根当且仅当 m = 2,4,p^a, 2*p^a这四种。其中p是奇素数。
• 求素数m原根的方法：对m-1进行质因数分解，即m-1 =
  p1^a1p2^a2p3^a3…pk^ak，若恒有g^(m-1)/pi != 1(mod m) 成立，则g就是m的原根。
  （对于合数求原根，只需把m-1换成phi(m)即可）。 聪明的读者已经发现了，phi(m) = m-1 当 m是素数。
• 原根的性质：如果一个数m有原根，其原根的个数为phi(phi(m))。

代码如下：

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e7+7;
int pri[maxn],ispri[maxn],divi[maxn];
int p = 2019;
int qpow(int a,int b)
{
	int ans = 1;
	while(b>0)
	{
		if(b&1) ans = ans*a%p;
		b>>=1;
		a = a*a%p;
	}
	return ans%p;
}
void getprime()
{
	memset(ispri,1,sizeof(ispri));
	ispri[0] = ispri[1] = 0;
	for(int i=2;i<maxn;i++)
	{
		if(ispri[i]) pri[++pri[0]] = i;
		for(int j=1;j<=pri[0] && i*pri[j]<maxn;j++)
		{
			ispri[i*pri[j]] = 0;
			if(i%pri[j]==0) break;
		}
	}
}
int findroot(int x)
{
	int tmp = x-1;
	for(int i=1;i<=pri[0];i++)
	{
		if(tmp%pri[i]==0)
		{
			divi[++divi[0]] = pri[i];
			while(tmp%pri[i]==0) tmp /= pri[i];
		}
	}
	for(int g=2;g<=x-1;g++)
	{
		bool found = false;
		for(int i=1;i<=divi[0];i++)
		{
			if(qpow(g,(x-1)/divi[i])==1) 
			{
				found = true;
				break;
			}
		}
		if(!found) return g;
	}
}
int main()
{
	getprime();
	printf("%d\n",findroot(2019));
	return 0;
}