使用扩展欧几里得算法对逆元求解

扩展欧几里得算法求逆元

解方程 ax + by = gcd(a, b)

在欧几里得算法中通过定理gcd(a, b) = gcd(b, a%b)，我们使用递归求得a与b的最大公约数，在递归边界当b = 0时，a = gcd(a, b)，此时显然有a * 1 + b * 0 = gcd。对于这样的一个式子，我们能否对它进行推广呢？

• 在递归的其中一步，通过相关定理（此处不叙述）可知 $a{x}_1+b{y}_1=gcd(a,b)$中解$x_1、y_1$一定存在
• 下一步计算$gcd(b,a\%b)$时，又有$b{x}_2+(a\%b)y_2=gcd(b,a\%b)$成立
• $而gcd(a,b)=gcd(b,a\%b)$
• 所以$a{x}_1+b{y}_1=b{x}_2+(a\%b)y_2$成立，这样我们就发现$x_1、x_2、y_1、y_2$之间存在关系
• $a\%b=a-(a/b)\ast b$成立（想一想为什么?）
• 代入得$a{x}_1+b{y}_1=b{x}_2+(a-(a/b)\ast b)y_2$
• 整理得$a{x}_1+b{y}_1=a{y}_2+b(x_2-(a/b)y_2)$
• 由于对于递归过程中任一步都满足，对比两式，可得递归表达式

$$\{\begin{array}{c}{x}_1=y_2\\ y_1=x_2-(a/b)y_2\end{array}$$

• 已知递归边界a * 1 + b * 0 = gcd，递归表达式，能够递归得到ax + by = gcd(a, b)的其中一个解

//求最大公约数的同时求ax + by = gcd(a, b)的解
//使用引用，exGcd函数结束时，x,y中就是所求解
int exGcd(int a, int b, int &x, int &y) {
    //递归边界
    if (b == 0) {
        x = 1; y = 0; return a;
    }

    //递归计算最大公约数gcd
    int gcd = exGcd(b, a % b, x, y);

    //递推公式，求解
    int temp = x;
    x = y;
    y = temp - a / b * y;

    return gcd;
}

乘法逆元的求解

• 由定义知，求a模m的逆元，就是求解同余式$ax\equiv 1(modm)$
• $ax\equiv 1(modm)$
• $(ax-1)\%m=0$
• 存在一个整数z满足$ax-1=mz$
• $ax-mz=1$
• 存在y = -z使得$ax+my=1$
• 所以求a模m的逆元，就是求解同余式$ax\equiv 1(modm)$，也就是求方程$ax+my=1$的解，并且在实际使用中，一般把x的最小正整数解称为a模m的逆元

//返回值是a模m的逆元
int inv(int a, int m) {
    int x;
    int y;
    int gcd = exGcd(a, m, x, y);    //此时得到的x是方程的一个解，但不一定是方程的最小正整数解，x可能为负
    return (x % m + m) % m;         //(x % m + m) % m 是方程最小正整数解，也就是a模m的逆元
}