若干类型的同余方程的解法

若干类型的同余方程的解法

所有的证明、推导都写在引用里了，不愿意看可以跳过.

大量参考了这篇文章，感谢这位大佬的教导.

所有的例题都在我的同余题单.

线性同余方程

解形如

$$ax\equiv c(\mathrm{mod}b)$$

的同余方程.

容易根据同余的定义转化为解形如

$$ax+by=c$$

的二元一次不定方程.

结论 1 该不定方程有整数解的充要条件为 $\gcd (a,b)\mid c$.（裴蜀定理）

证明：
必要性
若 $ax+by=c$ 有整数解 $(x_0,y_0)$，则由于 $\gcd (a,b)\mid a,\gcd (a,b)\mid b$，有 $\gcd (a,b)\mid (a{x}_0+b{y}_0)$，即 $\gcd (a,b)\mid c$.
充分性
若方程 $ax+by=\gcd (a,b)$ 有整数解 $(x_0,y_0)$，则方程 $ax+by=c$ 对应有整数解 $\left(x_0\cdot \frac{c}{\gcd (a,b)},y_0\cdot \frac{c}{\gcd (a,b)}\right)$
故下证 $ax+by=\gcd (a,b)$ 必有整数解.
首先由于 $\gcd (a,b)\mid a,gcd(a,b)\mid b$，则对于 $a,b$ 的任意线性组合，有 $\gcd (a,b)\mid (ax+by)$. 设 $s$ 为 $ax+by$ 的最小正值，且满足 $s=a{x}_0+b{y}_0$，则有 $\gcd (a,b)|s$.
但对于 $a\mathrm{mod}s=a-\lfloor \frac{a}{s}\rfloor s=a-\lfloor \frac{a}{s}\rfloor (a{x}_0+b{y}_0)=a(1-\lfloor \frac{a}{s}\rfloor x_0)+b(-\lfloor \frac{a}{s}\rfloor y_0)$ 也为 $a,b$ 的线性组合，则 $a\mathrm{mod}s\le 0$ 或 $a\mathrm{mod}s\ge s$. 故 $a\mathrm{mod}s=0$，即 $s\mid a$，同理 $s\mid b$，从而 $s\mid \gcd (a,b)$. 又 $\gcd (a,b)\mid s$，所以 $s=\gcd (a,b)$，故 $ax+by=\gcd (a,b)$ 有解 $(x_0,y_0)$.

结论 2 若 $ax+by=c,(a\perp b)$ 的一组特解为 $(x_0,y_0)$，则其所有的解为

$$\{\begin{array}{l}x=x_0+bt\\ y=y_0-at\end{array}(t\in \mathrm{Z})$$

证明：先证该组 $(x,y)$ 为原方程的解：
$ax+by=a(x_0+bt)+b(y_0-at)=a{x}_0+b{y}_0=1$；
再证这是所有的解. 对于 $ax+by=1$ 的任意一组解 $(x^{\prime },y^{\prime })$，有 $a{x}^{\prime }+b{y}^{\prime }=1$，故 $a(x^{\prime }-x_0)+b(y^{\prime }-y_0)=0$，故 $a\mid b(y^{\prime }-y_0)$
又 $a\perp b$，则 $a\mid (y^{\prime }-y_0)$，则 $y^{\prime }-y_0=ka$. 容易化为题给形式.

推论 若 $ax+by=c$ 的一组特解为 $(x_0,y_0)$，则其所有的解为

$$\{\begin{array}{l}x=x_0+\frac{b}{\gcd (a,b)}\cdot t\\ y=y_0-\frac{a}{\gcd (a,b)}\cdot t\end{array}(t\in \mathrm{Z})$$

解方程的过程

Step 1 对于方程 $ax+by=c$，其中 $\gcd (a,b)\mid c$，我们先找出 $ax+by=\gcd (a,b)$ 的一组特解.

由于 $ax+by=\gcd (a,b)=\gcd (b,a\mathrm{mod}b)$，若我们先解出了方程 $b{x}^{\prime }+(a\mathrm{mod}b)y^{\prime }=\gcd (b,a\mathrm{mod}b)$ 的一组特解 $(x_0^{\prime },y_0^{\prime })$（这是一个递归过程），我们想办法从 $(x_0^{\prime },y_0^{\prime })$ 推出 $(x,y)$ 的一组解.
由 $\gcd (a,b)=\gcd (b,a\mathrm{mod}b)$，则 $a{x}_0+b{y}_0=b{x}_0^{\prime }+(a\mathrm{mod}b)y^{\prime }$，即 $a{x}_0+b{y}_0=b{x}_0^{\prime }+(a-\lfloor a/b\rfloor b)y_0^{\prime }=a{y}_0^{\prime }+b(x_0^{\prime }-\lfloor a/b\rfloor y_0^{\prime })$.
所以可以令 $x_0=y_0^{\prime },y_0=x_0^{\prime }-\lfloor a/b\rfloor y_0^{\prime }$，得到方程 $ax+by=\gcd (a,b)$ 的一组解.

由上述过程我们容易得到算法：

void solve(int a, int b, int &x, int &y) {
   if(!b) x = 1, y = 0; //边界情况
   else {
      int xp, yp;
      solve(b, a % b, xp, yp);
      x = yp; y = xp - a / b * yp;
   }
}

或者简单写为

void solve(int a, int b, int &x, int &y) {
   if(!b) x = 1, y = 0;
   else solve(b, a % b, y, x), y -= a / b * x;
}

Step 2 我们得到 $ax+by=\gcd (a,b)$ 的特解 $(x_0,y_0)$，则 $ax+by=c$ 的解为 $\left(x_0\cdot \frac{c}{\gcd (a,b)},y_0\cdot \frac{c}{\gcd (a,b)}\right)$

Step 3 我们得到 $ax+by=c$ 的特解为 $(x_0,y_0)$（注意这里的定义和 Step 2 不一样），则其通解为

$$\{\begin{array}{l}x=x_0+\frac{b}{\gcd (a,b)}\cdot t\\ y=y_0-\frac{a}{\gcd (a,b)}\cdot t\end{array}(t\in \mathrm{Z})$$

一些细节：

1. 解方程之前需要用裴蜀定理判断是否有解.
2. 使用 solve 时保证 a, b 均非负.
3. 时间复杂度为 $O(\log V)$，$V$ 为值域上限.
4. 可以用来求逆元.

线性同余方程组

解形如

$$\{\begin{array}{l}x\equiv b_1(\mathrm{mod}{m}_1)\\ x\equiv b_2(\mathrm{mod}{m}_2)\\ ⋮\\ x\equiv b_n(\mathrm{mod}{m}_n)\end{array}$$

的方程组

特殊情况 （可使用中国剩余定理解决）

当 $m_1,m_2,\dots ,m_n$ 两两互质时，可以使用中国剩余定理解决.

首先：

结论1 对于正整数 $n,m$，且 $n\perp m$，若两个整数 $a,b$，满足 $0\le a<n$，$0\le b<m$，则存在且仅存在一个 $x=x_0\in [0,nm)\cap \mathrm{Z}$ 满足

$$\{\begin{array}{l}x\equiv a(\mathrm{mod}n)\\ x\equiv b(\mathrm{mod}m)\end{array}$$

，并且方程组的通解为 $x\equiv x_0(\mathrm{mod}nm)$
下面证明这个定理：

首先我们证明恰存在一个 $x_0$：
反证法，设有两个解 $x_1,x_2$ 满足 $0\le x_1,x_2<nm$，且 $x_1<x_2$，使得 $x_1,x_2$ 满足方程组.
令 $k=x_2-x_1>0$，容易证明 $n\mid k,m\mid k$，那么 可知 $k\ge \mathrm{lcm}(n,m)=nm$（因为 $n\perp m$），于是 $x_2=k+x_1\ge nm+x_1\ge nm$，与假设矛盾.

然后我们证明方程组的通解：
第一，对于 $x\equiv x_0(\mathrm{mod}nm)$，即 $x\equiv x_0\equiv a(\mathrm{mod}n)$，$x\equiv x_0\equiv b(\mathrm{mod}m)$，则 $x\equiv x_0(\mathrm{mod}nm)$ 确实为原方程的解.
第二，对于 $x\not\equiv x_0(\mathrm{mod}nm)$，由于 $n\perp m$，容易推得 $x\not\equiv x_0\equiv a(\mathrm{mod}n)$，或 $x\not\equiv x_0\equiv b(\mathrm{mod}m)$，不为原方程的解

推论 可归纳证明：对于上面的那个 $n$ 个线性同余方程构成的线性同余方程组，在 $[0,{\prod }_{i=1}^n{m}_i)$ 内有唯一解，且通解为 $x\equiv x_0(\mathrm{mod}{\prod }_{i=1}^n{m}_i)$，其中 $x_0$ 为特解.

问题到了：如何构造一个特解？

对于方程

$$\{\begin{array}{l}x\equiv b_1(\mathrm{mod}{m}_1)\\ x\equiv b_2(\mathrm{mod}{m}_2)\\ ⋮\\ x\equiv b_n(\mathrm{mod}{m}_n)\end{array}$$

其中 $m_i$ 两两互质，

我们定义 $M_i=\frac{{\prod }_{k=1}^n{m}_k}{m_i}$，$k_i\equiv b_i{M}_i^{-1}(\mathrm{mod}{m}_i)$，那么特解为 $x_0={\sum }_{i=1}^n{k}_i{M}_i$.

神乎其技

简单证明一下：
我们考虑第 $k$ 个式子 $x\equiv b_k(\mathrm{mod}{m}_k)$ 是否满足.
对于 $x_0={\sum }_{i=1}^n{k}_i{M}_i$ 来说，
当 $i\ne k$ 时，有 $M_i\equiv 0(\mathrm{mod}{m}_k)$，则 $k_i{M}_i$ 贡献为 0.
当 $i=k$ 时，有 $k_i{M}_i\equiv b_i{M}_i^{-1}{M}_i\equiv b_i(\mathrm{mod}{m}_i)$.
于是可以得到 $x_0\equiv {\sum }_{i=1}^n{k}_i{M}_i\equiv b_i(\mathrm{mod}{m}_i)$，证毕.

时间复杂度 $O(n\log V)$，$V$ 为值域上限.

一般情况（方程合并法）

对于

$$\{\begin{array}{l}x\equiv b_1(\mathrm{mod}{m}_1)\\ x\equiv b_2(\mathrm{mod}{m}_2)\\ ⋮\\ x\equiv b_n(\mathrm{mod}{m}_n)\end{array}$$

并不保证 $m_i$ 互质，如何是好？我们考虑简化问题，

$$\{\begin{array}{l}x\equiv a(\mathrm{mod}n)\\ x\equiv b(\mathrm{mod}m)\end{array}$$

其中 $n,m$ 不一定互质，

容易得到

$$\{\begin{array}{l}x=a+pn\\ x=b+qm\end{array}$$

其中 $p,q\in \mathrm{Z}$.

于是

$$\begin{array}{rl}a+pn=& b+qm\\ pn-qm=& b-a\end{array}$$
上面 $p,q$ 是未知量，而其他均为已知量，故解线性同余方程.

两种情况：

1. 若 $\gcd (n,m)\nmid (b-a)$，方程无解，原不等式无解.

2. 若 $\gcd (n,m)\mid (b-a)$，

   则方程有特解 $(p_0,q_0)$ 满足 $p_{0}n-q_{0}m=b-a$，
   通解为 $p=p_0+t\cdot \frac{m}{\gcd (n,m)},q=q_0+t\cdot \frac{n}{\gcd (n,m)},t\in \mathrm{Z}$.
   所以 $x=a+pn=a+(p_0+t\cdot \frac{m}{\gcd (n,m)})n=a+p_{0}n+t\cdot \mathrm{lcm}(n,m)$

   故只要求出 $p_{0}n-q_{0}m=b-a$ 的一组特解 $(p_0,q_0)$，方程合并后得到 $x\equiv a+p_{0}n(\mathrm{mod}\mathrm{lcm}(n,m))$.

两个如此， $n$ 个也是如此.

时间复杂度 $O(n\log V)$.