计算机算法基础-贝祖公式(基于欧几里得算法)

写在前面

贝祖公式的推导是建立在欧几里得算法的基础之上的，关于欧几里得算法的推导可以看我之前发的这个帖子：这里

为什么需要证明贝祖公式

对于互质的性质：若 整数 a 与 b 互质, 则存在 整数 x, y, 使得 xa + yb = 1 (贝祖公式)

公式内容

两个数的最大公约数，可以用两个数的整数倍相加获得，即
$$\forall a,b\in Z,gcd(a,b)=m，则\exists x,y\in Z,使得a\ast x+b\ast y=m$$

证明

要证：$$\forall a,b\in Z,gcd(a,b)=m，则\exists x,y\in Z,使得a\ast x+b\ast y=m$$
即证：$$关于x,y的方程：k1\ast x+k2\ast y=gcd(k1,k2)有整数解$$
关于 $x,y$ 的方程 $k1\ast x+k2\ast y=gcd(k1,k2)$, 由于这是一个一元一次方程，因此一定有解，需要证明的是它一定有整数解

不妨设$\{\begin{array}{ccc}& k1=a& \\ & k2=b& \end{array}$时，$x,y$的解是$x1,y1即a\ast x1+b\ast y1=m$

再设$\{\begin{array}{ccc}& k1=b& \\ & k2=a\mathrm{m}\mathrm{o}\mathrm{d}b& \end{array}$时，$x,y$的解是$x2,y2即b\ast x2+(a\mathrm{m}\mathrm{o}\mathrm{d}b)\ast y2=m$（根据欧几里得算法$gcd(a,b)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b$）
$\because a\ast x1+b\ast y1=b\ast x2+(a\mathrm{m}\mathrm{o}\mathrm{d}b)\ast y2$
$又\because a=(a\mathrm{m}\mathrm{o}\mathrm{d}b)+b\ast \lfloor a/b\rfloor$
$\therefore a\mathrm{m}\mathrm{o}\mathrm{d}b=a-\lfloor a/b\rfloor$
$\therefore a\ast x1+b\ast y1=b\ast x2+(a-\lfloor a/b\rfloor )\ast y2=a\ast y2+b\ast (x2-\lfloor a/b\rfloor \ast y2)$
$\therefore a\ast (x1-y2)+b\ast (y1-x2+\lfloor a/b\rfloor \ast y2)=0$
$则\{\begin{array}{cc}x1-y1=0& \\ y1-x2+\lfloor a/b\rfloor \ast y2=0& \end{array}解之：\{\begin{array}{cc}x1=y1& \\ y1=x2-\lfloor a/b\rfloor \ast y2& \end{array}$

$即对于该方程，a与b的解，可以用b,a\mathrm{m}\mathrm{o}\mathrm{d}b的解表示出来$
$无限递归下去，a,b,gcd(a,b)的解，可以用m,0,gcd(m,0)表示出来$（根据欧立德算法, m既是 a,b最大公约数，递归结果）

不妨设$\{\begin{array}{cc}k1=m& \\ k2=0& \end{array}$时，$x,y的解为x_n,y_n，m\ast x_n+0\ast y_n=gcd(m,0)=m$
$解之：x_n=1,y_n=\forall R$
回溯递归即 $\{\begin{array}{cc}{x}_{n-1}=y_n& \\ y_{n-1}=x_n-\lfloor a/b\rfloor \ast y_n& \end{array}$也一定是整数，无限向上递归，

最后可以确定 $\exists x,y\in Z,使得a\ast x+b\ast y=m$