离散对数(BSGS)及(exBSGS)

最新推荐文章于 2021-04-23 16:03:12 发布

CaptainHarryChen

最新推荐文章于 2021-04-23 16:03:12 发布

阅读量1k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： OI知识总结计算机数学

本文链接：https://blog.youkuaiyun.com/can919/article/details/80786283

OI知识总结同时被 2 个专栏收录

37 篇文章

订阅专栏

计算机数学

29 篇文章

订阅专栏

本文详细介绍了解同余方程ax≡b (mod P)的两种BSGS算法：普通BSGS算法适用于a与P互质的情况，扩展BSGS算法适用于a与P不互质的情况。文章给出了算法的具体步骤和实现代码。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

普通BSGS

问题

已知 $a, b, P$ ，且a与P互质，求解同余方程 $P)a^x \equiv b \space (mod \space P)$

算法推导

设 $⌉m=\left \lceil \sqrt P \space \right \rceil$
设 $x = i m + j$ ，即 $i=⌊xm⌋i=\left \lfloor \frac x m \right \rfloor$ ， $\space mod \space m$
所以满足 $\leq m$ ， $\leq m$
得 $P)a^{im+j} \equiv b \space (mod \space P)$
$P)a^j a^{im} \equiv b \space (mod \space P)$
$P)a^j \equiv a^{-im}b \space (mod \space P)$

算法流程

先枚举j，将得到的 $a^j$ 存入hash表；
再枚举i，计算 $a^{-im}b$ ，是否存在hash表中，如果有，则找到了解，输出即可。

为保证答案最小，必须保证先枚举j，再枚举i，且保证 $a^j$ 在hash表中j从小到大排列

代码

long long work(long long a,long long b,long long P)
{
	long long m,v,e=1,i;
	m=ceil(sqrt(P+0.5));
	
	v=pow_mod(a,m,P);
	v=inv(v);//v=a^(-im)
	
	Hash::clear();
	Hash::add(1,0);
	for(i=1;i<m;i++)
	{
		e=(e*a)%P;
		if(Hash::get(e)==-1)
			Hash::add(e,i);
	}
	for(i=0;i<m;i++)
	{
		int t=Hash::get(b);
		if(t!=-1)
			return i*m+t;
		b=(b*v)%P;
	}
	return -1;
}

模板题

POJ2417

扩展BSGS

问题

已知 $a, b, P$ ，a与P 不一定互质，求解同余方程 $P)a^x \equiv b \space (mod \space P)$

算法推导

$d = g c d (a, P)$
如果 $\space mod \space d ≠0$
显然无解

则可以得
$Pd)\frac {a^x} d \equiv \frac b d \space \left( mod \space \frac P d \right)$
$Pd)\frac {a} d \times a^{x-1} \equiv \frac b d \space \left( mod \space \frac P d \right)$

然后继续找 $d_2=gcd(a,P/d)$
$Pd×d2)\frac {a^2} {d \times {d_2}} \times a^{x-2} \equiv \frac b {d \times {d_2}} \space \left( mod \space \frac P {d \times {d_2}} \right)$

如此反复，令 $D={d_1}{d_2}{d_3}...{d_k}$
得
$PD)\frac {a^k} D \times a^{x-k} \equiv \frac b D \space \left( mod \space \frac P D \right)$

前k次暴力求解，后面 $PD\frac P D$ 与 $a^{x-k}$ 互质，然后用普通BSGS就可以做了

代码

#include<cstdio>
#include<cmath>
#include<cstring>
const long long HASHSIZE=31643,MAXNODE=31650;

namespace Hash
{
	long long id[MAXNODE],val[MAXNODE],nxt[MAXNODE];
	long long head[HASHSIZE],ncnt;
	void Clear()
	{
		ncnt=0;
		memset(head,-1,sizeof head);
	}
	void Insert(long long i,long long v)
	{
		id[ncnt]=i;
		val[ncnt]=v;
		nxt[ncnt]=head[i%HASHSIZE];
		head[i%HASHSIZE]=ncnt++;
	}
	long long Get(long long i)
	{
		for(long long u=head[i%HASHSIZE];u!=-1;u=nxt[u])
			if(id[u]==i)
				return val[u];
		return -1;
	}
}

long long gcd(long long a,long long b)
{
	if(b==0)
		return a;
	return gcd(b,a%b);
}
void exgcd(long long a,long long b,long long &x,long long &y)
{
	if(b==0)
	{
		x=1;y=0;
		return;
	}
	long long tx,ty;
	exgcd(b,a%b,tx,ty);
	x=ty;
	y=tx-(a/b)*ty;
}
long long PowMod(long long a,long long b,long long P)
{
	long long ret=1;
	while(b)
	{
		if(b&1)
			ret=(1LL*ret*a)%P;
		b>>=1;
		a=(1LL*a*a)%P;
	}
	return ret;
}
long long inv(long long a,long long P)
{
	long long x,y;
	exgcd(a,P,x,y);
	x=(x%P+P)%P;
	return x;
}

long long exBSGS(long long A,long long B,long long P)//A^X=B (mod P)
{
	long long k=0,d,a2=1;
	B%=P;
	A%=P;
	if(B==1)
		return 0;
	if(A==B)
		return 1;
	for(;;)//前k次，找公因数d
	{
		d=gcd(A,P);
		if(B%d!=0)
			return -1;//无解
		if(d==1)
			break;//已经互质，退出
		B/=d;P/=d;
		a2=(1LL*a2*A/d)%P;//记录(a^k)/(D)
		k++;
		if(a2==B)
			return k;
	}
	Hash::Clear();
	long long m=ceil(sqrt(P));
	long long s=B%P,iA=inv(A,P),Am=PowMod(A,m,P);
	//s=(B/D)*(a^-1)
	for(long long j=0;j<=m;j++)
	{
		if(Hash::Get(s)==-1)
			Hash::Insert(s,j);
		s=(1LL*s*iA)%P;
	}
	//s=(a^k)/D*(A^im)
	s=a2;
	for(long long i=0;i<=m;i++)
	{
		long long j=Hash::Get(s);
		if(j!=-1)
			return i*m+j+k;
		s=(1LL*s*Am)%P;
	}
	return -1;
}

int main()
{
	long long A,B,P,ans;
	for(;;)
	{
		scanf("%lld%lld%lld",&A,&P,&B);
		if(A==0&&P==0&&B==0)
			break;
		ans=exBSGS(A,B,P);
		if(ans==-1)
			puts("No Solution");
		else
			printf("%lld\n",ans);
	}
	return 0;
}