C++入门小案例_018：最大公约数

最大公约数

基本概念

最大公因数，也称最大公约数、最大公因子，指两个或多个整数共有约数中最大的一个。a，b的最大公约数记为（a，b），同样的，a，b，c的最大公约数记为（a，b，c），多个整数的最大公约数也有同样的记号。求最大公约数有多种方法，常见的有质因数分解法、短除法、辗转相除法、更相减损法。与最大公约数相对应的概念是最小公倍数，a，b的最小公倍数记为[a，b]。

求解方法

质因数分解法

这种方法不适用于编程方法

• 概念：

  把每个数分别分解质因数，再把各数中的全部公有质因数提取出来连乘，所得的积就是这几个数的最大公约数。

  例如：求24和60的最大公约数，先分解质因数，得24=2×2×2×3，60=2×2×3×5，24与60的全部公有的质因数是2、2、3，它们的积是2×2×3=12，所以，（24，60）=12。

短除法

这种方法不适用于编程方法

• 概念：

  短除法求最大公约数，先用这几个数的公约数连续去除，一直除到所有的商互质为止，然后把所有的除数连乘起来，所得的积就是这几个数的最大公约数。
  
  短除法的本质就是质因数分解法，只是将质因数分解用短除符号来进行。

---

无论是短除法，还是分解质因数法，在质因数较大时，都会觉得困难。这时就需要用新的方法。

辗转相除法

• 概念：

  辗转相除法是求两个自然数的最大公约数的一种方法，也叫欧几里德算法。

  例如，求（319，377）：
  ∵ 319÷377=0（余319）
  ∴（319，377）=（377，319）；
  ∵ 377÷319=1（余58）
  ∴（377，319）=（319，58）；
  ∵ 319÷58=5（余29）
  ∴ （319，58）=（58，29）；
  ∵ 58÷29=2（余0）
  ∴ （58，29）= 29；
  ∴ （319，377）=29。
  可以写成右边的格式。
  用辗转相除法求几个数的最大公约数，可以先求出其中任意两个数的最大公约数，再求这个最大公约数与第三个数的最大公约数，依次求下去，直到最后一个数为止。最后所得的那个最大公约数，就是所有这些数的最大公约数。

• 个人解读

以因质分解的角度去理解辗转相除：
1.有两个整数 m 、n （m < n）;
2.如果m<n，辗转相除第一步相当于一个初始化过程，换成 n、m；也就是说 (m,n) == (n,m)；此时大数在前，小数在后；
3.继续辗转相除， n / m = i ， n % m = j ；也就是说 n = m*i + j， (n,m) == (m,j)；
4.理解上一步后面就能理解了：

在 n > m 的情况下:
如果 n 能被 m 整除，那么 m 就是最大公约数；
如果 n 不能被 m 整除，那么 n = m*i + j ；如果找到一个最大数，既能被 m 整除 又能被 j 整除，这个数 ret 就是最大公约数；假设 m = ret * x , j = ret * y , n = (ret * x) * i  + (ret * y) = ret * (x*i + y) , 所以 n = ret * (x*i + y) , m =  ret * x , (n, m) = (m, j) = ret;
如果 m % j !=0 ,那就重复上一步，直到可以整除为止。

更相减损法

出自《九章算术》的一种求最大公约数的算法，它原本是为约分而设计的，但它适用于任何需要求最大公约数的场合。
《九章算术》是中国古代的数学专著，其中的“更相减损术”可以用来求两个数的最大公约数，即“可半者半之，不可半者，副置分母、子之数，以少减多，更相减损，求其等也。以等数约之。”

翻译成现代语言如下：
第一步：任意给定两个正整数；判断它们是否都是偶数。若是，则用2约简；若不是则执行第二步。
第二步：以较大的数减较小的数，接着把所得的差与较小的数比较，并以大数减小数。继续这个操作，直到所得的减数和差相等为止。
则第一步中约掉的若干个2与第二步中等数的乘积就是所求的最大公约数。
其中所说的“等数”，就是最大公约数。求“等数”的办法是“更相减损”法。所以更相减损法也叫等值算法。

例1．用更相减损术求98与63的最大公约数。
解：由于63不是偶数，把98和63以大数减小数，并辗转相减：
98-63=35
63-35=28
35-28=7
28-7=21
21-7=14
14-7=7
所以，98和63的最大公约数等于7。

• 个人解读

  ① 奇 - 奇 = 偶，偶 - 奇 = 奇；所以第一次除2以后，后面都不会出现 偶 - 偶 的情况；

  ② 任意给定两个正整数 m, n (m < n) ，假设两个数的最大公约数为 ret，
  如果 m 就是 最大公约数，n 就能被 m 一直减到 0；
  如果不能，假设 n - m = j ，找到m、j的最大公约数即可，即 (n,m) = (m,j)；

  ③在代码实现工程中，并不需要判断两个数的奇偶性，直接进行辗转相减即可。

完整代码

/*
encoding: GB2312(936)
*/

#include <iostream>
using namespace std;

void findCommonDivisor_Jx()
{
	int div_a, div_b;
	cout << "请输入两个数：" << endl;
	cin >> div_a >> div_b;
	int temp = 0;
	temp = div_a > div_b ? div_b : div_a;
	for (int i = temp; i >= 1; i--)
	{
		if (div_a % i == 0 && div_b % i == 0)
		{
			cout << "最大公约数为：" << i << endl;
			break;
		}
	}
}

void findCommonDivisor_Runoob()
{
	int n1, n2;

	cout << "输入两个整数：" << endl;
	cin >> n1 >> n2;

	while (n1 != n2)
	{
		if (n1 > n2)
			n1 -= n2;
		else
			n2 -= n1;
	}

	cout << "最大公约数为：" << n1 << endl;
}

int main()
{
	//findCommonDivisor_Jx();
	findCommonDivisor_Runoob();
	system("pause");
	return 0;
}

更多经典案例

• https://github.com/Jacob-xyb/Cpp_Note