洛谷-4195 【模板】exBSGS/Spoj3105 Mod

题目描述
已知数a,p,b，求满足a^x≡b(mod p)的最小自然数x。
输入格式
每个测试文件中最多包含100组测试数据。
每组数据中，每行包含3个正整数a,p,b。
当a=p=b=0时，表示测试数据读入完全。
输出格式
对于每组数据，输出一行。
如果无解，输出“No Solution”（不含引号），否则输出最小自然数解。

输入输出样例
输入 #1
5 58 33
2 4 3
0 0 0

输出 #1 复制
9
No Solution

说明/提示
100%的数据，a,p,b≤1e9。

#include <cstdio>
#include <cmath>
#include <unordered_map>
using namespace std;
#define ll long long
unordered_map <ll,int> Hash;
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
ll exbsgs(ll A,ll B,ll p){
    if(B==1) return 0;
    ll ct=0,d,k=1;
    while((d=gcd(A,p))^1){
        if(B%d) return -1;
        B/=d,p/=d,++ct;
        k=k*(A/d)%p;
        if(k==B) return ct;
    }
    ll t=sqrt(p)+1,kt=1;
    Hash.clear();
    for(int i=0;i<t;i++){
        Hash[kt*B%p]=i;
        kt=kt*A%p;
    }
    k=kt*k%p;
    for(int i=1;i<=t;i++){
        if(Hash.find(k)!=Hash.end()) return i*t-Hash[k]+ct;
        k=kt*k%p;
    }
    return -1;
}
int main(){
    ll a,p,b;
    while(scanf("%lld%lld%lld",&a,&p,&b),a&&p&&b){
        ll ans=exbsgs(a,b,p);
        if(~ans) printf("%lld\n",ans);
        else puts("No Solution");
    }
    return 0;
}