C++欧几里得算法求最大公约数和最小公倍数

定义

最大公约数即为 Greatest Common Divisor,常缩写为 gcd。
一组整数的公约数,是指同时是这组数中每一个数的约数的数。
一组整数的最大公约数,是指所有公约数里面最大的一个。
那么如何求最大公约数呢?我们先考虑两个数的情况。

欧几里得算法

过程

如果我们已知两个数 ab,如何求出二者的最大公约数呢?
不妨设a>b
我们发现如果 b 是 a 的约数,那么 b 就是二者的最大公约数。 下面讨论不能整除的情况,即a=bq+r,
其中r<b
我们通过证明可以得到gcd(a,b)=gcd(b,amodb),过程如下:
a=bk+c,显然有 c=amodb。设 da, db,则c=abk, cd=adbdk
由右边的式子可知cd 为整数,即 dc,所以对于 a,b 的公约数,它也会是 b,amodb 的公约数。
反过来也需要证明:
db, d(amodb),我们还是可以像之前一样得到以下式子
amodbd=adbdk, amodbd+bdk=ad
因为左边式子显然为整数,所以ad 也为整数,即 d \mid a,所以 b,a\bmod b 的公约数也是 a,b 的公约数。
既然两式公约数都是相同的,那么最大公约数也会相同。
所以得到式子gcd(a,b)=gcd(b,amodb)
既然得到了 gcd(a,b)=gcd(b,r),这里两个数的大小是不会增大的,那么我们也就得到了关于两个数的最大公约数的一个递归求法。

实现

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

最小公倍数

int gcd(int a, int b) {
    if(b == 0) return a;
    return gcd(b, a % b);
}
int lcm = a * b / gcd(a, b);

© 版权声明
THE END
喜欢就支持一下吧
点赞0

Warning: mysqli_query(): (HY000/3): Error writing file '/tmp/MY5BHd9R' (Errcode: 28 - No space left on device) in /www/wwwroot/583.cn/wp-includes/class-wpdb.php on line 2345
admin的头像-五八三
评论 抢沙发
头像
欢迎您留下宝贵的见解!
提交
头像

昵称

图形验证码
取消
昵称代码图片