【C语言】辗转相除法和更相减损术求最大公约数与最小公倍数

1. 基本概念

1.1 约数

约数，又称因数。在实数域上，整数a除以整数b(b≠0) 除得的商正好是整数而没有余数，我们就说a能被b整除，或b能整除a。a称为b的倍数，b称为a的约数。

1.2 最大公约数

最大公约数又叫最大公因数，也称最大公约数、最大公因子，指两个或多个整数共有约数中最大的一个。

1.3 最小公倍数

两个或多个整数公有的倍数叫做它们的公倍数，其中除0以外最小的一个公倍数就叫做这几个整数的最小公倍数。若已知两个实数m和n，且g代表m和n的最大公约数，则最大公约数与最小公倍数的关系式如下：
最小公倍数=(m*n)/g

2. 辗转相除法

辗转相除法又叫欧几里得算法，它用于计算两个非负整数a，b的最大公约数。假如需要求 1997 和 615 两个正整数的最大公约数,用欧几里得算法，是这样进行的：
1997 / 615 = 3 (余 152)
615 / 152 = 4(余7)
152 / 7 = 21(余5)
7 / 5 = 1 (余2)
5 / 2 = 2 (余1)
2 / 1 = 2 (余0)
至此，最大公约数为1

3. 更相减损术

更相减损术是出自《九章算术》的一种求最大公约数的算法，它原本是为约分而设计的，但它适用于任何需要求最大公约数的场合，其原文如下：
“可半者半之，不可半者，副置分母、子之数，以少减多，更相减损，求其等也。以等数约之。”
（如果需要对分数进行约分，那么）可以折半的话，就折半（也就是用2来约分）。如果不可以折半的话，那么就比较分母和分子的大小，用大数减去小数，互相减来减去，一直到减数与差相等为止，用这个相等的数字来约分。

4. 代码实现

根据以上算法描述，使用C语言实现算法如下：

#include <stdio.h>
#include <math.h>
//辗转相除法求最大公约数
void test1(int m, int n) 
{
	//先选出最大的数
	int max = fmax(m,n);
	int min = fmin(m,n);
	int temp1 = max;//存储相除取余的临时变量
	while (temp1 != 0)
	{
		temp1 = max % min;
		max = min;
		min = temp1;
	}
	printf("%d和%d的最大公约数是%d\n", m, n, max);
	printf("%d和%d的最小公倍数是%d\n", m, n, (m * n) / max);
}

//更相减损术求最大公约数
void test2(int m, int n)
{
	int max = fmax(m, n);
	int min = fmin(m, n);
	int temp2 = max;//存储相减的临时变量
	while (temp2 != 0)
	{
		temp2 = max - min;
		if (temp2 < min)
		{
			max = min;
			min = temp2;
		}
		else
		{
			max = temp2;
		}
	}
	printf("%d和%d的最大公约数是%d\n", m, n, max);
	printf("%d和%d的最小公倍数是%d\n", m, n, (m * n) / max);
}
int main()
{
	test1(35, 49);
	test2(35, 49);
	return 0;
}