啃了好几天。。。赶紧记一下
一、欧几里得
欧几里德算法又称辗转相除法,用于计算两个整数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。。