线性方程与最大公因数（扩展欧几里得算法）

1 线性方程与最大公因数

  已知两个整数$a$与$b$，我们考察由$ax+by$得到的所有正整数，其中$x$与$y$可为任何整数。（$x,y$可取负值）。例如下表为$a=42,b=30$的情况下$ax+by$可取的一些值。

我们来观察这个表格，首先我们会发现这个表格中的每一个数都能被$6$整除，原因也很显然，因为$42$和$30$都能被$6$整除，所以$42x+30y$自然也能被$6$整除，推广到一般情况，我们能够自然地得出下面这个结论：

$$形如ax+by的每个数都能被gcd(a,b)整除$$

其次我们再观察数据会发现$gcd(42,30)=6$也出现在表格里面，这揭示了下述结论：

$$形如ax+by的最小正整数等于gcd(a,b)$$

下面我们来证明这个结论。

2 证明

  由于每个数$ax+by$都能整除$gcd(a,b)$，$ax+by$的最小正整数值恰好是$gcd(a,b)$。所以我们只需要找出求解方程:
$$ax+by=gcd(a,b)$$
的整数解$x,y$的方法就可以了。下面便给出如何使用欧几里得算法来求解该方程。以$22x+60y=gcd(22,60)$为例。

  使用欧几里得算法计算最大公因数：
$$\begin{gathered} 60=2\ast 22+16 \\ 22=1\ast 16+6 \\ 16=2\ast 6+4 \\ 6=1\ast 4+2 \\ 4=2\ast 2+0 \end{gathered}$$

我们得出$gcd(22,60)=2$，接下来我们使用算法的中间商和余数来求解$22x+60y=2$。

  首先， 将第一个式子改写成：
$$16=a-2b$$
然后用这个等式替换第二个式子中的$16$得：

$$6=b-1\ast (a-2b)=-a+3b$$

用这两个等式替换第三个式子中的$6$和$16$得：
$$4=(a-2b)-2\ast (-a+3b)=3a-8b$$
最后将第四个式子的$4$和$6$做替换得：

$$-a+3b=3a-8b+2$$
移项得：
$$-4a+11b=2$$

因此我们便得到了$22x+60y=2$的解：
$$\begin{gathered} x=-4 \\ y=11 \end{gathered}$$

  这个方法的可行性在于，欧几里得算法的每一个步骤产生的新的余数都可以表示为$a$的倍数加$b$的倍数，即：

$$最新余数=a的倍数+b的倍数$$
我们最终得到的非零余数就是$gcd(a,b)$，因此方程的解也就自然而然地得出来了。

3 解的讨论

  现在我们知道方程$ax+by=gcd(a,b)$总有整数解$x,y$，并且可以通过欧几里得算法求得。下面来讨论一下方程有多少个解以及怎样表述所有解。

  首先考虑$a,b$互素的情况，即$gcd(a,b)=1$，假设$(x_1,y_1)$是方程$ax+by=1$的解，那么通过$x_1$加上$b$的倍数，$y_1$减去$a$的倍数我可以得到其他的解，即：
$$对任何整数k，我们得到新解(x_1+kb,y_1-ka)$$
计算验证也很简单：
$$a(x_1+kb)+b(y_1-ka)=a{x}_1+b{y}_1=1$$
另外，为了证明这种方法可以给出所有解，我们假设$(x_1,y_1),(x_2,y_2)$是方程$ax+by=1$的两个解，即：
$$\begin{gathered} a{x}_1+b{y}_1=1 \\ a{x}_2+b{y}_2=1 \end{gathered}$$
我们用$y_2$乘以第一个方程，用$y_1$乘以第二个方程，再相减得到：
$$a{x}_1{y}_2-a{x}_2{y}_1=y_2-y_1$$
类似地，我们用$x_2$乘第一个方程，用$x_1$乘第二个方程，在相减得到：
$$b{x}_2{y}_1-b{x}_1{y}_2=x_2-x_1$$
令$k=x_2{y}_1-x_1{y}_2$，则有：
$$\begin{gathered} x_2=x_1+kb \\ {y}_2=y_1-ka \end{gathered}$$
这样我们便证明了解$(x_2,y_2)$可以由$(x_1,y_1)$得到，因此这种方法可以找到$ax+by=1$的每一个解。

  那么如果$gcd(a,b)>1$呢？和简单，我们只需要在方程两端同时除以$g=gcd(a,b)$不就变成了上面这种情况了嘛。即：
$$\frac{a}{g}x+\frac{b}{g}y=1$$
然后我们直接套用前面的方法，就可以完成对$ax+by=g$的解的描述了。

4 线性方程定理

  设$a$与$b$是非零整数，$g=gcd(a,b)$。方程：
$$ax+by=g$$
总是有一个整数解$(x_1,y_1)$，它可以由欧几里得算法得到，则方程的每一个解可以由：
$$(x_1+k\frac{b}{g},y_1-k\frac{a}{g})$$
得到，其中$k$可为任意整数。

5 扩展欧几里得算法

  给定$ax+by=g$，扩展欧几里得算法可求得上述整数解。来看下面两个式子：
$$\begin{gathered} a{x}_1+b{y}_1=gcd(a,b) \\ b{x}_2+a\%b{y}_2=gcd(b,a\%b) \end{gathered}$$
如果我们能够知道$x_1,x_2,y_1,y_2$之间的关系，那么我们就可以递归求解。我们知道$a\%b=a-\lfloor \frac{a}{b}\rfloor \ast b$
，且$gcd(a,b)=gcd(b,a\%b)$，通过上式可求得：
$$\begin{gathered} x_1=y_2 \\ {y}_1=x_2-\lfloor \frac{a}{b}\rfloor \ast y_2 \end{gathered}$$
欧几里得算法的停止状态是$a=gcd,b=0$，放到这里的状态就是：
$$a\ast 1+b\ast 0=a$$
不管$b$的系数是多少，但$a$的系数一定得是$1$。现在我们可以令上式$x_2=1,y_2=0$然后一直往回推，就可以得到最初状态了。

def ext_gcd(a, b):
    if not b:
        return 1, 0, a
    x, y, g = ext_gcd(b, a % b)
    x, y = y, x - a//b * y
    return x, y, g

print(ext_gcd(60, 22))

输出结果为

(-4, 11, 2)

6 裴蜀定理（贝祖定理）

  在数论中，裴蜀定理是一个关于最大公约数的定理。裴蜀定理得名于法国数学家艾蒂安·裴蜀，说明了对任何整数$a,b$和它们的最大公约数$g$，关于未知数$x$和$y$的线性丢番图方程（称为裴蜀等式）：
$$ax+by=m$$
有解当且仅当$m$是$g$的倍数。裴蜀等式有解时必然有无穷多个整数解，每组解$x,y$都称为裴蜀数，可用扩展欧几里得算法求得。

  特别来说，方程 $ax+by=1$ 有解当且仅当整数$a$和$b$互素。

  裴蜀等式也可以用来给最大公约数定义：$g$其实就是最小的可以写成$ax+by$形式的正整数。这个定义的本质是整环中“理想”的概念。因此对于多项式整环也有相应的裴蜀定理。

7 参考资料

• 《数论概论》第四版 P24-P28
• 百度百科 裴蜀定理