Matrix[矩阵乘法]

传送门


原博客大佬太巨了 

 

 

 

 

 

 

 


#include<bits/stdc++.h>
#define N 15
#define LL long long 
using namespace std;
struct Matrix{
	LL a[N][N];
	Matrix(){memset(a,0,sizeof(a));}
};
int n,m; LL Mod=10000007;
Matrix mul(Matrix A,Matrix B){
	Matrix C;
	for(int i=1;i<=n+2;i++)
		for(int j=1;j<=n+2;j++)
			for(int k=1;k<=n+2;k++)
				C.a[i][j] = (C.a[i][j] + A.a[i][k]*B.a[k][j]) % Mod;
	return C;
}
Matrix power(Matrix X,int M){
	Matrix ans;
	for(int i=1;i<=n+2;i++)
		ans.a[i][i]=1;
	for(;M;M>>=1){
		if(M&1) ans = mul(X,ans);
		X = mul(X,X); 
	} return ans;
}
int main(){
	while(~scanf("%d%d",&n,&m)){
		Matrix A,B;  
		for(int i=1;i<=n;i++){
			scanf("%d",&A.a[i+1][1]);
		}A.a[1][1]=23 , A.a[n+2][1]=3;
		for(int i=1;i<=n+1;i++)
			B.a[i][1]=10;
		for(int i=1;i<=n+2;i++)
			B.a[i][n+2]=1;
		for(int i=1;i<=n+1;i++)
			for(int j=2;j<=i;j++)
				B.a[i][j]=1;
		B = power(B,m); A = mul(B,A);
		
		printf("%lld\n",A.a[n+1][1] % Mod);
	}return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

FSYo

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值