浅谈BSGS

用于求解形如(a^x≡bmod p)的最小非负整数解(x).

由欧拉定理(a^{phi(p)}≡1mod p)可以知道，我们找的解如果有解则一定在(phi(p))范围内，而最大的范围就是当(p)为质数时，等于(p-1).

一种暴力方法是枚举指数验证。由于(gcd(a,p)=1).则(a)在(mod p)意义下必有逆元。所以，我们考虑分解一下质数的表示形式。

知道最大范围不超过(p)，所以我们令(x=i*m-j),将(a^{-j})移项到左边，变成(a^{i*m}≡b*a^jmod p).这里(m=lceil{sqrt{p}} ceil).

我们枚举(b*a^j)，把它们(mod p)的值插入到( ext{Hash})表里面。我们再枚举(a^{i*m}),如果这个值在( ext{Hash})表里面出现过，则跳出循环，输出答案。

特判无解，一是没有找到解，另一个是无解，即(amod p=0,b!=0),或是(a=0)的时候。（就那几个特殊情况）

这就是(BSGS)的基本思路。下面是( ext{fast_power,Exgcd,BSGS})的三合一模板。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define int long long
int mul(int a,int b,int p){return (1ll*a*b)%p;}
int add(int a,int b,int p){return (1ll*a+b)%p;}
int qpow(int a,int b,int p){
	int res=1;
	while(b){
		if(b&1)res=mul(res,a,p);
		a=mul(a,a,p);b>>=1;
	}
	return res;
}
int Exgcd(int a,int b,int &x,int &y){
	if(!b){x=1,y=0;return a;}
	int res=Exgcd(b,a%b,x,y);
	int tmp=x;x=y;y=tmp-a/b*y;
	return res;
}
map<ll,ll>mp;
int T,opt;
signed main(){
	scanf("%lld%lld",&T,&opt);
	for(;T;--T){
		int y,z,p;
		scanf("%lld%lld%lld",&y,&z,&p);
		if(opt==1)printf("%lld
",qpow(y,z,p));
		else if(opt==2){
			int x,yy;
			int G=Exgcd(y,p,x,yy);
			if(z%G){puts("Orz, I cannot find x!");continue;}
			int tmp=p/G;while(x<0)x+=tmp;
			printf("%lld
",((x*z/G)%tmp+tmp)%tmp);
		}
		else{
			mp.clear();
			if(y%p==0&&z){puts("Orz, I cannot find x!");continue;}
			int m=ceil(sqrt(p)),f=qpow(y,m,p),now=z%p;
			for(int i=1;i<=m;++i)now=mul(now,y,p),mp[now]=i;
			now=1;int flag=1;
			for(int i=1;i<=m;++i){
				now=mul(now,f,p);
				if(mp[now]){
					int ans=add(mul(i,m,p),-mp[now],p);
					printf("%lld
",add(ans,p,p));
					flag=0;break;
				}
			}
			if(flag)puts("Orz, I cannot find x!");
		}
	}
	return 0;
} 

( ext{Hash})使用了( ext{map}).

时间复杂度(O(sqrt{phi(p)})),非扩展(BSGS)适用于(p)为质数情况。(gcd(p,a)=1).