拓展欧几里得算法其实也不难,主要是要去记公式、代码、必要的话连推导也记下来。
好了,不扯了,将正话。
概念: 要求解这样的方程:
推导: 对于这个我们其实就是利用辗转相除法,我们可以知道,我们辗转相除法的边界是a=d,b=0,(a和b为要求最大公约数的两个数,d为他们的公约数),此时我们可以知道a就是最大公约数,我们还可以知道,在这时,一定有个解是x=1,y=0,即a * 1+b * 0 =d,(a此时等于d,b等于0)。这个有什么用呢?
我们知道gcd(a,b)=gcd(b,a%b),如果我们可以推导出每一次的解x和y,与相除后的解x1和y1的关系我们就可以算出其中的一个接了,(x和y相当于是a和b的解,x1和y1是a已经变成了b,b变成了a%b时的解),我们下面推导一下,x、y和x1、y1之间有什么关系
这样我们就可以得到x、y和x1、y1之间的关系了,所以我们只要在回溯的时候计算一下就好了。
代码:
#include<iostream>
#include<cstdio>
using namespace std;
int x,y;
void gcd(int a,int b)
{
if(b==0)//边界
{
x=1;
y=0;
return ;
}
gcd(b,a%b);
int z=y; //因为x和y都会改变,所以我们开个变量储存
y=x-(a/b)*y;//直接代入我们的公式
x=z;
}
int main()
{
int a,b;
scanf("%d %d",&a,&b);//要求a*x+b*y=gcd(a,b)
gcd(a,b);
printf("%d %d",x,y);//输出
return 0;
}
希望大家看了这篇题解后可以明白拓展欧几里得,如果有不会的欢迎留言。