1.欧几里得算法

算法简介

欧几里得算法又名辗转相除法，迄今为止已知的最古老的算法，距今两千三百多年，用于快速计算两个数字的最大公约数，还可以用于快速求解 a*x + b*y = 1方程的一组整数解。

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

用较大数除以较小数(小数在前也可以，第一步会自动互换位置)，再用出现的余数（第一余数）去除除数，再用出现的余数（第二余数）去除第一余数，如此反复，直到最后余数是0为止，剩余的一个数为原本两个数的最大公约数。

/*************************************************************************
	> File Name: 1.gcd.c
	> Author: 陈杰
	> Mail: 15193162746@163.com
	> Created Time: 2021年04月15日 星期四 09时35分05秒
  > 欧几里得算法求两个正整数的最大公约数
************************************************************************/
#include<stdio.h>
int gcd(int a, int b) {
    return b ? gcd(b, a%b) : a;
}
int main(){
    int n, m;
    scanf("%d %d", &n, &m);
    printf("%d\n", gcd(n, m));
    return 0;
}

算法证明

1. 证明gcd(a, b)与gcd(b, a%b)等价

设：d为a、b的公约数，及d = gcd(a, b)；
  则有 a = k1 * d; b = k2 *d;

  c = a % b
  a = b + c

  k1*d = k2*d + c
  c = (k1-k2) * d;

  (a%b) = (k1 - k2) * d;
故此，d也是 a % b的公约数，则d 也是 gcd(b, a%b)的值

2. 证明这个公约数是最大公约数

设: gcd(x - ky, y) = d;
  y = nd;
  x - ky = md
  x = knd + md = d(kn + m);
重新表示a, b，则有 a = cd(kn + m), b = cdn;
则可得 gcd(a, b) >= cd，又因为gcd(a, b) = c ，所以d = 1;

欧几里得算法扩展

求解a*x + b*y = 1方程的一组整数解

1. 特殊情况：b = 0;
  {a = 1; x = 1; y任意}
设公约数为 d
d(ax + by) = 1;

f(1)成立 b = 0, x = 1且y = 0;

a*x0 + b*y0 = gcd(a,b);
a*x + b*y = 1;
> 下一项
b * x1 + a%b * y1 = 1;
b * x1 + (a- kb) * y1 = 1;
(x1-k*y1) * b + a * y1 = 1;

x = y1;
y = x1 - ky1;

/*************************************************************************
  > File Name: 1.gcd_ex.c
  > Author: 陈杰
  > Mail: 15193162746@163.com
  > Created Time: 2021年04月15日 星期四 11时40分13秒
  > 求 ax + by = 1；的一组整数解
************************************************************************/
#include <stdio.h>
int ex_gcd(int a, int b, int *x, int *y) {
    if(!b) {
        *x = 1;
        *y = 0;
        return a;
    }
    int ret = ex_gcd(b, a%b, y, x);
    *y -= a / b * (*x);
    return ret;
}
int main() {
    int a, b, x, y;
    scanf("%d %d", &a, &b);
    printf("ex_gcd(%d, %d) = %d\n", a, b, ex_gcd(a, b, &x, &y));
    printf("%d * %d + %d * %d = %d\n",a, x, b, y, a * x + b * y);
    return 0;
}