[数论]扩展欧几里得

导语

这是一个有关数论的知识，涉及最大公约数

正文

讲解

设 $a$ ， $b$ 是两个正整数，求 $x$ ， $y$ 满足
         ~~~~~~~~                   ~~~~~~~~         $\gcd (a,b)=ax+by(1)$
由 $b\ne 0\gcd (a,b)=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$ （辗转相除法）

设此时 $x^{\prime }$ ， $y^{\prime }$ 使得
                 ~~~~~~~~~~~~~~~~                 $b{x}^{\prime }+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y^{\prime }=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b)(2)$ (跟第二行差不多)
由 $(1)$ 和 $(2)$ 可得
                 ~~~~~~~~~~~~~~~~                 $b{x}^{\prime }+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y^{\prime }=ax+by$

注意到 $a\mathrm{m}\mathrm{o}\mathrm{d}b=a-\lfloor \frac{a}{b}\rfloor b$ （此处的括号为向下取整，$\lfloor \frac{a}{b}\rfloor$为 $a÷b$ 的商）
所以 $ax+by=b{x}^{\prime }+(a-\lfloor \frac{a}{b}\rfloor b)y^{\prime }$

注：我们注意到这个方程里面有个未知数却只有个方程，那么根据线性代数相关知识我们可以知道这东西是有无数组解的。注意到这个等式每一项只与a或b有关，所以我们只需要按照a和b分类一下这个东西就可以获得一组特解
—结尾dalao的博客

变形得
                 ~~~~~~~~~~~~~~~~                 $ax+by=a{y}^{\prime }+b(x^{\prime }-\lfloor \frac{a}{b}\rfloor y^{\prime })$

得一组解 $x=y^{\prime },y=x^{\prime }-\lfloor \frac{a}{b}\rfloor y^{\prime }$

可以发现这是一个递归的式子，当 $b=0$ 时是边界

因此，当 $b=0$ 时，显然 $x=1,y=0$ 时满足 $\gcd (a,b)=ax+by$

代码

#define pa pair <int ,int >
pa exgcd (int r1 ,int r2) {
 if (r2 == 0) return make_pair (1 ,0);
 pa tmp = exgcd (r2 ,r1 % r2);
 pa ans = make_pair (tmp.second ,tmp.first - (r1 / r2) * tmp.second);
 return ans;
}

结语

参考博客(写的真是好)

谢谢大家

——2020.9.30