欧几里得算法和扩展欧几里得算法——杨子曰数学

欧几里得算法和扩展欧几里得算法——杨子曰数学

超链接：数学合集

---

不说废话，咱们直接开始

欧几里得算法

一句话：
$$gcd(a,b)=gcd(b,amodb)$$
简单证明：

首先我们先来证一证这个：$gcd(a,b)=gcd(b,a-b)$

假设：$d=gcd(a,b),a=k_a\ast d,b=k_b\ast d（k_a和k_b互质）$

这样一来，我们会发现：
$$a-b=k_a\ast d-k_b\ast d=(k_a-k_b)\ast d$$
欧，有没有发现a-b和b还是都有公因子d

$a-b=(k_a-k_b)\ast d$

$b=k_b\ast d$

BUT，我们还要证明d是它们的最大公因子，也就是要证明，它们的系数$k_a-k_b$和$k_b$的互质（它们如果不互质，那么最大公因数就会变成$d\ast gcd(k_a-k_b,k_b)$）

我们用反证法：

假设$k_a-k_b$和$k_b$不互质

$$设:gcd(k_a-k_b,k_b)=\lambda （\lambda >1），k_a-k_b=\alpha \ast \lambda ，k_b=\beta \ast \lambda$$
于是
$$k_a=k_a-k_b+k_b=\alpha \ast \lambda +\beta \ast \lambda =(\alpha +\beta )\ast \lambda$$

有没有发现$k_a$和$k_b$有了公因数$\lambda$，别忘了上面还有一个条件叫$k_a和k_b互质$

于是得证！

然后我们就得到了：$gcd(a,b)=gcd(b,a-b)$

然后，我们把a-b看成等式左边新的b，然后一直拿a减，直到不能减了为止，其实就是:
$$gcd(a,b)=gcd(b,amodb)$$

得证！

---

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

---

扩展欧几里得算法

这个东西只有一个作用
$$求ax+by=gcd(a,b)的整数解$$

假设存在一组x’和y’使得：
$$b{x}^{\prime }+(amodb)y^{\prime }=gcd(b,amodb)$$

根据上面的欧几里得算法：
$$gcd(a,b)=gcd(b,amodb)$$

我们就得到了：
$$ax+by=b{x}^{\prime }+(amodb)y^{\prime }$$

幼儿园小盆友都知道：
$$amodb=a-\lfloor \frac{a}{b}\rfloor \ast b$$

然后把它带回去：
$$ax+by=b{x}^{\prime }+(a-\lfloor \frac{a}{b}\rfloor \ast b)y^{\prime }$$

把括号拆开：
$$ax+by=b{x}^{\prime }+a{y}^{\prime }-\lfloor \frac{a}{b}\rfloor \ast b{y}^{\prime }$$

再合并一下：
$$ax+by=a{y}^{\prime }+b(x^{\prime }-\lfloor \frac{a}{b}\rfloor \ast y^{\prime })$$

So，你有没有惊奇的发现一定有一组解长这样：
$$\{\begin{array}{rl}x& =y^{\prime }\\ y& =x^{\prime }-\lfloor \frac{a}{b}\rfloor \ast y^{\prime }\end{array}$$

也就是说我们只要求出x’和y’就可以求出x和y了

我们再来回头看一下它们的方程：
$$ax+by=gcd(a,b)$$

$$b{x}^{\prime }+(amodb)y^{\prime }=gcd(b,amodb)$$
有没有发现形式是一样滴，So，我们可以通过递归来解决它

那么临界情况是什么捏？

一个幼儿园小盆友都知道的事实：
$$gcd(a,0)=a$$
So，当b=0时：
$$ax+0=gcd(a,0)$$
然后就可以，美滋滋地让x=1，y是任意值了（一般是0，它会影响最后的解的大小）

OK，完事

void ex_gcd(int a,int b,int &x,int &y){
    if(!b) x=1,y=0;
    else ex_gcd(b,a%b,y,x),y-=(a/b)*x;
}

应用1：解方程ax+by=c

首先，在上面已经讲过了，扩展欧几里得可以求长这样的方程：
$$ax+by=gcd(a,b)$$
那有什么卵用！我从来没有在哪道题里见长成这样的方程

BUT，我们可以利用它来解长成这样的方程：
$$ax+by=c$$
让$d=gcd(a,b)$
这时，如果我们把这个方程同时乘上一个$\frac{d}{c}$

方程就变成了：
$$a\frac{d}{c}\ast x+b\frac{d}{c}\ast y=gcd(a,b)$$

我们就可以来解一个这样的方程：

$$a{x}^{\prime }+b{y}^{\prime }=gcd(a,b)$$

于是你就会发现，原方程的解就是：
$$\{\begin{array}{rl}x& =x^{\prime }\ast \frac{c}{d}\\ y& =x^{\prime }\ast \frac{c}{d}\end{array}$$

你一定会问了c如果不是d的倍数肿么办？

那就没有整数解呗！这有什么好问的？

应用2：解模线性方程：$ax\equiv b(modp)$

其实就是应用1，有没有发现$ax\equiv b(modp)$，就是：$ax+py=b$，然后再用上面的那个东东解出来就好了
解模线性方程，就可以做很多事情了，比如：求逆元，戳我

参考：
https://www.zybuluo.com/samzhang/note/541890
于HG机房