扩展欧几里德算法 递归法 递推法 手算法 原理及实现

扩展欧几里德算法 递归法 递推法 手算法 原理及实现

顾名思义, 扩展欧几里德算法 是在 欧几里德算法 基础上扩展的算法.

欧几里德算法 和 扩展欧几里德算法 在用途上的区别:

• 欧几里德算法(gcd): 即求两个整数的 最大公约数.
• 扩展欧几里德算法:
  • 用于求 乘法逆元.
  • 用于求 贝组等式 的一个解.

欧几里德算法

即 辗转相除法.

C 语言实现:

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

注意此算法的 终止条件 是 b == 0.

扩展欧几里德-递归法

设初始值: $a_1=m$, $b_1=n$, 且 $m>n$.

设 $a_1{x}_1+b_1{y}_1=\gcd (a_1,b_1)$, 则:

在 欧几里德算法 中, $a_{i+1}=b_i,b_{i+1}=a_i\mathrm{mod}{b}_i$, 所以:

• $a_2{x}_2+b_2{y}_2=b_1{x}_2+(a_1\mathrm{mod}{b}_1)y_2=\gcd (a_2,b_2)$.

因为 $a_1\mathrm{mod}{b}_1=a_1-(a_1/b_1)\cdot b_1$ (其中的 / 为整除). 得:

• $a_1{x}_1+b_1{y}_1=b_1{x}_2+(a_1-(a_1/b_1)b_1)y_2=a_1{y}_2+b_1(x_2-(a_1/b_1)y_2)$

所以得到 $x_1=y_2,y_1=x_2-(a_1/b_1)y_2$.

由以上推广得到:

推论 1: $\boxed{x_i=y_{i+1},y_i=x_{i+1}-(a_i/b_i)y_{i+1}}$

在 欧几里德算法 中, 迭代的终止条件是 b == 0, 并且此时 $a_i{x}_i+b_i{y}_i=\gcd (a_i,b_i)=\gcd (a_i,0)=a_i$, 显然此时的 $x_i=1,y_i=0$, 所以有:

推论 2: $\boxed{当b_i=0时,x_i=1,y_i=0.}$

以上标识符(比如 $a_1$, $b_2$)的下标数字是指在 欧几里德算法 的 gcd() 函数里的迭代次数.

为了更好的理解, 下面用实例进行说明, 假设 m = 1756, n = 551, $q$ 代表 $a/b$, / 表示整除:

(1) $a_1=1756b_1=551q_1=3$

(2) $a_2=551b_2=103q_2=5$

(3) $a_3=103b_3=36q_3=2$

(4) $a_4=36b_4=31q_4=1$

(5) $a_5=31b_5=5q_5=6$

(6) $a_6=5b_6=1q_6=5$

(7) $a_7=1b_7=0$ ($b_7$ 为 0, 终止, 所以 $x_7=1,y_7=0$)

因为我们要求的是 $x_1,y_1$ 于是乎:

• 递归最深处: 返回 $x_7=1,y_7=0$.

• 第 6 深: 返回

  $\begin{array}{rl}{x}_6& =y_7=0\end{array}$

  $\begin{array}{rl}{y}_6& =x_7-q_6{y}_7\\ & =1-0=1\end{array}$

• 第 5 深: 返回

  $\begin{array}{rl}{x}_5& =y_6=1\end{array}$

  $\begin{array}{rl}{y}_5& =x_6-q_5{y}_6\\ & =0-6\cdot 1=-6\end{array}$

• 第 4 深: 返回

  $\begin{array}{rl}{x}_4& =y_5=-6\end{array}$

  $\begin{array}{rl}{y}_4& =x_5-q_4{y}_5\\ & =1-1\cdot (-6)=7\end{array}$

• 第 3 深: 返回

  $\begin{array}{rl}{x}_3& =y_4=7\end{array}$

  $\begin{array}{rl}{y}_3& =x_4-q_3{y}_4\\ & =-6-2\cdot 7=-20\end{array}$

• 第 2 深: 返回

  $\begin{array}{rl}{x}_2& =y_3=-20\end{array}$

  $\begin{array}{rl}{y}_2& =x_3-q_2{y}_3\\ & =7-5\cdot (-20)=107\end{array}$

• 第 1 深: 返回

  $\begin{array}{rl}{x}_1& =y_2=107\end{array}$

  $\begin{array}{rl}{y}_1& =x_2-q_1{y}_2\\ & =-20-3\cdot 107=-341\end{array}$

即求得 $x=107,y=-341$, 代入 $1756x+551y=\gcd (1756,551)=1$ 中验算成立.

写出 C 语言代码为:

#include <stdio.h>

struct Solution {
  int x;
  int y;
};

struct Solution extended_euclid(int a, int b) {
  if (b == 0) {
    return Solution{1, 0};
  }
  struct Solution t = extended_euclid(b, a % b);
  int q = a / b;
  return Solution{t.y, t.x - q * t.y};
}

int main() {
    int m = 1756, n = 551;
    printf("for %dx + %dy = gcd(%d, %d)\n", m, n, m, n);
    struct Solution res = extended_euclid(m, n);
    printf("the solution is: x = %d, y = %d\n", res.x, res.y);
    return 0;
}

程序运行结果为:

for 1756x + 551y = gcd(1756, 551)
the solution is: x = 107, y = -341

扩展欧几里德-递推法

扩展欧几里德算法 在密码学中算是很常用的了, 古典密码中的 仿射密码 以及各种公钥密码好像都用的挺多, 递推法的一个好处是方便手算, 尤其是在期末考试 (ˉ▽ˉ；)…

跟递归法一样, 我们先设: $a_1=m$, $b_1=n$, 且 $m>n$.

然后以下式子成立是必然的:

$$\begin{array}{c}{a}_1=1\cdot a_1+0\cdot b_1\\ b_1=0\cdot a_1+1\cdot b_1\end{array}$$

我们根据上式 (1) 设 $x_{-1}=1,y_{-1}=0$.

根据上式 (2) 设 $x_0=0,y_0=1$.

设 $q$ 表示 a 和 b 商, $u$ 表示 $a÷b$ 的余数, 则:

• 第一轮:

  $$\begin{gathered} a_1=m=1\cdot a_1+0\cdot b_1 \\ {b}_1=n=0\cdot a_1+1\cdot b_1 \end{gathered}$$

  $q_1=a_1/b_1$ (/ 表示整除)

  $\begin{array}{rl}{u}_1& =a_1-q_1{b}_1(关键!将a_1,b_1代入)\\ & =(1\cdot a_1+0\cdot b_1)-q_1(0\cdot a_1+1\cdot b_1)\\ & =a_1+(-q_1)b_1\\ & =x_1{a}_1+y_1{b}_1\end{array}$

  得到: $x_1=1,y_1=-q_1$.

• 第二轮:

  $$\begin{gathered} a_2=b_1=x_0{a}_1+y_0{b}_1 \\ {b}_2=u_1=x_1{a}_1+y_1{b}_1 \end{gathered}$$

  $q_2=a_2/b_2$

  $\begin{array}{rl}{u}_2& =a_2-q_2{b}_2\\ & =(x_0{a}_1+y_0{b}_1)-q_2(x_1{a}_1+y_1{b}_1)\\ & =(x_0-q_2{x}_1)a_1+(y_0-q_2{y}_1)b_1\\ & =x_2{a}_1+y_2{b}_1\end{array}$

  得到: $x_2=x_0-q_2{x}_1,y_2=y_0-q_2{y}_1$

• 第三轮:

  $$\begin{gathered} a_3=b_2=x_1{a}_1+y_1{b}_1 \\ {b}_3=u_2=x_2{a}_1+y_2{b}_1 \end{gathered}$$

  $q_3=a_3/b_3$

  $\begin{array}{rl}{u}_3& =a_3-q_3{b}_3\\ & =(x_1{a}_1+y_1{b}_1)-q_2(x_2{a}_1+y_2{b}_1)\\ & =(x_1-q_3{x}_2)a_1+(y_1-q_3{y}_2)b_1\\ & =x_3{a}_1+y_3{b}_1\end{array}$

  得到: $x_3=x_1-q_3{x}_2,y_3=y_1-q_3{y}_2$

• …

• 第 i 轮:

  $$\begin{gathered} a_i=b_{i-1}=u_{i-2}=x_{i-2}{a}_1+y_{i-2}{b}_1 \\ {b}_i=u_{i-1}=x_{i-1}{a}_1+y_{i-1}{b}_1 \end{gathered}$$

  $q_i=a_i/b_i$

  $\begin{array}{rl}{u}_i& =a_i-q_i{b}_i\\ & =(x_{i-2}{a}_1+y_{i-2}{b}_1)-q_2(x_{i-1}{a}_1+y_{i-1}{b}_1)\\ & =(x_{i-2}-q_i{x}_{i-1})a_1+(y_{i-2}-q_i{y}_{i-1})b_1\\ & =x_i{a}_1+y_i{b}_1\end{array}$

  得到 推理: $\boxed{x_i=x_{i-2}-q_i{x}_{i-1},y_i=y_{i-2}-q_i{y}_{i-1}}$

关于终止条件:

因为我们要求解的是 $x{a}_1+y{b}_1=\gcd (a_1,b_1)$, 只有当 $u_{i+1}$ 为 0 时, 其第 i 轮的 $u_i$ 才等于 $\gcd (a_1,b_1)$. 所以我们判断终止有两种方法:

1. 判断 $u_{i+1}$ 为 0 时, 其解为 $x_i$ 和 $y_i$.
2. 判断 $u_i$ 为 0 时, 其解为 $x_{i-1}$ 和 $y_{i-1}$.

在手算时用第一种好一点. 在代码中我使用第二种.

实例讲解, 还是使用上面的例子:

假设 m = 1756, n = 551, $q$ 代表 $a/b$, / 表示整除:

• 第一轮:

  $a_1=1756b_1=551$

  $q_1=a_1/b_1=3$

  $u_1=a_1\mathrm{mod}{b}_1=103$

  $x_1=x_{-1}-q_1{x}_0=1-3\cdot 0=1$

  $y_1=y_{-1}-q_1{y}_0=0-3\cdot 1=-3$

  $u_2=b_1\mathrm{mod}{u}_1=36$ 不为 0, 继续.

• 第二轮:

  $a_2=b_1=551b_2=u_1=103$

  $q_2=a_2/b_2=5$

  $u_2=a_2\mathrm{mod}{b}_2=36$

  $x_2=x_0-q_2{x}_1=0-5\cdot 1=-5$

  $y_2=y_0-q_2{y}_1=1-5\cdot (-3)=16$

  $u_3=b_2\mathrm{mod}{u}_2=31$ 不为 0, 继续.

• 第三轮:

  $a_3=b_2=103b_3=u_2=36$

  $q_3=a_3/b_3=2$

  $u_3=a_3\mathrm{mod}{b}_3=31$

  $x_3=x_1-q_3{x}_2=1-2\cdot (-5)=11$

  $y_3=y_1-q_3{y}_2=-3-2\cdot 16=-35$

  $u_4=b_3\mathrm{mod}{u}_3=5$ 不为 0, 继续.

• 第四轮:

  $a_4=b_3=36b_4=u_3=31$

  $q_4=a_4/b_4=1$

  $u_4=a_4\mathrm{mod}{b}_4=5$

  $x_4=x_2-q_4{x}_3=-5-1\cdot 11=-16$

  $y_4=y_2-q_4{y}_3=16-1\cdot (-35)=51$

  $u_5=b_4\mathrm{mod}{u}_4=1$ 不为 0, 继续.

• 第五轮:

  $a_5=b_4=31b_5=u_4=5$

  $q_5=a_5/b_5=6$

  $u_5=a_5\mathrm{mod}{b}_5=1$

  $x_5=x_3-q_5{x}_4=11-6\cdot (-16)=107$

  $y_5=y_3-q_5{y}_4=-35-6\cdot 51=-341$

  $u_6=b_5\mathrm{mod}{u}_5=0$ 为 0, 终止!!!

即得到解 $x=x_5=107,y=y_5=-341$.

写出 C 语言代码为:

#include <stdio.h>

struct Solution
{
    int x;
    int y;
};

struct Solution extended_euclid2(int m, int n)
{
    // 初始
    int x_i_2 = 1; // x_{i - 2}
    int y_i_2 = 0; // y_{i - 2}
    int x_i_1 = 0; // x_{i - 1}
    int y_i_1 = 1; // y_{i - 1}

    int a = m;
    int b = n;
    int x, y;
    while (true)
    {
        int q = a / b; // q_i a_i b_i
        int u = a % b; // u_i
        // 当 u_i 为 0 则终止
        if (u == 0)
        {
            x = x_i_1; // 得出解 x = x_{i - 1}
            y = y_i_1; // 得出解 y = y_{i - 1}
            break;
        }
        x = x_i_2 - q * x_i_1; // x_i
        y = y_i_2 - q * y_i_1; // y_i
        x_i_2 = x_i_1;
        y_i_2 = y_i_1;
        x_i_1 = x;
        y_i_1 = y;
        a = b; // a_{i + 1}
        b = u; // b_{i + 1}
    }
    // 返回求得的 x y.
    return Solution{x, y};
}

int main()
{
    int m = 1756, n = 551;
    printf("for %dx + %dy = gcd(%d, %d)\n", m, n, m, n);
    struct Solution res = extended_euclid2(m, n);
    printf("the solution is: x = %d, y = %d\n", res.x, res.y);
    return 0;
}

程序运行结果为:

for 1756x + 551y = gcd(1756, 551)
the solution is: x = 107, y = -341

扩展欧几里德-手算法

其方法跟 扩展欧几里德-递推法 里的例子大同小异, 主要是我们密码学老师喜欢用这种方法手算:

建表, $u$ 代表余数, $q$ 代表商:

u	x	y	q
$a_1$	$x_{-1}$	$y_{-1}$
$b_1$	$x_0$	$y_0$	$q_1$
$u_1$	$x_1=x_{-1}-q_1{x}_0$	$y_1=y_{-1}-q_1{y}_0$	$q_2$
$u_2$	$x_2=x_0-q_2{x}_1$	$y_2=y_0-q_2{y}_1$	$q_3$
$u_3$	$x_3=x_1-q_3{x}_2$	$y_3=y_1-q_3{y}_2$	$q_4$
…
$u_i$	$x_i=x_{i-2}-q_i{x}_{i-1}$	$y_i=y_{i-2}-q_i{y}_{i-1}$	$q_{i+1}$

其 终止条件 是: 判断 $u_{i+1}$ 为 0, 解为 $x_i$ 和 $y_i$.

还是以 m = 1756, n = 551 为例, 建表:

u	x	y	q
1756	1	0
551	0	1	3
103	1	-3	5
36	-5	16	2
31	11	-35	1
5	-16	51	6
1	107	-347	0
0

嗯, 差不多就是这样╮(╯▽╰)╭.