扩展中国剩余定理

题意： 给定$n$个整数$a_1,a_2,...,a_n$和$n$个整数$m_1,m_2,...m_n$，$\forall i\in [1,n],x\equiv m_i(\mathrm{m}\mathrm{o}\mathrm{d}{a}_i)$。求存在的最小$x$
题解：
依次合并每个同余式。
下面以两个同余式为例：
$x\equiv m_1(\mathrm{m}\mathrm{o}\mathrm{d}{a}_1)$
$x\equiv m_2(\mathrm{m}\mathrm{o}\mathrm{d}{a}_2)$

等价于：
$x=a_1{k}_1+m_1$ ①
$x=a_2{k}_2+m_2$ ②

联立可得：
$a_1{k}_1+m_1=a_2{k}_2+m_2$

移项得：
$a_1{k}_1-a_2{k}_2=m_2-m_1$

可以看出来是$exgcd$的标准形式，故有解的条件是$gcd(a_1,a_2)\mid (m_2-m_1)$
目标为消去不定因素$k_1,k_2$

设$d=gcd(a1,a2)$
两边同除$d$得：
$\frac{a_1}{d}{k}_1-\frac{a_2}{d}{k}_2=\frac{m_2-m_1}{d}$

等价于：
$\frac{a_1}{d}{k}_1\equiv \frac{m_2-m_1}{d}(\mathrm{m}\mathrm{o}\mathrm{d}\frac{a_2}{d})$

到这里我们已经消除了$k_2$，继续消除$k_1$即可
对同余式两边同除$\frac{a_1}{d}$，得：
$k_1\equiv inv(\frac{a_1}{d},\frac{a_2}{d})\times \frac{m_2-m_1}{d}(\mathrm{m}\mathrm{o}\mathrm{d}\frac{a_2}{d})$
其中$inv(a,b)$表示$a$在模$b$意义下的逆元
将$k_1$用①式替换，得：
$\frac{x-m_1}{a_1}\equiv inv(\frac{a_1}{d},\frac{a_2}{d})\times \frac{m_2-m_1}{d}(\mathrm{m}\mathrm{o}\mathrm{d}\frac{a_2}{d})$
等价于：
$x\equiv inv(\frac{a_1}{d},\frac{a_2}{d})\times \frac{m_2-m_1}{d}\times a_1+m_1(\mathrm{m}\mathrm{o}\mathrm{d}\frac{a_1{a}_2}{d})$
至此，$x$已经变成了完全由已知的$a_1,a_2,m_1,m_2$所表达的同余式了。
形式与$x\equiv m(\mathrm{m}\mathrm{o}\mathrm{d}a)$相同
那么合并后的$m=inv(\frac{a_1}{d},\frac{a_2}{d})\times \frac{m_2-m_1}{d}\times a_1+m_1$，$a=lcm(a_1,a_2)$、
其中$inv(\frac{a_1}{d},\frac{a_2}{d})\times \frac{m_2-m_1}{d}$相当于求解$a_1{k}_1-a_2{k}_2=m_2-m_1$得到的最小非负$k_1$

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long ll;
int n;

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

ll gcd(ll a, ll b) {return b == 0 ? a : gcd(b, a % b);}
ll lcm(ll a, ll b) {return a / gcd(a, b) * b;}

int main()
{
    scanf("%d", &n);
    
    int flag = 1;
    ll a1, m1;
    for(int i = 1; i <= n; i++) {
        if(i == 1) scanf("%lld%lld", &a1, &m1);
        else {
            ll a2, m2;
            scanf("%lld%lld", &a2, &m2);
            ll k1, k2;
            ll d = exgcd(a1, a2, k1, k2); //求出两个mod的gcd
            //判断是否有解
            if((m2 - m1) % d) flag = 0;
            else {
                k1 *= (m2 - m1) / d; //扩大(m2-m1)/d倍
                ll t = a2 / d;//求出模:a2/d
                k1 = (k1 % t + t) % t;//求出最小的非负k1
                
                m1 = a1 * k1 + m1;//合并后的m
                a1 = abs(a1 / d * a2);//合并后的a
            }
        }
    }
    if(flag) printf("%lld\n", (m1 % a1 + a1) % a1);
    else puts("-1");
}