exBSGS学习笔记

Tags：数学

---

题目的话就做下洛谷的模板好了

// luogu-judger-enable-o2

#include<algorithm>

#include<iostream>

#include<cmath>

#include<map>

using namespace std;

int A,B,P;

map<int,int> H;

int ksm(int x,int k)

{

    int s=1;for(;k;k>>=1,x=1ll*x*x%P)

                if(k&1) s=1ll*s*x%P;return s;

}

int exBSGS(int A,int B,int P)//求A^t=B(mod P)的最小t

{

    int s=1,c=0,d,M=sqrt(P)+1;H.clear();

    for(d=__gcd(A,P);d!=1;d=__gcd(A,P))

    {

        if(B%d) return -1;

        c++;P/=d;B/=d;s=1ll*A/d*s%P;

        if(s==B) return c;

    }

    for(int i=0,t=B;i<M;i++,t=1ll*t*A%P) H[t]=i;

    for(int i=1,bs=ksm(A,M),t=1ll*s*bs%P;i<=M;i++,t=1ll*t*bs%P)

        if(H.count(t)) return i*M-H[t]+c;return -1;

}

int main()

{

    while(1)

    {

        cin>>A>>P>>B;

        if(!(A+B+P)) return 0;

        int ans=exBSGS(A,B,P);

        if(ans==-1) puts("No Solution");

        else cout<<ans<<endl;

    }

}