又见斐波那契
题目
这是一个加强版的斐波那契数列。
给定递推式
求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)
此题思路:
- F(i)= F(i-1)+ F(i-2)+ i^3 + i^2 + i + 1
- 根据递推式,确定目标,构造方阵
- 因为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;
}