基础数论part2 一些证明

GCD（辗转相除法）

最大公约数

求两个数a和b的最大公约数，考虑一般情况，设a > b

约定符号|，$A|B$表示A是B的约数

这里，假如一个数d，满足 $d|a$ ， 并且 $d|b$ ，（即：d为a与b的共同的约数

则 $d|（x\times a+y\times b）$ ,(x,y为任意整数) ，

把a拆开，a可以表示为$a=k\times b+r$ (这里的 r为$a\mathrm{m}\mathrm{o}\mathrm{d}b$， 即a除以b得到的余数)

$ k = \lfloor \frac{a}{b} \rfloor , 0 \leq r < b$

带入得 $r = a - \lfloor \frac{b}{a} \rfloor \times b $

满足前面$x\times a+y\times b$的形式，所以 $d|r$，即 $d|(a\mathrm{m}\mathrm{o}\mathrm{d}b)$

即，可以证明，任意d，即a与b任意的约数都同时是$a\mathrm{m}\mathrm{o}\mathrm{d}b$ 的约数,

所以求a和b的最大公约数可以转化为求b与$a\mathrm{m}\mathrm{o}\mathrm{d}b$的最大公约数

约定0和任意正整数的最大公约数为那个正整数本身，即

$gcd(a,0)=a$

---

裴蜀定理（贝祖定理）

扩展欧几里得定理（exgcd）

定义 ：设a，b是不全为0的整数，则存在整数x，y，使得$ax+by=gcd(a,b)$

---

（类比刚才的gcd的求法）

存在，$ax+by=gcd(a,b)$ ——（1）

也存在 $b{x}_1+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y_1=gcd(a,b)$ (因为这里的$gcd(a,b)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$)——（2）

（1）和（2）的等式右边都是$gcd(a,b)$，所以可以连等（

尝试把后面这个式子展开

可以得到
$$\begin{gathered} b{x}_1+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y_1 \\ =b{x}_1+(a-\lfloor \frac{a}{b}\rfloor \times b)y_1 \\ =b{x}_1+a{y}_1-\lfloor \frac{a}{b}\rfloor \times b{y}_1 \\ =a{y}_1+b(x_1-\lfloor \frac{a}{b}\rfloor \times y_1) \end{gathered}$$
观察（1）式，得到这里的
$$\{\begin{array}{l}x=y_1\\ y=x_1-\lfloor \frac{a}{b}\rfloor \times y_1\end{array}$$
可以注意到求$ax+by=gcd(a,b)$的x，y和求$b{x}_1+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y_1=gcd(a,b)$ 的$x_1,y_1$是可以互相转化的

而将前一个问题（（1）式的问题）转化为后一个问题（（2）式的问题） 可以将问题规模变小：

因为这里的$a>b$, $b>a\mathrm{m}\mathrm{o}\mathrm{d}b$）

可以证明他们是趋向于0的，

我们来考虑一下当$b=0$时候这个问题的解，

$a\times x+0\times y=gcd(a,0)=a$,

因为证明存在，我们可以任意取一组,比如
$$\{\begin{array}{l}x=1\\ y=0\end{array}$$
将它倒推回去，就可以得到原问题的解

---

费马小定理

定义 ： 若p为素数，gcd(a,p) = 1，则$a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

---

应用 ： 降幂 ：$a^b\equiv a^{b\mathrm{m}\mathrm{o}\mathrm{d}(p-1)}(\mathrm{m}\mathrm{o}\mathrm{d}p)$

---

逆元

定义 ： 如果一个线性同余方程$ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ ，则x称为a mod p的逆元，记作$a^{-1}$

---

求法 ：

• sol1 扩展欧几里得

$ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

可以转化为 $a\times x=1+p\times t$ （$t\in Z$）的形式

a，p互质，可以转化为求 $a\times x+p\times y=1$

---

• sol2 快速幂法

$ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ , 求x，费马 $a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

$ax\equiv a^{p-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$

$x\equiv a^{p-2}(\mathrm{m}\mathrm{o}\mathrm{d}p)$

再用快速幂

---

• sol3 线性求逆元 （公式 不要求掌握

n次计算算出1到n所有数对于p的逆元，复杂度是线性的

---

$1^{-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

把p拆开：

设 $p=k\times i+j$， $j<i,1<i<p$， 再放到mod p意义下 （这里i是1到 p-1 任意的整数）

先不要管为什么这样拆，先看下去）

$k\times i+j\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$

两边同时乘上$i^{-1},j^{-1}$

得到 $k{j}^{-1}+i^{-1}\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$

$i^{-1}\equiv -k{j}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 代入 $p=k\times i+j$

$i^{-1}\equiv -(\lfloor \frac{p}{i}\rfloor )(p\mathrm{m}\mathrm{o}\mathrm{d}i{)}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$

---

代码

将i对于p的逆元记作inv[i]，

$-(\lfloor \frac{p}{i}\rfloor )\equiv (p-(\lfloor \frac{p}{i}\rfloor ))(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 在取模的意义下是一样的，是将负数变成正数

inv[1] = 1;
for (int i = 2; i <= n; ++i)
    inv[i] = (long long)(p - p / i) * inv[p % i] % p;

因为 j （即$p\mathrm{m}\mathrm{o}\mathrm{d}i$）小于 i，我们这里又是按顺序求每个 i 的逆元，在求i的逆元用到 j 的时候，可以保证一定求过 j

---