表达整数的奇怪方式(中国剩余定理)

题目链接
题意：给定2n个整数a1,a2,…,an和m1,m2,…,mn,求一个最小的非负整数x，满足∀i∈[1,n],x≡mi(mod ai)。
输入格式
第1行包含整数n。
第2…n行：每i+1行包含两个整数ai和mi,数之间用空格隔开。
输出格式
输出最小非负整数x，如果x不存在，则输出-1。
如果存在x，则数据保证x一定在64位整数范围内。
数据范围
1≤ai≤2^31−1,
0≤mi<ai
1≤n≤25
输入样例：
2
8 7
11 9
输出样例：
31
思路：要求找到最小的一个x满足∀i∈[1,n],x % ai =mi,直接就有个中国剩余定理（前提条件模数ai必需两两互质）

• A = a1 * a2 * ……*an
• Ai = A / ai
• Ai^-1 表示 Mi模mi的逆元
• x = m1 * A1 *A1^-1 + m2 * A2 * A2^-1……mn * An * An^-1

比如把x带入x ≡ m1(mod a1)中证明：因为A1*A1^-1相当于1，而后面的A2,A3,An模上a1都是0，所有这个x对于每个表达式成立。
显然这个题并不满足互质的条件，就得从定理的根本原理来下手了。
转换一下：
x = k1 * a1 + m1
x = k2 * a2 + m2
k1 * a1 - k2 * a2 = m2 - m1
( a1,a2,m1,m2都已知，相当于求a1 * x + a2 * y = c；就是扩展欧几里得了啊。)
如果gcd(a1, a2) | m2 - m1，便有解
整理下就有:
x = a1 * k1 + m1 + k * gcd(a1, a2)
x = x0 + k * a (用x0和a来替代)，显然可以把n个式子两两结合最后变成一个 这个时候要求的x % a = x0.
代码实现：

#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; 
    scanf("%d", &n);
    ll a1, m1;
    cin >> a1 >> m1;
    bool has_answer = 1;
    for(int i = 0; i < n - 1; i ++){
        ll a2, m2;
        cin >> a2 >> m2;
        ll k1, k2;
        ll d = exgcd(a1, a2, k1, k2);
        if((m2 - m1) % d){
            has_answer = 0;
            break;
        }
        k1 *= (m2 - m1) / d;
        ll t = a2 / d;
        //避免k1很大，让它变成最小的一个正解
        k1 = (k1 % t + t) % t;
        m1 = a1 * k1 + m1;
        a1 = abs(a1 / d * a2);
    }
    if(has_answer){
        //取一个正的余数
        cout << (m1 % a1 + a1) % a1 << endl;
    }
    else cout << -1 << endl;
    return 0;
}