求最大公约数的方法

         转自我的新浪博客，搬家了！！！

一、更相减损术

   更相减损术是出自《九章算术》的一种求 最大公约数 的 算法。

  算法思想：

       《九章算术》是中国古代的数学专著，其中的“更相减损术”可以用来求两个数的最大公约数，即“可半者半之，不可半者，副置分母、子之数，以少减多，更相减损，求其等也。以等数约之。”

　　翻译成现代语言如下：

　　第一步：任意给定两个正整数；判断它们是否都是偶数。若是，则用2约简 (意思是如果都是偶数，则都除以2化简) ；若不是则执行第二步。

　　第二步：以较大的数减较小的数，接着把所得的差与较小的数比较，并以大数减小数。继续这个操作，直到所得的减数和差相等为止。

　　则第一步中约掉的若干个2与第二步中等数的乘积就是所求的最大公约数。

　　其中所说的“等数”，就是最大公约数。求“等数”的办法是“更相减损”法。

    还是举个实例吧，例如91和49,

91   49

49   42

42   7

35   7

28   7

21   7

14   7

7    7(此时减数和差相等)

    由于第一部没有约掉任何2，所以最大公约就是7了。同理：(102,78)--->(51,39)-->(39,12)--->(27,12)--->(15,12)--->(12,3)--->(9,3)--->(6,3)--->(3,3),第一步约掉2，最后减数和差都是3，所以最大公约数就是2*3=6。(在计算后我还以为错了呢，我以为是3，呵呵)

二、辗转相除法

     辗转相除法 , 又名欧几里德算法（ Euclidean algorithm ）乃求两个正整数之最大公因子的算法。它是已知最古老的算法 , 其可追溯至前 300 年。它首次出现于欧几里德的《几何原本》（第 VII 卷，命题 i 和 ii ）中，而在中国则可以追溯至东汉出现的《九章算术》。它并不需要把二数作质因子分解。

证明：

设两数为a、b(b＜a)，求它们最大公约数(a、b)的步骤如下：用b除a， 得a＝bq......r1(0≤r)。若r1=0，则(a，b)＝b；若r1≠0，则再用r1除b， 得b＝r1q......r2(0≤r2).若r2＝0，则(a，b)＝r1，若r2≠0，则继续 用r2除r1,……如此下去，直到能整除为止。其最后一个非零余数即为(a，b)。

辗转相除法是利用以下性质来确定两个正整数a和b的最大公因子的：

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

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

另一种写法是：

1.a÷b，令r为所得余数（0≤r＜b）。若r=0，算法结束；b即为答案。

2.互换：置a←b，b←r，并返回第一步。

实现代码：

//辗转相除法--纯循环法
 
int gcd(int a,int b)
{
    int res;
    while(b!=0) //while(b)
    {
        res = a%b;
        a = b;
        b = res;
    }
    return(a);
}

//辗转相除法---递归求解
 
int gcd(int a, int b)
{
    if(b == 0)
        return a;
    else
        return gcd(b, a%b);
}

//更相减损法
 
int gcd(int a, int b)
{
    whiel(a != b)
    {
        if(a > b)
            a -= b;
        else
            b -= a;
    }
    return a;
}

Ps:更相减损术的百度百科里竟然没有C语言的实现代码，后来我去编辑了。其实维基百科一直比百度百科好，但维基没有这词条。

pps：有些是看的某位网友的，但现在忘记是谁的了，对不起啦。