欧几里德算法求最大公约数

欧几里德算法求最大公约数

求最大公约数,貌似最近考得有点多,写篇博客记录下。如有错误或者不足,欢迎交流讨论。
1. 原理
2. 具体代码实现

原理

欧几里德算法,又称辗转相除法,用公式表示即为gcd(a,b)=gcd(b,a mod b),其中a>b,mod为取余运算。
证明(还是来自百度):
a可以表示成a = kb + r(a,b,k,r皆为正整数),则r = a mod b
假设d是a,b的一个公约数,记作d|a,d|b,即a和b都可以被d整除。
而r = a - kb,两边同时除以d,r/d=a/d-kb/d=m,等式左边可知m为整数,因此d|r
因此d也是(b,a mod b)的公约数
因此(a,b)和(b,a mod b)的公约数是一样的,其最大公约数也必然相等,得证。

代码实现

Java实现

package internetCorporationRecruiting;

import java.util.Scanner;
/** * 欧几里德算法 * @author 小锅巴 * @date 2016年3月29日下午8:44:40 * http://blog.csdn.net/xiaoguobaf */
public class gcd {
    //从gcd(a,b)=gcd(b,a mod b)可以看出这是个递归定义,先实现递归版本
    public static int gcd_1(int a, int b){
        if(a < b){
            a ^= b;
            b ^= a;
            a ^= b;
        }

        if(b == 0)
            return a;//递归结束判断条件,若a%b等于0,则最大公约就是a
        else
            return gcd_1(b, a%b);
    }

    //非递归版本
    public static int gcd_2(int a, int b){
        if(a < b){
            a ^= b;
            b ^= a;
            a ^= b;
        }
        while(b > 0){
            int temp = b;
            b = a % b;
            a = temp;
        }
        return a;
    }
    //单元测试
    public static void main(String[] args){
        System.out.println("请输入两个整数:");
        Scanner s = new Scanner(System.in);
        int a = s.nextInt();
        int b = s.nextInt();
        System.out.println(a+"和"+b+"的最大公约数是:"+gcd_1(a, b));
        System.out.println(a+"和"+b+"的最大公约数是:"+gcd_2(a, b));
    }
}

在eclipse中的执行情况,注意每次输入一个数要回车
欧几里德算法求最大公约数_第1张图片
C实现

#include<stdio.h>
//递归版本
int gcd_1(int a, int b){
    //确保a>=b
    if(a < b){
        a ^= b;
        b ^= a;
        a ^= b;
    }
    //设a=nb+r,r = a mod b,r=0时,显然b就是公约数
    if(b == 0)
        return a;
    else
        return gcd_1(b, a % b);
}

//非递归版本
int gcd_2(int a, int b){
    if(a < b){
        a ^= b;
        b ^= a;
        a ^= b;
    }

    while(b>0){
        int Rem = a % b;
        a = b;
        b = Rem;
    }
    return a;
}
void main(){
    int a, b;
    printf("请输入两个数:\n");
    scanf("%d,%d",&a,&b);
    printf("他们的最大公约数是:%d\n",gcd_1(a,b));
    printf("他们的最大公约数是:%d\n",gcd_2(a,b));
    return 0;
}

通过CodeBlocks运行

你可能感兴趣的:(求最大公约数)