C++数学-数论最大公约数和最小公倍数流食般投喂

理解算法流程，掌握核心知识~

最大公约数和最小公倍数：

        在小学，我们其实已经学习了。

        如果 a 除以 b 没有余数，那么 a 就是 b 的倍数， b 就是 a 的约束，记作 b | a。约数也称因数。

        📌b | a 的含义是：b 能整除 a。

        最大公约数 Greatest Common Divisor，常缩写为 gcd。

        最小公倍数 Least Common Multiple，常缩写为 lcm。

🐅对于两个数求解 gcd 和 lcm ，有如下性质：

        🐅对于两个数 a 和 b ，gcd(a, b) * lcm(a, b) = a * b。也就是最大公约数乘以最小公倍数等于两个数的乘积。

        🐅因此，由于我们找到了最大公约数和最小公倍数的关系，我们常常对求解最小公倍数的题，先求解最大公约数，然后通过两者之间的关系求解最小公倍数。   

🦁欧几里得算法：

        🦁欧几里得算法又称辗转相除法，可以求出两个整数的最大公约数。

📌算法流程：

        设 a > b：

        📌如果 b 是 a 的约数，那么 b 就是两者的最大公约数；

        📌如果 b 不是 a 的约数，那么 gcd(a，b) = gcd(b，a % b)

        因为 a % b 会不断减小，因此可以用递归进行求解。

gcd（8，14）

gcd（14，8）

gcd（8， 6）

gcd（6，2）

gcd（2，0）

LL gcd(LL a, LL b)
{
    if(!b) return a; // 如果 b 等于 0，说明 a 就是最⼤公约数
    return gcd(b, a % b);

    // return b == 0 ? a : gcd(b, a % b);
}

🦁时间复杂度：

求  gcd(a, b)  会遇到两种情况：

1、a < b， gcd(a, b) = gcd(b, a)

2、a > b， gcd(a, b) = gcd(b, a % b)

第二种情况至少让 a 折半，因此最多 log n 次。

第一种情况如上述表格，只会反转一次，因此时间复杂度为 O(log n)。

OJ题来源：洛谷

OJ题名：[信息与未来 2018] 最大公约数

OJ题归属：数学-数论【最大公约数与最小公倍数】

        解题算法：欧几里得算法

#include<iostream>

using namespace std;

int gcd(int a, int b)
{
	return b == 0 ? a : gcd(b, a % b);
}

int main()
{
	int x, y, z; cin >> x >> y >> z;

	cout << gcd(gcd(x, y), z);

	return 0;
}

📌此处穿插讲解一个算法，给下一个题做铺垫——秦久韶算法：

        📌秦久韶算法是将一元 n 次多项式求值问题转换成 n 个一次式的算法。

f(x) = 5x^4 + 4x^3 + 3x^2 + 2x + 1

不断提公因式                 = (5x^3 + 4x^2 + 3x + 2) * x + 1

= ((5x^2 + 4x + 3) * x + 2) * x + 1

= (((5x + 4) * x + 3) * x + 2) * x + 1

x = 6                               = (((5 * 6 + 4) * 6 + 3) * 6 + 2) * 6 + 1

如表所示：不断相乘再相加，一个 for 循环 O(n) 就可以求出 n 次多项式的值。

        单次是相乘再相加，只有乘法和加法，所以该算法可以处处取模~（整体可分阶段取模）

OJ题来源：洛谷

OJ题名：小红的gcd

OJ题归属：数学-数论【最大公约数与最小公倍数】

        解题算法：欧几里得算法  + 秦久韶算法

#include<iostream>
#include<string>

using namespace std;

string a;
int b;

// 欧几里得算法
int gcd(int a, int b)
{
    return b == 0 ? a : gcd(b, a % b);
}

// 秦久韶算法
int clac()
{
    long long t = 0;
    for (auto e : a)
    {
        t = t * 10 + e - '0'; // t * 10 的瞬间会临时变大，如果此时在 int 最大值附近，乘以2就会溢出
        t %= b;
    }
    
    return t;
}

int main()
{
    cin >> a >> b;
    
    // 此处为了 a 不溢出，直接就 gcd(b, a % b) 了
    cout << gcd(b, clac()) << endl;
    
    return 0;
}