拓展BSGS算法

扩展BSGS算法（ExBSGS）

ExBSGS算法用于解决离散对数问题，以解出关于x的同余方程的最小非负整数解：
$a^x\equiv b(\mathrm{mod}m)$

不保证$m\perp a$。

ExBSGS算法的做法是：

若$a/\perp m$，则同余方程等价于不定方程$a\cdot a^{x-1}+m\cdot y=b$，根据裴蜀定理，若$\gcd (a,m)|b$，有解，此时根据技巧，有$\frac{a}{g}\cdot a^{x-1}\equiv \frac{b}{g}\left(\mathrm{mod}\frac{m}{g}\right)$，否则无解。

然后此时$\frac{a}{g}\perp \frac{m}{g}$，有逆元，除过去：
$a^{x-1}\equiv \frac{b}{g}{\left(\frac{a}{g}\right)}^{-1}\left(\mathrm{mod}\frac{m}{g}\right)$

若$a$与$\frac{m}{g}$仍不互质，就再除。直到互质，就变成了BSGS问题。很轻松求解出$x$，最终答案就是$x+除了多少次$。

因此有递归函数：
$ExBSGS(a,b,m)=\{\begin{array}{c}a\perp m:BSGS(a,b,m)\\ \{\begin{array}{c}gcd(a,m)|b:ExBSGS\left(a,\frac{b}{g}\cdot {\left(\frac{a}{g}\right)}^{-1}\left(mod\frac{m}{g}\right),\frac{m}{g}\right)+1\\ else:无解\end{array}\end{array}$

还有一种做法是，直接把$\frac{a}{g}$留在左边，作为$BSGS$时$a^x$的系数，也是可以的。

时间复杂度$O(lo{g}^{2}n+\sqrt{n})=O(\sqrt{n})$

后记

于是皆大欢喜。