中国剩余定理【数论】

使用中国剩余定理来求解上面的“物不知数”问题，便可以理解《孙子歌诀》中的数字含义。这里的线性同余方程组是：

三个模数m₁3, m₂5, m₃7的乘积是M105，对应的M₁35, M₂21, M₃15. 而可以计算出相应的数论倒数：t₁2, t₂1, t₃1. 所以《孙子歌诀》中的70，21和15其实是这个“物不知数”问题的基础解：

而将原方程组中的余数相应地乘到这三个基础解上，再加起来，其和就是原方程组的解：

这个和是233，实际上原方程组的通解公式为：

《孙子算经》中实际上给出了最小正整数解，也就是k-2时的解：x23.

///n个mi互质
const LL maxn = 20;
LL a[maxn], m[maxn], n;
LL CRT(LL a[], LL m[], LL n)
{
    LL M = 1;
    for (int i = 0; i < n; i++) M *= m[i];
    LL ret = 0;
    for (int i = 0; i < n; i++)
    {
        LL x, y;
        LL tm = M / m[i];
        ex_gcd(tm, m[i], x, y);
        ret = (ret + tm * x * a[i]) % M;
    }
    return (ret + M) % M;
}

模不两两互质的同余式组可化为模两两互质的同余式组，再用孙子定理直接求解。

84=2²×3×7,160=2⁵×5,63=3²×7，由推广的孙子定理可得  与  同解。

///n个mi不互质
const LL maxn = 1000;
LL a[maxn], m[maxn], n;
LL CRT(LL a[], LL m[], LL n) {
    if (n == 1) {
        if (m[0] > a[0]) return a[0];
        else return -1;
    }
    LL x, y, d;
    for (int i = 1; i < n; i++) {
        if (m[i] <= a[i]) return -1;
        d = ex_gcd(m[0], m[i], x, y);
        if ((a[i] - a[0]) % d != 0) return -1;   //不能整除则无解
        LL t = m[i] / d;
        x = ((a[i] - a[0]) / d * x % t + t) % t; //第0个与第i个模线性方程的特解
        a[0] = x * m[0] + a[0];
        m[0] = m[0] * m[i] / d;
        a[0] = (a[0] % m[0] + m[0]) % m[0];
    }
    return a[0];
}