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=1;
    for(;b;base=base*base%md,b>>=1)   if(b&1) ans=ans*base%md;
    return ans;
}
LL logm(LL a,LL b,LL md)
{
    if(md==0)   return -1;
    if(a%md==0) return b?-1:(a==0);
    if(b==1)    return a?0:-1;
    //printf("a%lld %lld %lld\n",a,b,md);
    //以上为特判，记一下
    LL g,num=0,d=1;
    while((g=__gcd(a,md))!=1)
    {
        if(b%g) return -1;
        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-1)/sz+1,i,now,t=poww(a,sz,md);//sz块大小，sz1块数
    unordered_map<LL,LL> ma;
    for(i=0,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=1,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 -1;
}
LL a,md,b,ans;
int main()
{
    while(1){
    scanf("%lld%lld%lld",&a,&md,&b);
    if(a==0&&md==0&&b==0)   break;
    ans=logm(a,b,md);
    ans==-1?puts("No Solution"):printf("%lld\n",ans);
    }
    return 0;
}

压一下：

#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=1;
    for(;b;base=base*base%md,b>>=1)   if(b&1) ans=ans*base%md;
    return ans;
}
LL logm(LL a,LL b,LL md)
{
    if(md==0)   return -1;
    if(a%md==0) return b?-1:(a==0);
    if(b==1)    return a?0:-1;
    LL g,num=0,d=1;
    while((g=__gcd(a,md))!=1)
    {
        if(b%g) return -1;
        num++;b/=g;md/=g;d=d*(a/g)%md;
        if(b==d)    return num;
    }
    LL sz=sqrt(md+0.5),sz1=(md-1)/sz+1,i,now,t=poww(a,sz,md);
    unordered_map<LL,LL> ma;
    for(i=0,now=b%md;i<sz;i++,now=now*a%md)    ma[now]=i;
    for(i=1,now=d;i<=sz1;i++)
    {
        now=now*t%md;
        if(ma.count(now))   return i*sz-ma[now]+num;
    }
    return -1;
}
LL a,md,b,ans;
int main()
{
    while(1){
    scanf("%lld%lld%lld",&a,&md,&b);
    if(a==0&&md==0&&b==0)   break;
    ans=logm(a,b,md);
    ans==-1?puts("No Solution"):printf("%lld\n",ans);
    }
    return 0;
}