最大公约数的求法在大一时候就学过了,但是当时的算法是效率很低的。然而还有一个简单、高效,而且相当优美的算法——辗转相除法,也就是欧几里得算法
先了解一个定理——GCD递归定理:
对任意非负整数a和任意正整数b,gcd(a,b) = gcd(b,a mod b)
这个定理是这个算法的核心部分。
这个定理的证明就不证明了,记住即可。
所以这下来求gcd(a,b)就十分简单了,找到递归基:当b=0 , gcd(a,0)=a 。
那么算法就是:
#include<cstdio>
using namespace std;
int gcd(int a,int b)
{
return b==0 ? a : gcd(b , a % b);
}
int main()
{
int a,b;
scanf("%d%d",&a,&b);
printf("%d\n",gcd(a,b));
return 0;
}
那么求a和b的最小公倍数呢?
求出最大公因数就很简单了:a*b/gcd(a,b);