计算机基础-取模

引言.

在做算法题的过程中，经常需要取模操作，对结果模上1e9+7或者1e9+9或者998244353。但有时会发现题目A不了最后是取模的问题，为了在日后的题目完成过程中避免因取模而出错，因此整理一下取模相关的知识。

1. 取模运算和取余运算：

这一部分的内容介绍取模和取余的区别，是作为基础补充，到后续部分会用到，可先跳过该节，到时候反过来看

对于整数a和b，$a\%b$ 的运算方式如下（取余和取模一样）：

• 求整数商：$c=[a/b]$
• 计算模或者余数：$r=a-c\times b$

在整数商中有一个取整操作，取模和取余的不同在于：

运算	取整方式
取模	求整数商时, 取整向负无穷靠近, 即: $-1/3\to -1$
取余	求整数商时, 取整向0靠近, 即: $-1/3\to 0$

很容易发现, 当$a/b$为正数时, 取余和取模的结果是一样的, 只有当其为负数时, 两者才存在区别. 为了方便理解, 下面给出一个例子: 求$-10\%3$ 的模和余数.

运算	整数商	结果
取模	-4	2
取余	-3	-1

接下来要记住一个结论:

• C/C++, Java 的%为**取余**
• Python的%为**取模**

可以看到, 取余运算在a或者b为负数时得到的余数仍为负数, 而取模则能得到正数, 这一点在后续的做题过程中会意识到.

2. 为什么总是对这几个数取模?

按照规律一条条陈述, 结论不言而喻

• 整数类型(int, 64位操作系统上, 占4个字节)的范围是: $-2^{31}\sim 2^{31}-1$
• 1e9+7, 1e9+9, 998244353这三个数都小于$2^{30}$, 因而对于任意两个模过之后的数$a$和$b$:
  • $\text{a,b}<2^{30}$
  • $a+b<2^{31}$
  • $ab<2^{60}$, 在long long范围内不会溢出
  • 可以使用扩展欧几里德算法

3. 扩展欧几里德算法¹:

3.1 欧几里德算法(也叫辗转相除法):

gcd(a, b) = gcd(b, a mod b)

3.2 扩展欧几里德算法:

贝祖定理(也叫裴蜀定理)： 给定两个整数a, b, 必定存在整数x、y使得$ax+by=gcd(a,b)$ $\Rightarrow$ 如果ax+by=1有解，那么gcd(a,b)=1；

扩展欧几里德算法可以用来求模反元素（模逆元）的。

3.2.1 符号约定:

符号	含义
$gcd(\cdot )$	欧几里德算法
$[\cdot ]$	向0取整

3.2.2 算法思想:

$$\begin{gathered} ax+by=gcd(a,b)=gcd(b,a\%b)=gcd(b,a-[\frac{a}{b}]b) \\ 取某一x^{\prime },y^{\prime }得:ax+by=b{x}^{\prime }+(a-[\frac{a}{b}]b)y^{\prime }=gcd(a,b) \\ 故而:a(x-y^{\prime })+b(y-x^{\prime }+[\frac{a}{b}]y^{\prime })=0恒成立 \\ \Rightarrow \{\begin{array}{l}x=y^{\prime }\\ y=x^{\prime }-[\frac{a}{b}]y^{\prime }\end{array} \end{gathered}$$

类似于gcd($\cdot$), x’和y’的求法也是递归进行的, gcd($\cdot$)的终结条件是: b = 0, 此时:
$$\{\begin{array}{l}x=y^{\prime }=1\\ y=0\end{array}$$
得到递归基后, 可以将递归基作为新的x’, y’, 来求上一步的x, y

3.2.3 求模逆元:

对于$ab\%n=1$, b称为a的模逆元(反之亦然)

3.2.3.1 费马小定理:

假设p为质数, 那么$a^{(p-1)}\%p=1$.

推论: $a^{p-2}$为a%p的乘法逆元 ($a\times a^{p-2}\%p=a^{p-1}\%p=1$)

3.2.3.2 扩展欧几里德算法: [待完善, 还没完全理清]

$$\begin{gathered} ab\%p=1\Rightarrow ab\%p+kp\%p=1\text{ (kp\%p必然为0, 因此最外层不需要加\%)} \\ 若gcd(a,p)=1,则:ax+py=1 \end{gathered}$$

参考资料:

---

1. https://zhuanlan.zhihu.com/p/58241990 ↩︎