数学 - 扩展欧几里得算法

扩展欧几里得算法用于解决二元一次不定方程、一次同余方程、逆元等问题

二元一次不定方程

• 二元一次不定方程的定义：$$ax+by=c，(a,b,c\in Z且a,b\ne 0)\text{\hspace{1em}(1)}$$
• 如果$c=0$，则方程的整数解为：$x=-bt，y=at，(t\in Z)$
• 如果$c$为负数，则(1)式两边乘以$-1$，等号右边就变成正整数；因为所求的$x$，$y$可正可负，所以我们只用讨论$a$，$b$是正整数的情况。
• 设$d$为$a$和$b$的最大公因数，即 $d=gcd(a,b)，a=a_1\ast d，b=b_1\ast d，(a_1,b_1)=1$，则(1)式变为$$d\ast (a_{1}x+b_{1}x)=c\text{\hspace{1em}(2)}$$

结论：当 $d\nmid c$ 时，(1)式没有整数解；当 $d\mid c$ 时，(1)式一定有整数解，这可以通过贝祖定理（又称裴蜀定理）得到。

贝祖定理：如果$a$和$b$是整数，且 $gcd(a,b)=d$，那么对于任意的整数$x，y，ax+by$都一定是$d$的倍数；特别的，一定存在整数$x，y$，使$ax+by=d$成立。

• 贝祖定理的一个重要推论是：$a$和$b$互质的充分必要条件是存在整数$x，y$，使$ax+by=1$
• $ax+by$能够得到的最小结果就是它们的最小公倍数

我们通过几个例子，解释$x$和$y$的求解步骤，并在此基础上得出扩展欧几里得算法。

例1：$17x+5y=1$

$\because gcd(17,5)=1，\therefore 方程有解$

$(3\ast 5+2)x+5y=1$，整理得：
$5(3x+y)+2x=1$
用换元法，如果我们求出如下方程的解，
$5x^{{}^{\prime }}+2y^{{}^{\prime }}=1$，则：$x=y^{{}^{\prime }}，y=x^{{}^{\prime }}-3y^{{}^{\prime }}$

$(2\ast 2+1)x^{{}^{\prime }}+2y^{{}^{\prime }}=1$，整理得：
$2(2x^{{}^{\prime }}+y^{{}^{\prime }})+x^{{}^{\prime }}=1$
用换元法，如果我们求出如下方程的解，
$2x^{{}^{\prime \prime }}+y^{{}^{\prime \prime }}=1$，则：$x^{{}^{\prime }}=y^{{}^{\prime \prime }}，y^{{}^{\prime }}=x^{{}^{\prime \prime }}-2y^{{}^{\prime \prime }}$

$(2\ast 1+0)x^{{}^{\prime \prime }}+y^{{}^{\prime \prime }}=1$，整理得：
$1(2x^{{}^{\prime \prime }}+y^{{}^{\prime \prime }})+0x^{{}^{\prime \prime }}=1$
用换元法，如果我们求出如下方程的解，
$1x^{{}^{\prime \prime \prime }}+0y^{{}^{\prime \prime \prime }}=1$，则：$x^{{}^{\prime \prime }}=y^{{}^{\prime \prime \prime }}，y^{{}^{\prime \prime }}=x^{{}^{\prime \prime \prime }}-2y^{{}^{\prime \prime \prime }}$

显而易见，方程 $1x^{{}^{\prime \prime \prime }}+0y^{{}^{\prime \prime \prime }}=1$ 有解$x^{{}^{\prime \prime \prime }}=1，y^{{}^{\prime \prime \prime }}=0$。
$x^{{}^{\prime \prime \prime }}$的系数$1$，就是原方程中$17$和$5$的最大公约数

回溯可得，
$x^{{}^{\prime \prime }}=0，y^{{}^{\prime \prime }}=1$
$x^{{}^{\prime }}=1，y^{{}^{\prime }}=-2$
$x=-2，y=7$

例2：$27x+21y=3$

$\because gcd(27,21)=3，\therefore 方程有解$

$(1\ast 21+6)x+21y=3$，整理得：
$21(x+y)+6x=3$
用换元法，如果我们求出如下方程的解，
$21x^{{}^{\prime }}+6y^{{}^{\prime }}=3$，则：$x=y^{{}^{\prime }}，y=x^{{}^{\prime }}-y^{{}^{\prime }}$

$(3\ast 6+3)x^{{}^{\prime }}+6y^{{}^{\prime }}=3$，整理得：
$6(3x^{{}^{\prime }}+y^{{}^{\prime }})+3x^{{}^{\prime }}=3$
用换元法，如果我们求出如下方程的解，
$6x^{{}^{\prime \prime }}+3y^{{}^{\prime \prime }}=3$，则：$x^{{}^{\prime }}=y^{{}^{\prime \prime }}，y^{{}^{\prime }}=x^{{}^{\prime \prime }}-3y^{{}^{\prime \prime }}$

$(2\ast 3+0)x^{{}^{\prime \prime }}+3y^{{}^{\prime \prime }}=3$，整理得：
$3(2x^{{}^{\prime \prime }}+y^{{}^{\prime \prime }})+0x^{{}^{\prime \prime }}=3$
用换元法，如果我们求出如下方程的解，
$3x^{{}^{\prime \prime \prime }}+0y^{{}^{\prime \prime \prime }}=3$，则：$x^{{}^{\prime \prime }}=y^{{}^{\prime \prime \prime }}，y^{{}^{\prime \prime }}=x^{{}^{\prime \prime \prime }}-2y^{{}^{\prime \prime \prime }}$

显而易见，方程 $3x^{{}^{\prime \prime \prime }}+0y^{{}^{\prime \prime \prime }}=3$ 有解$x^{{}^{\prime \prime \prime }}=1，y^{{}^{\prime \prime \prime }}=0$。
$x^{{}^{\prime \prime \prime }}$的系数$3$，就是原方程中$27$和$21$的最大公约数

回溯可得，
$x^{{}^{\prime \prime }}=0，y^{{}^{\prime \prime }}=1$
$x^{{}^{\prime }}=1，y^{{}^{\prime }}=-3$
$x=-3，y=4$

对于方程：$ax+by=gcd(a,b)$

$(\lfloor a/b\rfloor \ast b+a\%b)x+by=gcd(a,b)$，整理得：
$b(\lfloor a/b\rfloor \ast x+y)+a\%b\ast x=gcd(a,b)$
用换元法，如果我们求出如下方程的解，
$b\ast x^{{}^{\prime }}+a\%b\ast y^{{}^{\prime }}=gcd(a,b)$，则：$x=y^{{}^{\prime }}，y=x^{{}^{\prime }}-\lfloor a/b\rfloor \ast y^{{}^{\prime }}$

不断递归，最后我们得到方程：
$gcd(a,b)\ast \tilde{x}+0\ast \tilde{y}=gcd(a,b)$，

可得该方程的一组解：$\tilde{x}=1，\tilde{y}=0$。
$\tilde{x}$的系数，就是原方程中$a$和$b$的最大公约数

不断回溯，最后得到原方程的解

扩展欧几里得算法的参考代码：

//二元一次不定方程：ax + by = gcd(a, b)
//函数返回 gcd(a, b)；通过引用传递，返回一组解 (x, y)
int exgcd(int a, int b, int &x, int &y) {	
    if (b == 0) {
        x = 1, y = 0;
        return a;
    }

    int gcd = exgcd(b, a % b, x, y);
    int t = x;
    x = y;
    y = t - a / b * y;
    return gcd;
}

总结：

对于一般情况，$ax+by=c$，($a,b,c$为正整数)，设 $d=gcd(a,b)$

• 如果$d\nmid c$，方程无解
• 如果$d\mid c$，
  • 先求出方程$ax+by=d$的一组解 $x_0,y_0$
  • 原方程的一组解为 $x^{{}^{\prime }}=x_0\ast c/d，y^{{}^{\prime }}=y_0\ast c/d$
  • 原方程的通解为 $x=x^{{}^{\prime }}+k\ast b/d，y=y^{{}^{\prime }}-k\ast a/d$，$k$为任意整数
  • 设$p=b/d$，$x$的最小非负整数解为 $(x^{{}^{\prime }}\%p+p)\%p$
  • 设$q=a/d$，$y$的最小非负整数解为 $(y^{{}^{\prime }}\%q+q)\%q$
  • 如果$a$或$b$是负数，则函数返回的最大公约数可能为负数，将$p$取绝对值，才能用以上方法得到$x$的最小非负整数解

求解一次同余方程 $ax\equiv b(modm)$

• 一次同余方程亦称线性同余方程，指未知数仅出现一次幂的同余方程。
• 该同余方程等价于 $ax=m\ast (-y)+b$，即$ax+my=b$
• 如果$gcd(a,m)\nmid b$，该方程无解
• 如果$gcd(a,m)\mid b$，用上述方法求解

求解乘法逆元

• 对于正整数$a,m$，如果有$ax\equiv 1(modm)$，则称$x$的最小正整数解为$a$模$m$的逆元，也叫做数论倒数，记作$a^{-1}$。$a$和$x$一定都和$m$互质，如果$a$和$m$不互质就不存在逆元。
• 乘法运算中$a\ast a^{-1}=1$，模运算中，$a\ast a^{-1}\equiv 1(modm)$

扩展欧几里得算法求逆元：

• 见求解一次同余方程

费马小定理求逆元：

• 如果模数$m$为质数，可利用费马小定理求逆元：$a^{m-1}\equiv 1(modm)$
• 因为$a\ast a^{m-2}\equiv 1(modm)$，所以$a$的逆元是$a^{m-2}\%m$，用快速幂计算
• 例如，$m$为7时，$3\ast 3^5\equiv 1(mod7)$，$3^5\%7$为$3$模$7$的逆元

欧拉函数求逆元：

• 欧拉定理：如果$a$和$m$互质，设$\varphi (m)$为小于$m$且与$m$互质的正整数的个数，则有$a^{\varphi (m)}\equiv 1(modm)$
• 因为$a\ast a^{\varphi (m)-1}\equiv 1(modm)$，所以$a$的逆元是$a^{\varphi (m)-1}\%m$，用快速幂计算
• $m$可以不是质数
• 例如，$m$为8时，$\varphi (8)=4$，$3\ast 3^3\equiv 1(mod8)$，$3^3\%8$为$3$模$8$的逆元

逆元的用处：

• 取余运算有如下性质
  • $(a+b)\%m=(a\%m+b\%m)\%m$
  • $(a-b)\%m=(a\%m-b\%m)\%m$
  • $(a\ast b)\%m=(a\%m\ast b\%m)\%m$
  • $a^b\%m=(a\%m{)}^b\%m$
• 但对除法不成立
  • $(a/b)\%m\ne (a\%m/b\%m)\%m$
  • 例如，$(8/2)\%7=4$，而 $(8\%7)/(2\%7)\%7=0$
• 如果$b\mid a$，则$(a/b)\%m=((a/b)\ast (b\ast b^{-1}))\%m=(a/b\ast b)\ast (b^{-1})\%m=(a\ast b^{-1})\%m$。这样我们就把除法取余运算转化成了乘法取余运算。
• 例如，求3的幂的和