exBSGS学习笔记

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;
    }
}