求两个或N个数的最大公约数(gcd)和最小公倍数(lcm)的较优算法

求两个或N个数的最大公约数(gcd)和最小公倍数(lcm)的较优算法

求两个或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));

}