可以证明gcd(a,b) = gcd(b,a % b)(a > b)
所以,最简单的欧几里得算法求gcd,这是最简单的O(logN)的方法。
int gcd(int a,int b)
{
return b == 0 ? a : gcd(b,a % b);
}
//这个玩意可以这么读:b == 0否?b == 0的话返回a,b != 0的话返回gcd(b,a%b)
//就是辗转相除,除到0了那么上一个除数就是最大公约数
对了,lcm和gcd毕竟是一家子。
lcm:
int lcm(int a,int b)
{
return a * b / gcd(a,b);
}
其实如果遇到需要用高精度的数,高精%(就是高精/和高精-(双高精))可能不那么好打,还有O(N)的(自己猜的)
更相减损法
int gcd(int a,int b)
{
if(a == b) return a;
if(a < b) swap(a,b);
return gcd(b,a - b);
}
只要一个高精减。
额要想再快一点可以二进制(不过好像要高精*/2和-)
int gcd(int m,int n)
{
if(m == n) return m;
if(m < n) return gcd(n,m);
if(m & 1 == 0) return (n & 1 == 0) ? 2 * gcd(m / 2,n / 2) : gcd(m / 2,n);
return (n & 1 == 0) ? gcd(m,n / 2) : gcd(n,m - n);
}
凸显一下二进制
int gcd(int m,int n)
{
if(m == n) return m;
if(m < n) return gcd(n,m);
if(m & 1 == 0) return (n & 1 == 0) ? gcd(m << 1,n << 1) >> 1: gcd(m << 1,n);
return (n & 1 == 0) ? gcd(m,n << 1) : gcd(n,m - n);
}