一 辗转相除法
辗转相除法, 又名欧几里得算法(Euclidean algorithm),目的是求出两个正整数的最大公约数。它是已知最古老的算法, 其可追溯至公元前300年前。这条算法基于一个定理:两个正整数a和b(a>b),它们的最大公约数等于a除以b的余数c和b之间的最大公约数。比如10和25,25除以10商2余5,那么10和25的最大公约数,等同于10和5的最大公约数。
public class LZN {
public static int getMaxCommonDivider(int numA, int numB) {
int result = 1;
if (numA > numB) {
result = gcd(numA, numB);
} else {
result = gcd(numB, numA);
}
return result;
}
private static int gcd(int a, int b) {
if (a % b == 0) {
return b;
} else {
return gcd(b, a % b);
}
}
public static void main(String[] args) {
System.out.println(getMaxCommonDivider(20, 30));
}
}
当两个整型数比较大时,取模运算的性能比较低,有什么好办法吗?
二 更相减损术
更相减损术, 出自于中国古代的《九章算术》,也是一种求最大公约数的算法。他的原理更加简单:两个正整数a和b(a>b),它们的最大公约数等于a-b的差值c和较小数b的最大公约数。比如10和25,25减去10的差是15,那么10和25的最大公约数,等同于10和15的最大公约数。
public class LZN {
public static int gcd(int a, int b) {
if (a == b) {
return a;
}
if (a > b) {
return gcd(b, a - b);
} else {
return gcd(b, b - a);
}
}
public static void main(String[] args) {
System.out.println(gcd(60, 30));
}
}
通过更相减损术,我们避免了大整数取模运算的性能问题,但是相减的方法的递归次数比辗转相除法多,这个又应该如何避免呢?
三 终极解法
众所周知,移位运算的性能非常快。对于给定的正整数a和b,不难得到如下的结论。其中gcb(a,b)的意思是a,b的最大公约数函数:
当a和b均为偶数,gcb(a,b) = 2*gcb(a/2, b/2) = 2*gcb(a>>1, b>>1)
当a为偶数,b为奇数,gcb(a,b) = gcb(a/2, b) = gcb(a>>1, b)
当a为奇数,b为偶数,gcb(a,b) = gcb(a, b/2) = gcb(a, b>>1)
当a和b均为奇数,利用更相减损术运算一次,gcb(a,b) = gcb(b, a-b), 此时a-b必然是偶数,又可以继续进行移位运算。
public class LZN {
public static int gcd(int a, int b) {
if (a == b) {
return a;
}
if (a < b) {
return gcd(b, a);
} else {
// 和1做按位与运算 判断奇偶
if ((a & 1) == 0 && (b & 1) == 0) {
return gcd(a >> 1, b >> 1) << 1;
} else if ((a & 1) == 1 && (b & 1) == 0) {
return gcd(a, b >> 1);
} else if ((a & 1) == 0 && (b & 1) == 1) {
return gcd(a >> 1, b);
} else {
return gcd(b, a - b);
}
}
}
public static void main(String[] args) {
System.out.println(gcd(5, 45));
}
}