大致意思是定义 f0=f1=1f_0=f_1=1f0=f1=1,fi=fi−1+fi−2f_i=f_{i-1}+f_{i-2}fi=fi−1+fi−2。
求Σi=0nfifn−i\Sigma^n_{i=0}f_if_{n-i}Σi=0nfifn−i。
令Fn=Σi=0nfifn−iF_n=\Sigma^n_{i=0}f_if_{n-i}Fn=Σi=0nfifn−i。
Fn−Fn−1=Σi=0nfifn−i−Σi=0n−1fifn−1−iF_n-F_{n-1}=\Sigma^n_{i=0}f_if_{n-i}-\Sigma^{n-1}_{i=0}f_if_{n-1-i}Fn−Fn−1=Σi=0nfifn−i−Σi=0n−1fifn−1−i
=f0fn+Σi=1nfifn−i−Σi=0n−1fifn−1−i=f_0f_n+\Sigma^{n}_{i=1}f_if_{n-i}-\Sigma^{n-1}_{i=0}f_if_{n-1-i}=f0fn+Σi=1nfifn−i−Σi=0n−1fifn−1−i
=f0fn+Σi=0n−1fifn−i−Σi=0n−1fifn−1−i=f_0f_n+\Sigma^{n-1}_{i=0}f_if_{n-i}-\Sigma^{n-1}_{i=0}f_if_{n-1-i}=f0fn+Σi=0n−1fifn−i−Σi=0n−1fifn−1−i
=fn+Σi=0n−1fi(fn−i−fn−1−i)=f_n+\Sigma^{n-1}_{i=0}f_i(f_{n-i}-f_{n-1-i})=fn+Σi=0n−1fi(fn−i−fn−1−i)
=fn+Σi=0n−1fi(fn−i−2)=f_n+\Sigma^{n-1}_{i=0}f_i(f_{n-i-2})=fn+Σi=0n−1fi(fn−i−2)
=fn+Fn−2=f_n+F_{n-2}=fn+Fn−2
所以,Fn=Fn−1+Fn−2+fnF_n=F_{n-1}+F_{n-2}+f_nFn=Fn−1+Fn−2+fn。
于是,Fn−1=Fn−2+Fn−3+fnF_{n-1}=F_{n-2}+F_{n-3}+f_nFn−1=Fn−2+Fn−3+fn,Fn−2=Fn−3+Fn−4+fnF_{n-2}=F_{n-3}+F_{n-4}+f_nFn−2=Fn−3+Fn−4+fn,二式三式相加并与一试作差,可消去 fnf_nfn。得到 FFF 的递推式:Fn=2Fn−1+Fn−2−2Fn−3−Fn−4F_n=2F_{n-1}+F_{n-2}-2F_{n-3}-F_{n-4}Fn=2Fn−1+Fn−2−2Fn−3−Fn−4,用矩阵加速求即可。
转移矩阵为:
{{2,1,0,0},
{1,0,1,0},
{-2,0,0,1},
{-1,0,0,0}}
初始化:F0=1,F1=2,F2=5,F3=10F_0=1, F_1=2, F_2=5, F_3=10F0=1,F1=2,F2=5,F3=10。