[数论]-----中国剩余定理（扩展中国剩余定理）

中国剩余定理

中国剩余定理（CRT）用于求形如：
$$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\\ \cdots \cdots \\ x\equiv a_k(mod{m}_k)\end{array}$$
这样一系列同余方程组的最小非负整数解x。

情况1：模数两两互质

模数$m_1,m_2...m_k$ 两两互质的情况，相对来说比较简单。
我们令：

1. $M$表示所有方程组的模的乘积。 $M={\prod }_{i=1}^k{m}_i$
2. $M_i$ 表示除第 i 个方程外，其余方程组模的乘积。 $M_i=\frac{M}{m_i}$
3. $t_i$ 为 $M_i$模$m_i$意义下的逆元。 $M_i{t}_i\equiv 1(mod{m}_i)$

其中$1\le i\le k$。

从这里我们就可以知道，为什么这种情况需要模数两两互质。因为$M_i=\frac{M}{m_i}$，$M$中包含$m_i$这个因子，当模数两两互质时，$M_i$也和$m_i$互质，这样$M_i$才有$m_i$意义下的逆元。

然后我们可以构造出一个解：$x={\sum }_{i=1}^k{a}_i{M}_i{t}_i$
由此，任意解$x_0$即为$x+k\ast M$，最小正整数解$x_{min}=x_0\%M$

代码：（洛谷P1495）

#include <cstdio>
#define ll long long int
using namespace std;
int N;
ll a[15], b[15], Mi[15], X, M = 1;
ll exgcd(ll a, ll b, ll &x, ll &y)
{
	if(b==0)
	{
		x = 1; y = 0;
		return a;
	}
	ll r = exgcd(b, a%b, x, y);
	ll temp = y;
	y = x-(a/b)*y;
	x = temp;
	return r;
}
int main()
{
	scanf("%d", &N);
	for(int i=1;i<=N;i++)
	{
		scanf("%d%d", &a[i], &b[i]);
		M *= a[i];
	}
	for(int i=1;i<=N;i++)
	{
		Mi[i] = M/a[i];
		ll x = 0, y = 0;
		exgcd(Mi[i], a[i], x, y);         //exgcd求逆元
		X += b[i]*Mi[i]*(x<0 ? x+a[i] : x);
	}
	printf("%lld\n", X%M);
	return 0;
}

扩展中国剩余定理

ex-CRT，其实就是考虑上例中，模数两两不一定互质的情况，不知道谁给它取了一个扩展中国剩余定理的名字。模数不互质，就不能用借助逆元的方法了，需要另一种（小学数学）方法。

假设已经求出了前 $k-1$ 个方程的解 $x_{k-1}$，设$M=LCM(a_i)$，即 $M$ 为前 $k-1$ 个除数的最小公倍数，则：
对于前 $k-1$ 个方程，都满足 $x_{k-1}+tM\equiv a_i(mod{b}_i)$ $(t\in Z)$
即：前 $k-1$ 个方程，通解为$x_{k-1}+tM$ $(t\in Z)$。

想求得第 $k$ 个方程的解，并且同时满足前$k-1$ 个方程，那么需要使第 $k$ 个方程的解，为前 $k-1$ 的方程的通解的同时，也满足第 $k$ 个方程。

设第$k$个方程的解 $x_k=x_{k-1}+tM$ $(t\in Z)$
将此解代入第 $k$ 个方程中，可得：
$x_{k-1}+tM\equiv a_k(mod{b}_k)$
即：$tM\equiv a_k-x_{k-1}(mod{b}_k)$
其中：$M,a_k,x_{k-1},b_k$ 都是已知的。

使用exgcd解出这个同余方程，得到 t 的值，然后反代回 $x_k=x_{k-1}+tM(t\in Z)$，就可以得到 $x_k$ 的值。进行 k 次上述操作，可以得到全部 k 个方程的解。

例题

模数互质模板：洛谷P1495 【模板】中国剩余定理(CRT).

模数不互质模板：洛谷P4777 【模板】扩展中国剩余定理（EXCRT）.