扩展BSGS算法

普通 $BSGS$ 有一个缺点，就是当 $\gcd (n,p)\ne 1$ 时，他不能用。

于是就出现了扩展 $BSGS$。

前置知识

普通 $BSGS$。

扩展欧几里得算法。

算法用途

用于求解 $\gcd (n,p)\ne 1$ 时，$n^x\equiv m(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 的方程

（除了 $x$ 都是常数）

算法复杂度

时间

？

空间

？

算法实现

我们先把 $n^x\equiv m(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 转化成

$$n^{x-1}\times n+kq=m$$

可以看成 $an+bq=m$，就是二元一次不定方程的，设 $d=\gcd (n,q)$

根据扩欧，$d|m$，否则无解。

所以，可以把原式变成 $$n^{x-1}\times \frac{n}{d}+k\frac{q}{d}=\frac{m}{d}$$
如果 $\gcd (n,\frac{q}{d})\ne 1$，就继续进行，如下

$$n^{x-2}\times \frac{n}{d}\times \frac{n}{d_2}+k\frac{q}{d\times d_2}=\frac{m}{d\times d_2}$$

一直到 $\gcd (n,q)=1$，设进行了 $a$ 次，中间所有最小公因数的乘积为 $b$。得到

$$n^{x-a}\times \frac{n^a}{b}+k\frac{q}{b}=\frac{m}{b}$$

为了方便表示，就可以转换回。

$$n^{x-a}\times \frac{n^a}{b}\equiv \frac{m}{b}(\mathrm{m}\mathrm{o}\mathrm{d}k\frac{q}{b})$$

然后用普通 $BSGS$ 搞定就行了。