扩展欧几里得以及求逆元算法

首先，我们先了解一下裴蜀定理：

设 a, b是不全为零的整数，则存在整数x,y , 使得ax + by = gcd(a, b) .

用于求解方程 ax+by=gcd(a,b)的解

当 b=0时

ax+by=a

故而 x=1,

y=0
当 b≠0时

因为
gcd(a,b)=gcd(b,a%b)
而
bx′+(a%b)y′=gcd(b,a%b)
bx′+(a−⌊a/b⌋∗b)y′=gcd(b,a%b)
ay′+b(x′−⌊a/b⌋∗y′)=gcd(b,a%b)=gcd(a,b)
故而

x=y′,        y=x′−⌊a/b⌋∗y′

#include <iostream>

using namespace std;

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;//迭代完成后原来的x变成了y，原来的y变成了x
    return d;//一系列迭代之后最终d还是最小公约数，但是x，y是求出的最小的满足ax+by=gcd(a，b)的结果
}

int main()
{
    int n;
    scanf("%d", &n);
    
    while(n--)
    {
        int a, b;
        scanf("%d%d", &a, &b);
        int x, y;
        exgcd(a, b, x, y);
        printf("%d %d\n", x, y);
    }
    return 0;
}

题目：线性同余方程

#include <iostream>

using namespace std;
typedef long long LL;

// a * x ≡ b (mod m)

// 变形为拓展欧几里得形式：a * x + b * y = gcd(a, b)
// 原式变为：   a * x = m * y + b  (注：mod m 为 b, 则相当于结果为 m 的倍数和 b 的和)
//              a * x - m * y = b
// 另y1 = -y得：a * x + m * y1 = b
// 根据拓展欧几里得定理(也可以说裴蜀定理)，只要 b 是 gcd(a, m)的倍数即有解！
// 另d = gcd(a, m), 我们得到的式子其实是：a * x + m * y1 = gcd(a, m) = d (注；上面的b其实就是d的倍数)
// 所以左右同乘 b / d 即可转化为：a * x * b / d + m * y1 * b / d = b * b / d = b
// 即最后答案为：res = x * d / b % m

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;
}

int main()
{
    int n;
    cin >> n;
    while(n--)
    {
        int a, b, x, y, m;
        cin >> a >> b >> m;
        int d = exgcd(a, m, x, y);
        //只要b是d的倍数就有解
        if(b % d)
            printf("impossible\n");
        else
            printf("%d\n", (LL) x * b / d % m);//模m，不然可能结果会爆int
    }
    return 0;
}