【数论】扩展欧几里得

欧几里得算法用于求解最大公因数
主要原理：

int gcd(int a,int b){return b?gcd(b,a%b):a;}

扩展欧几里得算法用于求解二元一次方程

首先了解一下裴蜀定理：$ax+by=c$有解当且仅当$gcd(a,b)|c$
那么我们就可以求解$ax+by=gcd(a,b)$

根据这个原理递归求解即可

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

根据这个求得的$x，y$只是一组可行解

如果$gcd(a,b)=1$,那么如果得到一组解$x_0,y_0$，通解就可以表示为$x=x_0+b\times t$,$y=y_0-a\times t$
如果需要求最小整数解x，我们可以令$t=\frac{　b　 }{(a　,　b)}$,那么$x=(x_0\%t+t)\%t$即可