gcd

可以证明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);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值