原创文章,转载请注明: 转载自 SunliyMonkey技术博客
本文链接地址: 状态转移思想解读:辗转相除(欧几里德)算法及扩展
URL: http://blog.csdn.net/sunliymonkey/article/details/48102067
1 算法实现
1.1 欧几里德算法(辗转相除法)
欧几里德算法,也被称为辗转相除法,其被用于求解两个数之间的最大公约数,常用 gcd(A,B) 进行表示。
它的算法实现十分容易,如下:
int gcd(int a,int b)
{
if(b==0)
return a;
else
return gcd(b,a%b);
}
1.2 扩展欧几里德算法
扩展欧几里德算法是用来求解这样一个式子:
A∗x+B∗y=gcd(A,B)
求解: x=? , y=?
这里的x,y不是唯一解,当得知一组可行解之后,其它的可行解都可以表示出来:
xk=x+(k∗Bgcd(A,B))
yk=y−(k∗Agcd(A,B))
k 为任意整数
在这里不阐述,x,y的值有何用处,如果你是acmer,应该发现有不少题,需要利用这个性质。
扩展欧几里德算法的实现也很简单,在原有欧几里德算法的基础上,少许额外的代码,即能实现:
int exGcd(int a,int b,int &x,int &y)
{
if(b==0)
{
x = 1;
y = 0;
return a;
// 此时a = gcd(A,B)
// a * 1 + b * 0 = gcd(A,B)
// 由于b=0,y可以任意取值
}
int d = exGcd(b, a%b, y, x); // 注意这里x,y的位置对调
y -= a/b * x;
return d; // 返回最大公约数
}
2 状态转移思想
看了上面的实现,相信你的第一感觉是十分简洁,确实也如此,对于一个不懂这两种算法的人,如果提醒你从“动态规划”,“状态转移”等角度思考这两个问题,相信会有不少人能够AC它们。
先来看看两个明显的结论:
如果A=B,gcd(A,B)=A或B 如果A=0