快乐地打牢基础（7）——同余

一、基本定义和定理

---

定义1： 给定正整数$m$,若用两个整数$a$和$b$所得余数相同，称$a$和$b$所得余数相同，称$a$和$b$对模$m$同余，记作$a\equiv b(mod)m$，并称该式子为同余式；否则称$a$和$b$对模$m$不同余。

定义2： 整数的集合被分为$m$个不同的集合，这些集合被称为模$m$剩余类（同余类）。每个同余类中的任意两个整数都是模$m$同余的。

定理1： $a\equiv b(modm)$，当且仅当$m|(a-b)$。

定理2： $a\equiv b(modm)$，当且仅当存在整数$k$，使得$a=b+km$

费马小定理： 若$p$是质数，则对于任意整数$a$，有$a^p\equiv a(modp)$。

欧拉定理： 若正整数$a，n$互质，则$a^{\phi (n)}\equiv 1(modn)$，其中$\phi (n)$为欧拉函数。

欧拉定理的推论：
若正整数$a，n$互质，则对于任意正整数$b$，有$a^b\equiv a^{bmod\phi (n)}(modn)$。

证明：
设$b=q\times \phi (n)+r$，其中$0\le r<\phi (n)$，即$r=bmod\phi (n)$。于是：$a^b\equiv a^{q\times \phi (n)+r}\equiv (a^{q\times \phi (n)}{)}^q\times a^r\equiv 1^q\times a^r\equiv a^r\equiv a^{bmod\phi (n)}(modn)$。

许多计数类的题目要求我们把答案对一个质数$p$取模后输出，面对$a+b,a-b,a\times b$这样的算，可以在计算前先把$a,b$对$p$取模。面对乘方算式，根据欧拉定义的推论，可以先把底数对 $p$取模、指数对$\phi (p)$取模，在计算乘方。

二、同余的性质

---

对于整数$a,b,c$和自然数$m,n$，对模$m$同余满足：

• 自反性：$a\equiv a(modm)$。
• 对称性：若$a\equiv b(modm)$，则$b\equiv a(modm)$。
• 传递性：若$a\equiv b(modm)$，$b\equiv c(modm)$，则$a\equiv c(modm)$。
• 同加性：若$a\equiv b(modm)$，则$a+c\equiv b+c(modm)$。
• 同乘性：若$a\equiv b(modm)$，则$ac\equiv bc(modm)$。
  一般情况，$a\equiv b(modm)$，$c\equiv d(modm)$，则有$ac\equiv bd(modm)$
• 同幂性：若$a\equiv b(modm)$，则$a^n\equiv b^n(modm)$。
• 若$amodp=x，amodq=x$，其中$p,q$互质，则$amod(pq)=x$
  证明：
  因为$amodp=x，amodq=x$，则一定存在整数$s,t$，使$a=sp+x，a=tq+x$，所以$sp=tq$，
  又因为$p,q$互质，肯定是$s$和$p$存在一个倍数的关系，则一定存在整数$r$，使$s=rq$，
  将$s=rq$带回到$a=sp+x$得到$a=rpq+x$，得出：$amod(pq)=x$。

三、扩展欧几里得算法

---

欧几里得算法： gcd(a,b) = gcd(b,a%b)。

定理1： 设$a,b$不全为0，则存在整数$x$和$y$，使得$ax+by=gcd(a,b)$。
证明：

ll exgcd(ll a,ll b,ll &x,ll &y){
   
	if(b == 0){
   
		x = 1;
		y = 0;
		return a;
	}
	ll t = exgcd(b,a%b,x,y);
	ll x0 = x,y0 = y;//x0 和 y0 是下一次 递归中 x 和 y 的解 
	x = y0;
	y = x0 - a / b * y0;
	return t;
}

注意一下此处递归边界的问题：
当到达递归边界的时候，b==0，a=gcd(a,b) 这时可以观察出来这个式子的一个解：a1+b0=gcd(a,b)，故x=1,y=0。

定理2： 对于不定方程$ax+by=c$，当且仅当gcd(a,b)|c时，方程有解。

扩展欧几里得算法的作用：求解不定方程 $ax+by=gcd(a,b)$ 的一组特解$（x^{{}^{\prime }},y^{{}^{\prime }}）$。
原理分析：
首先由定理2可知，当gcd(a,b)不整除c的时候是无解的。所以过程中gcd(a,b)都是可以整除c的
因为gcd(a,b)可以整除c,所以设
$$g=gcd(a,b)a^{{}^{\prime }}=a/g,b^{{}^{\prime }}=b/g,c^{{}^{\prime }}=c/g$$
可以得到式子：$a^{{}^{\prime }}x+b^{{}^{\prime }}y=c^{{}^{\prime }}$，这个式子中的$x,y$就是原式的解。

我们可以用扩展欧几里得exgcd($a^{{}^{\prime }},b^{{}^{\prime }}$)求出不定方程$a^{{}^{\prime }}{x}^{{}^{\prime }}+b^{{}^{\prime }}{y}^{{}^{\prime }}=1$的整数解。

因为$a^{{}^{\prime }},b^{{}^{\prime }}$是互质的，所以gcd($a^{{}^{\prime }},b^{{}^{\prime }}$) =1，故可以使用扩展欧几里得求出不定方程$a^{{}^{\prime }}{x}^{{}^{\prime }}+b^{{}^{\prime }}{y}^{{}^{\prime }}=1$的整数解。

那么$a^{{}^{\prime }}{c}^{{}^{\prime }}{x}^{{}^{\prime }}+b^{{}^{\prime }}{c}^{{}^{\prime }}{y}^{{}^{\prime }}=c^{{}^{\prime }}$

$a^{{}^{\prime }}g{c}^{{}^{}}$