题目
https://www.luogu.com.cn/problem/P6091
思路
首先引入阶的概念:
设 m > 1 m>1 m>1,且 g c d ( a , m ) = 1 gcd(a,m)=1 gcd(a,m)=1,那么使得 a r ≡ 1 ( m o d m ) a^r≡1(mod \ m) ar≡1(mod m) 成立的最小的正整数 r r r 称为 a a a 对模 m m m 的阶,记为 δ m ( a ) \delta_m(a) δm(a)
性质1:若 m > 1 m>1 m>1 并且 g c d ( a , m ) = 1 gcd(a,m)=1 gcd(a,m)=1 ,又满足 a n ≡ 1 ( m o d m ) a^n≡1(mod\ m) an≡1(mod m),那么 δ m ( a ) ∣ n \delta_m(a)|n δm(a)∣n
性质2:由欧拉定理和性质1得 δ m ( a ) ∣ ϕ ( m ) \delta_m(a)|\phi(m) δm(a)∣ϕ(m)
原根的定义:
设 m m m 为正整数, a a a 为整数,若 δ m ( a ) = ϕ ( m ) \delta_m(a)=\phi(m) δm(a)=ϕ(m) 则称 a a a 是 m m m 的一个原根
性质1:一个正整数有原根的充要条件: m = 2 、 4 、 p 、 2 p m=2、4、p、2p m=2、4、p、2p 其中 p p p 是奇质数
性质2:一个数的原根个数:
ϕ
(
ϕ
(
m
)
)
\phi(\phi(m))
ϕ(ϕ(m)) 且每个原根的形式为
g
p
g^p
gp 且
g
c
d
(
ϕ
(
m
)
,
p
)
=
1
gcd(\phi(m),p)=1
gcd(ϕ(m),p)=1
性质3:
g
,
g
2
,
g
3
…
…
g
ϕ
(
m
)
g,g^2,g^3……g^{\phi(m)}
g,g2,g3……gϕ(m) 是与
m
m
m 互质且小于的数的排列。
找原根
先将
ϕ
(
m
)
\phi(m)
ϕ(m) 分解质因数,枚举g
如果g是原根,那么
g
ϕ
(
m
)
=
1
g^{\phi(m)}=1
gϕ(m)=1 且恒满足
g
ϕ
(
m
)
p
i
!
=
1
g^{\frac{\phi(m)}{p_i}}!=1
gpiϕ(m)!=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");
}
}