hdu5015 233 Matrix

一列一列地推就可以了

#include <iostream>
#include <cstdio>
using namespace std;
typedef long long ll;
int n, m;
const int mod=10000007;
struct Matrix{
	int num[15][15];
	Matrix operator*(const Matrix &x)const{
		Matrix re;
		for(int i=0; i<=n+1; i++)
			for(int j=0; j<=n+1; j++){
				re.num[i][j] = 0;
				for(int k=0; k<=n+1; k++)
					re.num[i][j] = (re.num[i][j] + (ll)num[i][k]*x.num[k][j]) % mod;
			}
		return re;
	}
}dan, zhu, yua;
Matrix ksm(int b){
	Matrix re=dan;
	while(b){
		if(b&1)	re = re * zhu;
		zhu = zhu * zhu;
		b >>= 1;
	}
	return re;
}
int main(){
	for(int i=0; i<=11; i++)
		dan.num[i][i] = 1;
	while(scanf("%d %d", &n, &m)!=EOF){
		yua.num[0][0] = 23;
		for(int i=1; i<=n; i++)
			scanf("%d", &yua.num[0][i]);
		yua.num[0][n+1] = 1;
		for(int i=0; i<=n+1; i++)
			for(int j=0; j<=n+1; j++){
				if(j==n+1){
					if(i==n+1)	zhu.num[i][j] = 1;
					else	zhu.num[i][j] = 0;
				}
				else if(i==0)	zhu.num[i][j] = 10;
				else if(i==n+1)	zhu.num[i][j] = 3;
				else if(i<=j)	zhu.num[i][j] = 1;
				else	zhu.num[i][j] = 0;
			}
		cout<<(yua*ksm(m)).num[0][n]<<endl;
	}
	return 0;
}
posted @   poorpool  阅读(172)  评论(0编辑  收藏  举报
编辑推荐:
· 深入理解 Mybatis 分库分表执行原理
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
· 现代计算机视觉入门之:什么是图片特征编码
· .NET 9 new features-C#13新的锁类型和语义
阅读排行:
· 手把手教你在本地部署DeepSeek R1,搭建web-ui ,建议收藏!
· Spring AI + Ollama 实现 deepseek-r1 的API服务和调用
· 《HelloGitHub》第 106 期
· 数据库服务器 SQL Server 版本升级公告
· C#/.NET/.NET Core技术前沿周刊 | 第 23 期(2025年1.20-1.26)
点击右上角即可分享
微信分享提示