扩展欧几里得(逆元)

/*
    "数学的恐怖qwq"
    想了半天终于明白了, 这里尽量通俗的写出来
    扩展欧几里得算法有很多版本
    这里写两个, 选择喜欢的使用
    (注意本证明极其复杂, 请直接背模版和感性理解)
    
    扩欧可以解决两项未知解, 具体原理来自裴蜀定理
    裴蜀定理:设 a,b 是不全为零的整数, 则存在整数 x,y, 使得 ax+by=gcd(a,b).
    而扩欧就可以在求出gcd的过程中把其中的x和y求出来(一种方案)
    就可以得到式子ax + by = gcd(a, b);中所有的数值
    因为知道了gcd(a, b)就可以通过这样求出每个二元一次方程组的一个解
    例如已知a, b, c, 求满足ax + by = c中想x, y的一个解
    通过扩欧可以求出a*x1 + b*y1 = gcd(a, b);
    这个式子左右两边都乘上c / gcd(a, b)；
    就可以得到a*x2 + b*y2 = c; 这时候x2和y2我们知道, 而这就是ax + by = c的一个解
    通过这个还可以求逆元
    
    而具体实现
    普通欧几里得如下
    int gcd(int a, int b)
    {
        return b ? gcd(b, a % b) : a;
    }
    
    或者应该这样
    int gcd(int a, int b)
    {
        if (b == 0)
        {
            return a;
        }
        int d = gcd(b, a % b);
        return d;
    }
    
    我们可以多带两个值x, y, 因为递归下一层内的数组没法直接传上来, 这里就通过引用的方式, 
    取址传上来, 不是取值
    就变成
    int ecgcd(int a, int b, int &x, int &y)
    {
        if (b == 0)
        {
            return a;
        }
        int d = ecgcd(b, a % b, ?, ?); // 这两个问号下面会补上
        return d;
    }
    那俩问号暂且不管
    
    那么问题来了, 怎么求x, y, 
    首先你要知道, 只有一个等式是无法求出二元一次方程的唯一解的
    所以我们的解有无穷个, 而对于最后终止的那里, 也就是b == 0的情况(gcd(a, b), 注意对应)
    ax + by = d, 因为b == 0, 所以a = d(最后这里的a是d是最大公约数)
    所以x = 1, 这样就确定了x, 但是y是可以变化的, 这也是为什么有无数解的原因
    因此这里的y可以随便去一个, 通常选y = 0(为什么下面会说),
    
    // 请把所有的x和y看成未知数, 不要想位置关系, 这里只是为了方便
    
    那么终止的x, y选完了, 接着就是看看怎么把x和y带回去
    原式子是ax + by = d; (1)
    它的递归下一层就是b*y1 + a%b * x1 == d (2)
    a % b == a - a/b * b (3)
    由(1)(3)可得 (a%b + a/b * b)*x + by = d (4)
    化简得 a%b * x + b(y + a/b * x) = d (5)
    通过(2)和(5) 可以发现, x1 = x, y1 = y + a/b * b (6)
    y = y1 - a/b * x, x = x1; (7)
    
    而x1和y1就是我想通过x和y传回来的值, 这时候看看上面两个问号, 就有两种填法
    ecgcd(b, a % b, x, y) 和 ecgcd(b, a % b, y, x)
    
    如果是ecgcd(b, a % b, x, y)
    那么本次递归ax + by = d, 下次b * x2 + a%b * y2 = d; 注意这里和(2)不一样
    传回来的x = x2 = y1, y = y2 = x1; 注意这里x和y数值上是下一层递归里面传上来的
    根据 y = y1 - a/b * x, x = x1; (7) 可得
    本次递归里的y = 下个递归里的x2 - 本次的a/b * 本次的x, 本次的x = 下次的y2
    而在实现上就是int t = y; y = x - a / b * t; x = t;
    这里有点难懂, 要注意x和y数值上是下一层递归里面传上来的, 用的还是本次递归的x, y
    作为容器, 我们要利用x, y, 更改他们自己
    
    具体代码就是
    int exgcd(int a, int b, int &x, int &y)
    {
        if (b == 0)
        {
            x = 1;
            y = 0;
            return a;
        }
        int d = exgcd(b, a % b, x, y);
        int t = y;
        y = x - a / b * t;
        x = t;
        return d;
    }
    
    如果是如果是ecgcd(b, a % b, y, x)
    那么本次递归ax + by = d, 下次b * y3 + a%b * x3 = d; // 这里的x3 y3都是未知数, 没什么具体含义, 所有的关系就是下面的那一句话
    传回来的y = y3 = y1, x = x3 = x1;
    根据 y = y1 - a/b * x, x = x1; (7) 可得
    本次的y = 下次的y3 - 本次的a/b * 本次的x, 本次的x = 下次的x3
    实现就是 y = y - a / b * x, x = x ==> y -= a / b * x; 
    唉?这时候我们传上来的x恰好等于这次的x, x就不用动了
    y的话就相当于自减, 非常巧妙
    
    实现
    int exgcd(int a, int b, int &x, int &y)
    {
        if (b == 0)
        {
            x = 1;
            y = 0;
            return a;
        }
        int d = exgcd(b, a % b, y, x);
        y -= a / b * x;
        return d;
    }
    
    到这里就解释的差不多了, 注意=的关系, 

    b == 0时最好y = 0, 因为我们会把y传回去, 他会经历一系列变化
    取0的话可以防止出现int爆炸的情况, 不然你可以试试这题y = 100
    至此此题结束
    
*/
#include <iostream>
#include <algorithm>
#include <cstring>

using namespace std;

int n, m;

// int exgcd(int a, int b, int &x, int &y)
// {
//     if (b == 0)
//     {
//         x = 1;
//         y = 0;
//         return a;
//     }
//     int d = exgcd(b, a % b, x, y);
//     int t = y;
//     y = x - a / b * t;
//     x = t;
//     return d;
// }

int exgcd(int a, int b, int &x, int &y)
{
    if (b == 0)
    {
        x = 1;
        y = 0;
        return a;
    }
    int d = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}

int main()
{
    int T;
    cin >> T;
    while (T -- )
    {
        int x, y;
        scanf("%d%d", &n, &m);
        exgcd(n, m, x, y);
        printf("%d %d\n", x, y);
    }
    
    return 0;
}

求逆元

什么是逆元, 逆元就是形如 a * x ≡ 1 (mod p) 其中x就叫做a对于p的一个逆元
相对的a就是x关于p的一个逆元

知道了定义 a * x ≡ 1 (mod p) 这个式子我们可以把他展开
变成a * x + p * y = 1 (y是任意整数)
这个式子形如 ax + by = d
所以这个式子可以用exgcd进行求解
也就是求出里面的x

求逆元还有一种方法, 也是一种特殊的情况, 当且仅当p为质数且(a, p) = 1
这时候它满足费马小定理 ：如果p是一个质数，而整数a不是p的倍数，
也就是当且仅当p为质数且(a, p) = 1, 则有a ^ (p - 1) ≡ 1 (mod p)
这时候x也就是a的逆元就是a ^ (p - 2) ,
而这个东西可以通过快速幂求解
这就是逆元的特殊求法

欧拉定理 & 费马小定理 - OI Wiki (oi-wiki.org)

exgcd通解
欧几里得算法