CINTA作业二：GCD与EGCD

1.证明Bezout定理

证明 构造集合

S = { am − bn : m, n ∈ Z 且 am + bn ≥ 0 } .

显然，集合 S 非空，根据良序原理，取其中最小值 d = ar + bs 。然后证明两个属性： d

是 a 和 b 的公因子；如果存在 a 和 b 的公因子 d ′ ，则 d ′ | d 。

（1）证明d 是 a 和 b 的公因子：                                                                                                    根据除法定理，有a = q d + r ， 0 < = r < d ，a = qd+r，0<=r<da=qd+r，0<=r<d，q ∈ Z ，       则有r =a mod d= a − q d = a − q ( a x + b y ) = ( 1 − q x ) a + ( − q y ) b ∈ S ，                           又因为d是S中的最小值，所以r = 0，则a|d,同理b|d。                                                                   所以d为a和b的公因子。

(2) 如果存在 a 和 b 的公因子 d′，则 d′ | d:

设a=xd',b=yd';                                                                                                                                     则有d=ar+bs=xd'r+yd's=(xr+ys)d'，所以d'|d.

所以d = gcd(a, b)。

2.实现GCD算法的迭代版本

int gcd(int a, int b) {
	int temp;
	while (b) {
		temp = a % b;
		a = b;
		b = temp;
	}
	return a;
}

（3）实现EGCD算法。输入：a、b两个整数，输出：r、s、d三个整数，满足ar + bs =d。

#include<iostream>
using namespace std;
int mat[250][4];
void egcd(int a, int b, int* r, int* s, int* d) {
	int t;
	if (a < b) {
		t = a;
		a = b;
		b = t;
	}
	mat[0][0] = 1; mat[0][2] = a;
	mat[1][1] = 1; mat[1][2] = b;
	mat[2][0] = 1 - 0; mat[2][1] = 0 - a / b; mat[2][2] = a % b;
	t = 2;
	while (mat[t - 1][2] % mat[t][2]) {
		t++;
		mat[t][0] = mat[t - 2][0] - (mat[t - 2][2] / mat[t - 1][2]) * mat[t - 1][0];
		mat[t][1] = mat[t - 2][1] - (mat[t - 2][2] / mat[t - 1][2]) * mat[t - 1][1];
		mat[t][2] = mat[t - 2][2] % mat[t - 1][2];
	}
	*r = mat[t][0];
	*s = mat[t][1];
	*d = mat[t][2];
}
int main() {
	int a, b, r, s, d;
	cin >> a >> b;
	egcd(a, b, &r, &s, &d);
	cout << "r=" << r << " " << "s=" << s << " " << "d=" << d << endl;
	return 0;
}

4.批处理版本的GCD算法

#include<iostream>
using namespace std;
int gcd(int a, int b) {
	int temp;
	while (b) {
		temp = a % b;
		a = b;
		b = temp;
	}
	return a;
}
int bgcd(int a[], int len) {
	int temp = a[0];
	for (int i = 1; i < len; i++) {
		temp = gcd(temp, a[i]);
	}
	return temp;
}
int main() {
	int a[6] = { 56,48,16,64,800,8 };
	int len = sizeof(a) / sizeof(int);
	cout << "d=" << bgcd(a, len) << endl;
	return 0;
}