EXBSGS模板

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define mod 600000
using namespace std;
int cnt,tot;
int head[600005];
long long a,b,c,x,y;
struct Edge{
    int fr;
    int val;
    int nxt;
    long long to;
}edge[600005];
void init(){
    cnt=0;
    memset(head,-1,sizeof(head));
}
void addedge(int u,long long v,int w){
    edge[cnt].fr=u;
    edge[cnt].to=v;
    edge[cnt].val=w;
    edge[cnt].nxt=head[u];
    head[u]=cnt++;
}
long long exgcd(long long i,long long j,long long &x,long long &y){
    if(!j){
        x=1,y=0;
        return i;
    }
    long long tmp=exgcd(j,i%j,y,x);
    y-=(i/j)*x;
    return tmp;
}
long long BSGS(long long A,long long B,long long C){
    init();
    int m=sqrt(C);
    long long mul1=1;
    long long mul2=1;
    for(int i=1;i<=m;i++){
        mul1=(mul1*A)%C;
        long long ad=(mul1*B)%C;
        if(!ad)return -1;
        addedge(ad%mod,ad,i);
    }
    for(int i=1;i<=m+3;i++){
        mul2=(mul2*mul1)%C;
        for(int j=head[mul2%mod];j!=-1;j=edge[j].nxt){
            long long v=edge[j].to;
            if(v==mul2)return 1ll*i*m-1ll*edge[j].val;
        }
    }
    return -1;
}
void EXBSGS(long long A,long long B,long long C){
    tot=0;
    long long d=exgcd(A,C,x,y);
    long long mul=1;
    while(d!=1){
        if(B%d){
            if(B==1){
                printf("0
");
                return;
            }else{
                printf("No Solution
");
                return;
            }
        }
        B/=d,C/=d,tot++,mul*=(A/d);
        d=exgcd(A,C,x,y);
    }
    d=exgcd(mul,C,x,y);
    x=(x%C+C)%C;
    long long ans=BSGS(A,(B*x)%C,C);
    if(ans!=-1)printf("%lld
",ans+tot);
    else printf("No Solution
");
}
int main(){
    while(1){
        scanf("%lld%lld%lld",&a,&b,&c);
        if(!a&&!b&&!c)break;
        EXBSGS(a,c,b);
    }
    return 0;
}