最大公约数和最小公倍数

意义

例： 把一张长60厘米，宽40厘米的长方形纸板剪成边长是整数厘米数的小正方形，且无剩余，最少可以剪成多少块？

计算最大公约数(60 ，40） = 20
(60/20)*(40/20) = 6(块)

例：用长5CM，宽3CM的长方形硬纸片摆成一个正方形（中间无空隙），至少要用几个长方形硬纸片？

计算最小公倍数：[5,3] = 15
则需要(15/5)*(15/3) = 15(片)

计算方法

两个整数的最大公约数

辗转相除法
两个数a,b（a>b）的最大公约数
如果a%b != 0 则等于b和a%b的最大公约数。
如果a%b == 0 则等于b。

更相减损术
两个数a,b（a>b）的最大公约数
如果a-b != 0 则等于b和a-b的最大公约数。
如果a-b == 0 则等于b。

综合方法
两个数a,b的最大公约数为gcd(a，b)，则

• 如果a == b，则gcd(a,b) = b。
• gcd(a,b) = 2*gcd(a>>2,b>>2)，如果a，b均为偶数。
• gcd(a,b) = gcd(a>>1,b)，如果a为偶数，b为奇数。
• gcd(a,b) = gcd(a/2,b>>1)，如果a为奇数，b为偶数。
• gcd(a,b) = gcd(b,a-b)，如果a，b均为奇数。

多个整数的最大公约数

多个数的最大公约数先求前两个数的最大公约数，再求这个最大公约数与下一个数的最大公约数。这样依次下去，直到最后一个数为止。最后所得的一个最大公约数，就是所求的多个数的最大公约数。

两个数的最小公倍数

先计算出两个数的最大公约数，在两数相乘除以最大公约数即为最小公倍数。例如：12 * 16 / 4 = 48.

多个数的最小公倍数

先求先两个数的最小公倍数，再计算这个最小公倍数与下一个数的最小公倍数，依次迭代即可。
例如：(48，18)表示计算最大公约数。
12 16 18
前两个数的最小公倍数为48，
然后 48 * 18 / (48,18) = 48 * 18 / 6 = 48 * 3 = 144。

代码

注：只考虑正整数情况，未考虑负数和大数情况

#include <iostream>
#include <vector>
#include <cassert>
using namespace std;

class Solution {
public:
    /*
     * gcd1:辗转相除法，两个数的最大公约数等于两个数的余数和这两个数中较小值的最大公约数，缺点：求余数运算较慢
    */
    int gcd1(int a,int b)
    {
        int max = a > b ? a:b;
        int min = a > b ? b:a;
        if((max % min) != 0)
        {
            int tmp = max % min;
            return gcd1(min,tmp);
        }
        return min;
    }
    int gcd2(int a,int b)
    {
        if(a < b)
            return gcd2(b,a);
        if(a - b !=0)
        {
            return gcd2(b,a-b);
        }
        return b;
    }
    int gcd3(int a,int b)
    {
        if(a == b)
            return b;
        if(a < b)
            return gcd3(b,a);
        if((a&0x1) == 0 && (b & 0x1) == 0)
            return gcd3(a >> 1,b >> 1) << 1;
        else if((a&0x1) == 1 && (b & 0x1) == 0)
            return gcd3(a,b >> 1);
        else if((a&0x1) == 0 && (b & 0x1) == 1)
            return gcd3(a >> 1,b);
        else if((a&0x1) == 1 && (b & 0x1) == 1)
            return gcd3(b,a-b);
    }
    int gcd(vector<int> a)
    {
        assert(a.size() >=2);
        int tmp = gcd3(a[0],a[1]);
        for(auto c = a.begin()+2;c!=a.end();++c)
        {
            tmp = gcd3(tmp,*c);
        }
        return tmp;
    }

    int lcm1(int a,int b)
    {
        return a * b / (gcd3(a,b));
    }
    int lcm(vector<int> a)
    {
        assert(a.size() > 1);
        int tmp = lcm1(a[0],a[1]);
        for(auto c = a.begin()+2;c!=a.end();++c)
        {
            tmp = lcm1(tmp,*c);
        }
        return tmp;
    }
};

int main()
{
    Solution s;
    cout << s.gcd1(12,16) << endl;
    cout << s.gcd2(12,16) << endl;
    cout << s.gcd3(12,16) << endl;
    vector<int> a = {12,16,18};
    cout << s.gcd(a) << endl;
    cout << s.lcm1(12,16) << endl;
    cout << s.lcm(a) << endl;
}