算法笔记——拓展欧几里得定理

求 $a{x}_1+b{y}_1=\gcd (a,b)$

设 $b{x}_2+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y_2=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$

$\therefore \gcd (a,b)=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$

$\because a{x}_1+b{y}_1=b{x}_2+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y_2$

又 $\therefore a\mathrm{m}\mathrm{o}\mathrm{d}b=a-\lfloor \frac{a}{b}\rfloor \times b$

$\because a{x}_1+b{y}_1=b{x}_2+(a-\lfloor \frac{a}{b}\rfloor \times b)y_2$

我们利用乘法分配律，把这个式子展开

$b{x}_2+(a-\lfloor \frac{a}{b}\rfloor \times b)y_2=b{x}_2+a{y}_2-\lfloor \frac{a}{b}\rfloor \times b{y}_2$

然后，我们把 $b$ 提出来

$b{x}_2+a{y}_2-\lfloor \frac{a}{b}\rfloor \times b{y}_2=a{y}_2-b(x_2-\lfloor \frac{a}{b}\rfloor \times y_2)$

然后，我们发现 $x_1=y_2,y_1=x_2-\lfloor \frac{a}{b}\rfloor \times y_2$

这样，代码就好写了。

int exgcd(int a,int b,int &x,int &y){
    if(!b){
        x=1;
        y=0;
        return a;
    }
    int d=exgcd(b,a%b,x,y);
    int t=x;
    x=y;
    y=t-a/b*y;
    return d;
}