最大公约数和最小公倍数
一、最大公约数
1.1 穷举法
两个数 a 和 b ( a > b ), 从较小者 b 开始向下递减,如果出现能被两个数整除的数,则为最大公约数。但是效率很低,例如10000和10001求最大公约数,需要遍历10000次才能求出结果
public static int getGreatestCommonDivisor(int a, int b) {
int smallNum = a < b ? a : b;
int bigNum = a > b ? a : b;
for (int i = smallNum / 2; i > 1; i--) {
if (smallNum % i == 0 && bigNum % i == 0) {
return i;
}
}
return 1;
}
1.2 辗转相除法
两个数 a 和 b ( a > b ),a 对 b 取模,模为 c
若 c 为零,则 c 为最大公约数;
若 c 不为零, 则 a 和 b 的最大公约数为 c 和 b 的最大公约数
但是当两个整数较大的时候,模运算的效率会比较的低
public static int getGreatestCommonDivisor1(int a, int b) {
int smallNum = a < b ? a : b;
int bigNum = a > b ? a : b;
if (bigNum % smallNum == 0) {
return smallNum;
}
return getGreatestCommonDivisor1(bigNum % smallNum, smallNum);
}
1.3 更相减损法
两个数 a 和 b ( a > b ),a 减 b, 差为 c
若 c 为零,则 c 为最大公约数
若 c 不为零,则 a 和 b 的最大公约数为 c 和 b的最大公约数
但是,求差运算的次数远大于取模运算的次数,例如10000和1求最大公约数,需要递归9999次
public static int getGreatestCommonDivisor2(int a, int b) {
if (a == b) {
return a;
}
int smallNum = a < b ? a : b;
int bigNum = a > b ? a : b;
return getGreatestCommonDivisor2(bigNum - smallNum, smallNum);
}
1.4 最优解
辗转相除法和更相减损法相结合,使用移位运算进行求解
当a为偶数,b为偶数时,gcd(a,b)= 2 * gcd(a / 2, b / 2)= 2 * gcd(a >> 1, b >> 1)
当a为偶数,b为奇数时,gcd(a,b)= gcd(a / 2, b)= gcd(a >> 1, b)
当a为奇数,b为偶数时,gcd(a,b)= gcd(a, b / 2)= gcd(a, b >> 1)
当a为奇数,b为奇数时,先利用更相减损法运算一次,gcd(a,b)= gcd(a -b, b ),此时 a -b 一定是偶数,再次进行移位运算
public static int gcd(int a, int b) {
if (a == b) {
return a;
} else if ((a & 1) == 0 && (b & 1) == 0) {
// a 为偶数,b 也为偶数
return gcd(a >> 1, b >> 1) << 1;
} else if ((a & 1) == 0 && (b & 1) != 0) {
// a 为偶数,b 为奇数
return gcd(a >> 1, b);
} else if ((a & 1) != 0 && (b & 1) == 0) {
// a 为奇数,b 为偶数
return gcd(a, b >> 1);
} else {
int smallNum = a < b ? a : b;
int bigNum = a > b ? a : b;
return gcd(bigNum - smallNum, smallNum);
}
}
二、最小公倍数
有两个整数 a 和 b
最小公倍数 = a * b / 最大公约数
public static int getLowestCommonMultiple(int a, int b) {
return a * b / gcd(a, b);
}