一:欧几里得算法
1,欧几里德算法又称为辗转相除法,主要用于计算两个整数a,b的最大公约数。
2,原理:
//递归写法
int gcd(int a,int b)
{
if(b==0)
return a;
return gcd(b,a%b);
}
//非递归的
int gcd(int a,int b)
{
while(b)
{
int t=a;
a=b;
b=t%b;
}
return a;
}
二:欧几里得扩展原理:
1,对于不完全为 0 的非负整数 a,b,d=gcd(a,b)表示 a,b 的最大公约数,必然存在整数对 x,y ,使得d=ax+by成立。欧几里得扩展原理就是求x,y的算法。
2,推倒:
(1) 若 d=0 很显然 x=1 y=0 ;使得 ax+by=d=a。
(2) 若 d!=0
由上面的欧几里得算法可以得到 gcd(a,b)=gcd(b,a%b); 于是我们假设存在两对实数(x1,y1) 与 (x2,y2) ,使得 :
x1*a+y1*b=gcd(a,b) , x2*b+y2*(a%b)=gcd(b,a%b) 成立。
所以有 x1*a+y1*b = x2*b+y2*(a%b) ,而 a%b=(a-(a/b)*b)代入原式中得到 :
x1*a+y1*b = x2*b+y2*((a-(a/b)*b))=y2*a +(x2- (a/b)*y2)*b 所以:
可以得到递推公式:
int ex_gcd(int a,int b,int &x,int &y)
{
if(b==0)
{
x=1; y=0;
return a;
}
int d=ex_gcd(b,a%b,y,x);
y=y-(a/b)*x;
return d;
}