BSGS算法

问题：给定三个数a,b,p,保证a,p互质，求解$a^x\equiv b(modp)$
最小正整数$x$ $1<=a,b,p<=10^9$

BSGS（拔山盖世 大步小步算法）
根据欧拉定理：$a^{\phi (n)}\equiv 1$ $(mod$ $n)$
可以看出指数的同余的循环节
直接暴力时间复杂度$O(\phi (n))$
考虑优化
令$x=p^{\prime }q-r$ （同余常见搞法，设$a=pb+r$ 加减一样）
则$A^{p^{\prime }q}\equiv b\ast A^r$$(mod$ $n)$
首先在$O(p^{\prime })$内求出$a$至$a^{p^{\prime }-1}$每一个的值（$a,a^2,a^3...$）
并把这些值存入哈希表
然后枚举$q$，每次从哈希表中查询是否有满足原式的值
时间复杂度$O(max(p^{\prime },\phi (p)/p^{\prime }))$
显然最优是令$p^{\prime }=\sqrt{\phi (p)}$
总的时间复杂度$O(\sqrt{\phi (p)})$
That’s all.

$Code$

($p$为质数，$a,p$互质还可以求$\phi (p)$)
模板题：SDOI2011 计算器
TJOI2007 可爱的质数

#include <bits/stdc++.h>
using namespace std;
#define int long long
const int Mod=1e5+7;
const int MAXN=1<<20;
struct w{
	int to,nx,s;
	w(){s=-1;}
}Hash[MAXN+10];
int a[Mod+10],tot=0;

inline void add(int,int,int,int);
void Insert(int);
int Query(int);
void BSGS(int,int,int);

signed main(){
	//freopen ("std.in","r",stdin);
	//freopen ("std.out","w",stdout);
	int p,b,n;
	scanf("%lld%lld%lld",&p,&b,&n);
	BSGS(p,n,b);
	return 0;
}

inline void add(int x,int i,int val,int s){Hash[i].to=val,Hash[i].s=s,Hash[i].nx=a[x],a[x]=i;}
void Insert(int val,int s){
	int x=val%Mod;
	add(x,++tot,val,s);
}

int Query(int val){
	int x=a[val%Mod];
	while (Hash[x].to!=val && Hash[x].to)	x=Hash[x].nx;
	return Hash[x].s;
}

void BSGS(int p,int n,int a){
	int m=sqrt(p)+1,t=1,tt;
	for (register int i=0;i<m;++i){
	 	Insert((t*n)%p,i);
	 	t=(t*a)%p;
	}
	tt=t;
	for (register int i=1;i<=m;++i){
		int k=Query(t);
		t=(t*tt)%p;
		if (k==-1)	continue;
		printf("%lld\n",i*m-k);
		return;
	}
	printf("no solution\n");
}