The greatest common divisor gcd(最大公约数)

文章某些部分有借鉴大神处,共勉之。

说到最大公约数的算法,最熟悉的还是辗转相除法,又名欧几里德算法(Euclidean algorithm)。在求最大公数的时候,辗转相除法是利用以下性质来确定两个正整数 和 的最大公因子的:

1. 若 是 ÷ 的余数, 则gcd(a,b) = gcd(b,r)

2. a 和其倍数之最大公因子为 a

另一种写法是:

1. a ÷ b,令r为所得余数(0rb),若 r = 0,算法结束;即为答案。

2. 互换:置 abbr,并返回第一步。

第一种写法是利用递归实现,第二种可以更清楚的看出其辗转过程。实现代码如下:

int gcd(int m,int n)   //递归 
{
    if(n==0)   
      return  m;
    else 
      return  gcd(n,m%n);
}
//调用:  n=gcd(a,b);

int gcd(int m,int n)   //递归 
{
    if(n==0)   
      return  m;
    else 
      return  gcd(n,m%n);
}
//调用:  n=gcd(a,b);

今天第一次了解另一种求最大公约数算法——更相减损法(等值算法),其实个人认为也可理解为辗转相减法。。。

此处借用百科解释一下其原理:《九章算术》是中国古代的数学专著,其中的“更相减损术”可以用来求两个数的最大公约数,即“可半者半之,不可半者,副置分母、子之数,以少减多,更相减损,求其等也。以等数约之。”
翻译成现代语言如下: 
第一步:任意给定两个正整数,判断它们是否都是偶数。若是,用2约简;若不是执行第二步;
第二步:以较大的数减较小的数,接着把所得的差与较小的数比较,并以大数减小数。循环执行这个操作,直到所得的减数和差相等为止。
则第一步中约掉的若干个2与第二步中等数的乘积就是所求的最大公约数。
其中所说的“等数”,就是最大公约数。求“等数”的办法是“更相减损”法。

代码如下:

int gcd(int m,int n)   //更相减损法 
{
    while(m!=n)     //辗转相减
    {
      if(m>n)
        m-=n;
       else
         n-=m; 
    } 
    return m;
} 
实现代码如下,由于本人属于小水一名,简单的东西容易忘记,代码中附加了几种变量交换方式 ,共勉之。

#include 
#include 
using namespace std;

int gcd(int m,int n)   //递归 
{
    if(n==0)   
      return  m;
    else 
      return  gcd(n,m%n);
}

/*
int gcd(int m,int n)   //辗转相除法 
{
    int p;
    while(n)
    {
      p=m%n;
      m=n;
      n=p;
    }
    return m;
}

int gcd(int m,int n)   //更相减损法 
{
    while(m!=n)
    {
      if(m>n)
        m-=n;
       else
         n-=m; 
    } 
    return m;
} 
 */

void swap(int &m,int &n)   //引用 ,交换变量 (无返回值)
{
     int temp;
     temp=m;
     m=n;
     n=temp;     
 } 

/*void swap(int *m,int *n)  //指针,交换变量 
{
     int temp;
     temp=*m;
     *m=*n;
     *n=temp;
 } 
 调用:swap(&m,&n); 
 
 int swap(int &m,int &n)    //不用辅助变量,交换 (此两种方法指针通用)
 {
     m=m+n;
     n=m-n;
     m=m-n; 
 } 
 调用:swap(m,n);
  
 int swap(int &m,int &n)
 {
     m=m^n;
     n=n^m;
     m=m^n;
 } 
 调用:swap(m,n);
 */
 
int main()
{
     int a,b;
     
     while(1)
     { 
     printf("please enter the two positive integers: ");
     scanf("%d%d",&a,&b);
     if(a


你可能感兴趣的:(算法,gcd,算法)