求最大公约数和最小公倍数算法

时间:2020-12-04 00:29:45

一、求最大公约数:欧几里德算法

欧几里德算法又称辗转相除法,用于计算两个正整数a,b的最大公约数

其计算原理依赖于下面的定理:

定理:两个整数的最大公约数等于其中较小的那个数和两数相除余数的最大公约数。最大公约数缩写为gcd。

gcd(a, b) = gcd(b, a mod b)(前提条件是假设a > b 且 r = a mod b, r 不为0)


c++代码:

//两个数的最大公约数--欧几里德算法
int gcd(int a, int b){
//总是将较小的数放在b中
if(a < b){
swap(a, b);
}
if(b == 0){
return a;
}
else{
return gcd(b, a%b);
}
}


二、求最小公倍数

在已经算出整数a、b的最大公约数的基础上,我们可以通过下面的公式来求出它们的最小公倍数

lcm(a, b) = (a * b)/gcd(a, b)

//求a、b的最小公倍数
int lcm(int a, int b){
return a*b/gcd(a, b);
}


三、求n个数的最大公约数

思路:把n个数保存为一个数组,参数为数组的指针和数组的大小(需要计算的数的个数),然后先求出gcd(a[0], a[1]),然后将所求的gcd与数组的下一个元素作为gcd的参数继续求gcd,这样就产生一个递归的求ngcd的算法


c++代码:

int ngcd(int *a, int n)

{

if (n == 1) return *a;

return gcd(a[n-1], ngcd(a, n-1));

}


四、求n个数的最小公倍数

思路:算法过程和n个数的最大公约数求法类似,求出头两个的最小公倍数,再讲其和下一个元素求最小公倍数直到数组末尾,这样就产生一个地柜的求nlcm的算法


c++代码:

int nlcm(int *a, int n)

{

if (n == 1)

return *a;

else

return lcm(a[n-1], nlcm(a, n-1));

}