BZOJ_5015_[Snoi2017]礼物_矩阵乘法

#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
typedef long long ll;
int m;
ll n,mod=1000000007,c[15][15];
struct Mat {
	ll v[14][14];
	Mat(){memset(v,0,sizeof(v));}
	Mat operator*(const Mat &x)const {
		Mat re;int i,j,k;
		for(i=1;i<=m;i++) {
			for(j=1;j<=m;j++) {
				for(k=1;k<=m;k++) {
					(re.v[i][j]+=v[i][k]*x.v[k][j]%mod)%=mod;
				}
			}
		}
		return re;
	}
};
void print(Mat x) {
	int i,j;
	for(i=1;i<=m;i++) {
		for(j=1;j<=m;j++) {
			printf("%lld ",x.v[i][j]);
		}
		puts("");
	}
}
Mat qp(Mat x,ll y) {
	Mat I;
	int i;
	for(i=1;i<=m;i++) I.v[i][i]=1;
	while(y) {
		if(y&1ll) I=I*x;
		x=x*x;
		y>>=1ll;
	}
	return I;
}
Mat build() {
	Mat x;
	int i,j;
	x.v[1][1]=2; x.v[2][1]=1;
	for(i=2;i<=m;i++) {
		int k=m-i+1;
		for(j=0;j<=k;j++) {
			x.v[j+i][i]=c[k-1][j];
		}
	}
	/*for(i=1;i<=m;i++) {
		for(j=1;j<=m;j++) {
			printf("%lld ",x.v[i][j]);
		}
		puts("");
	}*/
	return x;
}
void init() {
	int i,j;
	for(i=0;i<=m;i++) c[i][0]=c[i][i]=1;
	for(i=1;i<=m;i++) {
		for(j=1;j<=m;j++) {
			c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
		}
	}
}
int main() {
	scanf("%lld%d",&n,&m); m+=2;
	init();
	Mat A;
	int i;
	for(i=2;i<=m;i++) A.v[1][i]=1;
	Mat x=build();
	Mat T=A*qp(x,n-1);
	//print(T);
	printf("%lld\n",((T*x).v[1][1]-T.v[1][1]+mod)%mod);
}