poj3243 Clever Y（EXBSGS+Hash表+欧拉函数）

扩展BSGS，经典的BSGS只能求解p为质数时的离散对数。

那么如果p不是质数呢？其实我们发现p是不是质数不重要，重要的是(a,p)==1.

$a^t\equiv b\mod p$
如果(a,p)!=1,我们设为k,a=Ak,p=Pk,b=Bk(显然当b不是k的倍数时无解)
那么我们有${Ak}^t\equiv {Bk}\mod {Pk}$
可以把k除掉，得到$Aa^{t-1}\equiv {B}\mod {P}$
一直下去，直到(a,p)==1，每次b/=k,p/=k,++cnt，d*=(a/k)

那么最后就是$da^{t-cnt}\equiv b\mod p$的形式，(a,p)==1.
我们令$t-cnt=im+j$，类似BSGS同样地去做即可。最后答案就是im+j+cnt
注意到这样我们就枚举不到0~cnt-1的答案了，我们可以提前枚举了。
(cnt是log mod级别的。)
求逆元可以欧拉定理也可以扩展欧几里得。

然后这题比较毒，有可能b>=mod，此时不一定是no solution 而是要先b%=mod

改了半天，原来是输出的s是大写的S。。。气到没脾气【再见】

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100100
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(S==T) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int mod,m,h[N],num=0;
struct Hash_table{
    int next,key,val;
}data[N];
inline void ins(int key,int val){
    int x=key%N;
    for(int i=h[x];i;i=data[i].next){
        if(data[i].key==key) return;
    }data[++num].next=h[x];h[x]=num;data[num].key=key;data[num].val=val;
}
inline int hs(int key){
    int x=key%N;
    for(int i=h[x];i;i=data[i].next){
        if(data[i].key==key) return data[i].val;
    }return -1;
}
inline int gcd(int x,int y){return y?gcd(y,x%y):x;}
inline int ksm(ll x,int k){
    ll res=1;for(;k;k>>=1,x=x*x%mod) if(k&1) res=res*x%mod;return res;
}
inline int phi(int x){
    int res=x,xx=x;
    for(int i=2;i*i<=x&&xx!=1;++i){
        if(xx%i) continue;
        res-=res/i;while(xx%i==0) xx/=i;
    }if(xx!=1) res-=res/xx;return res;
}
inline int EXBSGS(int a,int b){
    a%=mod;b%=mod;if(a==0) return b==0?0:-1;
    ll x=1,d=1;int cnt=0;
    for(int i=0;i<=30;++i){
        if(x==b) return i;x=x*a%mod;
    }while(1){
        int g=gcd(a,mod);if(g==1) break;
        if(b%g) return -1;b/=g;mod/=g;++cnt;d=d*(a/g)%mod;
    }m=ceil(sqrt(mod));memset(h,0,sizeof(h));num=0;
    int e=ksm(a,phi(mod)-1);x=b;
    for(int i=0;i<m;++i){
        ins(x,i);x=x*e%mod;
    }e=ksm(a,m);x=d;
    for(int i=0;i<=m;++i){
        int j=hs(x);if(j!=-1) return i*m+j+cnt;x=x*e%mod;
    }return -1;
}
int main(){
//  freopen("a.in","r",stdin);
    while(1){
        int a=read();mod=read();int b=read();
        if(!a&&!mod&&!b) break;
        int res=EXBSGS(a,b);
        if(res==-1) puts("No Solution");
        else printf("%d\n",res);
    }return 0;
}