【数学】中国剩余定理

给定 $n$ 个正整数 $a_i$ 和 $n$ 个非负整数 $b_i$，求解关于 $x$ 的线性同余方程组：
$$\{\begin{array}{l}x\equiv b_1(\mathrm{m}\mathrm{o}\mathrm{d}{a}_1)\\ x\equiv b_2(\mathrm{m}\mathrm{o}\mathrm{d}{a}_2)\\ \cdots \\ x\equiv b_n(\mathrm{m}\mathrm{o}\mathrm{d}{a}_n)\end{array}$$
其中 $\forall i,j\in [1,n],i\ne j$ 满足 $a_i,a_j$ 互质。

首先令 $m=\prod _{i=1}^n{a}_i$。

如果有一个解 $x_1$，那么通解就是 $x=x_1+tm$（$t$ 为整数）。

所以我们只需要求出最小解即可，那么这个最小解一定 $<m$。

中国剩余定理（$\mathrm{C}\mathrm{h}\mathrm{i}\mathrm{n}\mathrm{e}\mathrm{s}\mathrm{e}\mathrm{R}\mathrm{e}\mathrm{m}\mathrm{a}\mathrm{i}\mathrm{n}\mathrm{d}\mathrm{e}\mathrm{r}\mathrm{T}\mathrm{h}\mathrm{e}\mathrm{o}\mathrm{r}\mathrm{e}\mathrm{m},\mathrm{C}\mathrm{R}\mathrm{T}$）给出了构造解的方法：

$x=(\sum _{i=1}^n{b}_i{m}_i{M}_i)\mathrm{m}\mathrm{o}\mathrm{d}m$

其中 $m_i=\frac{m}{a_i}$，$M_i$ 为 $m_i$ 在模 $a_i$ 意义下的逆元。

则 $\forall i,j\in [1,n],i\ne j$：
$$\begin{array}{rl}{b}_j{m}_j{M}_j& \equiv b_j\frac{m}{a_j}{M}_j\\ & \equiv b_j\frac{a_1{a}_2\cdots a_i\cdots a_n}{a_j}{M}_j\\ & \equiv b_j\cdot 0\cdot M_j\\ & \equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}{a}_i)\end{array}$$
而 $\forall i\in [1,n]$：
$$\begin{array}{rl}{b}_i{m}_i{M}_i& \equiv b_i(m_i{M}_i)\\ & \equiv b_I(m_i{m}_i^{-1})\\ & \equiv b_i\cdot 1\\ & \equiv b_i(\mathrm{m}\mathrm{o}\mathrm{d}{a}_i)\end{array}$$
所以 $\forall i\in [1,n]$：
$$\begin{array}{rl}\sum _{i=1}^n{b}_i{m}_i{M}_i& \equiv \sum _{j=1,j\ne i}^n{b}_j{m}_j{M}_j+b_i{m}_i{M}_i\\ & \equiv 0+b_i\\ & \equiv b_i(\mathrm{m}\mathrm{o}\mathrm{d}{a}_i)\end{array}$$
这就是一组满足条件的解了。

接下来证明在模 $m$ 的意义下只有一个解。

假设有 $x\equiv y(\mathrm{m}\mathrm{o}\mathrm{d}m)$，那么 $\forall i\in [1,n]$，有

$$\begin{array}{rl}x& \equiv x\mathrm{m}\mathrm{o}\mathrm{d}m\\ & \equiv y(\mathrm{m}\mathrm{o}\mathrm{d}{a}_i)\end{array}$$

所以 $x$ 和 $y$ 是同一个解。

时间复杂度为 $\mathcal{O}(n\log n)$。

P1495 【模板】中国剩余定理(CRT)/曹冲养猪

$\text{Code}$

//18 = 9 + 9 = 18.
#include <iostream>
#include <cstdio>
typedef long long ll;
using namespace std;

const int MAXN = 15;

ll x, y;

void exgcd(ll a, ll b)
{
	if (!b)
	{
		x = 1, y = 0;
		return;
	}
	exgcd(b, a % b);
	ll tmp = x;
	x = y;
	y = tmp - a / b * y;
}

ll inv(ll a, ll p)
{
	exgcd(a, p);
	x = (x % p + p) % p;
	return x;
}

ll qmul(ll a, ll b, ll p)
{
	ll base = a, ans = 0;
	while (b)
	{
		if (b & 1)
		{
			ans = (ans + base) % p;
		}
		base = (base + base) % p;
		b >>= 1;
	}
	return ans;
}

int a[MAXN], b[MAXN];

ll CRT(int n)
{
	ll m = 1;
	for (int i = 1; i <= n; i++)
	{
		m *= a[i];
	}
	ll ans = 0;
	for (int i = 1; i <= n; i++)
	{
		ll mi = m / a[i];
		ll Mi = inv(mi, a[i]);
		ans = (ans + qmul(b[i], qmul(mi, Mi, m), m)) % m;
	}
	return ans;
}

int main()
{
	int n;
	scanf("%d", &n);
	for (int i = 1; i <= n; i++)
	{
		scanf("%d%d", a + i, b + i);
	}
	printf("%lld\n", CRT(n));
	return 0;
}

双倍经验：P3868 [TJOI2009] 猜数字