算法笔记之欧几里得算法&扩欧算法&乘法逆元

欧几里得算法：

①欧几里得算法（辗转相除法）简介：

欧几里得算法也称为辗转相除法，常用于求解最大公约数，一般用gcd(a,b)来表示a和b的最大公约数。若a、b均为正整数，则gcd(a,b) = gcd(b,a%b)。

②欧几里得算法的代码实现：

递归式：

int gcd(int a,int b){
	if(b == 0) 
		return a;//递归边界
	else 
		return gcd(b,a % b);
}

简洁写法：

int gcd(int a,int b){
	return !b ? a: gcd(b,a % b);
}

扩展欧几里得算法（扩欧）：

①扩展欧几里得算法（扩欧）简介：

扩展欧几里得算法主要用来解决这样的问题：给定两个非零整数a和b,求一组整数解(x,y)，使得ax+ by = gcd(a,b)成立，其中gcd(a,b)表示a和b的最大公约数。

①在欧几里得算法，将gcd(a, b)转化为求解gcd(b, a%b)，当b变为0时返回a，此时的a = gcd。即欧几里得算法结束的时候变量a中存放的是gcd，变量b中存放的是0，因此此时有a1 + b0= gcd成立，此时有x=1、y=0成立。

②利用欧几里得算法的过程来计算x和y。当递归边界成立时x=1、y=0，需反推出初始的x和y。当计算gcd(a, b)时，有ax1 +by1=gcd成立，当计算gcd(b,a%b)时，又有bx2 +(a%b)y2=gcd成立。因此ax|1+by1=bx2+(a%b)y2成立。因此ax1 + by1 = ay2 + b(x2 - (a/b)y2)。因此可得递推公式：

③因此只需要在达到递归边界、不断退出的过程中根据g递推公式计算x和y,就可以得到一组解啦！！！

②扩展欧几里得算法（扩欧）代码实现：

int exgcd(int a,int b,int &x,int &y){
	if(b == 0){
		x = 1;
		y = 0;
		return a;
	}
	int g = exgcd(b,a%b,x,y);  //递归计算exgcd(b,a%b)
	int t =x;  //存放x的值
	x = y;  //更新x1=y2
	y = t-a/b*y;  //更新y1=x2-a/b*y2
	return g; //g为gcd
}

乘法逆元：

①乘法逆元简介：

设a、b、m是整数, m>1,且有ab ≡ 1(mod m)成立，那么就说a和b互为模m的逆元，一般也记作a ≡ 1/b(mod m)或b ≡ 1/a(mod m)。如果两个整数的乘积模m后等于1,就称互为m的逆元。

②乘法逆元应用------ (b/a)%m 的计算：

对乘法来说（ba)%m = (b%m)(a%m)%m成立，但对于除法来说， (b/a)%m = ((b%m)/(a%m))%m不成立，(b/a)%m = (b%m)(a%m)%m也不成立。这时就需要使用逆元来计算(b/a)%m。通过找到a模m的逆元x，即有(b/a)%m = (bx)%m成立，因此就可以把除法取模转换成乘法取模。

③求解同余式ax ≡ 1(mod m)：

一般把x的最小正整数解称为a模m的逆元。同余式ax ≡ 1(modm)是否有解取决于1%gcd(a,m)是否为0,也就是gcd(a,m)是否为1：

①如果gcd(a,m)≠1，则同余式ax ≡ 1(mod m)无解，a不存在模m的逆元。

②如果gcd(a,m)=1, 则同余式ax ≡ 1(mod m)在(0,m)上有唯一解。(因为gcd(a,m)=1，因此ax+my=1=gcd(a,m)可直接使用扩展欧几里得算法解出x之后即用(x%m+m)%m得到(0,m)范围内的解，即所需逆元。)

④使用扩展欧几里得算法求解逆元：

int inverse(int a,int m){
	int x,y;
	int g = exgcd(a,m,x,y);  //求解ax+my=1
	return (x%m+m)%m; //a模m的逆元为(x%m+m)%m
}

⑤费马小定理：

如果m为素数，且a不是m的倍数，则可使用费马小定理来得到逆元。

费马小定理：设m是素数，a是任意整数且a ≠ 0(mod m)，则a^(m-1) ≡ 1(mod m)。

使用费马小定理来推导逆元：由a^(m-1) ≡ 1(modm)可知a*a^(m-2) ≡ 1(mod m)，直接由逆元的定义可知a^(m-2)%m即a模m的逆元。

当gcd(a,m)≠1时， 扩展欧几里得算法和费马小定理都会失效，此时a模m的逆元从概念上来说不存在，但(b/a)%m仍然是有值的。

在a和m有可能不互素的情况下，可使用公式(b/ a)%m = (b%(am))/a来计算(b/a)%m的值，需要注意的是a和m的乘积可能会太大而导致溢出。因此一般尽量使用扩展欧几里得算法或者费马小定理求解，条件不成立时再采用这种办法。