中国剩余定理

求解一元线性同余方程组：

这里要求mi两两互质。
令 M = m1 * m2 * m3 *···* mk，中国剩余定理说明：假设整数m₁,m₂, … ,m_n两两互质，则对任意的整数：a₁,a₂, … ,a_n，方程组 (S) 有解，并且通解可以用如下方式构造得到：
设 M = m₁ * m₂ * ···* m_n 是整数m₁,m₂, … ,m_n的乘积，并设 Mi = M/m_i， 是除了m_i以外的n - 1个整数的乘积。
设 ti = M_i^-1 为 Mi 模 m_i 的数论倒数( t_i 为 M_i 模 mi 意义下的逆元) M_it_i
方程组 (S) 的通解形式为x = a₁t₁M₁ + a₂t₂M₂ + ··· + a_nt_nM_n + kM = kM + ∑a_it_iM_i
在模 M 的意义下，方程组 (S) 只有一个解：x = (∑a_it_iM_i) mod M

void Exgcd(int A,int B,int &x,int &y){
    if(B == 0){
        x = 1;
        y = 0;
    }
    else {
        Exgcd(B,A%B,y,x);
        y -= (A/B)*x;
    }
}

int CRT(int a[],int m[],int n){
    int M = 1,Mi,ans = 0;
    int x,y;
    for(int i=0;i<n;i++)
        M *= m[i];  //mi两两互质，即求的lcm

    for(int i=0;i<n;i++){
        Mi = M / m[i];
        Exgcd(Mi,m[i],x,y);   //Exgcd求逆元
        ans = (ans + a[i] * Mi * x) % M;
    }
    return (ans+M) % M;
}

---

mi不互质时，就是将各个式子两两合并，再利用拓展欧几里得求 x 的过程，这里是在网上找到的图片——推导过程。

// 不互质的中国剩余定理 http://poj.org/problem?id=2891
// #include <bits/stdc++.h>
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int N = 1e5+10;
typedef long long ll;

void Exgcd(ll a,ll b,ll &gcd,ll &x,ll &y){
	if(b == 0){
		x = 1;
		y = 0;
		gcd = a;
	}
	else {
		Exgcd(b,a%b,gcd,y,x);
		y -= (a/b)*x;
	}
}

int main(){
	ll n,a1,a2,r1,r2;
	while(~scanf("%lld",&n)){
		scanf("%lld%lld",&a1,&r1);
		int flag = 1;
		for(int i=1;i<n;i++){
			scanf("%lld%lld",&a2,&r2);
			ll x,y,gcd;
			Exgcd(a1,a2,gcd,x,y);
			ll c = r2 - r1;
			if(c % gcd != 0){
				flag = 0;
				continue;
			}

			x = x * (c/gcd);
			ll t = a2/gcd;
			if(t < 0)	t = -t;
			x = (x%t + t)%t;
			// ll X = a1 * x + r1;
			// ll Lcm = a1*a2/gcd;
			// a1 = Lcm;
			// r1 = X % Lcm;

	        r1 += x*a1;
	        a1 *= a2/gcd;
	        r1 = (r1%a1+a1)%a1;
		}
		printf("%lld\n",flag ? r1 : -1);
	}

    return 0;
}

学习博客：
https://blog.csdn.net/qq_29980371/article/details/71053219
https://www.cnblogs.com/MashiroSky/p/5918158.html