exBSGS板子

/*bzoj2480*/
#include <map>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
typedef long long LL;
std::map<LL,LL> mmp;
inline LL Pow(LL x,LL y,LL P){
  LL ret=;
  while(y){
    if(y&)ret=ret*x%P;
    x=x*x%P,y>>=;
  }
  return ret;
}
inline LL gcd(LL x,LL y){
  return x==?y:gcd(y%x,x);
}
inline LL BSGS(LL A,LL B,LL C,LL D){
  int len=ceil(sqrt(C+0.5));
  LL temp=B,sum=D;
  mmp.clear();
  for(int i=;i<=len;++i)mmp[temp=temp*A%C]=i;
  temp=Pow(A,len,C);
  for(int i=;i<=len;++i){
    sum=sum*temp%C;
    if(mmp.count(sum))return i*len-mmp[sum];
  }
  return -;
}
inline LL exBSGS(LL A,LL B,LL C){
/*A>=0(这里涉及到0^0是否有意义(是否为1),我们假设他是1好了) B>=0 C>0(%0和/0是一个效果啊)*/
/*假设我们求的解是最小自然数(无解输出-1)*/
  if(C==)return ;
  A%=C,B%=C;
  if(B==)return ;
/*判掉了解为0的情况*/
  if(A==){
    if(B==)return ;
    return -;
  }
/*现在A>0 B>=0 B!=1 C>1*/
  LL D=;int cnt=;
  for(LL g=gcd(A,C);g!=;g=gcd(A,C)){
    if(B%g!=)return -;
    C/=g,B/=g,D=D*(A/g)%C,++cnt;
    if(D==B)return cnt;
  }
/*得到的C一定不为1*/
  A%=C;
/*现在我们得到了A的cnt次方与C除去gcd后的结果(分别是D和C),并且现在A与C互质了*/
  LL ret=BSGS(A,B,C,D);
  return ret==-?-:ret+cnt;
}
int main(){
  LL a,p,b,ans;
  while(true){
    scanf("%lld%lld%lld",&a,&p,&b);
    if(p==)return ;
    ans=exBSGS(a,b,p);
    if(ans==-)puts("No Solution");
    else printf("%lld\n",ans);
  }
}