「扩展欧几里得」学习笔记

$\mathtt{date}:\mathrm{2023.1.15}$

如果觉得写得好就给个赞吧
博客没有赞，作者两行泪

---

为了介绍扩展欧几里德，我们先介绍一下贝祖定理：

即如果 $a,b$ 是整数，那么一定存在整数 $x,y$ 使得 $ax+by=gcd(a,b)$

可以得到：如果 $ax+by=c$ 有解，那么 $gcd(a,b)|c$ 。

可以采用构造性的数学归纳法证明。

• 当 $b=0$ 时，$gcd(a,b)=a$。此时 $x=1,y=0$。
• 当 $b$ 不为 $0$ 时,根据欧几里得定理可得,即 $ax+by=b{x}^{\prime }+(amodb)y^{\prime }=b{x}^{\prime }+(a-b\times \lfloor \frac{a}{b}\rfloor )y^{\prime }$
• 移项得 $ax+by=b{x}^{\prime }+(amodb)y^{\prime }=a{y}^{\prime }+b(x^{\prime }-\lfloor \frac{a}{b}\rfloor )y^{\prime }$
• 于是就有：
  $$\{\begin{array}{l}x=y^{\prime }\\ y=x^{\prime }-\lfloor \frac{a}{b}\rfloor y^{\prime }\end{array}$$

我们可以在递归的最内层（即 $b=0$ 时）给 $x,y$ 赋值，也就是 $x=1,y=0$，这样就可以把每一层的 $x,y$ 都推算出来。

于是就有了代码：

void exgcd(int a, int b, int &x, int &y) {
    if(b==0) {
        x=1; y=0; // 赋值
        return;
    }
    exgcd(b, a%b, x, y); // 先递归到最里层
    int tmp=x;
    x=y;
    y=tmp-(a/b)*y;
}