阶,原根,指标和N次剩余

阶,原根,指标和$N$次剩余

参考博客

阶

阶的定义

设$a$,$p$是互质整数,存在$x$满足$a^x\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$,最小的满足条件的正整数$x$称为$a$模$p$的阶,记为$Or{d}_p(a)$

性质

性质一

设$a$,$p$是互质整数,存在$x^{\prime }$满足$a^{x^{\prime }}\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$,则$Or{d}_p(a)|x^{\prime }$

证明

反证法

设$a$,$p$是互质整数,存在整数$x^{\prime }$满足$a^{x^{\prime }}\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$且$Or{d}_p(a)\nmid x^{\prime }$

设$x\equiv x^{\prime }\mathrm{m}\mathrm{o}\mathrm{d}Or{d}_p(a)$且$1\le x<Or{d}_p(a)$

显然$a^x\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$,与$Ord$的定义矛盾

性质二

设$a$,$p$是互质整数,$Or{d}_p(a)|\phi (p)$

证明

$a^{\phi (p)}\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$

原根

原根的定义

设$a$,$p$是互质整数,$Or{d}_p(a)=\phi (p)$,则称$a$是模$p$的一个原根

性质

性质一

只有$2,4,p^k,2p^k$有原根($p$为奇素数)

性质二

一个数$p$若存在原根,则有$\phi (\phi (p))$个原根

证明

假设已经找到了一个原根$g$,那么对于所有$\gcd (x,\phi (p))=1$的$x$,$g^x$也是$p$的原根

性质三

若$a$是$p$的原根,$a^0,a^1...a^{\phi (p)-1}$两两不同,即构成了模$p$的简化剩余系

如何求原根

先求最小正原根

从小到大枚举$g$并检验$g$是否是$p$的原根,枚举次数是$p^{0.25}$级别的

如何检验?

满足$\forall i<varphi(p)$,$g^i\ne 1\mathrm{m}\mathrm{o}\mathrm{d}p$,且$g^{\phi (p)}\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$

根据阶的性质,只需要检验$\phi (p)$的所有真因数即可,考虑检验所有真因数的倍数

设$\phi (p)={\prod }_i{q}_i^{c_i}$

$\forall i$,检验$\frac{\phi (p)}{q_i}$即可

指标

定义

设$g$是$p$的原根,$a$,$p$是互质整数,存在唯一的整数$x$满足$g^x\equiv a\mathrm{m}\mathrm{o}\mathrm{d}p$,称$x$为以$g$为底对模$p$的一个指标,记为$x=in{d}_{g}a$

性质

性质一

若$a\equiv b\mathrm{m}\mathrm{o}\mathrm{d}p$,则$in{d}_{g}a\equiv in{d}_{g}b\mathrm{m}\mathrm{o}\mathrm{d}\phi (p)$

性质二

$in{d}_g(ab)\equiv in{d}_g(a)+in{d}_g(b)\mathrm{m}\mathrm{o}\mathrm{d}\phi (p)$

性质三

$in{d}_g(a^k)\equiv k\times in{d}_g(a)\mathrm{m}\mathrm{o}\mathrm{d}\phi (p)$

如何求指标

就是求离散对数,直接$BSGS$

BSGS

求$x$满足$A^x\equiv B\mathrm{m}\mathrm{o}\mathrm{d}p$,满足$\gcd (A,p)=1$

显然$x\in [0,p)$

设$M=\lceil \sqrt{(}p)\rceil$

设$x=i\times M-c,c\in [0,M)$

$$A^{i\times M-c}\equiv B\mathrm{m}\mathrm{o}\mathrm{d}p$$
$$A^{i\times M}\equiv B\times A^c\mathrm{m}\mathrm{o}\mathrm{d}p$$
枚举$c$,把$B\times A^c$插进一个map,最后枚举$i$求解即可

$N$次剩余

求所有$x$满足$x^A\equiv B\mathrm{m}\mathrm{o}\mathrm{d}C$

分解问题

考虑将$C$质因数分解

子问题为求所有$x$满足$x^A\equiv B\mathrm{m}\mathrm{o}\mathrm{d}{p}_i^{{\alpha }_i}$

每个方程可能有多种解,每个方程中选一个解,这样每组选择方法唯一对应一种答案,用$ExCrt/Crt$合并即可

处理每个$x^A\equiv B\mathrm{m}\mathrm{o}\mathrm{d}{p}^{\alpha }$

若$\gcd (B,p^{\alpha })=1$

设$M=p^{\alpha }$,$g$是$M$的原根,$y=in{d}_{M}x,z=in{d}_{M}B$
$$A\times y\equiv z\mathrm{m}\mathrm{o}\mathrm{d}\phi (M)$$
若$z\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}\phi (M)$
,设$qwq=\gcd (A,\phi (M))$

$$y\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}\frac{\phi (M)}{qwq}$$
$y\in [0,\phi (M))$,共有$qwq$组解

否则,设$qwq=\gcd (A,z,\phi (M))$
$$y\equiv \frac{z}{qwq}\times (\frac{A}{qwq}{)}^{-1}\mathrm{m}\mathrm{o}\mathrm{d}\frac{\phi (M)}{qwq}$$
$$x\equiv g^y\mathrm{m}\mathrm{o}\mathrm{d}M$$
$y\in [0,\phi (M))$,有$qwq$组解

若$\gcd (B,p^{\alpha })=p^q(1\le q<\alpha )$

$x^A\equiv B\mathrm{m}\mathrm{o}\mathrm{d}{p}^{\alpha }$

设$B=z\times p^q$,$x=y\times p^a$

$$(y\times p^a{)}^A\equiv z\times p^q\mathrm{m}\mathrm{o}\mathrm{d}{p}^{\alpha }$$
$$y^A\times p^{aA}\equiv z\times p^q\mathrm{m}\mathrm{o}\mathrm{d}{p}^{\alpha }$$
显然,$A|q$,反之无解

求解$y^A\equiv z\mathrm{m}\mathrm{o}\mathrm{d}{p}^{\alpha -q}$,$x=y\times p^{\frac{q}{A}}$即可

每组$y^A\equiv z\mathrm{m}\mathrm{o}\mathrm{d}{p}^{\alpha -\frac{q}{A}}$的解唯一对应一组$x^A\equiv B\mathrm{m}\mathrm{o}\mathrm{d}{p}^{\alpha }$的解

所以每个$y$分别对应$p^{q-\frac{q}{A}}$个$x$

若$B\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}{p}^{\alpha }$

$$x\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}{p}^{\lceil \frac{\alpha }{A}\rceil }$$

$x\in [0,p^{\alpha })$,共有$p^{\alpha -\lceil \frac{\alpha }{A}\rceil }$组解

若$p$为$2$

$x^A\equiv B\mathrm{m}\mathrm{o}\mathrm{d}{2}^{\alpha }$

求出$y^A\equiv B\mathrm{m}\mathrm{o}\mathrm{d}{2}^{\alpha -1}$,$x=y$或$x=y+2^{\alpha -1}$

如何用$ExCrt$合并答案?

以上选出的答案都为同余方程的形式,考虑如何合并两个同余方程

$$\{\begin{array}{r}x\equiv c_1\mathrm{m}\mathrm{o}\mathrm{d}{p}_1\\ x\equiv c_2\mathrm{m}\mathrm{o}\mathrm{d}{p}_2\end{array}$$
$$c_1+p_1\times k_1=c_2+p_2\times k_2$$
$$p_1\times k_1-p_2\times k_2\equiv c_2-c_1$$
设$qwq=\gcd (p_1,p_2)$,若$qwq\nmid (c_2-c_1)$,无解

否则,令$P_1=\frac{p_1}{qwq},P_2=\frac{p_2}{qwq},C=\frac{c_2-c_1}{qwq}$

求$k_1$,$k_1\equiv \frac{C}{P_1}\mathrm{m}\mathrm{o}\mathrm{d}{P}_2$

$$x_0=c_1+k_1\times p_1$$
$$=c_1+(k_1+\frac{p_2}{qwq}{k}_2)\times p_1$$
$$x\equiv k_1\times p_1+c_1\mathrm{m}\mathrm{o}\mathrm{d}\frac{p_1{p}_2}{\gcd (p_1,p_2)}$$

代码

#include<bits/stdc++.h>
using namespace std;
#define Ri register
template<typename T>inline T read(Ri T& t)
{Ri T f=1;Ri char ch=getchar();t=0;
while(ch<'0'||ch>'9'){if(ch=='-')f=-f;ch=getchar();}
while(ch>='0'&&ch<='9')t=t*10+ch-'0',ch=getchar();t*=f;return t;}
template <typename T,typename...Args>
inline void read(T&t,Args&...args)
{read(t);read(args...);}
#define ll long long
inline ll power(Ri ll x,Ri ll k,Ri ll p=1e18)
{
	ll re=1;
	for(;k;k>>=1,x=x*x%p)if(k&1)re=re*x%p;
	return re; 
}
inline int BSGS(ll b,ll n,ll p)//B^i=n/mod p
{
	if(n==1)return 0;
	map<int,int> mp;
	int m=ceil(sqrt(p))+1;
	ll t=n;
	for(int i=0;i<=m;i++,t=t*b%p)mp[t]=i;
	t=power(b,m,p);
	for(ll i=1,s=t;i<=m;i++,s=s*t%p)
		if(mp.count(s))
			return i*m-mp[s];
}
inline ll gcd(ll a,ll b)
{
	if(b==0)return a;
	return gcd(b,a%b); 
}
int pri[1000005],pri_top;
bool check(int n,int g,int varphi)
{
	if(power(g,varphi,n)!=1)return 0;
	for(int i=1;i<=pri_top;i++)
		if(power(g,varphi/pri[i],n)==1)
			return 0;
	return 1;
}
int get_g(int n,int varphi)
{
	int phi=varphi;
	pri_top=0;
	for(int i=2;i*i<=phi;i++)
	{
		if(phi%i)continue;
		pri[++pri_top]=i;
		while(phi%i==0)phi/=i;
	}
	if(phi!=1)pri[++pri_top]=phi;
//		printf("varphi=%d,phi:",varphi[n]);for(int i=1;i<=pri_top;i++)printf("%d ",pri[i]);printf("\n");
	int g=1;while(check(n,g,varphi)==0)g++;
	return g;
}
#define Pii pair<int,int>
#define Vii vector<pair<int,int> >
int Exgcd(int a,int b,ll & x,ll & y)
{
	if(b==0)return x=1,y=0,a;
	int re=Exgcd(b,a%b,x,y);
	ll qwq=x;
	x=y;y=qwq-(a/b)*y;
	return re;
}
int Inv(int a,int p)
{
	ll x,y;
	Exgcd(a,p,x,y);
	return (x%p+p)%p;
}
Vii SOLVE(int A,int B,int P,int Alpha)//x^A=B mod P^Alpha
{
	int M=power(P,Alpha);B%=M;	
//	printf("SOLVE:x^%d=%d mod %d^%d\n",A,B,P,Alpha);
	if(P==2)
	{
		if(Alpha==0)return {make_pair(0,1)};
		Vii pre=SOLVE(A,B,2,Alpha-1);
		Vii re;
		for(auto i:pre)
		{
			ll qwq=i.first;
			if(power(qwq,A,M)==B)re.push_back(make_pair(qwq,M));
			qwq=(qwq+M/2)%M;
			if(power(qwq,A,M)==B)re.push_back(make_pair(qwq,M));
		}
		return re;
	}
	Vii re;int var_M=M/P*(P-1);
	if(B==0)
	{
		int qwq=power(P,ceil((double)Alpha/A));
		for(int i=0;i*qwq<M;i++)
			re.push_back(make_pair(i*qwq,M));
		return re;
	}

	if(gcd(B,M)==1)
	{
		int g=get_g(M,var_M);
//		printf("g=%d var_M=%d\n",g,var_M);
		int z=BSGS(g,B,M);
		if(z==0)
		{
			int qwq=gcd(A,var_M);
			for(int i=0;i<var_M;i+=var_M/qwq)
				re.push_back(make_pair(power(g,i,M),M));
			return re; 
		}
		int qwq=gcd(A,gcd(z,var_M));
//		printf("%d^%d=%d mod %d\n",g,z,B,M);
		if(gcd(A/qwq,var_M/qwq)!=1)return {};
		int y=(ll)z/qwq*Inv(A/qwq,var_M/qwq)%(var_M/qwq);
		for(int i=0;i<qwq;i++)
//			printf("%d^%d\n",g,y+var_M/qwq*i),
			re.push_back(make_pair(power(g,y+var_M/qwq*i,M),M));
	}
	else
	{
		int gc=gcd(B,M);
		int q=1,qaq=P;
		while(qaq!=gc)
			qaq=qaq*P,q++;
		if(q%A)return {};
		int z=B/qaq;
		Vii pre=SOLVE(A,z,P,Alpha-q);
		int qwq=power(P,q/A,M);
		int qAq=power(P,q-q/A);
		int QaQ=power(P,Alpha-q);
		for(auto i:pre)
			for(ll j=0;j<qAq;j++)
				re.push_back(make_pair((i.first+j*QaQ%M)*qwq%M,M));
	}
//	printf("SOLVE:x^%d=%d mod %d^%d\n",A,B,P,Alpha);
//	for(auto i:re)printf("%lld mod %lld\n",i.first,i.second);
	return re;
}
Vii solution[1005];
int solution_top;
vector<long long> ans;
bool without_ans;
void DFS(int now,ll m,ll p)
{
	if(without_ans)return;
	if(now>solution_top)return ans.push_back(m),void();
	for(auto i:solution[now])
	{
		ll qwq=gcd(p,i.second);
		ll C=(i.first-m%i.second+i.second)%i.second;
		if(C%qwq)return without_ans=1,ans={},void();
		ll P1=p/qwq,P2=i.second/qwq,c=C/qwq;
		ll k1=c*Inv(P1,P2)%P2;
		ll p2=p*i.second/gcd(p,i.second);
		DFS(now+1,(k1*p+m)%p2,p2);
	}
}
void solve(int A,int B,int C)//x^A=b mod C
{	
	int M=C;solution_top=0;
	for(int i=2;i*i<=M;i++)
		if(M%i==0)
		{
			int sum=0;
			while(M%i==0)
				M/=i,sum++;
			solution[++solution_top]=SOLVE(A,B,i,sum);
//			printf("Devided Question:x^%d=%d mod %d^%d\n",A,B,i,sum);
//			for(auto i:solution[solution_top])printf("%lld mod %lld\n",i.first,i.second);
		}
	if(M!=1)solution[++solution_top]=SOLVE(A,B,M,1);
//	printf("Devided Question:x^%d=%d mod %d^%d\n",A,B,M,1);
//	for(auto i:solution[solution_top])printf("%lld mod %lld\n",i.first,i.second);
	ans.clear();	
	without_ans=0;
	DFS(1,0,1);
	sort(ans.begin(),ans.end());
	printf("%d\n",ans.size());
	if(ans.size()==0)return;
	for(auto i:ans)printf("%lld ",i);
	printf("\n");
}
int main()
{
	int T;
	read(T);
	while(T--)
	{
		int A,B,C;
		read(A,C,B);
		solve(A,B,C);
	}
	return 0
}