数论：拓展欧几里得求解二元一次不定方程

拓展欧几里得算法求解二元一次不定方程

本节内容，主要围绕 为什么、怎么做、怎么写 这三个关系来讲解。

为什么？

为什么拓展欧几里得算法可以求二元一次不定方程的解？

首先明确：对于任意的$ax+by=c$都可以转化为$ax+by=gcd(a,b)$的求解。

所以只需要关注：$ax+by=gcd(a,b)$的求解过程即可。

而对于$ax+by=gcd(a,b)$的求解，也只是在$a^{\prime }x+b^{\prime }y=1$的基础上进行的拓展。

所以本质就是用$exgcd$求解$a^{\prime }x+b^{\prime }y=1$

由裴蜀定理可知，若$gcd(a,b)=d$，则一定存在一组整数解$(x,y)$满足$ax+by=d$​。

由欧几里得算法：$gcd(a,b)=gcd(b,amodb)$

所以，原问题求$(ax+by=d)$的通解，等价于求：$bx+(amodb)y=d$的通解。

而求$bx+(amodb)y=d$的通解，等价于求$(amodb)x+bmod(amodb)y=d$的通解…。

很眼熟对吧？这类似于欧几里得算法的求解过程，这个步骤是有限的，我们只要不断的递归下去，一定有尽头。

尽头就是：$b=0$​

当我们求到尽头的时候，原方程转化为：$a^{\prime }x=d$ 有整数解，解为：$x=\frac{d}{a^{\prime }}$​

一个约定俗成的事是可以直接令$x=1,y=0$​

我们一般关注的不是特解，而是通解/最小正整数解，所以特解是什么都无所谓。

怎么做？

我们发现，即使我们不知道$ax+by=c$的解，我们仍然可以通过上述方法不断迭代，直到迭代到尽头。此时，我们就可以求出

考虑任意两层$(x,y)$之间的转化：

设当前层为：$(x_1,y_1)$，下一层为：$(x_2,y_2)$

联立：
$$\{\begin{array}{l}a{x}_1+b{y}_1=d\\ b{x}_2+[amodb]y_2=d\end{array}$$
中间状态：
$$\{\begin{array}{l}a{x}_1+b{y}_1=d1\\ b{x}_2+(a-[\frac{a}{b}]\ast b)y_2=d\to a{y}_2+b(x_2-[\frac{a}{b}]y_2)=d\end{array}$$
系数对齐：
$$\{\begin{array}{l}a{x}_1+b{y}_1=d\\ a{y}_2+b(x_2-[\frac{a}{b}]y_2)=d\end{array}$$

得到了相邻两层之间的整数解$(x,y)$的转移方程：
$$\{\begin{array}{l}{x}_1=y_2\\ y_1=x_2-[\frac{a}{b}]y_2\end{array}$$
先直接地给出代码：

int exgcd (int a, int b, int &x, int &y) {
        if (!b) {x = 1; y = 0; return a;}
        int d = exgcd (b, a % b, x, y);
        int x2 = x, y2 = y;
        x = y2, y = x2 - a / b * y2; 
        return d;
}

怎么写？

对于第i层，$x_i=y_{i+1}$，因为算法的特性，我们一开始并不知道$x,y$的值，所以我们不能对其进行任

何运算，而是需要递归到底部对$x,y$进行初始化，并在回溯过程中才能更新每一层的$x,y$。（要完成这样的操作，我们需要传入$x,y$​的引用。）

而回溯的时候，$x,y$ 的状态都是上一层的状态：$(x_2,y_2)$表示上一层的$(x,y)$

int x2 = x, y2 = y;
x = y2, y = x2 - a / b * y2; 

这并不难记，一旦忘记了相邻两层之间的$x,y$的传递关系，利用欧几里得算法推导即可。
$$\{\begin{array}{l}a{x}_1+b{y}_1=d\\ b{x}_2+[amodb]y_2=d\end{array}$$
正常我们直接抄板子就行了：

int exgcd (int a, int b, int &x, int &y) {
        if (!b) {x = 1; y = 0; return a;}
        int d = exgcd (b, a % b, y, x);
        y -= a / b * x;
        return d;
}

值得注意的是，我们求出来的$(x_1,y_1)$是关于方程：$ax+by=1$ 的一个解。

如何去求通解？

其实通过数学直觉，我们不难想出，如果$x$增加$\Delta x$ , 相当于方程左边增加了$a\Delta x$ ,我们要想方程右边不变的话，设让$y$减去$\Delta y$

$a(x+\Delta x)+b(y-\Delta y)=1$

可以得到等式：$a\Delta x=b\Delta y$​

$\therefore \frac{\Delta x}{\Delta y}=\frac{b}{a}$

接下来我们将右边的分式约分为最简式子即可：$\frac{b}{a}=\frac{\frac{b}{gcd(a,b)}}{\frac{a}{gcd(a,b)}}$

即除以它们的最大公约数。

通解：
$$\{\begin{array}{l}x=x_1+\frac{b}{gcd(a,b)}\times k\\ y=y_1-\frac{a}{gcd(a,b)}\times k\end{array}$$