中国剩余定理（孙子定理）

中国剩余定理

\quad 在《孙子算经》中有这样一个问题：“今有物不知其数，三三数之剩二（除以$3$余$2$），五五数之剩三（除以$5$余$3$），七七数之 剩二（除以$7$余$2$），问物几何？”这个问题称为“孙子问题”，该问题的一般解法国际上称为“中国剩余定理”。

具体解法分三步：

1. 找出三个数：从$3$和$5$的公倍数中找出被$7$除余$1$的最小数$15$，从$3$和$7$的公倍数中找出被$5$除余$1$的最小数$21$，最后从$5$和$7$的公倍数中找出除$3$余$1$的最小数$70$。
2. 用$15$乘以$2$（$2$为最终结果除以$7$的余数），用$21$乘以$3$（$3$为最终结果除以$5$的余数），同理，用$70$乘以$2$（$2$为最终结果除以$3$的余数），然后把三个乘积相加$15\ast 2+21\ast 3+70\ast 2$得到和$233$。
3. 用$233$除以$3，5，7$三个数的最小公倍数$105$，得到余数$23$，即$233$。这个余数$23$就是符合条件的最小数。

\quad 就这么简单。我们在感叹神奇的同时不禁想知道古人是如何想到这个方法的，有什么基本的数学依据吗？

\quad 我们将“孙子问题”拆分成几个简单的小问题，从零开始，试图揣测古人是如何推导出这个解法的。

\quad 首先，我们假设$n_1$是满足除以$3$余$2$的一个数，比如$2，5，8$等等，也就是满足$3\ast k+2（k>=0)$的一个任意数。同样，我们假设$n_2$是满足除以$5$余$3$的一个数，$n_3$是满足除以$7$余$2$的一个数。

\quad 有了前面的假设，我们先从$n_1$这个角度出发，已知$n_1$满足除以$3$余$2$，能不能使得$n_1+n_2$的和仍然满足除以$3$余$2$？进而使得$n_1+n_2+n_3$的和仍然满足除以$3$余$2$？

\quad 这就牵涉到一个最基本数学定理，如果有$a\%b=c$,则有$(a+k\ast b)\%b=c$($k$为非零整数)，换句话说，如果一个除法运算的余数为$c$，那么被除数与$k$倍的除数相加（或相减）的和（差）再与除数相除，余数不变。这个是很好证明的。

\quad 以此定理为依据，如果$n_2$是$3$的倍数，$n_1+n_2$就依然满足除以$3$余$2$。同理，如果$n_3$也是$3$的倍数，那么$n_1+n_2+n_3$的和就满足除以$3$余$2$。这是从$n_1$的角度考虑的，再从$n_2,n_3$的角度出发，我们可推导出以下三点：

• 为使$n_1+n_2+n_3$的和满足除以$3$余$2$，$n_2$和$n_3$必须是$3$的倍数。
• 为使$n_1+n_2+n_3$的和满足除以$5$余$3$，$n_1$和$n_3$必须是$5$的倍数。
• 为使$n_1+n_2+n_3$的和满足除以$7$余$2$，$n_1$和$n_2$必须是$7$的倍数。

\quad 因此，为使$n_1+n_2+n_3$的和作为“孙子问题”的一个最终解，需满足：

1. $n_1$除以$3$余$2$，且是$5$和$7$的公倍数。
2. $n_2$除以$5$余$3$，且是$3$和$5$的公倍数。
3. $n_3$除以$7$余$2$，且是$3$和$5$的公倍数。

\quad 所以，孙子问题解法的本质是从$5$和$7$的公倍数中找一个除以$3$余$2$的数$n_1$，从$3$和$7$的公倍数中找一个除以$5$余$3$的数$n_2$，从$3$和$5$的公倍数中找一个除以$7$余$2$的数$n_3$，再将三个数相加得到解。在求$n_1,n_2,n_3$时又用了一个小技巧，以$n_1$为例，并非从$5$和$7$的公倍数中直接找一个除以$3$余$2$的数，而是先找一个除以$3$余$1$的数，再乘以$2$。也就是先求出$5$和$7$的公倍数模$3$下的逆元，再用逆元去乘余数。

\quad 最后，我们还要清楚一点，$n_1+n_2+n_3$只是问题的一个解，并不是最小的解。如何得到最小解？我们只需要从中最大限度的减掉掉$3，5，7$的公倍数$105$即可。道理就是前面讲过的定理“如果$a\%b=c$，则有$(a-k\ast b)\%b=c”$。所以$（n_1+n_2+n_3）\%105$就是最终的最小解。

这样就得到了中国剩余定理的公式：

设正整数$m_1,m_2,,,,,m_k$两两互素，则同余方程组

那么再看提到的一元线性同余方程：(先做几个设定：设$N_i$为能够被$m_1,m_2,......,m_{i-1},m_{i+1},.....m_k$整除，但是除以$m_i$正好余1)

$X=N_1\ast a_1+N_2\ast a_2+.......+N_n\ast a_n$

就是我们要求的一个解，解集为$X\%LCM(m_1,m_2,....,m_n)$。

剩下的问题就变成了如何求解$ N_{1}, N_{2} … N_{n}$，我们继续向下看：

假设$m=LCM(m_1,m_2,....,m_n)$ $x^{\prime },y^{\prime }$为任意整数。

因为$N_i$的性质，$N_i$可以表示为：$N_i==m/m_i\ast x^{\prime }==m_i\ast y^{\prime }+1==>m/m_i\ast x^{\prime }+(-m_i)\ast y^{\prime }==1.$

推到现在有没有感觉很熟悉，对的！这个就是扩展欧几里德：

对于$gcd(a,b)=d$, 存在$a\ast x+b\ast y==gcd(a,b)$;

存在$gcd(-m_i,m/m_i)==1$,

对于$m/m_i\ast x^{\prime }+(-m_i)\ast y^{\prime }==1,$

套用一下扩展欧几里德求出$x^{\prime },$ 就可以求解出$N_i$。有了$N_i$求出$X$就是分分钟的事情！

Code：

//中国剩余定理模板
ll china(int a[],int b[],int n)//a[]为除数，b[]为余数
{
    int M=1,y,x=0;
    for(int i=0; i<n; ++i) //算出它们累乘的结果
        M*=a[i];
    for(int i=0; i<n; ++i)
    {
        int w=M/a[i];
        int tx=0;
        int t=exgcd(w,a[i],tx,y);//计算逆元
        x=(x+w*(b[i]/t)*x)%M;
    }
    return (x+M)%M;
}