【模板】原根

题目

https://www.luogu.com.cn/problem/P6091

思路

首先引入阶的概念：

设 $m>1$，且 $gcd(a,m)=1$，那么使得 $a^r\equiv 1(modm)$ 成立的最小的正整数 $r$ 称为 $a$ 对模 $m$ 的阶，记为 ${\delta }_m(a)$

性质1：若 $m>1$ 并且 $gcd(a,m)=1$ ，又满足 $a^n\equiv 1(modm)$，那么 ${\delta }_m(a)|n$

性质2：由欧拉定理和性质1得 ${\delta }_m(a)|\varphi (m)$

原根的定义：

设 $m$ 为正整数，$a$ 为整数，若 ${\delta }_m(a)=\varphi (m)$ 则称 $a$ 是 $m$ 的一个原根

性质1：一个正整数有原根的充要条件：$m=2、4、p、2p$ 其中 $p$ 是奇质数

性质2：一个数的原根个数：$\varphi (\varphi (m))$ 且每个原根的形式为 $g^p$ 且 $gcd(\varphi (m),p)=1$
性质3：$g,g^2,g^3\dots \dots g^{\varphi (m)}$ 是与 $m$ 互质且小于的数的排列。

找原根

先将 $\varphi (m)$ 分解质因数，枚举g
如果g是原根，那么 $g^{\varphi (m)}=1$ 且恒满足 $g^{\frac{\varphi (m)}{p_i}}!=1$

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1e6+77;
int phi[N],p[N],b[N],rt[N],f[N],cnt,tot,a[N];
int gcd(int x,int y)
{
	return y==0?x:gcd(y,x%y);
}
void init(int n)
{
	phi[1]=1;
	for(int i=2; i<=n; i++)
	{
		if(!b[i]) p[++tot]=i,phi[i]=i-1;
		for(int j=1; j<=tot&&i*p[j]<=n; j++)
		{
			b[i*p[j]]=1;
			if(i%p[j]==0)
			{
				phi[i*p[j]]=phi[i]*p[j]; break;
			}
			phi[i*p[j]]=phi[i]*(p[j]-1);
		}
	}
//	for(int i=1; i<=10; i++) printf("\n%d ",phi[i]);
	rt[2]=rt[4]=1;
	for(int i=2; i<=tot; i++)
	{
		for(int j=1; (1ll*j*p[i])<=n; j*=p[i]) rt[j*p[i]]=1;
		for(int j=2; (1ll*j*p[i])<=n; j*=p[i]) rt[j*p[i]]=1;
	}
}
void getrt(int n,int x)
{
	int p=1;
	for(int i=1; i<=phi[n]; i++)
	{
		p=p*x%n;
		if(gcd(i,phi[n])==1) a[++a[0]]=p;
	}
}
ll power(ll x,ll t,ll mod)
{
	ll b=1;
	while(t)
	{
		if(t&1) b=b*x%mod; x=x*x%mod; t>>=1;
	}
	return b;
}
bool check(int x,int n)
{
	if(power(x,phi[n],n)!=1) return 0;
	for(int i=1; i<=cnt; i++) 
		if(power(x,phi[n]/f[i],n)==1) return 0;
	return 1;
}
int T,n,d;
int main()
{
	scanf("%d",&T);
	init(1000005);
	while(T--)
	{
		scanf("%d%d",&n,&d);
		if(rt[n])
		{
			cnt=0; int t=phi[n];
			for(int i=2; i*i<=t; i++) if(t%i==0)
			{
				f[++cnt]=i;
				while(t%i==0) t/=i;
			}
			if(t>1) f[++cnt]=t; int yjy=0;
			for(int i=1; i<n; i++) if(check(i,n))
			{
				yjy=i; break;
			}
			a[0]=0;
			getrt(n,yjy);
			printf("%d\n",a[0]);
			sort(a+1,a+a[0]+1);
			for(int i=d; i<=a[0]; i+=d) printf("%d ",a[i]);
		}
		else printf("0\n");
		printf("\n");
	}
}