NOI大纲——普及组——辗转相除法

辗转相除法

辗转相除法（Euclidean Algorithm）是一种用于计算两个整数的最大公约数（Greatest Common Divisor, GCD）的算法。这种方法基于一个重要的性质：两个数的最大公约数等于其中较小的数与两数之差的最大公约数。进一步简化，可以用模运算来替代减法操作。

算法原理

假设我们要计算两个非负整数 (a) 和 (b) 的最大公约数，且 (a \geq b)。辗转相除法的步骤如下：

1. 如果 (b = 0)，则 (a) 就是 (a) 和 (b) 的最大公约数，算法结束。
2. 否则，将 (a) 取模 (b) 得到余数 (r)，用 (b) 和 (r) 递归地继续这个过程。

用公式表示就是：
[ \text{GCD}(a, b) = \text{GCD}(b, a % b) ]

这个过程会一直递归下去，直到余数为 0。最终的非零除数就是两个数的最大公约数。

C++ 实现

下面是辗转相除法的递归和迭代实现的 C++ 代码：

递归实现

#include <iostream 
using namespace std;

int gcd_recursive(int a, int b) {
    if (b == 0)
        return a;
    return gcd_recursive(b, a % b);
}

int main() {
    int a, b;
    cout << "Enter two numbers: ";
    cin >> a >> b;
    cout << "GCD of " << a << " and " << b << " is " << gcd_recursive(a, b) << endl;
    return 0;
}

迭代实现

#include <iostream>
using namespace std;

int gcd_iterative(int a, int b) {
    while (b != 0) {
        int temp = b;
        b = a % b;
        a = temp;
    }
    return a;
}

int main() {
    int a, b;
    cout << "Enter two numbers: ";
    cin >> a >> b;
    cout << "GCD of " << a << " and " << b << " is " << gcd_iterative(a, b) << endl;
    return 0;
}

运行示例

假设输入两个数 (48) 和 (18)：

• 递归过程：

  1. (\text{GCD}(48, 18) = \text{GCD}(18, 48 % 18) = \text{GCD}(18, 12))
  2. (\text{GCD}(18, 12) = \text{GCD}(12, 18 % 12) = \text{GCD}(12, 6))
  3. (\text{GCD}(12, 6) = \text{GCD}(6, 12 % 6) = \text{GCD}(6, 0))
  4. (\text{GCD}(6, 0) = 6)

• 迭代过程：

  1. 初始 (a = 48), (b = 18)
  2. (b \neq 0)：(a % b = 48 % 18 = 12)，交换 (a = 18), (b = 12)
  3. (b \neq 0)：(a % b = 18 % 12 = 6)，交换 (a = 12), (b = 6)
  4. (b \neq 0)：(a % b = 12 % 6 = 0)，交换 (a = 6), (b = 0)
  5. (b = 0)，停止迭代，结果为 (a = 6)

无论是递归还是迭代，结果都是 6。

时间复杂度

辗转相除法的时间复杂度是 (O(\log \min(a, b)))，因为每次操作后，较大数会大致减半。这种高效性使其在计算最大公约数时非常实用。