扩展欧几里德算法

设$gcd(a,b)=1,求ax+by=1,a>b$，找出满足条件的$x,y$。
递归反推思路
$$\begin{array}{rl}1& =a{x}_0+b{y}_0\\ & =b{x}_1+a\%b{x}_2\\ & =...\\ & =1x_n+0y_n\end{array}$$
$$\begin{gathered} a{x}_0+b{y}_0=b{x}_1+a\%b{y}_1\to \\ a{x}_0+b{y}_0=b{x}_1+(a-a/b\ast b)y_1\to \\ a{x}_0+b{y}_0=a{y}_1+b(x_1-a/b{y}_1)\to \\ {x}_0=y_1,y_0=x_1-a/b{y}_1 \end{gathered}$$

迭代思路
$$\begin{gathered} a=b{q}_1+r_1 \\ b=r_1{q}_2+r_2 \\ {r}_1=r_2{q}_3+r_3 \\ {r}_2=r_3{q}_4+r_4 \\ ... \\ {r}_{k-1}=r_k{q}_{k+1}+r_{k+1} \end{gathered}$$
上文的求解过程中，$r_i$最终会缩减到$\to gcd(a,b)\to 0$。
令$r_i=a{\lambda }_i+b{\beta }_i$。注意$gcd(a,b)=1$，故左式一定可以成立。
则有
$$r_{i+1}=r_{i-1}-q_{i+1}{r}_i$$
因而有
$$\begin{array}{rl}a{\lambda }_{i+1}+b{\beta }_{i+1}& =a{\lambda }_{i-1}+b{\beta }_{i-1}-q_{i+1}(a{\lambda }_i+b{\beta }_i)\\ & =a({\lambda }_{i-1}-q_{i+1}{\lambda }_i)+b({\beta }_{i-1}-q_{i+1}{\beta }_i)\end{array}$$
故有${\lambda }_{i+1}=({\lambda }_{i-1}-q_{i+1}{\lambda }_i),{\beta }_{i+1}=b({\beta }_{i-1}-q_{i+1}{\beta }_i)$。
注意到上面的式子仅开头两个格式有点不统一，令$a=r_{-1},b=r_0$。则形式统一。其中$q_{i+1}=r_{i-1}/r_i$

代码