浅谈斐波那契数列——从递推到矩阵乘法

说在前面

相信大家都已经知道这个中外著名的费波纳切数列了吧，关于费波那契数列有很多有趣的性质，但我们这里不讲，在这里我们只是利用斐波那契数列来引出另一个神奇的东西，矩阵乘法，递推在这里是起一个对比与铺垫的作用，（没有对比就不知道矩阵乘法有多快）。

斐波那契数列

定义一个数列

F(n)=⎧⎩⎨0,1，Fn−1+Fn−2,n=0 n=1n≥2

递推

递推（我当大家都会了），但是这样的做法时间复杂度是 O(n) ，的，当n很大时，这个算法就会T。

通项公式

斐波那契数列有一个通项公式，

F(n)=15√×⎡⎣(1+5√2)n−(1−5√2))n⎤⎦

,直接用快速幂计算两个log，但是如果要取模怎么办？求一个非整数的逆元？通项公式似乎很难。

进入正题——矩阵乘法

知识预备：
①设

A=[1−10321]

B=⎡⎣⎢321110⎤⎦⎥

AB=[(1×3+0×2+2×1)(−1×3+3×2+1×1)(1×1+0×1+2×0)(−1×1+3×1+1×0)]=[5412]

②设

A=[abcd]

B=[ef]

AB=[ae+cfbe+df]

有了这些东西，接下来就可以做斐波那契了，把相邻两项写在一个2*1的矩阵中，

[FnFn−1]=[Fn−1+Fn−2Fn−1]=[Fn−1×1+Fn−2×1Fn−1×1+Fn−2×0]=[1110]×[Fn−1Fn−2]=[1110]n−1×[F1F0]=[1110]n−1×[10]

那么问题的本质，就是如何解决计算这一块

[1110]n−1

算完之后，取矩阵第一行第一列的元素，计算这个其实又是快速幂，说些什么呢？谢谢蒙格马利吧。

代码

快速幂：

function f(a,b:extended;n:int64):int64;
var
        t,y:extended;
begin
        t:=1;
        y:=a;
        while b<>0 do
        begin
                if(b and 1)=1 then
                t:=t*y mod n  ;
                y:=y*y mod n;
                b:=b shr 1;
        end;
        exit(t);
end;

用快速幂求出斐波那契的第n项。

type
        matrix=array[0..1,0..1] of int64;
const
        mo=trunc(1e+8+7);
var
    a,ans:array[0..1] of int64;
    b,c,d:matrix;
    i,j,k:longint;
        n:int64;
procedure work(var a,b,c:matrix);//两个2*2的矩阵相乘
var
        i,j,k:longint;
begin
        for i:=0 to 1 do
        begin
                for j:=0 to 1 do
                begin
                        for k:=0 to 1 do
                        begin
                                c[i,j]:=(c[i,j]+a[i,k]*b[k,j] mod mo) mod mo;
                        end;
                end;
        end;
end;
procedure fastpower(x:int64);
begin
        c[0,0]:=1;
    c[0,1]:=1;
    c[1,0]:=1;
        c[1,1]:=0;
    b[0,0]:=1;//开始将B赋值为单位矩阵，相当于实数中1在乘法中的作用
    b[1,1]:=1;
    while x<>0 do
    begin
        if (x and 1)=1 then
        begin
            fillchar(d,sizeof(d),0);
                        work(b,c,d);
                        b:=d;
        end;
        fillchar(d,sizeof(d),0);
                work(c,c,d);
                c:=d;
                x:=x div 2;
    end;
end;
begin
        readln(n);
        if n=0 then
                writeln(0);
        if (n=1)or(n=2) then
                writeln(1);
    fastpower(n-1);
    a[0]:=1;
    a[1]:=0;
    for i:=0 to 1 do
    begin
        for j:=0 to 1 do
        begin
            ans[i]:=(ans[i]+a[j]*b[j][i] mod mo) mod mo;
        end;
    end;
    writeln(ans[0] mod mo);
end.

可以发现其实两者实质上是一样的，只是在计算的时候，一个用矩阵计算，另一个直接用一个变量计算，（矩阵其实也就是一堆变量）。

the end

由于我的水平有限，难免有些会写错的，希望大家多多包容，批评指正，谢谢。