ACM中一些常用数学公式定理的简单推理及证明

数学蒟蒻表示这也太难了(sad)

---

更相减损术

$gcd(a,b)=gcd(a,b-a)$

$证明：$

若有$gcd(a,b)=d$，则有$a=k_{1}d,b=k_{2}d$，且$k_1$与$k_2$一定互质，即$gcd(k_1,k_2)=1$。
设$gcd(a,b-a)=xd$，也就是$gcd(k_{1}d,(k_2-k_1)d)=xd$，化简为$gcd(k_1,k_2-k_1)=x$。
同样的，可以设$k_1=t_{1}x,k_2-k_1=t_{2}x$。
有$k_1=t_{1}x,k_2=(t_1+t_2)x$。

因为前面有$gcd(k_1,k_2)=1$，所以此处$x$只能等于$1$。

所以，有$gcd(a,b)=gcd(a,b-a)$。

并且，该定理可以推广到多个数，即$gcd(a,b,c,...)=gcd(a,b-a,c-b,...)$。
简单证明：
$$\begin{gathered} gcd(a,b,c)=gcd(a,gcd(b,c)) \\ =gcd(a,b,c-b) \\ =gcd(gcd(a,b),c-b) \\ =gcd(a,b-a,c-b) \end{gathered}$$。

---

辗转相除法（欧几里得算法gcd）

$gcd(a,b)=gcd(b,amodb)$

$证明:$
设$a=kb+r$，则$gcd(b,amodb)=gcd(b,r)$
若$gcd(a,b)=d$,则有$d|a$且$d|b$，而$r=a-kb$，所以一定也有$d|r$。
所以$gcd(b,r)=d$

证得:$gcd(a,b)=gcd(b,amodb)$

---

扩展欧几里得算法（exgcd）

请看这里。

---

快速幂

快速幂是用来快速的求得$a^b$。

$简单推理:$
将b按二进制权展开，有
$b=k_n{2}^n+k_{n-1}{2}^{n-1}+...+k_1{2}^1+k_0{2}^0$
$k$的取值要么是$0$,要么是$1$。
则$$a^b=a^{k_n{2}^n+k_{n-1}{2}^{n-1}+...+k_1{2}^1+k_0{2}^0}$$
因为$a^{2^i}\ast a^{2^i}=a^{2^{i+1}}$，由这个性质可以快速的计算$a^b$。

具体操作是：
以$3^{45}$为例 \2,954,312,706,550,833,698,643
$45=(101101{)}_2$
从二进制的低位到高位计算。
初始时答案$ans=1$，目前的次方数为$p=3^{2^0}$
如果二进制第0位是1，让$ans=ans\ast p$。
如果是$0$就不需要乘。
每操作完一位后，让$p=p\ast p$，此时$p=3^{2^1}$，然后进行下一位的判断。
.
最后得到$ans=3^{2^0}\ast 3^{2^2}\ast 3^{2^3}\ast 3^{2^5}$。
如果一个一个乘，我们需要乘$45$次，而快速幂只需要$lo{g}_{2}45$次。

---

逆元

请看这里。

---

线性同余方程

线性同余方程就是给出一个这种形式的方程，求出x:
$$ax\equiv b(modm)$$

等价于求:

$ax=my+b$
$ax-my=b$（因为y是一个未知数，可以把y前面的负号合并到y里面，不影响x的值）
也就是:$ax+my=b$

由exgcd知，一定有$a{x}_0+m{y}_0=d$，（$d$是$gcd(a,m)$）。

$a{x}_0+m{y}_0=d$

两边同时乘一个$\frac{b}{d}$:
$\frac{b}{d}\ast a{x}_0+m{y}_0=\frac{b}{d}\ast d$

也就是:
$a(\frac{b}{d}{x}_0)+m(\frac{b}{d}{y}_0)=b$

得到一组特解为:
$\{\begin{array}{l}x=\frac{b}{d}{x}_0\\ y=\frac{b}{d}{y}_0\end{array}$

通解为:

设有一组特解$x0,y0$
$\{\begin{array}{l}x=x_0+k\frac{m}{d}\\ y=y_0-k\frac{a}{d}\end{array}$

注意：这里的$x_0,y_0$不是上面的$x_0,y_0$。

---

约数和质数

1、一个数的约数个数和约数之和的求法。

2、线性筛法简单证明：
线性筛法的思想是对于每一个数，只用它的最小的质因子去筛它。

例，求$1$~$n$内的所有质数。

具体执行是:

首先，按顺序枚举$2$—>$n$的每一个数。
枚举到$i$的时候，如果$i$没有被$1$~$i-1$中的任何一个质数筛掉，说明$i$本身就是一个质数，那 么将$i$加入质数集合中,如果被筛掉了(说明不是质数)，就不需要加。

接下来，在质数集合中从小到大枚举现有的所有质数。
枚举到$p_j$时，令$x=p_j\ast i$。
如果$i\%p_j!=0$，说明$p_j$一定是$p_j\ast i$的最小质因子，将它筛去(因为它不是质数)。
否则，$i\%p_j==0$，说明$p_j$一定是$i$的最小质因子(因为$p_j$是从小到大枚举的，当第一次出现$i\%p_j==0$的情况时，一定成立)，也是$p_j\ast i$的最小质因子，同样把它筛掉，然后停止枚举(也就是，当出现了$i\%p_j==0$的情况后，就不在往下枚举了)。

我们可以证明，对于任意一个数$x$，设$x$的最小质因子是$p$。
当上面的$i$枚举到$\frac{x}{p}$的时候，$x$一定会在这个时候被筛去。所以不会存在某一个数被漏掉的情况。

---

欧拉函数

请看这里。

---

费马小定理与欧拉定理

费马小定理：
如果$p$是一个质数，而整数$a$与$p$互质，则有$$a^{p-1}\equiv 1（modp）$$

欧拉定理：
对于任意互质的$a$和$n$，有$$a^{\phi (n)}\equiv 1(modn)$$
其中$\phi (n)$是$n$的欧拉函数。

(费马小定理是欧拉定理的一个特例，即$n$为质数的时候)。

$证明：$
设所有与$n$互质的数的集合$S=${$x_1,x_2,...,x_{\phi (n)}$}，共$\phi (n)$个。
在设另一个集合$Z=${$y_1,y_2,...,y_{\phi (n)}$}，其中，$y_i=a{x}_i(modn)$。

则有：在模n的意义下，集合$S=Z$。

$Why?$
因为$a$与$n$也是互质的，所以对于每一个$y_i$，都有$y_i$与n互质，且它们都是互不相同的。而与$n$互质的数就$\phi (n)$个，集合Z正好也有$\phi (n)$个与$n$互质的数。所以，这$\phi (n)$个数就是那$\phi (n)$个数，集合$S$与集合$Z$表示的其实都是同一个集合。

所以，有
$$\begin{gathered} y_1\ast y_2\ast ...\ast y_{\phi (n)}\equiv x_1\ast x_2\ast ..\ast x_{\phi (n)}(modn) \\ {a}^{\phi (n)}\ast x_1\ast x_2\ast ..\ast x_{\phi (n)}\equiv x_1\ast x_2\ast ..\ast x_{\phi (n)}(modn) \end{gathered}$$
化简一下有：
$$a^{\phi (n)}\equiv 1(modn)$$。

---

中国剩余定理与扩展中国剩余定理

代码展示请看这里。

中国剩余定理是指：

解这样一个同余方程组，其中$m_1,m_2,...,m_n$两两互质。
$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\\ .....\\ x\equiv a_n(mod{m}_n)\end{array}$

做法：

令$M=m_1\ast m_2\ast ...\ast m_n$。
在模$M$的意义下，方程只有一个解。
再令$M_i=\frac{M}{m_i}$，对于每一个$M_i$，求出$M_i^{-1}$（指$M_i$模$m_i$的逆元）。
(如果上面的模数不是两两互质的，会存在有的$M_i$ 的逆元求不出来)

则可以构造出：
$$x=a_1{M}_1{M}_1^{-1}+a_2{M}_2{M}_2^{-1}+...+a_n{M}_n{M}_n^{-1}(modM)$$

简单证明：
检查该解是否满足上面的$n$个线性同余方程。
对于第一个方程：$x\equiv a_1(mod{m}_1)$。
因为$M_1^{-1}是M_1在模m_1下的逆元$。所以，在模$m_1$的情况下，$M_1{M}_1^{-1}=1$。所以解$x$的第一项$a_1{M}_1{M}_1^{-1}mod{m}_1=a_1$，而除此之外，其他的任意一个$M_i(i!=1)$，$M_i$都是$m_1$的倍数，它们在模$m_1$的情况下都是等于$0$的。所以整个式子模$m_1$是等于$a_1$的。
同样的，对于其他方程也可以得出类似的结论。
所以，该式子是正确的。

---

中国剩余定理要求$m_1,m_2,...,m_n$两两互质，如果不互质怎么办那？这里就要用到扩展中国剩余定理了。

鸽。

扩展中国剩余定理的做法完全不同于上面，它是采用了两两合并的思想。
对于$n$个线性同余方程组，它每次合并两个，合并$n-1$次之后就剩下了一个方程，之后在解就十分简单了。

来看一看它是怎么样合并的。

$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\\ .....\\ x\equiv a_n(mod{m}_n)\end{array}$

我们只看前两个方程：
$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\end{array}$

对于这两个方程，改写一下方程的表示，可以这样写：
$$x=k_1{m}_1+a_1\text{\hspace{1em}(1)}$$
$$x=k_2{m}_2+a_2\text{\hspace{1em}(2)}$$

于是有：
$k_1{m}_1+a_1=k_2{m}_2+a_2$
整理一下：
$k_1{m}_1-k_2{m}_2=a_2-a_1$（其中，$k_1$和$k_2$是未知量）
可以用exgcd先求出一组$k_1$和$k_2$的解，记为 $k^{{}^{\prime }}$ 和 $k^{{}^{\prime \prime }}$。
这里要判断一下，如果$gcd(m_1,m_2)$不能整除$(a_2-a_1)$方程是无解的。
那么$k_1$和$k_2$的通解形式可以表示为：

$\{\begin{array}{l}{k}_1=k^{{}^{\prime }}+k\frac{m_2}{d}\\ k_2=k^{{}^{\prime \prime }}+k\frac{m_1}{d}\end{array}$
（$d$是$gcd(m_1,m_2)$，$k$是任意整数）。

那么仅对于$(1)$方程和$(2)$方程来说，$x$的通解就是：
$x=(k^{{}^{\prime }}+k\frac{m_2}{d})m_1+a_1$，（将$k_1$的通解带入到$(1)$方程中）

化简一下有:

$x=m_1{k}^{{}^{\prime }}+k\frac{m_1{m}_2}{d}+a_1$
我们可以发现，$\frac{m_1{m}_2}{d}$就是$m_1$和$m_2$的最小公倍数，记为$[m_1,m_2]$。

于是有

$x=k[m_1,m_2]+m_1{k}^{{}^{\prime }}+a_1$，其中$[m_1,m_2]$和$m_1{k}^{{}^{\prime }}+a_1$都是定值。分别记为$m=[m_1,m_2]，a=m_1{k}^{{}^{\prime }}+a_1$。

那么最终可以化为：
$$x=km+a\text{\hspace{1em}(3)}$$

然后观察一下$(3)$式与$(1)$式和$(2)$式，可以惊奇的发现，它们好像是一样的？！

最终，两个方程就合并成了一个方程，我们只用方程$(3)$就可以代表方程$(1)$和方程$(2)$。

合并完之后，$n$个方程变成了$n-1$个方程，然后继续合并，直到方程数为$1$。

---

高斯消元解线性方程组

请看这里。

---

组合数

请看这里。

---

简单博弈论

请看这里。

---

整除分块

请看这里。

---

莫比乌斯反演

请看这里。