CINTA作业二、GCD与EGCD

1、给出Bezout定理的完整证明。

由题目条件，a和b为非零整数，设d=gcd(a,b)
则可得d | a，d | b，设存在任意整数m和n
易得d | am，d | bn
由命题1.1的证明可知，
d | (am + bn)
即有am+bn = kd
则d= a*(m/k)+b*(n/d)
由m和n为任意整数
则必存在整数 r=m\k，整数s=n\d
可得d=ar+bs
故得证

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

//输入：两个整数a和b，且a>b
//输出：a和b的最大公因子
int gcd(int a,int b)
{
	int temp;
	while (b != 0)
	{
		temp = a % b;
		a = b;
		b = temp;
	}
	return a;
}

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

//输入：两个整数a和b，且a>b
//输出：a和b的最大公因子d，以及满足d=a*r+b*s的系数r和s
int* egcd(int a, int b)
{
	int* arr = new int[3];
	int r0, r1, s0, s1;
	r0 = 1, r1 = 0, s0 = 0, s1 = 1;
	int temp, temp_r, temp_s;
	while (b != 0)
	{
		int q = a / b;
		temp = a % b;
		a = b;
		b = temp;
		temp_r = r0 - q * r1, r0 = r1, r1 = temp_r;
		temp_s = s0 - q * s1, s0 = s1, s1 = temp_s;
	}
	arr[0] = r0, arr[1] = s0, arr[2] = a;
	return arr;
}

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

//输入：一个整型数组和它长度
//输出：数组中所有元素的最大公因子
int arr_gcd(int* arr,int arr_len)
{
	if (arr_len == 0) return 0;
	else if (arr_len == 1) return *arr;
	else if (arr_len == 2) return gcd(arr[0],arr[1]);
	else
	{
		int arr_temp = gcd(arr[0], arr[1]);
		for (int i = 2; i < arr_len; i++)
			arr_temp = gcd(arr_temp, arr[i]);
		return arr_temp;
	}
}