大家一定都知道欧几里得算法gcd(a,b),它能够在O(lg(ab))的时间内求出整数a和b的最大公约数,是一个非常优秀的算法。这么一个优秀的算法,能不能在不增加时间复杂度的情况下,让它能处理更多的信息,以便解决更多的问题呢?今天,我们来介绍一下
扩展欧几里得算法。
算法解决对象
扩展欧几里得算法可以用来求形如ax+by=GCD(a,b)的方程的
整数解。(当然,a,b是整数)
算法证明
首先,一定存在解符合题意。(我也不知道为什么?)
其次,因为GCD(a,b)=GCD(b,a%b),所以x*a+y*b=GCD(a,b)=GCD(b,a%b)=p*b+q*a%b=p*b+q*(a-a/b*b)=q*a+(p-a/b*q)*b. 说成人话就是:
假设方程 b*x+a%b*y=GCD(b,a%b) 的解为x0,y0;那么原方程的解就是 x=y0,y=x0-(a/b)*y0.
基本情况(Base Case):当b=0时,x=1,y=0.
算法实现
#include<bits/stdc++.h>
using namespace std;
int extended_gcd(int a, int b, int &x, int &y){
int res,x0,y0;
if (!b){
x=1;
y=0;
return a;
}
res=extended_gcd(b,a%b,x0,y0);
x=y0; y=x0-a/b*y0;
return res;
}
int main(){
int a,b,x,y,z;
scanf("%d%d", &a, &b);
z=extended_gcd(a,b,x,y);
printf("%d %d %d\n",z,x,y);
return 0;
}
/*
模板名称:扩展欧几里得
模板作用:求方程ax+by=GCD(a,b)的一组整数解
*/