矩阵快速幂-又见斐波那契 （2021-8-17）

又见斐波那契

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

求F(n)的值，由于这个值可能太大，请对109+7取模。
输入
第一行是一个整数T(1 ≤ T ≤ 1000)，表示样例的个数。
以后每个样例一行，是一个整数n(1 ≤ n ≤ 1018)。
输出
每个样例输出一行，一个整数，表示F(n) mod 1000000007。
示例1
输入
4
1
2
3
100
输出
1
16
57
558616258
关于构造方阵的小技巧：
（1）递推式有几项，矩阵有几行，但是对于中间缺项的，要补上。
（2）配凑系数决定方阵
eg: F(n)=aF(n-1)+bF(n-3]+c
F(n) | a 0 b 1 | F(n-1)
F(n-1) |1 0 0 0 | F(n-2)
F(n-2) | 0 1 0 0 | F(n-3)
c | 0 0 0 1 | c
补上了F（n-2）

此题思路：

1. F（i）= F（i-1）+ F（i-2）+ i^3 + i^2 + i + 1
2. 根据递推式，确定目标，构造方阵
3. 因为F（1）=1，所以方阵为n-1次
   
   4.初始序列，i=2，代入，得{1，0，8，4，2，1}
   5.方阵的n-1的第一行*对应初始序列=F（n）
   6.矩阵快速幂模板（结构体版）

#include<iostream>
#include<string.h>
using namespace std;
#define int long long
const int mod=1e9+7;
const int maxn=105;
struct Mat{//1.模板记一下  2.找方阵竟然开始找错了
	int m[maxn][maxn];
};
Mat a,ans;
const int N=6;
Mat Mul(Mat a,Mat b){
	Mat c;
	for(int i=1;i<=N;i++){
		for(int j=1;j<=N;j++){
			c.m[i][j]=0;
		}
	}
	for(int i=1;i<=N;i++){
		for(int k=1;k<=N;k++){
			for(int j=1;j<=N;j++){
				c.m[i][j]=(c.m[i][j]+a.m[i][k]*b.m[k][j]%mod)%mod;
			}
		}
	}
	return c;
}
Mat power(Mat a,int n){
	memset(ans.m,0,sizeof(ans.m));
	for(int i=1;i<=N;i++){
		ans.m[i][i]=1;
	}
	while(n){
		if(n%2){
			ans=Mul(ans,a);
		}
		a=Mul(a,a);
		n/=2;
	}
	return ans;
}
void init(){
	a.m[1][1]=1;a.m[2][1]=1;a.m[3][1]=0;a.m[4][1]=0;a.m[5][1]=0;a.m[6][1]=0;
	a.m[1][2]=1;a.m[2][2]=0;a.m[3][2]=0;a.m[4][2]=0;a.m[5][2]=0;a.m[6][2]=0;
	a.m[1][3]=1;a.m[2][3]=0;a.m[3][3]=1;a.m[4][3]=0;a.m[5][3]=0;a.m[6][3]=0;
	a.m[1][4]=1;a.m[2][4]=0;a.m[3][4]=3;a.m[4][4]=1;a.m[5][4]=0;a.m[6][4]=0;
	a.m[1][5]=1;a.m[2][5]=0;a.m[3][5]=3;a.m[4][5]=2;a.m[5][5]=1;a.m[6][5]=0;
	a.m[1][6]=1;a.m[2][6]=0;a.m[3][6]=1;a.m[4][6]=1;a.m[5][6]=1;a.m[6][6]=1;
}
signed main(){
	int t;
	cin>>t;
	int n;
	init();
	int ceil[7]={0,1,0,8,4,2,1};
	int ret;
	while(t--){
		cin>>n;
		if(n==0){
			cout<<"0"<<endl;
		}
		else if(n==1){
			cout<<"1"<<endl;
		}
		else{
			ret=0;
			power(a,n-1);
//			for(int i=1;i<=N;i++){
//				for(int j=1;j<=N;j++){
//					cout<<ans.m[i][j]<<" "; 
//				} 
//				cout<<endl;
//			}
			for(int i=1;i<=N;i++){
				ret=(ret+(ans.m[1][i]*ceil[i]%mod))%mod;
			}	
			cout<<ret<<endl;
		}
		
		
	}
	
	return 0;
}