51NOD 1038：X^A Mod P——题解

http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1038

X^A mod P = B，其中P为质数。给出P和A B，求< P的所有X。

找了半天找到了一道模板题，死嗑了一晚上对着题解写完了。

道理应该都懂吧……不懂我也懒的再说一遍了，看其他人博客吧：https://blog.csdn.net/dreamzuora/article/details/52744666

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cctype>
#include<cstdio>
#include<vector>
#include<cmath>
using namespace std;
typedef long long ll;
const int N=;
const int M=;
struct HASH{
    int w,to,nxt;
}h[M];
int cnt,head[N];
inline void add(int x,int y){
    int t=x%N;
    for(int i=head[t];i;i=h[i].nxt){
        int v=h[i].to;
        if(v==x){
            h[i].w=y;//记大的
            return;
        }
    }
    h[++cnt].to=x;h[cnt].w=y;h[cnt].nxt=head[t];head[t]=cnt;
}
inline int query(int x){
    int t=x%N;
    for(int i=head[t];i;i=h[i].nxt){
        int v=h[i].to;
        if(v==x)return h[i].w;
    }
    return -;
}
int gcd(int a,int b){
    return (!b)?a:gcd(b,a%b);
}
int exgcd(int a,int b,int &x,int &y){
    if(!b){
        x=,y=;
        return a;
    }
    int ans=exgcd(b,a%b,y,x);
    y-=(ll)(a/b)*x;
    return ans;
}
int inv(int a,int c){
    int x,y;
    exgcd(a,c,x,y);
    return (x%c+c)%c;
}
//a^x=b(mod c);
int BSGS(int a,int b,int c){
    if(!a){
        if(!b)return ;
        return -;
    }
    int tot=,g,d=;
    while((g=gcd(a,c))!=){
        if(b%g)return -;
        ++tot;b/=g,c/=g;
        d=(ll)d*(a/g)%c;
    }
    b=(ll)b*inv(d,c)%c;
    cnt=;memset(head,,sizeof(head));
    int s=sqrt(c),p=;
    for(int i=;i<s;i++){
        if(p==b)return i+tot;
        add((ll)p*b%c,i);
        p=(ll)p*a%c;
    }
    int q=p;
    for(int i=s;i-s+<c;i+=s){
        int t=query(q);
        if(t!=-)return i-t+tot;
        q=(ll)q*p%c;
    }
    return -;
}
vector<int>v;
int qpow(int k,int n,int p){
    int res=;
    while(n){
        if(n&)res=(ll)res*k%p;
        k=(ll)k*k%p;
        n>>=;
    }
    return res;
}
bool pan(int g,int p){
    for(int i=;i<v.size();i++){
        if(qpow(g,(p-)/v[i],p)==)return ;
    }
    return ;
}
int primitive(int p){
    int res=p-;v.clear();
    for(int i=;i*i<=res;i++){
        if(res%i==){
            v.push_back(i);
            while(res%i==)res/=i;
        }
    }
    if(res!=)v.push_back(res);
    for(int i=;;i++){
        if(pan(i,p))return i;
    }
}
vector<int>X;
void solve(int p,int a,int b){
    X.clear();
    int g=primitive(p)%p;
    int t=BSGS(g,b,p);
    if(t==-)return;
    int A=a,B=p-,C=t,x,y;
    int d=exgcd(A,B,x,y);
    if(C%d)return;
    x=(ll)x*(C/d)%B;
    int delta=B/d;
    for(int i=;i<d;i++){
        x=((x+delta)%B+B)%B;
        X.push_back(qpow(g,x,p));
    }
    sort(X.begin(),X.end());
    X.erase(unique(X.begin(),X.end()),X.end());
}
int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        int p,a,b;
        scanf("%d%d%d",&p,&a,&b);
        solve(p,a,b);
        if(!X.size()){puts("No Solution");continue;}
        else{
            for(int i=;i<X.size();i++)printf("%d ",X[i]);
            puts("");
        }
    }
    return ;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+欢迎访问我的博客：http://www.cnblogs.com/luyouqi233/ +

+++++++++++++++++++++++++++++++++++++++++++