数论(一) 欧几里得算法与扩展欧几里得算法

1 欧几里得算法（Euclidean algorithm）

1.1 概述

欧几里得算法，又称辗转相除法，是用于求两个不全等于0的非负整数的最大公约数的算法。
一般用 $gcd(a,b)$ 表示 $a$ 与 $b$ 的最大公约数，其中 gcd 是 Greatest Common Divisor 的缩写。

欧几里得算法的核心思想可以归结为一个简短而又优雅的式子——$gcd(a,b)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$。(默认$a\ge b$)

通过这样的一个式子，我们可以把原问题的规模不断缩小，即一直递归直到到达 base case (即 $b=0$ )为止。当 $b=0$ 时显然有 $gcd(a,0)=a$，此时 base case 中的 $a$ 值即为原问题中的 $a$ 与 $b$ 的最大公约数。

下面证明为何 $gcd(a,b)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$ 成立。

1.2 正确性证明

设 $g_0=gcd(a,b)$，则 $g_0\mid a$，$g_0\mid b$。
设 $g_1=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$，则 $g_1\mid b$，$g_1\mid (a\mathrm{m}\mathrm{o}\mathrm{d}b)$。

因为 $a=\lfloor a/b\rfloor b+(a\mathrm{m}\mathrm{o}\mathrm{d}b)$，而 $g_1$ 同时整除 $b$ 和 $(a\mathrm{m}\mathrm{o}\mathrm{d}b)$，故有 $g_1\mid a$。
由 $g_1\mid a$ 和 $g_1\mid b$ 以及 $g_0$ 是 $a$ 和 $b$ 的最大公约数的事实，可得知 $g_1\mid g_0$。

因为 $a\mathrm{m}\mathrm{o}\mathrm{d}b=a-\lfloor a/b\rfloor b$，而 $g_0$ 同时整除 $a$ 和 $b$，故有 $g_0\mid (a\mathrm{m}\mathrm{o}\mathrm{d}b)$。
由 $g_0\mid b$ 和 $g_0\mid (a\mathrm{m}\mathrm{o}\mathrm{d}b)$ 以及 $g_1$ 是 $b$ 和 $(a\mathrm{m}\mathrm{o}\mathrm{d}b)$ 的最大公约数的事实，可得知 $g_0\mid g_1$。

因为 $g_1\mid g_0$ 及 $g_0\mid g_1$，故有 $g_0=g_1$。

证毕。

1.3 复杂度分析

随便 google 一下就能找到不少借助 Fibonacci 数列得出复杂度为 $\mathcal{O}(\log n)$ 的帖子。

这里介绍另一种稍微容易一点的分析方法。

设 $r_k$ 为进行了 $k$ 次取模后的余数，并用 $r_0$ 表示 $b$，则有：
$$\begin{gathered} a\mathrm{m}\mathrm{o}\mathrm{d}{r}_0=r_1 \\ {r}_0\mathrm{m}\mathrm{o}\mathrm{d}{r}_1=r_2 \\ ⋮ \\ {r}_k\mathrm{m}\mathrm{o}\mathrm{d}{r}_{k+1}=r_{k+2} \\ ⋮ \end{gathered}$$

观察若干步骤后会发现 —— $r_{k+2}\le 1/2r_k$，即递归过程每加深两层，问题规模必小于等于原来的一半。

下面给出证明。

1. 若 $r_{k+1}\ge 1/2r_k$，则有 $r_k-r_{k+1}\le 1/2r_k$，故 $r_{k+2}=r_k\mathrm{m}\mathrm{o}\mathrm{d}{r}_{k+1}\le r_k-r_{k+1}\le 1/2r_k$。
2. 若 $r_{k+1}<1/2r_k$，则有 $r_{k+2}=r_k\mathrm{m}\mathrm{o}\mathrm{d}{r}_{k+1}<r_{k+1}\le 1/2r_k$。

综上，必有 $r_{k+2}\le r_k$。

证毕。

故 Euclidean algorithm 的时间复杂度为 $\mathcal{O}(\log n)$。

1.4 代码

LL gcd(LL a, LL b)
{
   
    if (b > a) swap(a, b);
    if (b == 0) return a;
    return gcd(b, a % b);
}

2 扩展欧几里得算法（extended Euclidean algorithm）

2.1 概述

扩展欧几里得算法实际上就是在欧几里得算法的基础上给出了如何构造 $ax+by=gcd(a,b)$ 的解的方法，其中 $a$ 与 $b$ 为不全等于0的两个非负整数， $x$ 和 $y$ 为任意整数。

2.2 算法流程及相关证明

2.2.1 解的存在性证明

下面使用构造法来证明解的存在性

1. 当 $b=0$ 时，显然有 $gcd(a,b)=a$，故一组解为 $x=1,y=0$。

2. 当 $b\ne 0$ 时，设 $a{x}_0+b{y}_0=gcd(a,b)$，$b{x}_1+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y_1=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$。

   由 $gcd(a,b)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$ 可得 $ax$