欧几里德算法

欧几里德算法

本篇文章部分参考于 百度百科和丶Sj 博客

[https://blog.csdn.net/qq_40475868/article/details/80277705]

最大公约数(GCD)

算法原理

$$gcd(a,b)=gcd(b,a\%b)$$

证明

设$a>b$，$k=gcd(a,b)$表示$a$与$b$的最大公约数，有：
$$\{\begin{array}{l}a=ck\\ b=dk\end{array}(gcd(c,d)=1)$$设$p$为正整数，且满足：
$$r=a\%b=a-pb=(c-pd)k$$$c-pd与d$互质，若不然，则存在$m>1$，使得：
$$\{\begin{array}{l}c-pd=xm\\ d=ym\end{array}(gcd(x,y)=1)$$于是有：
$$\{\begin{array}{l}a=ck=(pd+xm)k=(py+x)mk\\ b=dk=ymk\end{array}$$此时$gcd(a,b)=mk$，与原题矛盾，故$c-pd与d$互质。
所以$r与b$的最大公约数也是$k$：
$$gcd(a,b)=gcd(b,r)=gcd(b,a\%b)=k$$

int gcd(int a,int b){
	return b ? gcd(b,a%b) : a;
}

最小公倍数（LCM）

证明

设$a,b$为整数，有：
$$\{\begin{array}{l}a=c\cdot gcd(a,b)\\ b=d\cdot gcd(a,b)\end{array}(gcd(c,d)=1)$$很明显最小公倍数为:$$lcm(a,b)=\frac{a}{gcd(a,b)}\cdot b$$

代码实现

int lcm(int a,int b){
	return a/gcd(a,b)*b;
}

求二元一次方程的整数特解

1. 关于$ax+by=1(gcd(a,b)=1)$方程总有整数解
   对方程$ax+by=k$两边同时取 a 模：
   $$by\equiv k(moda)$$对于 a 的一个完全剩余系 { y } ( y ∈ [ 1 , a ] ) \{y\}(y\in[1,a]) {y}(y∈[1,a])，由于 a 与 b 互质，由费马小定理中的引理2，知 { b y } ( y ∈ [ 1 , a ] ) \{by\}(y\in[1,a]) {by}(y∈[1,a])也是模 a 中的一个完全剩余系，故总存在一个 $y_1\in [1,a]s.t.$
   $$b{y}_1\equiv 1(moda)$$所以 $b{y}_1=ka+1$，所以：
   $$ax+by=a(-k)+b{y}_1=1$$即，$ax+by=1$总有整数解，且显然 $xy<0$
2. 关于$ax+by=1(gcd(a,b)=1)$的整数解的无限性
   若$(x_1,y_1)(x_2,y_2)$是$ax+by=1$的解，那么总有：
   $$a(x_2-x_1)=b(y_1-y_2)$$由于$gcd(a,b)=1$，所以：
   $$x_2-x_1=kb,y_1-y_2=ka(k是整数)$$所以：
   $$\{\begin{array}{l}{x}_2=x_1+kb\\ y_2=y_1-ka\end{array}$$由于 k 是可以任意取的，所以原方程的解有无限个。
3. 关于$ax+by=1(gcd(a,b)=1)$的整数特解
   设$r=a\%b=a-\frac{a}{b}$，有：
   $$ax+by=1\Rightarrow (r+\frac{a}{b}\cdot b)x+by=1\Rightarrow b(\frac{a}{b}\cdot x+y)+rx=1$$令
   $$\{\begin{array}{l}{x}^{\prime }=\frac{a}{b}\cdot x+y\\ y^{\prime }=x\end{array}$$有：
   $$\{\begin{array}{l}x=y^{\prime }\\ y=x^{\prime }-\frac{a}{b}{y}^{\prime }\end{array}$$关于 b 与 r 互质的事实已在上最大公约数算法原理中证实，于是可以发现这个问题被表达为另一个同类型的子问题，这可以写一个递归，结束条件为：
   $$r=0时\{\begin{array}{l}x=1\\ y=0\end{array}$$
   代码实现（顺带求gcd）

int gcd(int a,int b,int&x,int&y){
	if(!b){
		x=1;y=0;
		return a;
	}
	int ans = gcd(b,a%b,y,x);//a%b就是r,y和x换位是为了直接让x=y`而y=x`，之后再减一个a/b*x即可
	y-=a/b*x;
	return ans;
}

4. 其他二元一次方程及其整数解
   可以发现上面算法其实与方程等号右边是否为1无关，事实上，它同样适用于任意的：
   $$ax+by=gcd(a,b)$$两边同除$gcd(a,b)$即可变回原来分析的方程，而保持解$(x,y)$不变。
   而对于任意方程：
   $$ax+by=c$$其有整数解当且仅当$c\%gcd(a,b)=0$时（证明略），很明显，其解为：
   { x = k ⋅ x 0 y = k ⋅ y 0 { ( x 0 , y 0 ) ∣ a x 0 + b y 0 = g c d ( a , b ) } ( k = c g c d ( a , b ) ) \begin{cases}x=k\cdot{x_{0}}\\y=k\cdot{y_{0}}\end{cases}\{(x_{0},y_{0})|ax_{0}+by_{0}=gcd(a,b)\}(k=\frac{c}{gcd(a,b)}) {x=k⋅x0​y=k⋅y0​​{(x0​,y0​)∣ax0​+by0​=gcd(a,b)}(k=gcd(a,b)c​)
5. 求二元一次方程的最小非负整数解
   先对于$ax+by=1(gcd(a,b)=1)$，用扩展欧几里德算法得到的解为$(x_0,y_0)$，那么显然最小正整数解 $x$ 可以表示为：
   $$x=\{\begin{array}{l}{x}_0\%b+b,x\le 0\\ x_0\%b,x>0\end{array}=(x_0\%b+b)\%b$$类似的对于$ax+by=gcd(a,b)$：
   $$x=(x_0\%q+q)\%q(q=\frac{b}{gcd(a,b)})$$
   而对于$ax+by=c(c\%gcd(a,b)=0)$：
   $$x=(x_1\%q+q)\%q(q=\frac{b}{gcd(a,b)})(x_1=x_0\cdot \frac{c}{gcd(a,b)})$$相应的$y$带入原方程即可解得。