BSGS和EXBSGS

BSGS

全称

Baby Step Gaint Step（大步小步法）

内容

求解$A^x\equiv B(mod\ C)\quad(0\le x<C)$一类的方程（$C$为素数）。

$0\le x<C$的证明：

当$C$为素数时，根据欧拉定理得$A^{\varphi(C)}\equiv1(mod\ C)$
又$A^0\equiv1(mod\ C)$
$\therefore 0$到$\varphi(C)$是一个循环节。
$\therefore$如果$x$在$[0,C)$上无解，那么此方程也无解。

过程

以下默认$gcd(A,C)=1$。因为如果$gcd(A,C)\not=1$，由于$C$是素数，$A$就是$C$的倍数了，这种情况直接特判就好了。

令$m=\lceil\sqrt C\rceil,x=i\cdot m+j$，那么原式转化为$A^{i\cdot m+j}=B(mod\ C)$。这时我们只需要枚举$i$即可。时间复杂度$O(\sqrt C)$。
对于枚举出来的$i$，设$D=A^{i\cdot m},$原式转化为$D\cdot A^j\equiv B(mod\ C)$，用EXGCD就可以求出$D$。
求完$D$，接下来就是求$j$。我们可以事先$O(\sqrt C)$把$A^j$存到哈希表里，然后$O(1)$查询。

EXBSGS

内容

求解$A^x\equiv B(mod\ C)\quad(0\le x<C)$一类的方程（$C$不一定为素数）。

过程

令$d=gcd(A,C),A=a*d,B=b*d,C=c*d$

则$(a*d)^x\equiv b*d(mod\ c*d)$
$a*(a*d)^{x-1}\equiv b(mod\ c)$

然后我们把多出来的那个$a$乘到$D$里面（就是上面的那个$D$），就又循环了一遍。
最后原式就变成$D*A^{x-cnt}\equiv B(mod\ C)$，令$x=i*m+j+cnt$，然后和BSGS一样做就好了。
但是因为$i\ge0,j\ge0$，所以还要枚举$0$到$cnt-1$这一段。