扩展欧几里得详解

给出a、b、c，，给出$a\ast x+b\ast y=c$中，x、y满足一定条件的解。
这时我们就要用到扩展都几里德了。
首先要知道一个性质：若$a\ast x+b\ast y=c$有解，则有$c=k\ast gcd(a,b),k\in Z$，这一点可以利用唯一分解定理证得。
利用上面的性质，我们将等式两边同除于$k$得到：      $a\ast x$₁$+b\ast y$₁$=gcd(a,b)$  $(1)$
用欧几里得得到：                      $b\ast x$₂$+a\%b\ast y$₂$=gcd(b,a\%b)$
一直递归重复计算，终点是： $b\ast x$_n$+0\ast y$_n$=gcd(b,0)=b$
得一组解：$x_n=1$ ,$y_n=0$ ($y_n$可以回溯到前两组解证明确实是0)
接下来我们要利用回溯：
我们知道 ：$a\%b=a-[a/b]\ast b$
得到：$b\ast x$₂ + { a − [ a / b ] ∗ b } ∗ y +\{a-[a/b]*b\}*y +{a−[a/b]∗b}∗y₂$=gcd(b,a\%b)$
化简得: b ∗ { x b*\{x b∗{x₂$-[a/b]\ast y$₂$\}+a\ast y$₂$=gcd(b,a\%b)$ $(2)$
联立$(1)(2)$得$x_1=y_2$，$y_1=x_2-[a/b]\ast y_2$
这样我们就可以通过$x$_i-1,$y$_i-1得到：$x$_i,$y$_i
回溯到最终时，我们会得到$x$₁,$y$₁，但是别忘了这是我们通过对等式作了改变得到的，我们还需要将结果在$\ast k$。
另外，我们还有一直通解：
$x$_通=$x+b/gcd(a,b)$ ，$y$_通=$y-a/gcd(a,b)$

模板

typedef long long LL;
LL gcd(LL a,LL b){
	return b==0 ? a : gcd(b,a%b);
}

LL exGcd(LL &x,LL &y,LL a,LL b){
	if(b==0){
		x=1,y=0;
		return a;
	}
	LL GCD=exGcd(x,y,b,a%b);
	int t=x;
	x=y;
	y=t-a/b*y;
	return GCD;
}

再贴一道题 洛谷1082
题意：求$a\ast x-y\ast b=1$时,$x$的最小正整数解

ac代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
LL exgcd(LL &x,LL &y,LL a,LL b){
    if(b==0){
        x=1;y=0;
        return a;
    }
    LL GCD=exgcd(x,y,b,a%b);
    LL t=x;
    x=y;
    y=t-(a/b)*y;
    return GCD;
}
int main(){
    LL a,x,b,y,c;//a*x+b*y=c;
    cin>>a>>b;
    c=1;
    if(c%__gcd(a,b)==1) cout<<"no"<<endl;
    else {
        LL GCD=exgcd(x,y,a,b);
        int k=c/GCD;
        x*=k;
        y*=k;
        int tx=b/GCD,ty=-a/GCD;
        while(x<0) x+=tx,y+=ty;
        cout<<x<<endl;
    }
    return 0;
}