数论 欧几里德算法

欧几里得算法：欧几里德算法又称辗转相除法，是指用于计算两个正整数a，b的最大公约数。应用领域有数学和计算机两个方面。计算公式gcd(a,b) = gcd(b,a mod b)。

别称：辗转相除法

应用：计算两个正整数a，b的最大公约数

原理：gcd(a,b) = gcd(b,a mod b)      即，两个整数的最大公约数等于其中较小的那个数和两数相除余数的最大公约数。

注：gcd是最大公约数的简写，gcd(a,b)表示a和b的最大公约数，gcd(b,a mod b)表示b和a模b的最大公约数。

 

gcd函数的基本性质：gcd(a,b)=gcd(b,a)=gcd(-a,b)=gcd(|a|,|b|)

 

算法简介

欧几里德算法是用来求两个正整数最大公约数的算法。

eg：

假如需要求 1997 和 615 两个正整数的最大公约数,用欧几里德算法，是这样进行的：

当被加的数为 0 时，就得出了 1997 和 615 的最大公约数 1。

 

计算原理

其计算原理依赖于下面的定理：

定理：两个整数的最大公约数等于其中较小的那个数和两数相除余数的最大公约数。

最大公约数（Greatest Common Divisor）缩写为GCD。

gcd(a,b) = gcd(b,a mod b) (不妨设a>b 且r=a mod b ,r不为0)

 

公式证明     gcd(a,b) = gcd(b,a mod b)

gcd(a,b)=gcd(b,a mod b)

证明：a可以表示成a = kb + r，则r = a mod b

假设d是a,b的一个公约数，则有d|a, d|b，而r = a - kb，因此d|r，因此d是(b,a mod b)的公约数

假设d 是(b,a mod b)的公约数，则d | b , d |r ，但是a = kb +r，因此d也是(a,b)的公约数

因此(a,b)和(b,a mod b)的公约数是一样的，其最大公约数也必然相等，得证

 

程序设计

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

⒈ 若 r 是 a ÷ b 的余数，且r不为0， 则

gcd(a,b) = gcd(b,r)

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

另一种写法是：

⒈ 令r为a/b所得余数（0≤r<b）

若 r= 0，算法结束；b 即为答案。

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

算法版本

c语言版

/*
欧几里德算法：辗转求余
原理： gcd(a,b)=gcd(b,a mod b)
当b为0时，两数的最大公约数即为a
getchar()会接受前一个scanf的回车符
*/
#include<stdio.h>
unsigned int Gcd(unsigned int M,unsigned int N)
{
    unsigned int Rem;
    while(N > 0)
    {
        Rem = M % N;
        M = N;
        N = Rem;
    }
    return M;
}
int main(void)
{
    int a,b;
    scanf("%d %d",&a,&b);
    printf("the greatest common factor of %d and %d is ",a,b);
    printf("%d\n",Gcd(a,b));
    return 0;
}

 

递归形式：

 int gcd(int a,int b)
{
     if(b==0)
         return a;
     return 
         gcd(b,a%b);
}

简化：

int gcd(int a,int b)
{
    return b ? gcd(b,a%b) : a;
}

迭代形式：（非递归）

int Gcd(int a, int b)
{
     while(b != 0)
     {
     　　int r = b;
     　　b = a % b;
     　　a = r;
     }
     return a;
}

参考：

https://baike.baidu.com/item/%E6%AC%A7%E5%87%A0%E9%87%8C%E5%BE%B7%E7%AE%97%E6%B3%95/9002848?fr=aladdin

https://blog.csdn.net/zhj_fly/article/details/75667813

https://blog.csdn.net/w36680130/article/details/81259859