数论-扩展欧几里得

欧几里得算法

经典的欧几里得算法（辗转相除法） 即GCD算法

long long gcd(long long a,long long b)
{
    return b==0? a:gcd(b,a%b);//辗转相除求GCD（a,b）
}

扩展欧几里得算法

那么什么是扩展欧几里得算法呢？
简单来理解 在GCD的过程中
我们找到了一些不定方程的解

（裴蜀定理）

---

如此以来我们不难得到
ax+by=gcd(a,b)
//
b(x1)+(a%b)(y1)=gcd(b,a%b)
//
ax+by=b(x1)+(a%b)(y1)
//
a%b=a-floor(a/b)b
//
ax+by=b(x1)+(…)(y1)
//
ax+by=ay1+b((x1)-floor(a/b)(y1))

由上面的转换不难得到x和y的递推式

x=(y1)
y=(x1-floor(a/b)*y1)

base case a=1 b=0
所以有如下的exgcd

pair<int,int> exgcd(int a,int b)
{
    if(b==0) return make_pair(1,0);
    pair<int,int> tmp=exgcd(b,a%b);
    pair<int,int> ans=make_pair(tmp.second,tmp.first-a/b*tmp.second);
    return ans;
}

扩展欧几里得算法的应用

1. 求解不定方程
2. 求解乘法逆元
3. 求解线性同余方程