[数论+模板] 中国剩余定理+表达整数的奇怪方式(模板+数论难题)

1. 中国剩余定理证明+模板

204. 表达整数的奇怪方式

重点： 中国剩余定理、同余方程

这道题目很难。目前来讲。 貌似和中国剩余定理没啥用到…是个等价变形。

强烈推荐大佬题解

模板代码：

#include <iostream>
#include <algorithm>

using namespace std;

typedef long long LL;

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;
    
    LL x = 0, a1, m1;   // x作为当前是否无解的判断
    cin >> a1 >> m1;    // 读入第一个方程
    
    for (int i = 0; i < n - 1; ++i) {   // 读入后续n-1个方程
        LL a2, m2;
        cin >> a2 >> m2;
        
        LL k1, k2;  
        LL d = exgcd(a1, -a2, k1, k2);  // 扩欧计算k1, k2,在此a2正负无所谓,不影响最大公约数的计算
        
        if ((m2 - m1) % d) {            // 扩欧有解的充要条件
            x = -1;
            break;
        }
        
        k1 *= (m2 - m1) / d;            // 等式恒等变换,令m2-m1为等式右边项，扩欧是d为等式右边项
        k1 = (k1 % (a2/d) + a2/d) % (a2/d); // 令k1变成方程的最小正整数解
        
        x = k1 * a1 + m1;               // x的所有解,此时并不为最小非负整数
        
        LL a = abs(a1/d * a2);          // 临时变量a存一下变化a的值，要为正值，就是lcm(a1,a2)
        
        m1 = k1 * a1 + m1;              // 更新m1，其实就是x
        a1 = a;                         // a1 更新为a
    }

    if (x != -1) x = (x % a1 + a1) % a1;  // 最小正余数
    
    cout << x << endl;
    
    return 0;
}

// 理一下思路：
// 每次都是两个式子合并，关键是用exgcd计算得到k1,k2
// 但是我们在此仅使用k1即可，k2是不需要被用到的，但是没有k2是计算不出来k1的
// 所以每次都需要两个式子，这是相辅相成的
// 经过式子的等价变形，使得具有相同形式，那么就能重复合并操作
// 最后得到一个式子，且知道k,a,x0都已知，求x最小整数解也是水到渠成