中国剩余定理以及扩展中国剩余定理

中国剩余定理

设$m1,m2，\cdots ,m_n$是两两互质的整数，对于下式的一元线性同余方程组有整数解
$$\{\begin{array}{l}x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_1)\\ x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{m}_2)\\ ⋮\\ x\equiv a_n(\mathrm{m}\mathrm{o}\mathrm{d}{m}_n)\end{array}$$
解为：
$$x=\sum _{i=1}^n{a}_i\ast M_i\ast M_i^{-1}$$
其中：$m={\prod }_{i=1}^n{m}_i，M_i=\frac{m}{m_i}$，$M_i^{-1}$为$M_i$的逆元，可以使用扩展欧几里得求解

具体证明可参考百度百科：中国剩余定理

代码模板：

LL CRT(LL a[], LL m[], int n)
{
    LL M = 1, ans = 0;
    for(int i = 1; i <= n; ++i) M *= m[i];
    for(int i = 1; i <= n; ++i)
    {
        LL x, y;
        LL Mi = M / m[i];
        exgcd(Mi, m[i], x, y);
        x = ((x % m[i]) + m[i] ) % m[i];
        ans = (ans  + a[i] * Mi * x ) ;
    }
    return ans;
}

扩展中国剩余定理

原题链接（模板题）
同样的，给出一个一元线性同余方程组：
$$\{\begin{array}{l}x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_1)\\ x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{m}_2)\\ ⋮\\ x\equiv a_n(\mathrm{m}\mathrm{o}\mathrm{d}{m}_n)\end{array}$$
不保证$m_1,m_2,\cdots ,m_n$两两互质。

先考虑前两个式子拿出来：

$$\begin{gathered} x=k_1\ast a_1+m_1① \\ x=k_2\ast a_2+m_2 \end{gathered}$$
两式相减有：
$$k_1\ast a_1+k_2\ast (-a_2)=m_2-m_1②$$
这个式子两种情况：

• 设$d=gcd(a_1,-a_2)$，如果$d\nmid (m_2-m_1)$，则无解
• 否则可以用扩展欧几里得求出一组解${k_1}^{\prime }$和${k_2}^{\prime }$，设$d=gcd(a_1,-a_2),y=\frac{m_2-m_1}{d}$
  我们将${k_1}^{\prime }$和${k_2}^{\prime }$都扩大$y$倍后，能满足$②$式

接下来就是需要找到最小的解了。

我们有性质：
$$\begin{gathered} k_1=k_1+k\ast \frac{a_2}{d}③ \\ {k}_2=k_2+k\ast \frac{a_1}{d} \end{gathered}$$
其中：$k\in Z$，此时$k_1,k_2$仍然满足$②$式
证明如下：
我们将新的$k_1,k_2$带入②式有：
$$\begin{gathered} (k_1+k\ast \frac{a_2}{d})\ast a_1+(k_2+k\ast \frac{a_1}{d})\ast (-a_2) \\ =k_1\ast a_1+k\ast \frac{a_1\ast a_2}{d}+k_2\ast (-a_2)+k\ast \frac{a_1\ast (-a_2)}{d} \\ =k_1\ast a_1+k_2\ast (-a_2) \end{gathered}$$
证毕
我们需要一个最小的非负这整数解，我们对$③$式取模$abs(\frac{a_2}{d})$的到：$$k_1=k_1\mathrm{m}\mathrm{o}\mathrm{d}abs(\frac{a_2}{d})$$

• 取绝对值，是为了因为不知道$\frac{a_2}{d}$的正负性，同时使得$k_1$尽可能的小

最后我们把得到的$k_1$带回①式有：
$$\begin{gathered} x=(k_1+k\ast \frac{a_2}{d})\ast a_1+m_1 \\ x=k_1\ast a_1+m_1+k\ast \frac{a_1\ast a_2}{d} \\ {k}_i\ast a_i{m}_i \\ x=\overbrace{k\ast \frac{a_1\ast a_2}{d}}+\overbrace{k_1\ast a_1+m_1} \end{gathered}$$
这样：我们就把前两个式子合并为一个式子了，这样子取合并下面的式子。经过$n-1$次合并，就将$n$个式子合并为一个式子了。

扩展中国剩余定理模板：

LL exgcd(LL a, LL b, LL &x, LL &y)
{
    if(!b) {
        x = 1, y = 0;
        return a;
    }
    LL x1, y1, Gcd;
    Gcd = exgcd(b, a % b, x1, y1);
    x = y1;
    y = x1 - a / b * y1;
    return Gcd;
}
LL inline Mod(LL a, LL b)
{
    return ((a % b) + b) % b;
}
LL EXCRT(LL a[], LL m[], int n)
{
    LL a1 = a[1], m1 = m[1];
    for(int i = 2; i <= n; ++i)
    {
        LL a2 = a[i], m2 = m[i], k1, k2;
        LL d = exgcd(a1, -a2, k1, k2);
        if((m2 - m1) % d)
            return -1;
        LL y = (m2 - m1) / d;
        k1 = Mod(k1 * (m2 - m1) / d, abs(a2 / d));
        m1 = k1 * a1 + m1;
        a1 = abs(a2 / d * a1);
    }
    return m1;
}

另外：

AC代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 50;

LL a[N], m[N];

LL exgcd(LL a, LL b, LL &x, LL &y)
{
    if(!b) {
        x = 1, y = 0;
        return a;
    }
    LL x1, y1, Gcd;
    Gcd = exgcd(b, a % b, x1, y1);
    x = y1;
    y = x1 - a / b * y1;
    return Gcd;
}
LL inline Mod(LL a, LL b)
{
    return ((a % b) + b) % b;
}
LL EXCRT(LL a[], LL m[], int n)
{
    LL a1 = a[1], m1 = m[1];
    for(int i = 2; i <= n; ++i)
    {
        LL a2 = a[i], m2 = m[i], k1, k2;
        LL d = exgcd(a1, -a2, k1, k2);
        if((m2 - m1) % d)
            return -1;
        LL y = (m2 - m1) / d;
        k1 = Mod(k1 * (m2 - m1) / d, abs(a2 / d));
        m1 = k1 * a1 + m1;
        a1 = abs(a2 / d * a1);
    }
    return m1;
}

int main()
{
    int n;
    cin >> n;
    for(int i = 1; i <= n; ++i) cin >> a[i] >> m[i];
    
    cout << EXCRT(a, m, n) << endl;
}