【学术篇】SDOI2011 计算器

好一道三合一…(然而被我做成了四合一)
其实1 2 3是独立的OvO

然后就可以逐个分析了…

1

快速幂..就不说了..(我省选的时候有这么水的20pts部分分么←_←

2

两种做法(写在标题里面了)..

2.1(扩展欧几里得)

$xy\equiv z(mod\ p)$
很显然可以写成$xy=np+z$, 移项得$xy-np=z$
为保证$y,p$互质, 我们让两边同除$gcd(y,p)$, 也就是让$y,p,z$分别除$gcd(y,p)$
由扩展欧几里得, 如果$gcd(y,p)$不能整除$z$,那么方程无解…
处理完直接exgcd就行了.. 这样我们就能得到$xy-np=1$的解了.. 然后再乘个$z$后对p取模即可..
如果结果小于0呢? 我们可以$+z$,一直加到不再小于0
或者是可以用$p-(-p)\ mod\ b$这样的方式来处理…

2.2(费马小定理)

嘛,其实有一个更简单的角度OvO
$xy\equiv z(mod$ $p)$ ∴$x=y^{-1}*z$
所以我们求个逆元就完了OvO 由于$p$是质数, 而$y$一定不是$p$的倍数所以一定互质…
直接$x=z*y^{p-2}$就做完了…

3

BSGS板子…
BSGS算法 原名Baby steps Giant steps算法 又叫大步小步、北上广深、拔山盖世算法…
就是专门用来求对于给定$y,z$,满足$y^x\equiv z(mod\ p)$中的最小的x的方法…
其思路如下:
- 令$x=i*m-j$(这里写成减号方便移项), 则$y^{im-j}\equiv z(mod\ p)$
- 移项, $y^{im}\equiv z*y^j(mod\ p)$
然后我们枚举$i,j$即可…
这里的$m$我们取$\left \lceil \sqrt p \right \rceil$(可能是跟均值不等式有什么不可告人的关系吧??我也不知道)
之后先从$0\sim m$枚举$j$ ,将$z*y^j$存到表里…(开map大概就可以了OvO, 记得值相同就覆盖掉, 因为要求$im-j$的最小值, 所以要让$j$尽量大…)
然后从$1\sim m$枚举$i$ ,如果$y^{im}$在表里找到了相等的值, 那就找到了$x$, 而且$x$肯定是最小的..

这样就行了..具体实现可以看代码

代码

今天心情特别好, 所以写代码要压行… 然后这个三合一写了15行OvO
所以就出现了很难看的代码OvO…建议学习代码实现的还是去百度另寻高明吧OvO

#include <map>
#include <cmath>
#include <cstdio>
typedef long long LL;
std::map<LL,int> v;LL T,L,y,z,c;
#define NOANS   {puts("Orz, I cannot find x!");return;}
inline LL gn(LL a=0,char c=0){for(;c<48||c>57;c=getchar());for(;c>47&&c<58;c=getchar())a=a*10+c-48;return a;}
LL qpow(LL a,LL b,LL p,LL s=1){for(;b;b>>=1,a=a*a%p)if(b&1) s=s*a%p;return s;}
void B(){y%=c;z%=c;if(z&&!y)NOANS printf("%d\n",z*qpow(y,c-2,c)%c);} //这里用的费马(因为短)
void C(LL m=1,LL x=1,LL p=1){if(y%c==0)NOANS m=sqrt(c)+0.5;x=z;p=qpow(y,m,c);v.clear();
    for(LL j=0;j<=m;++j)v[x]=j,x=x*y%c;x=1; 
    for(LL i=1;i<=m;++i){x=x*p%c;if(v[x]){printf("%lld\n",i*m-v[x]);return;}}
NOANS} //BSGS
void WORK(){T=gn(),L=gn();while(T--){y=gn(),z=gn(),c=gn();if(L==1)printf("%d\n",qpow(y,z,c));if(L==2)B();if(L==3)C();}}
int main(){WORK();}
//LL gcd(LL a,LL b){return !b?a:gcd(b,a%b);}
//void exgcd(LL a,LL b,LL& x,LL& y){!b?x=1,y=0:(exgcd(b,a%b,y,x),y-=(a/b)*x);}
//void B(LL a=0,LL b=0,LL d=0,LL x=0){a=y;b=-c;d=gcd(a,b);if(z%d)NOANS
//  a/=d;b/=d;z/=d;exgcd(a,b,x,y);x=x*z%b;if(x<0)x=b-(-x)%b;printf("%d\n",x);
//} //被注释的几行是2.1(扩欧的实现)