【转】斐波那契数列取模（大数）分治算法

 

菲波那契数列如下：1，1，2，3，5，8，13，21，34......此中a[1] = 1， a[2] = 1， a[n]=a[n-1]+a[n-2]（n>=3）。对给定的下标n，求解a[n]％1997的值.

此中测试数据n是整数局限内。

这个题目，主如果用到很关键的一个数学常识，斐波那契数列的求法，可以转换为矩阵的连乘，矩阵的n此方算法又可以用分治的算法。

并且又有理论根据：（n*m）％c=[ （n％c）*（m％c） ]％c ; （n+m）％c=[ （n％c）+（m％c） ]％c ，所以过程中的成果可以随时取模，而不影响终极的成果

关于斐波那契数列的矩阵连乘求法如下：

我们知道我们若要简单策画f（n），有一种办法就是先保存
a=f（0），b=f（1），然后每次设：
a""=b b""=a+b

并用新的a""和b""来持续这一运算

若是大师熟悉哄骗“矩阵”这一对象的话，就知道，若是把a、b写成一个向量[a，b]，完成上述操纵相当于乘以矩阵
0 1
1 1
也就是说，若是我们请求第100个fibonacci数，只须要将矩阵
[0，1]乘上
0 1
1 1
的一百次方，再取出第一项

因为我们知道，矩阵运算满足连络律，一次次右乘那个矩阵完全可以用乘上那个矩阵的N次方庖代，更进一步，那个矩阵的N次方就是如许的情势：
f（n-1） f（n）
f（n） f（n+1）

而求矩阵的N次方，因为矩阵乘法满足连络律，所以我们可以用log（N）的算法求出——这个算法大师都邑么？
一个是二分，一个是基于二进制的求幂

二分的道理：请求矩阵的N次方A（N），设i=N/2若N％2==1， 则 A（N）=A（i）*A（i）*A（1）若N％2==0， 则 A（N）=A（i）*A（i）

基于二进制的道理：将N拆为二进制数，譬如13=1101那么 A^13= A^8 * A^4 * A^1 （这里^默示幂运算）

也就是说，由A^1开端，自乘获得A^2，然后自乘获得A^4，若是N对应位为1，则将这个成果乘到目标上去

如许的话，将所有乘法改为模乘，就可以获得一个较大Fibonacci数除以M的余数

若不消递归，其实类似