欧几里得与扩欧

欧几里得算法

$\gcd (a,b)=\gcd (b,a\%b)$

裴蜀定理

对于给定的正整数$a,b$方程$ax+by=c$有解的充要条件为c是$gcd(a,b)$的整数倍

扩展欧几里得

在求a、b最大公约数的同时，能找到整数x、y使得$ax+by=g\mathrm{cd}(a,b)$

给定$n$对整数$a_i,b_i$,对于每对数，求出一组$x_i,y_i$，使其满足$a_i\ast x_i+b_i\ast y_i=gcd\left(a_i,b_i\right)$

证明：因为$a\ast x+b\ast y=g\mathrm{cd}(a,b)$ 且$gcd(a,b)=gcd(b,a$%$b)$

$gcd(b,a\%b)$=$b\ast y+(amodb)x$

$=b\ast y+(a-a/b\ast b)x$

$=a\ast x+b\ast (y-a/b\ast x)$

#include<iostream>
#include<algorithm>
using namespace std;
int exgcd(int a,int b,int &x,int &y){
    if(!b){
        x = 1,y =0;
        return a;
    }
    int d = exgcd(b,a % b,y,x);
    y -= a / b * x;
    
    return d;
    
}
int main(){
    int t;cin >> t;
    while(t--){
        int a, b, x, y;
        cin >> a >> b;
        exgcd(a,b,x,y);
        cout << x << " " << y << endl;
    }
    
    return 0;
}

扩展欧几里得算法的几种用途

①求逆元

前提：$gcd(a,b)=1$

逆元的定义：若$a\ast x\equiv 1(modb)$且a与b互质，那么x为a的逆元

上式可以变形成为：$a\ast x=by+1$ ⇨$a\ast x+b\ast y=1$ 所以通过$exgcd(a,b,x,y)$求出的x即为a的乘法逆元

②求线性同余方程

给定$n$组数据$a_i,b_i,m_i$，对于每组数求出一个$x_i$，使其满足$a_i\ast x_i\equiv b_i(mod{m}_i)$

上式与求逆元一样变形得$a_i\ast x_i\equiv b_i+y\ast m_i$⇨$a_i\ast x_i+m_i\ast y=b$先求$ecgcd(a,m,x,y)$得到$d=gcd(a,m)$判断b是否为d倍数(裴蜀定理)如果是，此时x为$a_i\ast x_i+m_i\ast y=b$的解，所以$x\ast (b/d)$ % $m$为最终解 ，否则无解