【又见斐波那契——矩阵快速幂】

又见斐波那契

---

题目描述

这是一个加强版的斐波那契数列。
给定递推式

求F(n)的值，由于这个值可能太大，请对109+7取模。

---

输入描述:

第一行是一个整数T(1 ≤ T ≤ 1000)，表示样例的个数。
以后每个样例一行，是一个整数n(1 ≤ n ≤ 10^18)。

输出描述:

每个样例输出一行，一个整数，表示F(n) mod 1000000007。

示例

输入
4
1
2
3
100

输出
1
16
57
558616258

---

思路

构造矩阵
ans[6]=[1,0,8,4,2,1];
由题：f(i)=f(i-1)+f(i-2)+i^3 +i^2+i+1;
提出每一项的系数就得到了A[0]=[1,1,1,1,1,1]
考虑 f(i-1) = 0f(i)+1f(i-1)+0i^3 +0 * i^2+0 * i+0 * 1;
得 A [1]=[0,1,0,0,0,0]
(i+1)^3展开为 0 * f(i-1)+0 * f(i-2) +i^3+3i ^ 2+3*i+1,
得A[2]=[0,0,1,3,3,1]
以下就不写了
得到A=[ f( i ) , f(i-1) , (i+1)^3 , (i+1)^2 , i+1 , 1 ];
推出A后直接快速幂就行，然后就是记得取模

AC代码：

#include<iostream>
#include<cstring>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
struct node{
	ll num[6][6]={
       {1,1,1,1,1,1},
       {1,0,0,0,0,0},
       {0,0,1,3,3,1},
       {0,0,0,1,2,1},
       {0,0,0,0,1,1},
       {0,0,0,0,0,1}
    }; 
};
void init(node &temp){
	temp.num[0][0]=1;
	temp.num[1][0]=0;
	temp.num[2][0]=8;
	temp.num[3][0]=4;
	temp.num[4][0]=2;
	temp.num[5][0]=1;
}
node f(node a,node b){
	node temp;
	for(int i=0;i<6;++i)
		for(int j=0;j<6;++j){
			temp.num[i][j]=0;
			for(int k=0;k<6;++k)
				temp.num[i][j]+=(a.num[i][k]*b.num[k][j])%mod;
            temp.num[i][j]%=mod;
		}
	return temp;
}
int t;
ll n;
int main(){
	cin>>t;
	while(t--){
		cin>>n;
		if(n<=1){
			printf("%lld\n",n);
			continue;
		}
		node A,ans;
		init(ans);
		--n;
		while(n){
			if(n&1)
				ans=f(A,ans);
			A=f(A,A);
			n>>=1;
		}
		f(A,ans);
		printf("%lld\n",ans.num[0][0]%mod);
	}
} 

---