【数论】同余（四）：一元线性同余方程组（两两相消、中国剩余定理）

同余问题共7part，我的博客链接：

• 基本概念与性质
• 逆元：概念、求解方法与推导
• 线性同余方程
  • 一元线性同余方程
  • 一元线性同余方程组
  • 多元线性同余方程
• 高次同余方程：BSGS算法（大小步算法、拔山盖世算法）
• 快速幂与矩阵快速幂

一元线性同余方程组

形如：
$$\{\begin{array}{ll}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}$$

方法一：两两相消

先看两个方程：
$$\{\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)\end{array}$$
转化为
$$\{\begin{array}{l}x+m_1{k}_1=a_1\\ x+m_2{k}_2=a_2\end{array}\text{\hspace{1em}(1)}$$
则得到$$a_1-m_1{k}_1=a_2-m_2{k}_2$$ 即$$m_1{k}_1-m_2{k}_2=a_1-a_2$$ 则有解的等价条件为$$\gcd (m_1,m_2)\mid (a_1-a_2)$$ 用上一节的一元线性同余方程解法解出$k_1$ 的一个特解$k^{\prime }$ ，于是通解$$k_1=k^{\prime }+\frac{c{m}_2}{\gcd (m_1,m_2)},c\in N$$ 代入$(1)$ 式可得，$$x=a_1-m_1{k}^{\prime }+\frac{c{m}_1{m}_2}{\gcd (m_1,m_2)}$$ 于是转换为了新的同余方程：$$x\equiv a_1-m_1{k}^{\prime }(\mathrm{m}\mathrm{o}\mathrm{d}\text{lcm}(m_1,m_2))$$ 依次合并即可解得最终答案。

// 求解 x==a[] mod m[]，数组从0开始编号，共n个式子
int mod_equations(int a[], int m[], int n)
{
    int a1 = a[0], a2, m1 = m[0], m2, k1, k, d;
    for (int i = 1; i < n; i++)
    {
        a2 = a[i], m2 = m[i];
        int k2, tmp;
        exgcd(m1, m2, d, k1, k2);
        if ((a2 - a1) % d)
            return -1;
        tmp = m2 / d;
        k1 = (k1 * (a2 - a1) / d % tmp + tmp) % tmp; // 特解k'，保证最小正数
        a1 += m1 * k1;
        m1 *= m2 / d;
        a1 = (a1 + m1) % m1;
    }
    return a1; // 通解是a1+cm1,c是整数
}

方法二：中国剩余定理

条件：原方程组$m_1,m_2\cdots m_n$ 互素

定理：有模$M=m_1\cdot m_2\cdots m_n$ 的唯一解

推导：令$M_i=M/m_i$ ，由互素知$\gcd (M_i,m_i)=1$ ，又$M_i\cdot M_i^{-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$ ，则通解为$$cM+\sum _{i=1}^n{a}_i{M}_i{M}_i^{-1}$$

// 中国剩余定理解一元线性同余方程组，数组从0标号
ll CRT(int a[], int m[], int n)
{
    ll M = 1ll, ans = 0;
    for (int i = 0; i < n; i++)
        M *= m[i];
    for (int i = 0; i < n; i++)
        ans = (ans + a[i] * M % M * inv(M) % M) % M;
    return ans;
}