CINTA作业二：GCD与EGCD

1、给出Bezout定理的完整证明。
构造集合 S={am+bn:m、n ∈ 且am+bn  0} ，显然集合S非空，由良序原则，
 取其中最小值d = ar+bs。
要证明d是a和b的公因子，即证d|a、d|b。
假设d不能整除a，由除法定理有：

, 
r只能为0，即d|a。
同理可证d|b，故d是a和b的公因子。
取a和b的任意公因子d’，有d’|a、d’|b。
因为
所以d’|d，故d是a和b的最大公因子。
所以d=gcd(a,b)=ar+bs。

2、实现GCD算法的迭代版本。

#include<iostream>
using namespace std;

int gcd(int a, int b)
{
	if (a == 0 || b == 0) return -1;
	
	while(b != 0)
	{
		a = a % b;
		int tmp = a;
		a = b;
		b = tmp;
	}
	return a;
}

int main()
{
	int a, b,n;
	cout << "请输入测试次数：" << endl;
	cin >> n;
	while (n)
	{
		cout << "请输入a和b：" << endl;
		cin >> a >> b;
		if (gcd(a, b) == -1)
			cout << "输入错误!" << endl;
		else
		{
			cout << "a和b的最大公因数为：" << gcd(a, b) << endl;
		}
		n--;
	}
	return 0;
}

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

#include<iostream>
using namespace std;

int* egcd(int a, int b)
{
	if (a < b)
	{
		int temp = a;
		a = b;
		b = temp;
	}
	int arr[3] = {0};
	int temp_r, temp_b, temp_ab, n, r0=1, r1=0, b0=0, b1=1;//n用来计算要减的倍数
	while (b != 0)
	{
		n = a / b;
		//实现a和b交换
		a %= b;
		temp_ab = a;
		a = b;
		b = temp_ab;
		//实现系数交换
		r0 -= n * r1, b0 -= n * b1;
		temp_r = r0, temp_b = b0;
		r0 = r1, b0 = b1;
		r1 = temp_r, b1 = temp_b;
	}
	arr[0] = r0, arr[1] = b0, arr[2] = a;//此时输出的a即为d

	return arr;
}

int main()
{
	int a, b;
	cin >> a >> b;

	cout << "系数r：" <<egcd(a,b)[0]<< " " << "系数s：" << egcd(a, b)[1]<< " " 
		<< "最大公因子d：" << egcd(a, b)[2]<< endl;
	return 0;
}

4、实现一种批处理版本的GCD算法，即，给定一个整数数组，输出其中所有整数的最大公因子。输入：一个整数数组a；输出：一个整数d，是a数组中所有整数的最大公因子。

#include<iostream>
using namespace std;

int mul_gcd(int arr[],int n)
{
	int d;
	if (arr[0] < arr[1])
	{
		int tmp = arr[0];
		arr[0] = arr[1];
		arr[1] = tmp;
	}
	while (arr[1] != 0)
	{
		arr[0] = arr[0] % arr[1];
		int tmp = arr[0];
		arr[0] = arr[1];
		arr[1] = tmp;
	}
	d = arr[0];
	for (int i = 2; i < n; i++)
	{
		if (d < arr[i])
		{
			int tmp = arr[i];
			arr[i] = d;
			d = tmp;
			
		}
		while (arr[i] != 0)
		{
			d = d % arr[i];
			int tmp = d;
			d = arr[i];
			arr[i] = tmp;
		}

	}
	return d;
}

int main()
{
	int n;
	cout << "请输入数组的元素个数：" << endl;
	cin >> n;
	int* p = new int[n];
for (int i = 0; i < n; i++)
	{
		cin >> p[i];
	}
	cout << "最大公因数为：" << mul_gcd(p, n) << endl;
	return 0;
}