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