[AcWing]204. 表达整数的奇怪方式（C++实现）求解线性同余方程组的最小非负整数x（扩展中国剩余定理）模板题

1. 题目

2. 读题（需要重点注意的东西）

思路：

中国剩余定理
中国剩余定理是为了解如下线性同余方程组：
以上只为了解释中国剩余定理，本题的主要思路见下方。

---

本题的思路（扩展中国剩余定理）
核心思想：将线性同余方程组按如下方式两两合并即可！

3. 解法

---------------------------------------------------解法---------------------------------------------------

#include <iostream>
#include <algorithm>

using namespace std;

typedef long long LL;


// 扩展欧几里得算法求最大公约数和k1、k2
LL exgcd(LL a, LL b, LL &x, LL &y)
{
    if (!b)
    {
        x = 1, y = 0;
        return a;
    }

    LL d = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}


int main()
{
    int n;
    cin >> n;
    
    // 现有方程m1、a1
    LL x = 0, a1, m1;
    cin >> a1 >> m1;
    
    // -----核心：用现有的m1和a1合并m和a--------
    for (int i = 0; i < n - 1; i ++ )
    {
        LL a2, m2;
        cin >> a2 >> m2;
        
        // 算出一组k1、k2
        LL k1, k2;
        // d是a1和a2的最大公约数
        LL d = exgcd(a1, a2, k1, k2);
        
        // 如果无解
        if ((m2 - m1) % d)
        {
            x = -1; // 无解，输出-1
            break;
        }
        
        k1 *= (m2 - m1) / d;
        // 避免溢出，将k1变成最小的正整数解
        k1 = (k1 % (a2/d) + a2/d) % (a2/d);
        
        // 
        x = k1 * a1 + m1;
        // 求a1和a2的最小公倍数
        LL a = abs(a1 / d * a2);
        m1 = k1 * a1 + m1; // 将k1 * a1 + m1赋值给m1,用于下一次合并
        a1 = a; // 将a1和a2的最小公倍数再赋值给a1,用于下一次合并
    }
    
    // 求出最小正整数解
    if (x != -1) x = (m1 % a1 + a1) % a1;

    cout << x << endl;

    return 0;
}

---

4. 可能有帮助的前置习题

• [AcWing]877. 扩展欧几里得算法（C++实现）扩展欧几里得算法模板题

5. 所用到的数据结构与算法思想

• 中国剩余定理
• 扩展欧几里得算法

6. 总结

求解线性同余方程组的最小非负整数x（扩展中国剩余定理）模板题，理解公式并背下代码。