扩展欧几里得算法

Bézout’s定理

$\forall a,b\in Z$，$\exists x,y\in Z^{+}$，使得
$$ax+by=gcd(a,b)\text{\hspace{1em}(0)}$$

扩展欧几里得算法求该定理的解

已知
$$\begin{array}{r}gcd(a,b)=gcd(b,a\%b)\end{array}\text{\hspace{1em}(1)}$$
根据$B\acute{e}zou{t}^{\prime }s$定理，$\exists x^{{}^{\prime }},y^{{}^{\prime }}\in Z^{+}$，使得：

$$\begin{array}{r}b{x}^{{}^{\prime }}+(a\%b)y^{{}^{\prime }}=gcd(b,a\%b)\end{array}\text{\hspace{1em}(2)}$$
又

$$\begin{array}{r}a\%b=a-b\ast \lfloor a/b\rfloor \end{array}\text{\hspace{1em}(3)}$$
将$(3)$带入$(2)$，得：

$$\begin{array}{rl}b{x}^{{}^{\prime }}+(a-b\ast \lfloor a/b\rfloor )y^{{}^{\prime }}& =gcd(b,a\%b)\\ a{y}^{{}^{\prime }}+b(x^{{}^{\prime }}-\lfloor a/b\rfloor y^{{}^{\prime }})& =gcd(b,a\%b)\end{array}\text{\hspace{1em}(4)}$$
令$x=y^{{}^{\prime }}$，$y=x^{{}^{\prime }}-\lfloor a/b\rfloor y^{{}^{\prime }}$，将$(1)$代入$(4)$，整理得：

$$\begin{array}{r}ax+by=gcd(b,a\%b)=gcd(a,b)\end{array}\text{\hspace{1em}(5)}$$
$(5)$说明，知道$(2)$的解就可以推出定理的解。公式(0)到公式(2)明显是欧几里得算法求解最大公约数的过程。公式(3)，(4)，(5)是说明：公式(0)的解与公式(2)的解的关系。将上述过程归纳成代码形式：

void exgcd(int a, int b, int& x, int& y) {	//公式(0)
    if (b == 0) {
        x = 1;
        y = 0;
        return a;
    }
    exgcd(b, a % b, x, y);	//公式(2)
    int tp = x;
    x = y;
    y = tp - a / b * y;
}

参考

• https://www.cnblogs.com/fusiwei/p/11775503.html