浅谈用矩阵求斐波那契数列

导语

两个矩阵的乘法仅当第一个矩阵$A$
的列数和第二个矩阵$B$
的行数相等时才能定义。如
$A$
是$m\ast n$矩阵，$B$是$n\ast p$的矩阵，他们的乘积
$C$=$(c_{i,j})$,它的任意一个元素值为：

$c_{i,j}=a_{i,1}{b}_{1,j}+a_{i,2}{b}_{2,j}+...+a_{i,n}{b}_{n,j}=\sum _{n=1}^n{a}_{i,r}{b}_{r,j}$

求斐波那契数列数列有多种方法，比如用递推公式，或是用生成函数求出其通项公式，可惜这两种方法在遇到数据规模较大时会超时，下面介绍一种方法我们将用到矩阵，用矩阵快速幂的方法可以求得斐波那契数列的第i项。

证明

我们首先证明这样一个结论

$gcd(F_n,F_m)=F_{gcd(n,m)}$

证明1

设$F_n=a,F_{n+1}=b$

所以 $F_{n+2}=a+b,F_{n+3}=a+2b,F_{n+4}=2a+3b,F_{(}n+5)=3a+5b$

得到

$F_m=F_{m-n-1}\ast a+F_{m-n}\ast b=F_{m-n-1}\ast F_n+F_{m-n}\ast F_{n+1}$

假设$m>n$,需求$gcd(F_n,F_m)$

又因为$F_n|F_n\ast F_{m-n-1}$

所以

$gcd(F_n,F_m)⟺gcd(F_n,F_{m-n}{F}_{n+1})$

这里再证明一个结论

$gcd(F_n,F_{n+1})=1$

证明2

$gcd(F_n,F_{n+1})=gcd(F_n,F_{n+1}-F_n)=gcd(F_n,F_{n-1})=gcd(F_{n-2},F_{n-1})=...=gcd(F_1,F_2=1)$

故

$gcd(F_n,F_{n+1})=1$

回到证明1

因为$gcd(F_n,F_{n+1})=1$

所以

$gcd(F_n,F_{m-n}{F}_{n+1})=gcd(F_n,F_{m-n})$

$gcd(F_n,F_m)=gcd(F_n,F_{m-n})$

$gcd(F_n,F_{m-n})=gcd(F_n,F_{m-2n})$

…

$gcd(F_n,F_m)=gcd(F_n,F_{m-(m/n)\ast n})$

然后$m1=n,n1=m-(m/n)n$

再来

直到出现$F_0$,此时$F_{n1}$即我们要求的答案，$n1==gcd(n,m)$

证毕

证明2

对于斐波那契数列，我们可以这样理解

$\left[\begin{array}{cc}{F}_{n+1}& F_n\\ F_n& F_{n-1}\end{array}\right]$

即

$\left[\begin{array}{cc}{F}_n+F_{n-1}& F_n+F_{n-1}\cdot 0\\ F_n+F_{n-1}\cdot 0& F_{n-1}+F_n\cdot 0\end{array}\right]$

$\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]\times \left[\begin{array}{cc}{F}_n& F_{n-1}\\ F_{n-1}& F_{n-2}\end{array}\right]$

如此持续操作，直到出现

$\left[\begin{array}{cc}{F}_2& F_1\\ F_1& F_0\end{array}\right]$

即

$\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]$

此间有$n$个$\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]$相乘，即$\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]{}^n$

因此用矩阵乘法求快速幂即可

证毕