【文文殿下】ExBSGS

无需逆元版本：

#include<cstdio>
#include<cassert>
#include<cmath>
#include<map>
typedef long long ll;
ll gcd(ll a,ll b) {
	return b?gcd(b,a%b):a;
}
ll qpow(ll a,ll b,ll p) {
	ll ret = 1;
	while(b) {
		if(b&1) {
			ret=ret*a%p;
		}
		a=a*a%p;
		b>>=1;
	}
	return ret;
}
ll exgcd(ll a,ll b,ll &x,ll &y) {
	if(b==0) {
		x=1;y=0;
		return a;
	}
	ll d = exgcd(b,a%b,x,y);
	ll tmp = x;
	x = y;
	y=tmp-a/b*y;
	return d;
}
inline ll bsgs(ll a,ll b,ll p) {
	std::map<ll,ll> M;
	if(p==1) return 0;
	a%=p;b%=p;
	ll d,cnt=0,q=1;
	while((d=gcd(a,p))!=1) {
		if(b==1) return cnt;
		if(b%d) return -1;
		b/=d;
		p/=d;
		++cnt;
		q=a/d*q%p;
	}
	ll lmt = ceil(sqrt(p));
	ll tmp = b%p;
	for(int i = 0;i<lmt;++i,tmp=tmp*a%p) {
		M[tmp]=i;
	}
	tmp = qpow(a,lmt,p);
	for(int i = 1;i<=lmt+1;++i) {
		q=q*tmp%p;
		if(M.count(q)) {
			return i*lmt-M[q]+cnt;
		}
	}
	return -1;
}
int main() {
	ll a,b,p;
	scanf("%lld%lld%lld",&a,&p,&b);
	while(a||b||p) {
		ll ans = bsgs(a,b,p);
		if(ans==-1) puts("No Solution");
		else printf("%lld
",ans);
		scanf("%lld%lld%lld",&a,&p,&b);
	}
	return 0;
}