Modular Arithmetic 模算术

Modular Arithmetic 模算术

我们都见过时钟，在时钟上有12个刻度。假设某天晚上时针指向6的位置，那么它表示的是晚上六点钟，可是到了第二天早上，当时针再次指向6时，它表示的又是早上六点钟。

在上述例子中，我们实际上是用同一个刻度来表示两个不同的“六点钟”，这就是模算术的基本思路–用同一个刻度表示不同的值（也可以用轮回来理解）。

模算术的定义
令$a$，$b$和$n$为整数且$n>0$，那么当且仅当$n$可以整除$(a-b)$时，$a$和$b$是模$n$同余的，写作：

$$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}n)⟺n|(a-b)$$

这里的 $n|(a-b)$ 表示 $kn=(a-b),k$为整数。

什么是同余 Congruence Modulo
当$kn=(a-b)$时，$a$和$b$可写为：
$$\begin{array}{rl}a& =pn+r\\ b& =qn+r\end{array}$$
$r$是他们共同的余数，$n$为同余的模数（也就是除数）。

例如：$7(=1\times 6+1)$和$13(=2\times 6+1)$对模数6同余。 类似的数还有19，25，31，$\dots$，$(n\times 6+1)$。将这些数画在数轴上，我们可以发现在每一个长度为6的区间上，都只有一个这样的数存在，并且他们都出现在各自区间内“相同的位置”上。因此我们可以在不同区间内用同一个刻度表示不同的值。

由对于模$n$同余的所有整数组成的集合称为同余类。

同余类 Congruence Class

给定任意整数$b>0$，我们可以将数轴上的整数划分到长度为$b$的区间内。对于$a\in [0,b-1]$，集合 { a + k b } \{a + kb\} {a+kb}的所有数对于模$b$同余，余数为$a$，构成一个同余类，记作：
[ a ] b = { a + k b ∣ k ∈ Z } [a]_b = \{a+kb | k \in \mathbb{Z} \} [a]b​={a+kb∣k∈Z}

模算术的规则

• 加法：如果有$a\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}n)$和$b\equiv d(\mathrm{m}\mathrm{o}\mathrm{d}n)$，那么$a+b\equiv c+d(\mathrm{m}\mathrm{o}\mathrm{d}n)$

• 乘法：如果有$a\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}n)$和$b\equiv d(\mathrm{m}\mathrm{o}\mathrm{d}n)$，那么$ab\equiv cd(\mathrm{m}\mathrm{o}\mathrm{d}n)$

• $\forall a,a\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}n)$

• 如果有$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}n)$，那么$b\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}n)$

• 如果有$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}n)$和$b\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}n)$，那么$a\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}n)$

• 如果有$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}n)$和$n\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}m)$，那么$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$

• 如果有$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}n)$，那么$am\equiv bm(\mathrm{m}\mathrm{o}\mathrm{d}n)$和$a^m\equiv b^m(\mathrm{m}\mathrm{o}\mathrm{d}n)$

模的倒数 Multiplicative Inverse
对于$a\not\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}n)$，若有整数$a^{\prime }$满足$a{a}^{\prime }\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$，则称$a^{\prime }$为$a$的倒数。例如：$2\times 5\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}3)$，所以$a^{\prime }=5$是整数$2(\mathrm{m}\mathrm{o}\mathrm{d}3)$的倒数（反之亦然）。

并不是所有的整数对于任意模数都有倒数。对于$a(\mathrm{m}\mathrm{o}\mathrm{d}n)$，如果整数$a$和$n$不互质，那么$a$就没有对于模数$n$的倒数。

中国余数定理

古时有一位将军统领着一支1500人的军队，一场大战过后，大约400~500人死亡。为了统计剩下的人数，将军让士兵们站成方阵。当剩下的人3个站成一行时，多余2个人。当剩下的人5个站成一行时，多余4个人。当剩下的人7个站成一行时，多余6个人。那么这支军队总共还有多少士兵？

这个问题实际上就是求解同余方程组
$$\{\begin{array}{r}{c}_{1}x\equiv d_1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_1)\\ c_{2}x\equiv d_2(\mathrm{m}\mathrm{o}\mathrm{d}{m}_2)\\ \dots \\ c_{k}x\equiv d_k(\mathrm{m}\mathrm{o}\mathrm{d}{m}_k)\end{array}$$

先来思考如何求解 $ax\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}n)$，

1. 当$\mathrm{g}\mathrm{c}\mathrm{d}(a,n)=1$时：$a$和$n$互质，因此存在$a^{\prime }$使得$a{a}^{\prime }\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$。我们对式子两边同乘$a^{\prime }$，得
   $$x\equiv a^{\prime }b(\mathrm{m}\mathrm{o}\mathrm{d}n)$$
2. 当 $\mathrm{g}\mathrm{c}\mathrm{d}(a,n)=c>1$时：
   2.1. 如果$c$可以整除$b$：对于$ax\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}n)$可得
   $$\begin{array}{rl}& ax\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}n)\\ ⟺& ax=b+nk\\ ⟺& a_{1}cx=b_{1}c+n_{1}ck\\ ⟺& a_{1}x=b_{1}x+n_{1}k\end{array}$$
   此时$a_1,n_1$互质，因此可以归为$a,n$互质一类。
   2.2. 如果$c$不能整除$b$：
   $$\begin{array}{rl}& ax\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}n)\\ ⟺& ax=b+nk\\ ⟺& a_{1}cx=b+n_{1}ck\\ ⟺& b=c(a_{1}x-n_{1}k)\end{array}$$
   与假设矛盾，无解。

经过化简（当然，也存在无法化简的情况）以后，方程组可写为
$$\{\begin{array}{r}x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{n}_1)\\ x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{n}_2)\\ \dots \\ x\equiv a_k(\mathrm{m}\mathrm{o}\mathrm{d}{n}_k)\end{array}$$
假设$n_1,n_2,\dots ,n_k$之间两两互质，则该方程组对于模$N=n_1{n}_2\dots n_k$有唯一解，解为
$$x\equiv (N_1{N}_1^{\prime }{a}_1+N_2{N}_2^{\prime }{a}_2+\cdots +N_k{N}_k^{\prime }{a}_k)(\mathrm{m}\mathrm{o}\mathrm{d}N)$$
其中$N_i=\frac{N}{n_i}$，$N_i^{\prime }$为$N_i$对于模$n_i$的倒数。

取模运算
在计算机中，给定整数$a$，$b$，取模运算就是求$a$除以$b$的余数：$a$%$b$
关于取余运算的实现逻辑可参考汇编除法原理