参考图书:《程序设计算法基础》梁冰 第10章初等数论 10.2最大公约数和扩展欧几里得算法
最大公约数(Greatest Common Divisor,GCD)为几个整数的公约数中最大的一个,一般求解最大公约数的算法有列举法、质因数分解法和欧几里得算法等。
其中,效率较高且在算法竞赛中最常用欧几里得算法求最大公约数。
欧几里得算法
欧几里得算法(Euclidean algorithm)又称辗转相除法,用于计算两个正整数a、b的最大公约数。
定理:gcd(a,b) = gcd(b, a mod b) (b != 0)
写成分段函数的形式:
则欧几里得算法的递归求解方式如下:
//求两个非负整数的最大公约数
int GCD(int a, int b)
{
return b ? GCD(b, a % b) : a;
}
非递归形式:
int GCD(int a, int b)
{
while(b != 0)
{
swap(a, b);
b %= a;
}
return a;
}
在求多个整数的最大公约数的情况下,先求出前两个数的最大公约数,在用得到的最大公约数与第三个数进行最大公约数的运算,即可得到前三个数的最大公约数,用此方法不断迭代下去即可得到这些数的最大公约数。
两个整数的最小公倍数(Least Common Multiple,LCM)的计算可以利用最大公约数的结果,证明如下:
程序如下:
//求解两个非负整数的最小公倍数
int LCM(int a, int b)
{
return a / gcb(a,b) * b;//为防止a*b超过整数类型范围,先将a的最大公约数约掉
}
扩展欧几里得算法
扩展欧几里得算法(Extended Euclidean algorithm)是用来咋已知a,b的条件下求解一组x,y使得(根据数论中的相关定理,解一定存在),其中a、b、x、y都为整数。将这样的方程称为线性同余方程。
求解代码:
//求解使ax + by = gcd(a,b)成立的x,y的解,返回 a 和 b 的最大公约数
int EX_GCD(int a, int b, int &x, int &y)
{
if(b == 0)
{
x=1,y=0;
return a;
}
int g = EX_GCD(b, a % b, y, x);
y -= a/b * x;
return g;
}