欧几里德 与 扩展欧几里得 学习笔记

最近学了学扩展欧几里得，总结一下

欧几里得算法

欧几里得算法主要用来求a,b的最大公约数，又称为gcd

代码

int Gcd(int a,int b)
{
    if(b==0) return a;
    else return Gcd(b,a%b); 
} 

就是用递归来求解，当$b==0$时退出，证明也非常易懂

证明

就是证明$gcd(a,b)=gcd(b,a$ $mod$ $b)$
设$d=gcd(a,b),即d|a,d|b$
$∵$$a$ $mod$ $b=a-b*\lfloor a/b\rfloor）,d|a,d|b$
$∴d|（a-b*\lfloor a/b\rfloor）$
$又∵d|b$
$∴d|gcd(b,$ $a$ $mod$ $b)$

应用

应用就是求$a,b$的最大公约数。

扩展欧几里得算法

扩展欧几里得算法就是求一组$(x,y)$使得$ax+by=gcd(a,b)$,又称为$ex$_$gcd.$

代码

我们就是在求$gcd$的过程中对$x,y$进行计算。
当$b=0$时我们会发现$x=1,y=0$时，满足gcd的计算
然后我们再倒着推回去，证明下面；

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

证明

我们就是来证明递归的时侯$x,y$的计算式
首先我们设$ax+by=gcd(a,b)$
又因为$gcd(a,b)=gcd(b,a$ $mod$ $b)=ax+by$
把$gcd(b,a$ $mod$ b$)$代入得
$bx^1+(a-b*\lfloor a/b \rfloor )y^1=gcd(b,a$ $mod$ $b)$
去括号得$ay^1+b(x^1-y^1*\lfloor a/b \rfloor)=gcd(a,b)$
证完了

通解

我们使用扩展欧几里得算法求出来的解$x,y$只是$ax+by=gcd(a,b)$其中的一组可行解，有时候题目会要求我们求出一些有限制的解这时候我们就要根据这一组解来求出满足限制的解
通解：
我们设$x,y$是方程的解
$x^1=x+t*\frac{lcm(a,b)}{a},t\in Z$
$y^1=y-t*\frac{lcm(a,b)}{b},t\in Z$
$lcm(a,b)$表示$a,b$的最小公倍数 $lcm(a,b)=\frac {ab}{gcd(a,b)}$

应用

1. 求解不定方程$ax+by=c$

   这个看着和上面那一个差不多一样，我们先设$d=gcd(a,b)$
   我们把方程两遍同时$/gcd(a,b)$就可以得到
   $a^1x+b^1y=c^1$
   这时候我们发现因为$x,y$都是整数，所以方程有解得条件是$d|c$
   然后再同时$/c^1$得到
   $a^1\frac {x}{c^1}+b^1\frac{y}{c^1}=1$
   我们用Ex_gcd求出一组解再乘回去就可以了.

2. 求解逆元
   我们用$inv(x)$来表示x在$mod$ $p$下的逆元
   扩展欧几里得还可以用来求解$ax\equiv 1 (\pmod p)$
   转移一下得：$ax-pk=1$然后我们就可以用扩展欧几里得算法来求出x的值

有不懂的或者有错误的可以留言，谢谢