数论：欧几里得算法与斐波那契数列

欧几里得算法的一些秘密

1、欧几里得算法的最坏运行时间——斐波那契数列

欧几里得算法的一般运行次数是取决于$a,b$的位数。

欧几里得算法：$gcd(a,b)=gcd(b,amodb)$

显然，最终答案是$gcd(r,0)$ ，此时递归返回$gdc=r$

欧几里得算法的运行时间，就取决与递归的次数$k$​。

递归过程：

$gcd(a,b)\to gcd(b,amodb)\to ...\to gcd(r,0)$

逆过来：

$gcd(r,0)\to gcd(k_1(r),r)\to gcd(k_2(k_{1}r)+r,k_{1}r)\to gcd(k_3(k_2{k}_{1}r+r)+k_{1}r,k_2{k}_{1}r+r)\to ...\to gcd(a,b)$

$gcd(k4(k_3{k}_2{k}_{1}r+(k_1+k_3)r)+(k_2{k}_{1}r+r),k_3{k}_2{k}_{1}r+(k_1+k_3)r)$

我们可以观察，每个$gcd(a,b)$中$a$的变化：

$r\to k_{1}r\to k_2{k}_{1}r+r\to k_3(k_2{k}_{1}r+r)+k_{1}r\to k4(k_3(k_2{k}_{1}r+r)+k_{1}r)+(k_2{k}_{1}r+r)\to ...$

$f(0)\to f(1)\to f(2)\to ...\to f(n)$

可以发现满足一个这样的性质：

$f(i)modf(i-1)=f(i-2)$

所以，$a$的变化可以写成：$f(i)=kf(i-1)+f(i-2)\ge f(i-1)+f(i-2)=Fibonacci(i)=F(i)$

所以简化递归过程：

$f_a(0)=r\to f_a(1)=k_1{f}_a(0)\to f_a(2)=k_2{f}_a(1)+f_a(0)\to ...$

$f_b(0)=0\to f_b(1)=f_a(0)\to f_b(2)=f_a(1)\to ...$

所以，我们令$k,r=1$​

可以得到一个数组$F_a:1,1,2,3,...,$（$a$的变化过程）​

就可以得到$a,b$递归$k$次的最小取值：$a=F_a(k+2),b=F_a(k+1)$。

$gcd(F(k+2),F(k+1))=gcd(F(k+1),F(k+2)modF(k+1))=gcd(F(k+1),F(k))$

注意：上述的递归$k$次的最小取值，仅在$a>b\ge 0$ 的情况下。

所以，我们就可以通过，$a,b$的取值，大概判断，调用的次数。