【bzoj2242】计算器

 #include<bits/stdc++.h>
 #define inf 1000000000
 using namespace std;
 typedef long long ll;
 int gcd(int a,int b){return b==?a:gcd(b,a%b);}
 void exgcd(int a,int b,int &x,int &y){
     if(b==){x=;y=;return;}
     exgcd(b,a%b,x,y);
     int t=x;x=y;y=t-a/b*y;
 }
 int solve1(ll y,int z,int p){
     y%=p;ll ans=;
     for(int i=z;i;i>>=,y=y*y%p)if(i&)ans=ans*y%p;
     return ans;
 }
 void solve2(ll y,int z,int p){
     p=-p;int t=gcd(y,p);
     if(z%t){puts("Orz, I cannot find x!");return;}
     y/=t;z/=t;p/=t;int a,b;exgcd(y,p,a,b);
     a=1LL*a*z%p;while(a<)a+=p;
     printf("%d\n",a);
 }
 map<int,int>mp;
 void solve3(int y,int z,int p){
     y%=p;if(!y&&!z){puts("-1");return;}
     if(!y){puts("Orz, I cannot find x!");return;}
     mp.clear();
     ll m=ceil(sqrt(p)),t=;
     mp[]=m+;
     for(ll i=;i<m;i++){t=t*y%p;if(!mp[t])mp[t]=i;}
     ll tmp=solve1(y,p-m-,p),ine=;
     for(ll k=;k<m;k++){
         int i=mp[z*ine%p];
         if(i){
             if(i==m+)i=;
             printf("%lld\n",k*m+i);return;
         }
         ine=ine*tmp%p;
     }
     puts("Orz, I cannot find x!");
 }
 inline int read(){
     int f=,x=;char ch;
     do{ch=getchar();if(ch=='-')f=-;}while(ch<''||ch>'');
     do{x=x*+ch-'';ch=getchar();}while(ch>=''&&ch<='');
     return f*x;
 }
 int main(){
     int T=read(),k=read();
     while(T--){
         int y=read(),z=read(),p=read();
         if(k==)printf("%d\n",solve1(y,z,p));
         else if(k==)solve2(y,z,p);
         else solve3(y,z,p);
     }
     return ;
 }

奇奇怪怪的bsgs裸题。

不知道啥是bsgs的可以百度。

别名 拔山盖世算法，百事公司算法等。