扩展欧几里得、同余方程、乘法逆元

扩展欧几里得、同余方程、乘法逆元

1.欧几里得

欧几里德算法又称辗转相除法，用于计算两个整数a,b的最大公约数 gcd(a,b)。基本算法：设 a = qb + r，其中a，b，q，r都是整数，则 gcd(a,b) = gcd(b,r)，即 gcd(a,b) = gcd(b,a%b)。

证明：
a = qb + r
如果 r = 0，那么 a 是 b 的倍数，此时显然 b 是 a 和 b 的最大公约数。
如果 r ≠ 0，任何整除 a 和 b 的数必定整除 a - qb = r，而且任何同时整除 b 和 r 的数必定整除 qb + r = a，所以 a 和 b 的公约数集合与 b 和r 的公约数集合是相同的。特别的，a 和 b 的最大公约数是相同的。

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

2.扩展欧几里得

贝祖定理：ax+by=m 有整数解当且仅当m为gcd(a,b)的倍数

扩展欧几里得算法就是在求a,b最大公约数的同时，求出贝祖等式中x,y的一组整数解

公式推导
$$\begin{gathered} 设k=m/gcd(a,b),则求出ax+by=gcd(a,b)的解后，再乘以k就得到原来的解 \\ 已知a{x}_1+b{y}_1=gcd(a,b),且gcd(a,b)=gcd(b,a\%b) \\ 则b{x}_2+(a-a/b\ast b)y_2=gcd(b,a\%b) \\ 结合两条式子，得到a{x}_1+b{y}_1=b{x}_2+(a-a/b\ast b)y_2 \\ 整理，即a{x}_1+b{y}_1=a{y}_2+b(x_2-a/b/\ast y_2) \\ 所以得到递推式x_1=y_2,y_1=x_2-a/b\ast y_2 \end{gathered}$$
递推到最后b=0，此时x=1，y=0就是一组整数解。

int exgcd(int a, int b, int &x, int &y)
{
    if(b == 0)
    {//推理1，终止条件
        x = 1;
        y = 0;
        return a;
    }
    int r = exgcd(b, a%b, x, y);
    //先得到更底层的x2,y2,再根据计算好的x2,y2计算x1，y1。
    //推理2，递推关系
    int t = y;
    y = x - (a/b) * y;
    x = t;
    return r;
}

通解：
$$\begin{gathered} x=x_0+(b/gcd)\ast t \\ y=y_0-(a/gcd)\ast t \end{gathered}$$

同余方程

$$形如ax\equiv c(modb),表示ax\%b=c\%b$$

如何与上述贝祖定理联系呢？由上述关系可设:
$$\begin{gathered} ax=k_{1}b+r \\ c=k_{2}b+r \\ 即ax=k_{1}b+(c-k_{2}b) \\ ax+b(k_2-k_1)=c \end{gathered}$$
所以在已知a，b，c的值情况下**(满足解题条件**)，可以求出x的整数解，通解

乘法逆元

$$ax\equiv 1(modb)$$

同余方程的特殊情况，此时称x为a的乘法逆元，最小正整数解