求两个或N个数的最大公约数(gcd)和最小公倍数(lcm)的较优算法
//两个数的最大公约数--欧几里得算法
int gcd(int a, int b)
{
if (a < b)
swap(a, b);
if (b == 0)
return a;
else
return gcd(b, a%b);
}
//n个数的最大公约数算法
//说明:
//把n个数保存为一个数组
//参数为数组的指针和数组的大小(需要计算的数的个数)
//然后先求出gcd(a[0],a[1]), 然后将所求的gcd与数组的下一个元素作为gcd的参数继续求gcd
//这样就产生一个递归的求ngcd的算法
int ngcd(int *a, int n)
{
if (n == 1) return *a;
return gcd(a[n-1], ngcd(a, n-1));
}
//两个数的最小公倍数(lcm)算法
//lcm(a, b) = a*b/gcd(a, b)
int lcm(int a, int b)
{
return a*b/gcd(a, b);
}
//n个数的最小公倍数算法
//算法过程和n个数的最大公约数求法类似
//求出头两个的最小公倍数,再将欺和大三个数求最小公倍数直到数组末尾
//这样产生一个递归的求nlcm的算法
int nlcm(int *a, int n)
{
if (n == 1)
return *a;
else
return lcm(a[n-1], nlcm(a, n-1));
}
posted on 2010-10-02 14:20
Vontroy 阅读(3541)
评论(1) 编辑 收藏 引用 所属分类:
数论