中国剩余定理

概述

  中国剩余定理(CRT)，又称孙子定理，其典型例题为著名的“韩信点兵”。

求解类型

中国剩余定理所能求解的问题如下：

给定同余方程组：
$$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\\ x\equiv a_3(mod{m}_3)\\ ⋮\\ x\equiv a_n(mod{m}_n)\end{array}$$
求解 $x$。

  在这$n$个式子中，所有的$m_i(i\in [1,n])$都关于其他的$m_j(j\in [1,n])$互质。
  对于该式，普通算法是难以快速求解的，所以一般会使用CRT或exCRT求解。

推导过程

  CRT的思路其实是利用构造的方法求出一个特解 $x_0$ ，再由特解推出其他解。

  首先，对于一个同余式$x\equiv a_i(mod{m}_i)$而言，它能用扩展欧几里得算法求解得到一组特解 $x_0$ ，且通解式为 $x=x_0+k\times m$ 其中 $m$ 为模数。所以当我们算出第一个同余式的解时，我们就可以利用通解式算出满足其他同余式的解。

  可以想到，假如设$$M=m_1\times m_2\times m_3\times \cdots \times m_n$$$$M_i=M÷m_i$$

  那么，对于任意的 $M_i$ 而言，它一定只模 $m_i$ 不为0，而对于其它的 $m_j(j\ne i)$ 而言，其模数一定为0。所以，我们可以求常出 $M_i$ 关于 $m_i$ 的逆元 $in{v}_i$ 。这样的话，对任意一组同余式而言一定存在一组特解 $x_0$ 可以表示为 $x_0=a_i\times M_i\times in{v}_i$ 。

  同时，由于 $M_i\mathrm{m}\mathrm{o}\mathrm{d}{m}_{i-1}=0$ ，所以 $a_i\times M_i\times t_i$ 一定是 $m_{i-1}$ 的倍数，所以如果 $x_{i-1}$ 是第 $i-1$ 个同余式的解，那么 $x_{i-1}+a_i\times M_i\times t_i$一定是第$i-1$个同余式的解。

  由于$M_i=M_{i+1}\times m_{i+1}÷m_i$，所以$m_{i+1}$是$M_i$的倍数，如果从第一个同余式推出来的话，那么$x_{i-1}+a_i\times M_i\times t_i$一定是第$i$个式子的解。所以我们可以推导出公式：
$$x=\sum _{i=1}^n{a}_i{M}_i{t}_i$$
  $x$即为同余式解。

  对于此代码，我们可以用扩展欧几里得算法求逆元，这样的话实际代码难度不大。

code

int A[Size] , m[Size] , M[Size] , M1 , x , y , n , ans;
void exgcd(int n , int m)//扩展欧几里得
{
	if(m == 0)
	{
		x = 1;
		y = 0;
		return ;
	}
	exgcd(m , n % m);
	int z = x;
	x = y;
	y = z - y * (n / m);
	return ;
}

int CRT()
{
	for(int i = 0 ; i < n ; i ++)
	{
		M1 *= m[i];
	}
	for(int i = 0 ; i < n ; i ++)
	{
		M[i] = M1 / m[i];
		exgcd(M[i] , m[i]);
		ans += A[i] * x * M[i] % M1;
		ans %= M1;
	}
	return (ans + M1) % M1;//求最小非负整数解
}

通解

  我们可以发现，因为$M={\prod }_{i=1}^n{m}_i$。所以，当且仅当特解$x_0$乘上$M$的倍数时，对同余式没有影响。由此可得通解式$x=x_0+k\times M(k\in Z)$