状态转移思想解读：辗转相除（欧几里德）算法及扩展

原创文章，转载请注明: 转载自 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不是唯一解，当得知一组可行解之后，其它的可行解都可以表示出来：
　　　　　　　
　　　　　　　　$x^k = x + \left( k*\frac{B}{gcd(A, B)} \right)$
　　　　　　　　$y^k = y - \left( k*\frac{A}{gcd(A, B)} \right)$
　　　　　　　　$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