扩展欧几里得算法 - exgcd

学exgcd的时候没好好听课，几乎每次遇到都忘记。
于是打算写篇博客。

扩展欧几里得算法，就是欧几里得算法的扩展。

欧几里得算法，就是 gcd，共产党 ，用来求最大公约数的。

还是一样，首先搞明白他是干啥的。

我们有一个等式，$a,b$ 为常数：

$$ax+by=\gcd (a,b)$$

求这个等式的一组可行整数解。

由于我们知道，$\gcd (a,0)=a$（不要问我为什么，因为 $\gcd (a,0)$ 就等于 $a$）

所以当 $b=0$ 时，有 $ax=a$，显然 $x=1$。

由于可以随便取一组可行解，为了让之后的计算简便，我们让 $y=0$。

我们知道 $\gcd (a,b)=\gcd (b,a\mathrm{mod}b)$

所以
$$\begin{array}{rl}ax+by& =\gcd (a,b)\\ & =\gcd (b,a\mathrm{mod}b)\\ & =b{x}_0+(a\mathrm{mod}b)y_0\\ & =b{x}_0+(a-\lfloor \frac{a}{b}\rfloor \times b)y_0\\ & =b{x}_0+a{y}_0-\lfloor \frac{a}{b}\rfloor \times b{y}_0\\ & =a{y}_0+b(x_0-\lfloor \frac{a}{b}\rfloor \times y_0)\end{array}$$
这样我们可以得到：
$$x=y_0,y=x_0-\lfloor \frac{a}{b}\rfloor \times y_0$$

根据这个，我们可以从 $b=0$ 的时候的解反过来推出任意 $ax+by=\gcd (x,y)$ 的解

具体操作可以用普通 $\gcd$ 的函数，通过回溯的方法求解。