204. 表达整数的奇怪方式

输入样例：
2
8 7
11 9

输出样例：
31

---

想法：
   核心：让程序去执行合并的操作
    先用扩展欧几里得算法求得k，然后代入到原式内，此时扩展欧几里得算法求得的结果与原式成整数倍，所以求得的解需要 * (m2-m1)/d

x (mod a₁) = m₁ (mod a₁)
x (mod a₂) = m₂ (mod a₂)
x (mod a₃) = m₃ (mod a₃)
·
·
·
x (mod a_n) = m_n (mod a_n)

将等式两两合并成一个式子，直到最后只剩下一个式子
合并过程：
x mod a₁ = m₁           =>           x = a₁k₁ + m₁ ①
x mod a₂ = m₂          =>            x = a₂k₂ + m₂ ②
合并消去x ： k₁ * a₁ - k₂ * a₂ = m₂ - m₁ ③
此时将 k₁ k₂当作x,y ;     a₁ a₂ 当作a,b;
转换为式子 ： ax+by = gcd(a,b);
有解条件为： m₂-m₁ = gcd(a,b) 的倍数 d= gcd(a,b)
解出k:
         k₁= k₁+k(a₂/d) = k₁%（a₂/d）
         k₂= k₂+k(a₁/d) = k₂%（a₁/d）
带入③式恒成立
将 k1= k1+k(a2/d) 代入①式
x = a₁k₁ + m₁ = (k₁+k(a₂/d))*a₁ + m₁
化简后：
x = a₁k₁ + m₁ +k (a₂ * a₁ /d) ④
目标： 求出m和a 获取到新的等式 x = m + ka
再化简 （ 等式右侧的值均已知 ）
将 m= a₁k₁ + m₁ ； a=(a₂ * a₁/d)
x = m + ka ⑤
再将⑤与下一个方程合并求解

最后 x = m+ka;
因为 m = a_nk_n + m_n 和 a = | a₂ * a₁/d |
此时 x = m+ka; 减去一个数的整数倍等于mod上这个数
所以结果为 x = (m%a + a)%a 防止是负数，所以加上a再 %a
求解得出x 的值

代码

#include <cstdio>
#include <iostream>
using namespace std;
typedef long long LL;
int n;
LL exgcd(LL a, LL b, LL &x, LL &y){
    if(b == 0){
        x = 1, y = 0;
        return a;
    }

    LL d = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}
LL inline mod(LL a, LL b){
    return ((a % b) + b) % b;
}
int main(){
    scanf("%d", &n);
    LL a1, m1;
    scanf("%lld%lld", &a1, &m1);
    for(int i = 1; i < n; i++){
        LL a2, m2, k1, k2;
        scanf("%lld%lld", &a2, &m2);
        LL d = exgcd(a1, -a2, k1, k2);
        if((m2 - m1) % d){ puts("-1"); return 0; }
        k1 = mod(k1 * (m2 - m1) / d, abs(a2 / d));
        m1 = k1 * a1 + m1;
        a1 = abs(a1 / d * a2);
    }
    printf("%lld\n", m1);
    return 0;
}