中国剩余定理笔记

1.定义

孙子定理是中国古代求解一次同余式组（见同余）的方法。是数论中一个重要定理。又称中国余数定理。一元线性同余方程组问题最早可见于中国南北朝时期（公元5世纪）的数学著作《孙子算经》卷下第二十六题，叫做“物不知数”问题，原文如下：

有物不知其数，三三数之剩二，五五数之剩三，七七数之剩二。问物几何？即，一个整数除以三余二，除以五余三，除以七余二，求这个整数。《孙子算经》中首次提到了同余方程组问题，以及以上具体问题的解法，因此在中文数学文献中也会将中国剩余定理称为孙子定理。

翻译过来就是有一个数除以三余二，除以五余三，除以七余二，问这个数是多少。

《孙子算经》解这道题目的“术文”和答案是：
“三三数之剩二，置一百四十；
五五数之剩三，置六十三；
七七数之剩二，置三十。
并之，得二百三十三，以二百十减之，即得。”“答曰：二十三。”
这些话是什么意思呢？用通俗的话来说，就是：
先求被3除余2，并能同时被5、7整除的数，这样的数最小是35；
再求被5除余3，并能同时被3、7整除的数，这样的数最小是63； 然后求被7除余2，并能同时被3、5整除的数，这样的数最小是30。
于是，由35+63+30=128，得到的128就是一个所要求得的数。但这个数并不是最小的。
再用求得的“128”减去或者加上3、5、7的最小公倍数“105”的倍数，就得到许许多多这样的数：
{23，128，233，338，443，…}
从而可知，23、128、233、338、443、…都是这一道题目的解，而其中最小的解是23。

2.证明过程

首先知道的就是（a*k）%b=a%b+a%b+…+a%b=c+c+…+c=kc（k>0）,也就是说，如果一个除法的余数为c，那么被除数的k倍与除数相除的余数为kc。

我们应该先知道两个知识点：

1. 如果 a%b=c , 则有 (a+kb)%b=c (k为非零整数)。
2. 如果 a%b=c，那么 (a*k)%b=kc (k为大于零的整数)。证明在上面

因此它的解法总结一下就是，下面并不是直接求取模等于需要的数字，而是先算出来取模3， 5， 7等于1的数字，然后在乘以需要的数字。

1. 求出三和五的一个公倍数并且除以7得1，设为n₁。
2. 求出三和七的一个公倍数并且除以5得1，设为n₂。
3. 求出五和七的一个公倍数并且除以3得1，设为n₃。
4. 将n₁乘以2，n₂乘以3，n₃乘以2，在相加得到n₄。由上面的两个定理我们知道n₄一定是满足条件的值。
5. 这个时候取得的n₄，不一定是最小的解，我们还要模3， 5，7的最小公倍数105。

3.数学证明

用现代数学的语言来说明的话，中国剩余定理给出了以下的一元线性同余方程组：

有解的判定条件，并用构造法给出了在有解情况下解的具体形式。
中国剩余定理说明：假设整数m₁,m₂, … ,m_n两两互质，则对任意的整数：a₁,a₂, … ,a_n，方程组有解，并且通解可以用如下方式构造得到：

因此我们只需要计算出M_i和逆元t_i即可。逆元求法参考我这篇博客：https://blog.csdn.net/qq_43227036/article/details/99406800

对于数学的方法，和咱们上面的过程其实是想对应的，所有数乘积为M，去除掉第i个后得到M_i在求得其与a_i的逆元t_i，保证了t_i*M_i对a_i取余后的1，然后在乘以a_i，就得到的是取余后得a_i的数。将所有的数加起来，最后对M取余

4.实现代码

void exgcd(int a1,int b,int &x,int &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return ;
    }
    exgcd(b,a1%b,x,y);
    int t=x;
    x=y;
    y=t-(a1/b)*y;
}
int CRT(int a[],int m[],int n)
{
    int M=1,ans=0,t,x,y;
    for(int i=0; i<n; i++)
    {
        M*=m[i];///M为除数乘积
    }
    for(int i=0; i<n; i++)
    {
        t=M/m[i];///除了mi以外的n-1个整数乘积
        exgcd(t,m[i],x,y);///求逆元，由扩展欧几里得转换成t*ti+m[i]*y=1来求ti
        ans=(ans+a[i]*x*t)%M;
    }
    //防止结果为负值，为了防止爆int，应该判断一下，而不是直接写为下面的代码
    return (ans+M)%M;
}