初等数论【咕咕咕】

$\text{LATEX}$ 太那啥了啊（又悲又恼

如有错误或 $\text{LATEX}$ 不规范的地方，请在评论区指出

1. 整除

1.1. 定义

设 $a,b\in \mathbb{Z},a\ne 0$ ，如果 $\exists q$ ，使 $aq=b$ ，则 $b$ 能被 $a$ 整除，记作 $a\mid b$ ，否则， $b$ 不能被 $a$ 整除，记作 $a\nmid b$

1.2. 性质

1.2.1. 性质1

若 $a\mid b,b\mid c$ ，则 $a\mid c$

证明：

$$\begin{array}{rl}\because & a\mid b,b\mid c\\ & \text{令}ax=b,by=c(x,y\in \mathbb{Z},x\ne 0,y\ne 0)\\ \therefore & axy=c\\ \because & x,y\in \mathbb{Z}\\ \therefore & xy\in \mathbb{Z}\\ \therefore & a\mid c\end{array}$$

1.2.2. 性质2

若 $a\mid b,a\mid c$ ，则有 $\forall x,y\in \mathbb{Z}$ ，满足 $a\mid (bx+cy)$

证明：

$$\begin{array}{rl}\because & a\mid b,a\mid c\\ & \text{令}a{x}_1=b,a{y}_1=c(x_1,y_1\in \mathbb{Z},x_1\ne 0,y_1\ne 0)\\ & \begin{array}{rl}bx+cy& =a{x}_{1}x+a{y}_{1}y\\ & =a(x_{1}x+y_{1}y)\end{array}\\ \because & a\mid a(x_{1}x+y_{1}y)\\ \therefore & a\mid (bx+cy)\end{array}$$

1.2.3. 性质3

若 $a\mid b$ ，则 $am\mid bm(m\ne 0)$

证明：

$$\begin{array}{rl}\because & a\mid b\\ & \text{令}ax=b(x\in \mathbb{Z},x\ne 0)\\ \therefore & amx=bm\\ \therefore & am\mid bm\end{array}$$

1.2.4. 性质4

若 $\forall x,y\in \mathbb{Z}$ 满足 $ax+by=1$ ，且 $a\mid n,b\mid n$ ，则 $ab\mid n$

证明：

$$\begin{array}{rl}\because & a\mid n,b\mid n\\ & \text{令}a{x}_1=n,b{y}_1=n(x_1,y_1\in \mathbb{Z},x_1,y_1\ne 0)\\ & \begin{array}{rl}ax+by& =1\\ \frac{x}{b}+\frac{y}{a}& =\frac{1}{ab}\\ \frac{xn}{b}+\frac{yn}{a}& =\frac{n}{ab}\\ \frac{b{y}_{1}x}{b}+\frac{a{x}_{1}y}{a}& =\frac{n}{ab}\\ x{y}_1+x_{1}y& =\frac{n}{ab}\end{array}\\ \because & x,y,x_1,y_1\in \mathbb{Z}\\ \therefore & x{y}_1+x_{1}y\in \mathbb{Z}\\ \therefore & \frac{n}{ab}\in \mathbb{Z}\\ \therefore & ab\mid n\end{array}$$

1.2.5. 性质5

若 $b=dq+c,q\in \mathbb{Z}$ ，那么 $d\mid b$ 的充要条件为 $d\mid c$

证明1（若 $d\mid b$ ，则 $d\mid c$）：

$$\begin{array}{rl}\because & d\mid b\\ & \text{令}dk=b(k\in \mathbb{Z},k\ne 0)\\ \therefore & b=dq+c⟺dk=dq+c\\ \therefore & d\mid c\end{array}$$

证明2（若 $d\mid c$ ，则 $d\mid b$）：

$$\begin{array}{rl}\because & d\mid c\\ & \text{令}dk=c(k\in \mathbb{Z},k\ne 0)\\ \therefore & b=dq+c⟺b=d(p+k)\\ \therefore & d\mid b\end{array}$$

2. 模运算

2.1. 定义

定义：有 $\forall a,b\in \mathbb{Z},b\ne 0$ ，求 $a÷b$ 的余数，称 $a$ 模 $b$ ，记作 $a\mathrm{mod}b$

2.2. 性质

2.2.1. 性质1

$(a+b)\mathrm{mod}c=(a\mathrm{mod}c+b\mathrm{mod}c)\mathrm{mod}c$

2.2.2. 性质2

$(a-b)\mathrm{mod}c=(a\mathrm{mod}c-b\mathrm{mod}c)\mathrm{mod}c$

2.2.3. 性质3

$(ab)\mathrm{mod}c=[(a\mathrm{mod}c)\times (b\mathrm{mod}c)]\mathrm{mod}c$

证明：

$$\begin{array}{rl}\text{设}& a\mathrm{mod}c=m_a,b\mathrm{mod}c=m_b(m_a,m_b\in \mathbb{Z})\\ \because & a\mathrm{mod}c=m_a,b\mathrm{mod}c=m_b\\ \text{令}& a=m_a+x{s}_a,b=m_b+c{s}_b(s_a,s_b\in \mathbb{Z})\\ & \begin{array}{rl}(ab)\mathrm{mod}c& =[(m_a+c{s}_a)(m_b+c{s}_b)]\mathrm{mod}c\\ & =(m_a{m}_b+c{m}_a{s}_b+c{m}_b{s}_a+c^2{s}_a{s}_b)\mathrm{mod}c\\ & =(m_a{m}_b)\mathrm{mod}c\\ & =[(a\mathrm{mod}c)\times (b\mathrm{mod}c)]\mathrm{mod}c\end{array}\end{array}$$

2.2.4. 性质4

属于性质三的一种扩展

$(a^b)\mathrm{mod}c=(a\mathrm{mod}c{)}^b\mathrm{mod}c$

2.2.5. 性质5

如果 $a\mathrm{mod}b=c,d\ne 0,d\in \mathbb{Z}$ ，则 $(ad)\mathrm{mod}(bd)=cd$

证明：

$$\begin{array}{rl}\because & a\mathrm{mod}b=c\\ & \text{设}bx+c=a(x\in \mathbb{Z})\\ \therefore & bx+c=a⟺bdx+cd=ad\\ \therefore & (ad)\mathrm{mod}(bd)=cd\end{array}$$

2.2.6. 性质6

如果 $a\mathrm{mod}b=c,d\mid a,d\mid b$ ，则 $(\frac{a}{d})\mathrm{mod}(\frac{b}{d})=\frac{c}{d}$

证明：

$$\begin{array}{rl}\because & a\mathrm{mod}b=c\\ & \text{设}bx+c=a(x\in \mathbb{Z})\\ \therefore & bx+c=a⟺\frac{b}{d}x+\frac{c}{d}=\frac{a}{d}\\ \therefore & (\frac{a}{d})\mathrm{mod}(\frac{b}{d})=\frac{c}{d}\end{array}$$

3. 同余

3.1. 定义

设 $m$ 是给定的 ${\mathbb{Z}}^{+}$ ，若满足 $m\mid (a-b)$ ，则称 $a$ 与 $b$ 对模 $m$ 同余，记作 $a\equiv b(\mathrm{mod}m)$ ；若满足 $m\nmid (a-b)$ ，则称 $a$ 与 $b$ 对模 $m$ 不同余

如果将定义反过来，如果 $a\equiv b(\mathrm{mod}m)$ ，则 $m\mid (a-b)$（即 $a\equiv b(\mathrm{mod}m)⟺m\mid (a-b)$）

证明1（已知 $m\mid (a-b)$）：

$$\begin{array}{rl}\because & m\mid (a-b)\\ & \text{令}mk=a-b(k\in \mathbb{Z},k\ne 0)\\ \therefore & mk=a-b⟺mk+b=a⟺b\mathrm{mod}m=a\mathrm{mod}m\\ \therefore & a\equiv b(\mathrm{mod}m)\end{array}$$

证明2（已知 $a\equiv b(\mathrm{mod}m)$）：

$$\begin{array}{rl}\because & a\equiv b(\mathrm{mod}m)\\ \therefore & a\mathrm{mod}m=b\mathrm{mod}m\\ & \text{令}mx+k=a,my+k=b(x,y,k\in \mathbb{Z})\\ & \begin{array}{rl}a-b& =mx+k-my-k\\ & =m(x-y)\end{array}\\ \because & m\mid m(x-y)\\ \therefore & m\mid (a-b)\end{array}$$

3.2. 性质

3.2.1. 性质1

$a\equiv a(\mathrm{mod}m)$

3.2.2. 性质2

如果 $a\equiv b(\mathrm{mod}m)$ ，那么 $b\equiv a(\mathrm{mod}m)$

3.2.3. 性质3

如果 $a\equiv b(\mathrm{mod}m)$ ， $b\equiv c(\mathrm{mod}m)$ ，那么 $a\equiv c(\mathrm{mod}m)$

3.2.4. 性质4

如果 $a\equiv b(\mathrm{mod}m)$ ，那么 $(a+c)\equiv (b+c)(\mathrm{mod}m)$

证明：

$$\begin{array}{rl}\because & a\equiv b(\mathrm{mod}m)\\ & \text{令}mk=a-b(k\in \mathtt{Z})\\ & \begin{array}{rl}mk& =a-b\\ & =a+c-b-c\\ & =(a+c)-(b+c)\end{array}\\ \therefore & (a+c)\equiv (b+c)(\mathrm{mod}m)\end{array}$$

3.2.5. 性质5

如果 $a\equiv b(\mathrm{mod}m)$ ，那么 $(a-c)\equiv (b-c)(\mathrm{mod}m)$

证明方法类似于性质3

3.2.6. 性质6

如果 $a\equiv b(\mathrm{mod}m)$ ，那么 $(ac)\equiv (bc)(\mathrm{mod}m)$

证明：

$$\begin{array}{rl}\because & a\equiv b(\mathrm{mod}m)\\ & \text{令}mk=a-b(k\in \mathtt{Z})\\ \therefore & mk=a-b⟺mkc=ac-bc\\ \therefore & (ac)\equiv (bc)(\mathrm{mod}m)\end{array}$$

3.2.7. 性质7

如果 $a\equiv b(\mathrm{mod}m),c\mid a,c\mid b,(c,m)=1$ ，那么 $(\frac{a}{c})\equiv (\frac{b}{c})(\mathrm{mod}m)$

多说一句： $(a,b)=\gcd (a,b)$

证明：

$$\begin{array}{rlr}\because & a\equiv b(\mathrm{mod}m)& \\ & \text{令}mk=a-b(k\in \mathtt{Z})& \\ \therefore & mk=a-b⟺m\frac{k}{c}& =\frac{a}{c}-\frac{b}{c}\\ \because & c\mid a,c\mid b& \\ \therefore & \frac{a}{c},\frac{b}{c}\in \mathbb{Z}& \\ \therefore & \frac{mk}{c}\in \mathbb{Z}& \\ \because & (c,m)=1& \\ \therefore & \frac{m}{c}\notin \mathbb{Z}& \\ \therefore & \frac{k}{c}\in \mathbb{Z}& \\ \therefore & \frac{a}{c}\equiv \frac{b}{c}(\mathrm{mod}m)& \end{array}$$

3.2.8. 性质8

如果 $a\equiv b(\mathrm{mod}m)$ ，那么 $(a^c)\equiv (b^c)(\mathrm{mod}m)$

证明：

$$\begin{array}{rl}\because & a\equiv b(\mathrm{mod}m)\\ & \text{令}mk=a-b(k\in \mathtt{Z})\\ & \\ & \begin{array}{rl}mk& =a-b\\ mk\times \sum _{i=0}^{c-1}{a}^i{b}^{c-i-1}& =(a-b)\times \sum _{i=0}^{c-1}{a}^i{b}^{c-i-1}\\ mk\times \sum _{i=0}^{c-1}{a}^i{b}^{c-i-1}& =a^c-b^c\end{array}\\ & \\ \therefore & (a^c)\equiv (b^c)(\mathrm{mod}m)\end{array}$$

另注： $a^n-b^n=(a-b)\times \sum _{i=0}^{n-1}{a}^i{b}^{n-i-1}$

3.2.9. 性质9

如果 $a\equiv b(\mathrm{mod}m),c\equiv d(\mathrm{mod}m)$ ，那么 $(a+c)\equiv (b+d)(\mathrm{mod}m)$

证明：

$$\begin{array}{rl}\because & a\equiv b(\mathrm{mod}m),c\equiv d(\mathrm{mod}m)\\ & \text{令}mx=a-b,my=c-d(x,y\in \mathtt{Z})\\ & \begin{array}{rl}mx+my& =a-b+c-d\\ m(x+y)& =a+c-(b+d)\end{array}\\ \therefore & (a+c)\equiv (b+d)(\mathrm{mod}m)\end{array}$$

3.2.10. 性质10

如果 $a\equiv b(\mathrm{mod}m),c\equiv d(\mathrm{mod}m)$ ，那么 $(ac)\equiv (bd)(\mathrm{mod}m)$

证明：

$$\begin{array}{rl}\because & a\equiv b(\mathrm{mod}m),c\equiv d(\mathrm{mod}m)\\ & \text{令}mx=a-b,my=c-d(x,y\in \mathtt{Z})\\ \because & mx=a-b,my=c-d\\ \therefore & mxc=ac-bc,myb=bc-cd\\ \therefore & mxc+myb=ac-bc+bc-cd⟺m(xc+yb)=ac-bd\\ \therefore & (ac)\equiv (bd)(\mathrm{mod}m)\end{array}$$

3.2.11. 性质11

如果 $a\mathrm{mod}p=x,a\mathrm{mod}q=x,(p,q)=1$ ，那么 $a\mathrm{mod}(pq)=x$

证明：

$$\begin{array}{rl}\because & a\mathrm{mod}p=x,a\mathrm{mod}q=x\\ & \text{令}p{k}_1+x=a,q{k}_2+x=a(k_1,k_2\in \mathbb{Z})\\ \therefore & p{k}_1+x=q{k}_2+x\\ \therefore & \frac{p{k}_1}{q}=k_2\\ \because & k_2\in \mathbb{Z}\\ \therefore & \frac{p{k}_1}{q}\in \mathbb{Z}\\ \because & (p,q)=1\\ \therefore & \frac{p}{q}\notin \mathbb{Z}\\ \therefore & \frac{k_1}{q}\in \mathbb{Z}\\ \therefore & q{k}_3=k_1\\ \therefore & p{k}_1+x=a⟺pq{k}_3+x=a\\ \therefore & a\mathrm{mod}(pq)=x\end{array}$$

4. 排列组合

咕咕咕

5. 筛法

5.1. 暴力筛

暴力筛，顾名思义，非常暴力的筛法（？

暴力筛是所有筛法的祖先，需要好好掌握

我们知道，一个数的倍数是一个合数，那么，我们可以进行一个遍历，将每一个数的倍数全部标记一遍，剩下的未被标记的就是质数了

#include<cmath>
#include<cstdio>
bool flag[200000005];
int n,ans,maxn;			//ans 用于统计质数数量, maxn 用于记录最大质数
int main(){
	scanf("%d",&n);
	for(int i=2;i<=sqrt(n);i++){			//只需枚举到 sqrt(n) 即可
											//原理很简单，假设 n=pq(p<q) ，那么，我们一定会在 p 之前的 q 枚举，而 p 的最大值又是 sqrt(n) ，所以枚举到 sqrt(n) 即可
		if(!flag[i]){			//未被标记，说明是质数
			ans++,maxn=i;			//更新数量和最大值
		}
		for(int j=i+i;j<=n;j+=i){			//将其倍数设为合数
			flag[j]=1;
		}
	}
	for(int i=sqrt(n)+1;i<=n;i++){			//因为遍历到 sqrt(n) 就结束了，后面的素数还要再来一个循环统计一次
		if(!flag[i]){
			ans++;
			maxn=i;
		}
	}
	printf("%d %d",ans,maxn);
	return 0;
}

然而，这显然不是最优筛法

5.2. 埃式筛

我们发现，假设一个合数 $n$ 的其中一个质因子为 $k$ ，那么 $n=kx$ ， $n$ 与其倍数必然已经在遍历到 $k$ 的时候标记了，就不用再筛一遍了

所以，我们可以只用筛素数的倍数即可

#include<cmath>
#include<cstdio>
bool flag[200000005];
int n,ans,maxn;
int main(){
	scanf("%d",&n);
	for(int i=2;i<=sqrt(n);i++){
		if(!flag[i]){
			ans++;
			maxn=i;
		}
		if(flag[i]){			//当遍历到素数时，直接跳过，不进行遍历，剩下的和暴力筛一模一样
			continue;
		}
		for(int j=i+i;j<=n;j+=i){
			flag[j]=1;
		}
	}
	for(int i=sqrt(n)+1;i<=n;i++){
		if(!flag[i]){
			ans++;
			maxn=i;
		}
	}
	printf("%d %d",ans,maxn);
	return 0;
}

埃式筛已经够快了，但是，我们还有一个更快的解决方法

5.3. 欧拉筛

在埃式筛中，某些数依旧会被重复标记两次，比如 $12$ ，既会被 $3$ 标记，也会被 $2$ 标记

这时，我们尽量让每个数只被标记一次

而欧拉筛的原理是让每一个数被它的最小质因数标记

这里先上代码

#include<cstdio>
bool flag[200000005];
int n,ans,maxn;
int num[1500005],len;			//num 数组用于统计素数数组
int main(){
	scanf("%d",&n);
	for(int i=2;i<=n;i++){
		if(!flag[i]){			//遍历到素数
			num[++len]=i,ans++,maxn=i;
		}
		for(int j=1;j<=len&&num[j]*i<=n;j++){			//枚举可能的最小质因数
			flag[num[j]*i]=1;			//标记
			if(i%num[j]==0){			//后文会有详细解释
				break;
			}
		}
	}
	printf("%d %d",ans,maxn);
	return 0;
}

这里解释一下这一段代码：

if(i%num[j]==0){
	break;
}

$$\begin{array}{rl}& \text{假设}num[j]=k\\ \because & k\mid i\\ \therefore & i=kx(x\in {\mathbb{Z}}^{+})\\ & \text{当枚举到}num[j+1]=k_1\text{时}\\ & i{k}_1=kx{k}_1\\ \because & k<k_1\\ \therefore & k_1\text{不是最小质因数}\end{array}$$

显然，这样已经违背了欧拉筛的原理，所以退出

6. gcd/lcm

6.1. 定义

公约数定义：一般设有 $k$ 个数 $a_1,a_2,a_3,\cdots ,a_k\in {\mathbb{Z}}^{+}$ ，如果 $\exists d\in {\mathbb{Z}}^{+}$ ，使 $d\mid a_1,d\mid a_2,d\mid a_3,\cdots ,d\mid a_k$ ，那么 $d$ 为 $a_1,a_2,a_3,\cdots ,a_k$ 的公约数， $d$ 是不唯一的

最大公约数：在 $a_1,a_2,a_3,\cdots ,a_k$ 的公约数集合中，最大的那个公约数称为最大公约数，如果是 $a,b$ 两个元素的最大公约数，记为 $\gcd (a,b)$ ，简写为 $(a,b)$

公因数与最小公因数的定义与公约数与最大公约数的定义类似，注意一下符号即可：

$a,b$ 两个元素的最小公因数记为 $\mathrm{lcm}(a,b)$ ，简写为 $[a,b]$

6.2. 一个定理

$\gcd (a,b)\times \mathrm{lcm}(a,b)=ab$

证明：

将 a , b 进行质因子分解，设其质因子集合并集为 { p 1 , p 2 , ⋯   , p n } 设 a = p 1 k 1 p 1 k 2 ⋯ p n k n , b = p 1 j 1 p 1 j 2 ⋯ p n j n ( j i , k i ∈ Z , j i ≥ 0 , k i ≥ 0 ) ∵   gcd ⁡ ( a , b ) = p 1 min ⁡ ( k 1 , j 1 ) p 2 min ⁡ ( k 2 , j 2 ) ⋯ p n min ⁡ ( k n , j n ) lcm ⁡ ( a , b ) = p 1 max ⁡ ( k 1 , j 1 ) p 2 max ⁡ ( k 2 , j 2 ) ⋯ p n max ⁡ ( k n , j n ) min ⁡ ( a , b ) + lcm ⁡ ( a , b ) = a + b   ∴ gcd ⁡ ( a , b ) × lcm ⁡ ( a , b ) = p 1 min ⁡ ( k 1 , j 1 ) p 2 min ⁡ ( k 2 , j 2 ) ⋯ p n min ⁡ ( k n , j n ) p 1 max ⁡ ( k 1 , j 1 ) p 2 max ⁡ ( k 2 , j 2 ) ⋯ p n max ⁡ ( k n , j n ) = p 1 min ⁡ ( k 1 , j 1 ) + max ⁡ ( k 1 , j 1 ) p 2 min ⁡ ( k 2 , j 2 ) + max ⁡ ( k 2 , j 2 ) ⋯ p n min ⁡ ( k n , j n ) + max ⁡ ( k n , j n ) = p 1 k 1 + j 1 p 2 k 2 + j 2 ⋯ p n k n + j n = p 1 k 1 p 1 k 2 ⋯ p n k n p 1 j 1 p 1 j 2 ⋯ p n j n = a b \begin{aligned}&\text{将}a,b\text{进行质因子分解，设其质因子集合并集为}\{p_1,p_2,\cdots,p_n\} \\&\text{设}a=p_1^{k_1}p_1^{k_2}\cdots p_n^{k_n},b=p_1^{j_1}p_1^{j_2}\cdots p_n^{j_n}(j_i,k_i\in\mathbb{Z},j_i\ge0,k_i\ge0)\\&\because\ \gcd(a,b)=p_1^{\min(k_1,j_1)}p_2^{\min(k_2,j_2)}\cdots p_n^{\min(k_n,j_n)}\\&\operatorname{lcm}(a,b)=p_1^{\max(k_1,j_1)}p_2^{\max(k_2,j_2)}\cdots p_n^{\max(k_n,j_n)}\\&\min(a,b)+\operatorname{lcm}(a,b)=a+b\\&\ \begin{aligned}\therefore\gcd(a,b)\times\operatorname{lcm}(a,b)&=p_1^{\min(k_1,j_1)}p_2^{\min(k_2,j_2)}\cdots p_n^{\min(k_n,j_n)}p_1^{\max(k_1,j_1)}p_2^{\max(k_2,j_2)}\cdots p_n^{\max(k_n,j_n)}\\&=p_1^{\min(k_1,j_1)+\max(k_1,j_1)}p_2^{\min(k_2,j_2)+\max(k_2,j_2)}\cdots p_n^{\min(k_n,j_n)+\max(k_n,j_n)}\\&=p_1^{k_1+j_1}p_2^{k_2+j_2}\cdots p_n^{k_n+j_n}\\&=p_1^{k_1}p_1^{k_2}\cdots p_n^{k_n}p_1^{j_1}p_1^{j_2}\cdots p_n^{j_n}\\&=ab\end{aligned}\end{aligned} ​将a,b进行质因子分解，设其质因子集合并集为{p1​,p2​,⋯,pn​}设a=p1k1​​p1k2​​⋯pnkn​​,b=p1j1​​p1j2​​⋯pnjn​​(ji​,ki​∈Z,ji​≥0,ki​≥0)∵ gcd(a,b)=p1min(k1​,j1​)​p2min(k2​,j2​)​⋯pnmin(kn​,jn​)​lcm(a,b)=p1max(k1​,j1​)​p2max(k2​,j2​)​⋯pnmax(kn​,jn​)​min(a,b)+lcm(a,b)=a+b ∴gcd(a,b)×lcm(a,b)​=p1min(k1​,j1​)​p2min(k2​,j2​)​⋯pnmin(kn​,jn​)​p1max(k1​,j1​)​p2max(k2​,j2​)​⋯pnmax(kn​,jn​)​=p1min(k1​,j1​)+max(k1​,j1​)​p2min(k2​,j2​)+max(k2​,j2​)​⋯pnmin(kn​,jn​)+max(kn​,jn​)​=p1k1​+j1​​p2k2​+j2​​⋯pnkn​+jn​​=p1k1​​p1k2​​⋯pnkn​​p1j1​​p1j2​​⋯pnjn​​=ab​​

6.3. 欧几里得算法

一种非常常见的求最大公约数的算法：

$\gcd (a,b)=\gcd (b,a\mathrm{mod}b)$

代码形式：

#include<cstdio>
int gcd(int a,int b){
	return !b?a:gcd(b,a%b);
}
int main(){
	int n,m;
	scanf("%d%d",&n,&m);
	printf("%d",gcd(n,m));
	return 0;
}
//非常简单，应该不用注释了吧

证明一下：

$$\begin{array}{rl}& \text{设}\gcd (a,b)=p\\ \therefore & a=kp,b=gp(\gcd (k,g)=1)\\ \therefore & a\mathrm{mod}b=(kp)\mathrm{mod}(gp)=(k\mathrm{mod}g)p\\ & \begin{array}{rl}\gcd (b,a\mathrm{mod}b)& =\gcd (gp,(k\mathrm{mod}g)p)\\ & =p\\ & =\gcd (a,b)\end{array}\end{array}$$

7. 裴蜀定理

这是一个非常常见的定理，相当重要

7.1. 基本定理

对于不定方程 $ax+by=m(a,b,x,y,m\in \mathbb{Z})$ ，其有解的充分必要条件为 $\gcd (a,b)\mid m$

证明1（已知 $ax+by=m$ 有解）：

$$\begin{array}{rl}\because & ax+by=m\text{有解}\\ \therefore & \exists x,y\in \mathbb{Z},\text{满足}ax+by=m\\ \because & \gcd (a,b)\mid a,\gcd (a,b)\mid b\\ \therefore & \gcd (a,b)\mid (ax+by)(\text{整除性质2})\\ \therefore & \gcd (a,b)\mid m\end{array}$$

证明2（已知 $\gcd (a,b)\mid m$）：

咕咕咕

OI 中常见的有关裴蜀定理的出题形式（？

对于 $\forall x,y\in \mathbb{Z}$ ，函数 $f(a,b)=ax+by$ 的最小 ${\mathbb{Z}}^{+}$ 的取值为 $\gcd (a,b)$

证明：

$$\begin{array}{rl}& \text{设}m=ax+by\\ & \text{由裴蜀定理得：}\gcd (a,b)k=m(k\in {\mathbb{Z}}^{+})\\ \because & k\text{的最小值为}1,\gcd (a,b)\text{为定值}\\ \therefore & m\text{的最小值为}\gcd (a,b)\end{array}$$

7.2. 引理及推论

7.2.1. 引理1

上台社死后我发现，我换一个表述方式似乎就没什么大问题了QAQ

如果 $a,b\in {\mathbb{Z}}^{+},(a,b)=1$ ，则 $\nexists$ 小于 $b$ 的 $k\in {\mathbb{Z}}^{+}$ ，满足 $b\mid ak$

证明：

$$\begin{array}{rl}& \text{假设存在小于}b\text{的}k\in \mathbb{Z},\text{满足}b\mid ak\\ \therefore & ak\text{中包含}b\text{的全部质因子}\\ \because & (a,b)=1\\ \therefore & a\text{不可能包含}b\text{的任何质因子}\\ \therefore & k\text{中包含}b\text{的全部质因子}\\ \therefore & k=bx(x\in {\mathbb{Z}}^{+})\\ \because & b>k\\ \therefore & \text{矛盾，假设不成立}\end{array}$$

7.2.2. 推论

如果 $a,b\in {\mathbb{Z}}^{+},(a,b)=1$ ，则 $0,a,2a,\cdots ,(b-1)a$ 这些数模 $b$ 的余数互不相同

证明：

$$\begin{array}{rl}& \text{假设}ai,aj\text{模}b\text{的余数相同}(1<j<i<b)\text{,即}ai\mathrm{mod}b=aj\mathrm{mod}b\\ \therefore & ai\mathrm{mod}b-aj\mathrm{mod}b=0⟺(ai-aj)\mathrm{mod}b=0⟺a(i-j)\mathrm{mod}b=0\\ \therefore & b\mid a(i-j)(i-j<b)\\ \because & a,b\in {\mathbb{Z}}^{+},(a,b)=1\\ \therefore & \nexists \text{小于}b\text{的}k\in {\mathbb{Z}}^{+}\text{，满足}b\mid ak\\ \therefore & \text{矛盾，假设不成立}\end{array}$$

7.2.3. 引理2

如果 $(a,b)=1$ ，则 $\exists k\in \mathbb{Z}$ ，满足 $ak\mathrm{mod}b=1$

由推论可以很容易的得到该结论

8. 扩展欧几里得

扩展欧几里得（简称扩欧），是一种算法，一般封装为函数，其函数名为 exgcd

可以在已知 $a,b\in \mathbb{Z}$ 的情况下，求不定方程 $ax+by=\gcd (a,b)$ 的一组整数解，通常情况下用扩欧求得的解满足 $xy<0$

证明：

$$\begin{array}{rl}& \text{设}a{x}_1+b{y}_1=\gcd (a,b),b{x}_2+(a\mathrm{mod}b)y_2=\gcd (b,a\mathrm{mod}b)\\ & \text{由欧几里得算法得}\gcd (a,b)=\gcd (b,a\mathrm{mod}b)\\ \therefore & a{x}_1+b{y}_1=b{x}_2+(a\mathrm{mod}b)y_2\\ \because & a\mathrm{mod}b=a-b\lfloor \frac{a}{b}\rfloor \\ & \begin{array}{rl}a{x}_1+b{y}_1& =b{x}_2+(a\mathrm{mod}b)y_2\\ & =b{x}_2+(a-b\lfloor \frac{a}{b}\rfloor )y_2\\ & =b{x}_2+a{y}_2-b\lfloor \frac{a}{b}\rfloor y_2\\ & =a{y}_2+b(x_2-b\lfloor \frac{a}{b}\rfloor )\end{array}\\ \therefore & \{\begin{array}{l}{x}_1=y_2\\ y_1=x_2-b\lfloor \frac{a}{b}\rfloor \end{array}\text{是方程}a{x}_1+b{y}_1=\gcd (a,b)\text{的一组解，且该组解由}(x_2,y_2)\text{得来}\\ & \text{当第}n\text{组方程关于}\gcd (a,0)\text{时，有}a{x}_n+b{y}_n=a\\ \therefore & \{\begin{array}{l}{x}_n=1\\ y_n=\forall k\in \mathbb{Z}\end{array}\text{是方程}a{x}_n+b{y}_n=a\text{的一组解}\\ \therefore & \text{可以由解}(x_n,y_n)\text{得到解}(x_1,y_1)\end{array}$$

我们拿这道题来举个例子

8.1. 小练习

我们来处理一下原方程：

$$\begin{array}{rl}ax& \equiv 1(\mathrm{mod}b)\\ by+1& =ax(y\in \mathbb{Z})\\ ax-by& =1\end{array}$$

这样一来，直接套扩欧，再处理一下取值，基本上就没问题了

不 开 long long 见 祖 宗

#include<cstdio>
long long int a,b,c,x,y;
void exgcd(long long int a,long long int b){			//扩欧模板
	if(!b){
		x=1,y=0;
		return;
	}
	exgcd(b,a%b);
	long long int t=x;
	x=y;
	y=t-a/b*y;
}
int main(){
	scanf("%lld%lld",&a,&b);
	exgcd(a,b);
	if(x<0){			//处理结果，保证是最小的正整数
		x=x+(-x/b+1)*b;
	}else{
		x=x-x/b*b;
	}
	printf("%lld",x);			//输出
	return 0;
}
//水

9. 费马小定理

如果 $p$ 为质数，且 $a\mathrm{mod}p\ne 0$ ，则 $a^{p-1}\mathrm{mod}p=1$

证明：

$$\begin{array}{rl}\because & a(2a)(3a)\cdots [(p-1)a]\mathrm{mod}p=a^{p-1}\times (p-1)!\mathrm{mod}p\\ \text{又}\because & p\text{为质数，且}a\mathrm{mod}p\ne 0\\ \therefore & \text{由裴蜀定理推论得：}a\mathrm{mod}p,(2a)\mathrm{mod}p,(3a)\mathrm{mod}p,\cdots ,[(p-1)a]\mathrm{mod}p\text{互不相等，且在区间}(0,p)\text{内}\\ \therefore & a(2a)(3a)\cdots [(p-1)a]\mathrm{mod}p=(p-1)!\mathrm{mod}p\\ \therefore & a^{p-1}(p-1)!\equiv (p-1)!(\mathrm{mod}p)\\ \because & p\text{为质数}\\ \therefore & ((p-1)!,p)=1\\ \therefore & a^{p-1}\equiv 1(\mathrm{mod}p)\end{array}$$