扩展欧几里得算法及应用

欧几里得算法

即利用辗转相除法求最大公因数

ll gcd(ll a,ll b){
	return b==0?a:gcd(b,a%b);
}

裴蜀定理(贝祖定理)

若$a,b$是整数，那么一定存在整数$x,y$，使$ax+by=gcd(a,b)$成立

重要推论：$a$和$b$互质的充要条件是存在整数$x,y$，有$ax+by=1$

扩展欧几里得算法

扩展欧几里得算法又称$exgcd$，就是利用欧几里得算法，求出贝祖定理$ax+by=gcd(a,b)$的整数解

证明：

$\because gcd(a,b)=gcd(b,a\%b),gcd(b,a\%b)=bx+(a\%b)y$

$\therefore a_1{x}_1+b_1{y}_1=gcd(a_1,b_1)=gcd(b_1,a_1\%b_1)=b_1{x}_2+(a_1\%b_1)y_2$

令$a_2=b_1,b_2=a_1\%b_1$

$\therefore a_1{x}_2+b_2{y}_2=gcd(a_2,b_2)=gcd(b_2,a_2\%b_2))=b_2{x}_3+(a_2\%b_2)y_3$

同理…

最后可化简为$gcd(a_n,0)=a_n{x}_n+0\ast y_n=a_n$

得到最后一组解为$x_n=1,y_n=0$

$\because a\%b=a-\lfloor \frac{a}{b}\rfloor \ast b$，又因为第一次我们化简的等式中有$a_1{x}_1+b_1{y}_1=b_1{x}_2+(a_1\%b_1)y_2$

$\therefore a_1{x}_1+b_1{y}_1=b_1{x}_2+(a_1-a_1/b_1\ast b_1)y_2=a_1{y}_2+b_1(x_2-a_1/b_1\ast y_2)$

$\therefore x_1=y_2,y_1=x_2-a_1/b_1\ast y_2$

然后由最后一组解回溯，即可求得$x_1,y_1$

下面的代码中是简洁版本：我们都是由下一层得到上一层的，当我们由最后一层返回倒数第二层的时候，我们把$y$传给上一层的$x$，这样由上面公式来看是正确结果，而上一层的$y$怎么求呢，上一层的$y$等于下一层的$x$减去$a/b$乘以下一层的$y$。而下一层的$y$我们已经赋值给了当前层的$x$，下一层的$x$我们已经赋值给了当前层的$y$，因此仔细思考，$y-=(a/b)\ast x$;就可以得到当前层的$y$，再一层层向上传递即可

代码之所以返回$a,b$的最大公因数是因为我们可能用到，这样就不需要再求一遍$gcd$了

ll exgcd(ll a,ll b,ll &x,ll &y){
    if(!b){
        x=1,y=0;
        return a;
    }
    ll gcd=exgcd(b,a%b,y,x);
    y-=(a/b)*x;
    return gcd;     //返回的是a,b的最大公因数
}

应用

1.求解二元一次不定方程

已知$ax+by=c$，设$r=gcd(a,b)$，则

• 若$c\%r!=0$，该方程组无整数解

• 若$c\%r==0$，我们可以先求出$ax+by=r$对应的一组整数解$(x_1,y_1)$，然后再乘以$c$并除以$r$即可得到$ax+by=c$的最小整数解$x_0=x_1\ast c/r,y_0=y_1\ast c/r$，那么其通解为$x=x_0+k\ast (b/r),y=y_0-k\ast (a/r)$，$k$取任意整数

2.求逆元

若$a\ast x$在模$p$意义下的结果为$1$，那么称$x$为$a$在模$p$意义下的逆元，$x$又写作$a^{-1}$记作$ax\equiv 1(modp)$

那么我们可以转化为求解$ax+py=1$，由于最后的求得的$x$可能是负数，因此我们再进行处理$(x\%p+p)\%p$

int exgcd(int a,int b,int &x,int &y){
    if(!b){
        x=1;y=0;
        return a;
    }
    int gcd=exgcd(b,a%b,y,x);
    y-=(a/b)*x;
    return gcd;     //返回的是a,b的最大公因数
}

int inv(int a,int p){
	int x,y;
	exgcd(a,p,x,y);
    return (x%p+p)%p;
}

3.求解线性同余方程

已知线性同余方程$ax\equiv b(modp)$ ,则该同余方程可以化为求解二元一次不定方程$ ax+py=b$，设$r=gcd(a,p)$

显然当且仅当$b\%r==0$时方程组有解

我们先求$ax+py=r$的解$x_1$，接着可以得出$ax\equiv b(modp)$的一组解为$x_0=x_1\ast b/r$，设解之间的最小间隔为$s=p/r$，且可以证明该方程组一共有$r$个解

因此$ax\equiv b(modp)$的最小整数解为$(x_0\%s+s)\%s,ax\equiv b(modp)$的通解为$x=x_0+k\ast (p/r)(k\in [0,r-1])$

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

int exgcd(int a,int b,int &x,int &y){
    if(!b){
        x=1;y=0;
        return a;
    }
    int gcd=exgcd(b,a%b,y,x);
    y-=(a/b)*x;
    return gcd;     //返回的是a,b的最大公因数
}

int solve(int a,int b,int p){
	int x,y,r,x0;
    r=exgcd(a,p,x,y);
    int s=p/r;
    x=x*b/r;
    x0=(x%s+s)%s; //最小整数解
    return x0;
}