2021-03-09 欧几里得算法

欧几里得算法的个人理解

欧几里得算法可以快速地求得两个数的最大公约数，其原理大致如下：

假设有两个数a和b，则**a和b的公约数**与**b和a % b的公约数**相等。

原理证明：
设a和b的一个公约数为c，r=a%b。
我们可以令r=a-kb（k是一个整数），再对该式两边同时除以c，由于a和b都一定可以整除c，所以式子的右边一定是一个整数，因此r/c也一定是整数，所以c也是a%b和b的公约数。

• 为什么欧几里得算法可以算出最大公约数？
  上面证明的原理只能说明a和b的公约数与b和a % b的公约数相等，而要利用欧几里得算法实现快速算出最大公约数，需要我们在设计程序时采用递归的思想。
  假设一个求最大公约数的函数gcd( )，它的函数声明为int gcd(int a, int b)；，我们可以把a和b作为参数传入，然后在函数体中通过递归的设计不断地把b和a%b作为参数传入，传入的参数a和b都会不断减小（为便于理解，先假设a>b），并且a的减小速度比b的减小速度要慢，所以参数b会首先减小为0（b的值是a和b的反复求余，最后一定会等于0，如果理解不了，可以随便挑两个数算算，并不难算），而这就是递归的结束条件。

• 代码如下

#include <stdio.h>

//求最大公约数地函数
int gcd(int a, int b) {  //a和b的大小关系无要求
	if (!b) return a;//递归的结束条件：b=0
	return gcd(b, a % b); 
}

int main() {
    int a, b;
    while (~scanf("%d%d", &a, &b)) {
        printf("gcd(%d, %d) = %d\n", a, b, gcd(a, b));
    }
    return 0;
}

注：传入的参数a小于参数b也是可以的，因为如果a<b，则a%b=a，之后传入函数时，就会是a’=b，b’=a%b=a,两个参数正好交换了位置，使得a’>b’，可以正常地进行求余运算。

参考文章：
https://cloud.tencent.com/developer/article/1433267