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;
}
拙いものと思えども、
その手に握る其れこそが、
いつか幻想を生んでいく。
分类:
数学——矩阵相关
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 深入理解 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)