数论 欧几里得和扩展欧几里得

啃了好几天。。。赶紧记一下

一、欧几里得

欧几里德算法又称辗转相除法，用于计算两个整数a,b的最大公约数。记这个最大公约数c=gcd（a，b）

算法就是a%b=c，如果b不等于0，就a=b，b=c，然后一直继续

辗转相除法（摘自百度百科)
辗转相除法：辗转相除法是求两个自然数的最大公约数的一种方法，也叫欧几里德算法。
例如，求（319，377）：
∵ 319÷377=0（余319）
∴（319，377）=（377，319）；
∵ 377÷319=1（余58）
∴（377，319）=（319，58）；
∵ 319÷58=5（余29）
∴ （319，58）=（58，29）；
∵ 58÷29=2（余0）
∴ （58，29）= 29；
∴ （319，377）=29。

1.1代码

// 这是循环的
while(b!=0)
{
	c=a%b;
	a=b;
	b=c;
}

// 这是递归的
#include<iostream>
using namespace std;

int gcd(int a,int b)
{
	if(!b) return a;
	return gcd(b,a%b);
}
int main(){
	int a,b;
	cin>>a>>b;
	int c=gcd(a,b);
	cout<<"a和b的公约数是"<<c<<endl;
	return 0;
} 

1.2证明

只看了第一个证明，第二个证这个公约数是最大的还没看emmm，之后再补咯。

好了欧几里得就差不多先这样吧。。继续啃扩展

二、扩展欧几里得

除了计算a、b两个整数的最大公约数，此算法还能找到整数x、y（其中一个很可能是负数）。通常谈到最大公因子时, 我们都会提到一个非常基本的事实: 给予二整数 a 与 b, 必存在有整数 x 与 y 使得ax + by = gcd(a,b)——贝祖定理

有两个数a,b，对它们进行辗转相除法，可得它们的最大公约数——这是众所周知的。然后，收集辗转相除法中产生的式子，倒回去，可以得到ax+by=gcd(a,b)的整数解。

贝祖定理：若存在a、b是整数，则必存在整数x、y，满足ax+by=gcd(a,b)。

参考这儿的

证明：
　　　　当 b=0 时，gcd(a,b)=a，此时 x=1 , y=0
　　　　当 b!=0 时，
　　　　　　设 ax1+by1=gcd(a,b)=gcd(b,a%b)=bx2+(a%b)y2
　　　　　　又因 a%b=a-a/b*b
　　　　　　则 ax1+by1=bx2+(a-a/b*b)y2
　　　　　　　　  ax1+by1=bx2+ay2-a/b*by2
　　　　　　　　　　 =ay2+bx2-b*a/b*y2
　　　　　　　　　　 =ay2+b(x2-a/b*y2)
　　　　　　解得 x1=y2 , y1=x2-a/b*y2
　　　　因为当 b=0 时存在 x , y 为最后一组解
　　　　而每一组的解可根据后一组得到
　　　　所以第一组的解 x , y 必然存在
证毕。

（重点）解得 x1=y2 , y1=x2-a/b*y2
既然最后一组x和y有了，就可以借助欧几里得算法进行回溯，然后求得最开始的x和y。

康康代码应该能懂吧？？！

#include<iostream>
using namespace std;

int exgcd(int a,int b,in &x,int &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	// 重点是这里，首先这个最后一组exgcd求得的x=1和y=0会返回到这里
	int r=exgcd(b,a%b,x,y);
	//然后将返回的x（其实就是上面证明里面的x2）保存到t
	int t=x;
	
	x=y;	//然后x1=y（就是证明里面的y2）
	y=t-a/b*y;	//后面就看得懂了吧，就是y1=x2-a/b*y2
	//一直回溯到递归开始的地方，就可以求出初始的x和y了，并且返回最大公约数r
	return r;
}

int main(){
	int a,b,x,y;
	cin>>a>>b;
	int c=exgcd(a,b,x,y);
	cout<<"a和b的最大公约数"<<c<<endl;
	cout<<"x="<<x<<endl;
	cout<<"y="<<y<<endl;
	return 0;
}

应用

参考大神
扩展欧几里德算法的应用主要有以下三方面：

（1）求解不定方程；

（2）求解模线性方程（线性同余方程）；

（3）求解模的逆元；

对于（1）用上面的代码就可以求出a和b的最大公约数c，以及c的线性组合，也就是不定方程ax+by=c的一组解。
后面的（2）（3）用到再补充吧。。。

好了看了这么多就为了一道青蛙约会题，emmmm。。