辗转相除法
辗转相除法(Euclidean Algorithm)是一种用于计算两个整数的最大公约数(Greatest Common Divisor, GCD)的算法。这种方法基于一个重要的性质:两个数的最大公约数等于其中较小的数与两数之差的最大公约数。进一步简化,可以用模运算来替代减法操作。
算法原理
假设我们要计算两个非负整数 (a) 和 (b) 的最大公约数,且 (a \geq b)。辗转相除法的步骤如下:
- 如果 (b = 0),则 (a) 就是 (a) 和 (b) 的最大公约数,算法结束。
- 否则,将 (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):
-
递归过程:
- (\text{GCD}(48, 18) = \text{GCD}(18, 48 % 18) = \text{GCD}(18, 12))
- (\text{GCD}(18, 12) = \text{GCD}(12, 18 % 12) = \text{GCD}(12, 6))
- (\text{GCD}(12, 6) = \text{GCD}(6, 12 % 6) = \text{GCD}(6, 0))
- (\text{GCD}(6, 0) = 6)
-
迭代过程:
- 初始 (a = 48), (b = 18)
- (b \neq 0):(a % b = 48 % 18 = 12),交换 (a = 18), (b = 12)
- (b \neq 0):(a % b = 18 % 12 = 6),交换 (a = 12), (b = 6)
- (b \neq 0):(a % b = 12 % 6 = 0),交换 (a = 6), (b = 0)
- (b = 0),停止迭代,结果为 (a = 6)
无论是递归还是迭代,结果都是 6。
时间复杂度
辗转相除法的时间复杂度是 (O(\log \min(a, b))),因为每次操作后,较大数会大致减半。这种高效性使其在计算最大公约数时非常实用。