spoj3105 MOD - Power Modulo Inverted（exbsgs）

裸的ex_bsgs，不懂ex_bsgs的可以看看这篇博客---->超级跳转

#include<bits/stdc++.h>
#define ll long long
using namespace std;
map<ll,ll> mp; 
ll gcd(ll x,ll y)
{
	return y==0?x:gcd(y,x%y);
}
ll powmod(ll x,ll y,ll mod)
{
	ll sum=1;
	while(y)
	{
		if(y&1)
		{
			sum=sum*x%mod;
		}
		x=x*x%mod;
		y>>=1;
	}
	return sum;
}
void ex_bsgs(ll a,ll b,ll mod)
{
     ll i,j;
     if(b==1)
     {
     	printf("0\n");
     	return;
     }
     ll k=0,x=1;
     while(true)
     {
     	ll d=gcd(a,mod);
     	if(d==1) break;
     	if(b%d)
     	{
     		printf("No Solution\n");
     		return;
     	}
     	mod/=d,b/=d,++k,x=x*a/d%mod;
     	if(x==b)
     	{
     		printf("%lld\n",k);
     		return;
     	}
     }
     mp.clear();
     ll m=sqrt(mod)+1,t,tt;
     for(i=0,t=b;i<m;i++,t=t*a%mod) mp[t]=i;
     for(i=1,tt=powmod(a,m,mod),t=x*tt%mod;i<=m;i++,t=t*tt%mod)
     {
     	j=mp.find(t)==mp.end()?-1:mp[t];
		 if(j>=0&&i*m-j+k>=0)
		 {
		 	printf("%lld\n",i*m-j+k);
		 	return;
		 } 
     }
} 
int main()
{
	ll i,j,x,z,k;
	while(true)
	{
		scanf("%lld%lld%lld",&x,&z,&k);
		if(k==0&&z==0&&x==0) break;
		ex_bsgs(x,k,z);
	}
	return 0;
}