# BZOJ 2242 BSGS

## Description

1、给定y,z,p,计算Y^Z Mod P 的值；
2、给定y,z,p，计算满足xy≡ Z ( mod P )的最小非负整数；
3、给定y,z,p，计算满足Y^x ≡ Z ( mod P)的最小非负整数。

输入包含多组数据。

【样例输入1】
3 1
2 1 3
2 2 3
2 3 3
【样例输入2】
3 2
2 1 3
2 2 3
2 3 3
【数据规模和约定】

## Sample Output

【样例输出1】
2
1
2
【样例输出2】
2
1
0

#include <iostream>
#include <cstring>
#include <algorithm>
#include <stdio.h>
#include <math.h>
#include <map>
#define ll long long

template<typename _t>
_t x=0,f=1;
char ch=getchar();
for(;ch>'9'||ch<'0';ch=getchar())if(ch=='-')f=-f;
for(;ch>='0'&&ch<='9';ch=getchar())x=x*10+(ch^48);
return x*f;
}

namespace work1{
inline void Main(){
for(;k;k>>=1,x=x*x%P)if(k&1)Ans = Ans * x % P;
printf("%lld\n",Ans);
}
}

namespace work2{
inline ll gcd(ll a,ll b){return b==0?a:gcd(b,a%b);}
inline ll exgcd(ll a,ll b,ll &x,ll &y){
if(!b){x=1;y=0;return a;}
ll ans = exgcd(b,a%b,x,y);
ll t = x;x=y;y=t-a/b*y;return ans;
}

inline void Main(){
ll Ans = gcd(a,P);
if(b%Ans!=0){printf("Orz, I cannot find x!\n");return;}
a/=Ans;b/=Ans;P/=Ans;
exgcd(a,P,x,y);
x = ((x*b)%P+P)%P;
printf("%lld\n",x);
}
}

namespace work3{
std::map<ll,ll>ma;
inline ll qpow(ll x,ll k,ll P){
ll Ans = 1;
for(;k;k>>=1,x=x*x%P)
if(k&1)Ans = Ans * x % P;
return Ans;
}

inline void Main(){
ma.clear();
if(y%P==0){
printf("Orz, I cannot find x!\n");
return;
}
int len = ceil(sqrt((double)P));
ll M_pow = qpow(y,len,P);
ll now = z%P;ma[z%P] = 0;
for(int i=1;i<=len;i++){
now = now * y % P;
ma[now]=i;
}
now = 1;
for(int i=1;i<=len;i++){
now = now * M_pow % P;
if(ma.count(now)){
printf("%d\n",i*len-ma[now]);return;
}
}
printf("Orz, I cannot find x!\n");
}
}

int main(){
register int T,type;
}