EXBSGS

http://210.33.19.103/problem/2183

参考：https://blog.csdn.net/frods/article/details/67639410（里面代码好像不太对）（不用求逆元的方法）

https://blog.csdn.net/zzkksunboy/article/details/73162229（要逆元的方法）

参考代码：https://blog.csdn.net/clove_unique/article/details/51227328

以下代码pojA不掉，没有unordered_map，map又T掉；要手写哈希表

https://www.luogu.org/problemnew/show/P3846

https://cn.vjudge.net/problem/POJ-3243

https://cn.vjudge.net/problem/POJ-2417

 #include<cstdio>
 #include<algorithm>
 #include<cmath>
 #include<tr1/unordered_map>
 using namespace std;
 using namespace tr1;
 typedef long long LL;
 LL poww(LL a,LL b,LL md)
 {
     LL base=a,ans=;
     for(;b;base=base*base%md,b>>=)   if(b&) ans=ans*base%md;
     return ans;
 }
 LL logm(LL a,LL b,LL md)
 {
     if(md==)   return -;
     if(a%md==) return b?-:(a==);
     if(b==)    return a?:-;
     //printf("a%lld %lld %lld\n",a,b,md);
     //以上为特判，记一下
     LL g,num=,d=;
     while((g=__gcd(a,md))!=)
     {
         if(b%g) return -;
         num++;b/=g;md/=g;d=d*(a/g)%md;
         if(b==d)    return num;//特殊技巧，记一下;num==0的情况前面已经判了
         //设gg为所有已经除过的g的积，原来b为B，此时的b实际相当于B/gg
         //d=a^num/gg,如果a^num==B,则d==b
     }
     //printf("b%lld %lld %lld\n",a,b,md);
     LL sz=sqrt(md+0.5),sz1=(md-)/sz+,i,now,t=poww(a,sz,md);//sz块大小，sz1块数
     unordered_map<LL,LL> ma;
     for(i=,now=b%md;i<sz;i++,now=now*a%md)//x=i*sz-j,因此j属于[0,sz-1]
     {
         ma[now]=i;//应该使得x最小，则j应越大越好，因此相同的保留最大的
         //printf("c%lld %lld\n",i,now);
     }
     for(i=,now=d;i<=sz1;i++)//因此i属于[1,sz1]
     {
         now=now*t%md;
         if(ma.count(now))   return i*sz-ma[now]+num;
         //printf("d%lld\n",i);
     }
     return -;
 }
 LL a,md,b,ans;
 int main()
 {
     while(){
     scanf("%lld%lld%lld",&a,&md,&b);
     if(a==&&md==&&b==)   break;
     ans=logm(a,b,md);
     ans==-?puts("No Solution"):printf("%lld\n",ans);
     }
     return ;
 }

压一下：

 #include<cstdio>
 #include<algorithm>
 #include<cmath>
 #include<tr1/unordered_map>
 using namespace std;
 using namespace tr1;
 typedef long long LL;
 LL poww(LL a,LL b,LL md)
 {
     LL base=a,ans=;
     for(;b;base=base*base%md,b>>=)   if(b&) ans=ans*base%md;
     return ans;
 }
 LL logm(LL a,LL b,LL md)
 {
     if(md==)   return -;
     if(a%md==) return b?-:(a==);
     if(b==)    return a?:-;
     LL g,num=,d=;
     while((g=__gcd(a,md))!=)
     {
         if(b%g) return -;
         num++;b/=g;md/=g;d=d*(a/g)%md;
         if(b==d)    return num;
     }
     LL sz=sqrt(md+0.5),sz1=(md-)/sz+,i,now,t=poww(a,sz,md);
     unordered_map<LL,LL> ma;
     for(i=,now=b%md;i<sz;i++,now=now*a%md)    ma[now]=i;
     for(i=,now=d;i<=sz1;i++)
     {
         now=now*t%md;
         if(ma.count(now))   return i*sz-ma[now]+num;
     }
     return -;
 }
 LL a,md,b,ans;
 int main()
 {
     while(){
     scanf("%lld%lld%lld",&a,&md,&b);
     if(a==&&md==&&b==)   break;
     ans=logm(a,b,md);
     ans==-?puts("No Solution"):printf("%lld\n",ans);
     }
     return ;
 }