【数论】同余（五）：多元线性同余方程

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

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

多元线性同余方程

形如：$a_1{x}_1+a_2{x}_2+\cdots +a_n{x}_n\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$

由$a_{1}x+a_{2}y=\gcd (a_1,a_2)$ 有解可知左式能表达$c\gcd (a_1,a_2),c\in N$ ；引入$a_3$ ，则$gcd(a_1,a_2)x+a_{3}y$ 可以表示$c\gcd (a_1,a_2,a_3),c\in N$ 。以此类推，$a_1{x}_1+a_2{x}_2+\cdots +a_n{x}_n$ 可以表达$c\gcd (a_1,a_2,\cdots ,a_n),c\in N$ ，因此有解的条件为$\gcd (a_1,a_2,\cdots ,a_n,m)\mid b$ 。

一个很好理解的推导：我们将$m$ 添加到$a$ 数列中，仍设数量为$n$ 。设$d=\gcd (a_1,a_2,\cdots ,a_n),d_1=\gcd (a_1,a_2,\cdots ,a_{n-1})$ ，于是有$\gcd (d_1,a_n)=d$ ，则原式转化为$a_n{x}_n+d_{1}k=b$ ，可以通过前面一元线性同余方程的解法解出$x_n,k$ ，则问题转化为$a_1{x}_1+a_2{x}_2+\cdots +a_{n-1}{x}_{n-1}=d_1\ast k$，仍采取相同的解法。

但$O(n^2)$ 求每次的$\gcd$ 显然不现实。设$d_i=\gcd (a_1,a_2,\cdots ,a_i)$ ，$x[1...n]$ 表示解（从$1$ 标号，$a[n+1]$ 是$m$ ）则可以求解$d_i{X}_i+a_{i+1}{y}_i=d_{i+1}$ 的$X_i,y_i$ 。
观察最后一个式子：$$d_n{X}_n+a_{n+1}{y}_n=d_{n+1}$$ 要求解答案，等式两边乘上$c=\frac{b}{d_{n+1}}$ ，所以$$\begin{gathered} x_{n+1}=y_{n}c \\ b-a_{n+1}{x}_{n+1}=d_n{X}_{n}c \end{gathered}$$ 则前一个式子$$d_{n-1}{X}_{n-1}+a_n{y}_{n-1}=d_n$$ 要求解答案，等式两边乘上$c^{\prime }=c{X}_n$ ，依次求解即可。

// 求解a1x1+a2x2+...+anxn==b mod m，a从1开始标号
int X[N], y[N];
bool mod_Xn_equation(int a[], int b, int m, int n, int x[])
{
    a[n + 1] = m;
    int d = a[1];
    for (int i = 2; i <= n + 1; i++)
        exgcd(d, a[i], X[i - 1], y[i - 1]);
    if (b % d)
        return false;
    int c = b / d;
    x[n + 1] = y[n] * c;
    y[0] = 1;
    for (int i = n; i >= 1; i--)
    {
        c = (c * X[i] % m + m) % m;
        x[i] = (y[i - 1] * c % m + m) % m;
    }
    return true;
}