欧几里得算法 2018-10-18

欧几里得算法就是辗转相除。

欧几里得算法和辗转相除法都是求两数最小公倍数的算法。
所以，欧几里得算法（辗转相除法）的重要代码如下：

int gcd(int x,int y){
	if(y==0)
		return x;
	else
		return gcd(y,x%y);
}

欧几里得算法：可以求出两个数的最大公因数
那么，是否可以满足不定方程ax+by=gcd(a,b)？？？
这就需要用到扩展的欧几里得算法，这个算法就可以求出上面这个不定方程的一组特例解。
所以，欧几里得算法和扩展欧几里得算法的最大区别就是欧几里得算法只能求出两数的最大公因数，而扩展欧几里得算法就可以顺带求出ax+by=gcd(a,b)的解:

 x = x0 + (b/gcd)*t
 y = y0 – (a/gcd)*t

所以，扩展欧几里得的关键代码如下:

int e_gcd(int a,int b,int &x,int &y){
	if(b==0){
		x=1;
		y=0;
		return a;
	}
	int ans=e_gcd(b,a%b,x,y);
	int temp=x;
	x=y;
	y=temp-a/b*y;
	return ans;
}

开始做题吧！！！http://lxg.yali.edu.cn/contest/12
T1:这就是一个裸扩展欧几里得算法
代码:

#include <bits/stdc++.h>
using namespace std;
int a,b,x,y,ans;
int e_gcd(int a,int b){
	if(b==0){
		x=1;
		y=0;
		return a;
	}
	int ans=e_gcd(b,a%b);
	int temp=x;
	x=y;
	y=temp-a/b*y;
	return ans;
}
int main(){
	cin>>a>>b;
	e_gcd(a,b);
	cout<<(x+b)%b;
	return 0;
}

T2:直接a*b-a-b就可以A了2333扩展欧几里得算法没有推出来······

#include <bits/stdc++.h>
using namespace std;
int main(){
	long long a,b;
	cin>>a>>b;
	cout<<a*b-a-b;
}