1.最大公约数(Greatest Common Divisor)
1.1定义
能够整除多个整数的最大正整数。而多个整数不能都为零。例如8和12的最大公因数为4。
——维基百科
1.2 求解a,b的最大公约数
记a,b的最大公约数是gcd(a,b)
1.2.1 穷举思想
- gcd(a,b)一定小于等于a,b中的较小数
- a == b时,gcd(a,b) == a
根据这个特点,可以从1遍历至min(a,b),查找同时整除a,b的最大数。
int gcd(int a,int b) {
// 穷举法实现
if(a == b) {
return a;
}
int min = a < b ? a :b;
int gcd = 0;
for(int i = 1; i <= min; i++) {
if(a % i == 0 && b % i == 0) {
gcd = i;
}
}
return gcd;
}
// 时间复杂度O(min(a,b))
1.2.2 辗转相除法(欧几里得算法)
学过密码学的朋友应该对这个算法有所了解。
原理
两个整数的最大公约数等于其中较小的数和两数相除余数的最大公约数。
首先假设a <= b
,否则交换a跟b
所以gcd(a,b) == gcd(a,b%a)
// 代码实现
int gcd(int a,int b) {
if(a == b) {
return a;
}
if(a > b) {
int temp = a;
a = b;
b = temp;
}
int mod = 0; // 记录余数
while(a != 0) {
mod = b % a;
b = a;
a = mod;
}
return a;
}
// 递归实现
int gcd(int a,int b) {
if( b == 0) {
return a; // 任何非零数与0的最大公约数是其本身
} else {
return gcd(b,a%b); // gcd(a,b) = gcd(b,a%b)
// 这里有一个有趣的点需要注意:
// 为什么没有判断a 和 b 的大小呢?
// a >= b时 gcd(a,b) = gcd(b,a%b) 符合逻辑
// a < b时 a%b == a也就是说gcd(b,a%b) == gcd(b,a)进行了a,b的交换
}
}
2. 最小公倍数
2.1 定义
若有一个数x,可以被另外两个数a、b整除,且x大于(或等于)a和b,则为a和b的公倍数。a和b的公倍数有无限个,而所有的公倍数中,最小的公倍数就叫做最小公倍数。 ——维基百科
2.2 求解 a,b的最小公倍数
思路
两个数的乘积等于其最大公约数与最小公倍数的乘积
证明过程如下:
int lcm(int a,int b) {
return a / gcd(a,b) * b; // a * b可能会溢出,所以先除掉最大公约数,然后再乘b
}