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;
}
}