中国剩余定理学习 && 拓展中国剩余定理

中国剩余定理：

仅供自己复习时查看一下大佬笔记，详细学习过程在大佬的博客。

学习连接：https://www.cnblogs.com/MashiroSky/p/5918158.html

拓展中国剩余定理：

中国剩余定理在求解同余方程组的时候，要求模数之间是互质的，那么假如模数没有说明一定是互质的，这时就要用到了我们的拓展中国剩余定理($excrt$)，拓展中国剩余定理其实和中国剩余定理联系不大，我觉得是因为它是用$exgcd$退出，并且与其联系紧密所以叫这个名字。
同余方程组：
$$\{\begin{array}{l}x\equiv r_1(mod{m}_1)\\ x\equiv r_2(mod{m}_2)\\ ..\\ x\equiv r_n(mod{m}_n)\end{array}$$
考虑最开始的两个式子$$\{\begin{array}{l}x\equiv r_1(mod{m}_1)\\ x\equiv r_2(mod{m}_2)\end{array}$$
即$x=r_1+k_1\ast m_1$，$x=r_2+k_2\ast m_2$，得$r_1+k_1\ast m_1=r_2+k_2\ast m_2$，移项的，$r_2-r_1=k_1\ast m_1+(-k_2)\ast m_2=k_1\ast m_1+k_2\ast m_2$，这样这个式子就变成类似$ax+by=c$的形式了，于是我们求
$K_1\ast m_1+K_2\ast m_2=gcd(m1,m2)$，得到$K_1$后两边同时乘上$\frac{(r2-r1)}{gcd(m_1,m_2)}$，可得$k_1=K_1\ast \frac{(r2-r1)}{gcd(m_1,m_2)}$，这样把$k_1$在转化为最小的正数解之后，就可以用得到一个特解
$x^{\ast }=k_1\ast m_1+r_1$，怎么通过这个特解去获得一个同余方程呢，观察前面的式子
$$\{\begin{array}{l}x=r_1+K_1\ast \frac{(r2-r1)}{gcd(m_1,m_2)}\ast m_1\\ x=r_2+K_2\ast \frac{(r2-r1)}{gcd(m_1,m_2)}\ast m_2\end{array}$$
可以看到任意两个解之间的差值一定是$\frac{m_1\ast m_2}{gcd(m1,m2)}$的倍数，也就是我们的特解 $x^{\ast }+t\ast lcm(m1,m2)=x(通解)$，到此也就能看出新的同余方程就是$x\equiv x^{\ast }(modlcm(m1,m2))$。
所以对于n组方程我们只需要不断两两进行合并，更新合并出的方程组的对应的值，合并n-1次就可以得到一个最终的同余方程，最后的答案即为我们求得。

一些具体实现细节，需要结合代码加深理解。

例题：POJ2891 模板题

#include <algorithm>
#include <cmath>
#include <cstring>
#include <iostream>

using namespace std;
#define int long long
#define pb emplace_back
#define MP make_pair
#define pii pair<int,int>
#define pll pair<ll,ll>
#define lson rt<<1
#define rson rt<<1|1
typedef long long ll;
typedef double db;
const int INF = 0x3f3f3f3f;
const db eps = 1e-6;
const int MAXN = 1e5 + 7;
int n;
int mod[MAXN],r[MAXN];//分别代表模数数组 和 余数数组
int exgcd(int a,int b,int &x,int &y) {
	if(b == 0) {
		x = 1; y = 0;
		return a;
	}
	else {
		int res = exgcd(b,a%b,x,y);
		ll t = x;
		x = y;
		y = t - a / b * y;
		return res;
	}
}
int excrt() {//就是一个不断合并并且更新新的参数的过程
	int MOD = mod[1],R = r[1],x,y;
	for(int i = 2;i <= n;i ++) {
		int g = exgcd(MOD,mod[i],x,y);
		if((r[i] - R) % g != 0) return -1;//同余方程组无解
		x = x * (r[i] - R) / g;
		x = (x % (mod[i] / g) + (mod[i] / g)) % (mod[i] / g);//得到最小的解
		R = x * MOD + R;
		MOD = MOD / g * mod[i];//合并之后的模数 应该是lcm
		R = R % MOD;//重新取一遍得到合并后的新的同余方程
	}
	return (R % MOD + MOD) % MOD;
}

signed main() {
	while(~scanf("%lld",&n)) {
		for(int i = 1;i <= n;i ++) {
			scanf("%lld%lld",&mod[i],&r[i]);
		}
		int ans = excrt();
		printf("%lld\n",ans);
	}
}