于是,我们在完成神奇的前两步之后,来到了这个神奇的地方——欧几里得算法和扩展欧几里得算法。
那么,这是用来干什么的算法呢?
算最大公约数(GCD)~~~
好吧,考虑到有一些同学可能还不知道这是怎样一种神奇的东西,那么我就把这个东西的定义放到下面来:
最大公因数,也称最大公约数、最大公因子,指两个或多个整数共有约数中最大的一个。a和b的最大公约数可以用gcd(a,b)表示
嗯嗯,就是这样子的
那么,我们的欧几里得算法又是怎么算这样一个神奇的东西的呢?
它有一个更加生动的名字叫做辗转相除法
,因为,在我们实际执行的时候就可以发现,它其实就是这样除来除去的~~~
那么,我们在做除法的过程中呢,选择的是带余除法,说的简单粗暴一点,就是计算的结果是商和余数。
现在,欧几里得算法的过程大概如下
设两数为a、b(a>b),求a和b最大公约数(a,b)的步骤如下:用a除以b,得a÷b=q……r1(0≤r1)。若r1=0,则(a,b)=b;若r1≠0,则再用b除以r1,得b÷r1=q……r2 (0≤r2).若r2=0,则(a,b)=r1,若r2≠0,则继续用r1除以r2,……如此下去,直到能整除为止。其最后一个为被除数的余数的除数即为(a, b)。
当然,我们也可以简化出这样一种神奇的式子来总结一下gcd(a,b)=gcd(b,a mod b)
嗯,这就是我们欧几里得算法的精髓(其实就这一句话,别的都是把你搞晕过去的)
证明如下:(借鉴了一点百度百科)
设两数为a、b(a>b),用gcd(a,b)表示a,b的最大公约数,r=a mod b 为a除以b的余数,k为a除以b的商,即a÷b=k.......r。
第一步:令c=gcd(a,b),则设a=mc,b=nc
第二步:根据前提可知r =a-kb=mc-knc=(m-kn)c
第三步:根据第二步结果可知c也是r的因数
第四步:可以断定m-kn与n互质(假设m-kn=xd,n=yd (d>1),则m=kn+xd=kyd+xd=(ky+x)d,则a=mc=(ky+x)cd,b=nc=ycd,则a与b的一个公约数cd>c,故c非a与b的最大公约数,与前面结论矛盾),因此c也是b与r的最大公约数。
从而可知gcd(b,r)=c,继而gcd(a,b)=gcd(b,r)。
证毕。
以上步骤的操作是建立在刚开始时r!=0的基础之上的。即m与n亦互质。
那么,我们的代码如何实现呢?
int gcd(int a,int b){
if(b==0)return a;//显然,这就是递归的边界
else return gcd(b,a%b);//这也很显然(其实那个else是可以删掉的)
}
这样,gcd函数的返回值就是a和b的最大公约数的值了。
那么,我们来说说扩展欧几里得是个什么神奇的东西
首先,扩展欧几里得算法是用来干什么的呢?
其实,是解一种神奇的方程的—>ax+by=gcd(a,b)
在这里,a和b是常数,我们需要求出x,y的值使之成立。啊………
那么,它的用处也很明显了:求解线性模方程和方程组。
(当然这只是大多数用处)
那么,代码如下
#include<bits/stdc++.h>
using namespace std;
int exgcd(int a,int b,int &x,int &y){
if(!b){
x=1;y=0;
return a;
}
int rec=exgcd(b,a%b,y,x);
y-=a/b*x;
return rec;
}
int main(){
int i,j,k,m,n,a,b,x,y;
cin>>a>>b;
cout<<exgcd(a,b,x,y)<<endl;
cout<<x<<" "<<y<<endl;
return 0;
}
显然,神奇的代码就在这里了,由于和朴素的欧几里得算法没有什么新的难的地方(当然朴素的也没有什么难的),证明过程和详细的注释就省略了~
好了,我们 LG的数学计划 第三步就这么愉快的结束了,拜拜~~
第四步 神奇的传送门(当然现在这是我自己的了)
:)
当然你也可以选择访问LYQ的文章传送到神奇的LYQ
大家下一步再见咯~