原根的求解及应用

niiick

已于 2022-10-23 11:54:45 修改

阅读量1.8k

点赞数

分类专栏：数学;数论文章标签：算法

于 2021-10-13 17:10:16 首次发布

本文链接：https://blog.csdn.net/niiick/article/details/120746931

版权

数学;数论专栏收录该内容

36 篇文章 3 订阅

订阅专栏

原根的定义

首先引入数论中阶的定义

设 $a\in Z, m\in N^+$ ，若 $g c d (a, m) = 1$ ，则称使得 $a^n\equiv 1(\mod m)$ 成立的最小正整数 $n$ 为 $a$ 对膜 $m$ 的阶

接下来是原根的定义

设 $a\in Z, m\in N^+$ ，若 $g c d (a, m) = 1$ ，且 $a$ 膜 $m$ 的阶为 $\varphi(m)$ ，则 $a$ 为膜 $m$ 的原根

原根的求解

定理1：若 $m\notin \{2,4\}\bigcup\{p^\alpha,2p^\alpha\}$ （其中p为奇素数，α为整数），则膜m不存在原根

由该定理，我们可以先求素数表，然后预处理出每个数是否有原根

定理2：膜m的原根数量级别为 $O(m^{0.25})$

定理3：若膜m存在原根，则其有 $\varphi(\varphi(m))$ 个原根

假设我们已求得膜 $m$ 的最小原根 $g$ ，那么膜 $m$ 的所有原根就是 $\{g^k \ |\ gcd(k,\varphi(m)=1)\}$
所以我们可以直接从小到大枚举每个数判断其是否为膜m的原根
先得到膜m的最小原根g，再求g的幂次方以求得所有原根

为判断g是否为膜m的原根，在引入一个定理

定理4：设 $a\in Z, m\in N^+$ ，若 $g c d (a, m) = 1$ 且 $a^k\equiv 1(\mod m)$ ，则 $k|\varphi(m)$

易知若 $g$ 是膜 $m$ 的原根，则对所有 $k<\varphi(m)$ 均有 $g^k \mod m \not= 1$
结合定理4，我们只需要检查 $\varphi(m)$ 的真因数d，若存在 $g^d\equiv 1(\mod m)$ 则g不是原根

进一步的，我们只需要对 $\varphi(m)$ 质因数分解，检查所有 $d=\frac{\varphi(m)}{p_i}$ 即可

洛谷P6091 【模板】原根

const int maxn=1000010;
lt prim[maxn],phi[maxn],vis[maxn],cnt;
int has[maxn];
lt fac[maxn],facN;
lt rt[maxn],rtN;

lt gcd(lt a,lt b){ return b==0?a:gcd(b,a%b);}

lt qpow(lt a,lt k,lt mod)
{
	lt res=1;
	while(k){
		if(k&1) res=(res*a)%mod;
		a=(a*a)%mod; k>>=1;
	}
	return res;
}

void qPhi(lt n)
{
	phi[1]=1;
	for(lt i=2;i<=n;++i)
	{
		if(!vis[i]){ prim[++cnt]=i; phi[i]=i-1;}
		for(int j=1;j<=cnt;++j)
		{
			if(i*prim[j]>n) break;
			vis[i*prim[j]]=1;
			if(i%prim[j]==0){ phi[i*prim[j]]=phi[i]*prim[j]; break;}
			else phi[i*prim[j]]=phi[i]*phi[prim[j]];
		}
	}
}

void preWork(int lim)
{
	has[2]=has[4]=1;
	for(int i=2;i<=cnt;++i)
	{
		for(lt j=prim[i]; j<=lim; j*=prim[i]) has[j]=1; // p^α
		for(lt j=prim[i]*2; j<=lim; j*=prim[i]) has[j]=1; // 2p^α
	}
}

void div(int x)
{
	facN=0;
	for(int i=2;i*i<=x;++i)
	{
		if(x%i==0)
		{
			while(x%i==0) x/=i;
			fac[++facN]=i;
		}
	}
	if(x>1) fac[++facN]=x;
}

bool check(int g,int n)
{
	if(qpow(g,phi[n],n)!=1) return false; // g^φ(n) ≡ 1 (mod n)
	
	for(int i=1;i<=facN;++i){
		if(qpow(g,phi[n]/fac[i],n)==1) return false;
	}
	
	return true;
}

int minRt(int n)
{
	div(phi[n]); // 对φ(n)质因数分解并将质因子存在fac[]里
	for(int i=1;i<=n;++i)
	if(check(i,n)) return i;
	
	return 0;
}

void allRt(int g,int n)
{
	rtN=0;
	lt t=1;
	for(int i=1;i<=phi[n];++i)
	{
		t=t*g%n;
		if(gcd(i,phi[n])==1) rt[++rtN]=t;
	}
}

int main()
{
	qPhi(maxn-10); // 质数表和phi
    preWork(maxn-10); // 预处理每个数是否有原根
	
	int T=read();
    while(T--)
    {
    	int n=read(),d=read();
    	if(!has[n])
		{
			printf("0\n\n");
    		continue;
		}
		
		lt g=minRt(n); // 求最小原根
		allRt(g,n); // 求所有原根
		
		printf("%d\n",rtN);

		sort(rt+1,rt+1+rtN);
		for(int i=d;i<=rtN;i+=d)
		printf("%d ",rt[i]);
		
		printf("\n");
	}
	return 0;
}