HDU 1575 Tr A(矩阵快速幂)
时间:2022-05-08
本文章向大家介绍HDU 1575 Tr A(矩阵快速幂),主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
Problem Description
A为一个方阵,则Tr A表示A的迹(就是主对角线上各项的和),现要求Tr(A^k)%9973。
Input
数据的第一行是一个T,表示有T组数据。 每组数据的第一行有n(2 <= n <= 10)和k(2 <= k < 10^9)两个数据。接下来有n行,每行有n个数据,每个数据的范围是[0,9],表示方阵A的内容。
Output
对应每组数据,输出Tr(A^k)%9973。
Sample Input
2 2 2 1 0 0 1 3 99999999 1 2 3 4 5 6 7 8 9
Sample Output
2 2686
Author
xhd
Source
HDU 2007-1 Programming Contest
Recommend
linle | We have carefully selected several similar problems for you: 1757 1588 2256 2604 2254
矩阵快速幂的裸题!。
我们需要新建一个矩阵,
满足:对角线为1,其余为0
然后跑一下矩阵快速幂就好
1 #include<cstdio>
2 #include<cstring>
3 using namespace std;
4 const int MAXN=101;
5 inline void read(int &n){char c='+';bool flag=0;n=0;
6 while(c<'0'||c>'9') c=='-'?flag=1,c=getchar():c=getchar();
7 while(c>='0'&&c<='9') n=n*10+c-48,c=getchar();flag==1?n=-n:n=n;}
8 struct matrix
9 {
10 int m[11][11];matrix(){memset(m,0,sizeof(m));}
11 };
12 matrix ma;
13 int limit;
14 const int mod=9973;
15 matrix mul(matrix a,matrix b)
16 {
17 matrix c;
18 for(int k=0;k<limit;k++)
19 for(int i=0;i<limit;i++)
20 for(int j=0;j<limit;j++)
21 c.m[i][j]=(c.m[i][j]+(a.m[i][k]*b.m[k][j]))%mod;
22 return c;
23 }
24 matrix fast_martix_pow(matrix ma,int p)
25 {
26 matrix bg;
27 for(int i=0;i<limit;i++)
28 for(int j=0;j<limit;j++)
29 bg.m[i][j]=(i==j);
30 while(p)
31 {
32 if(p&1) bg=mul(bg,ma);
33 ma=mul(ma,ma);
34 p>>=1;
35 }
36 return bg;
37 }
38 int main()
39 {
40 int T;read(T);
41 while(T--)
42 {
43 read(limit);int n;read(n);
44 for(int i=0;i<limit;i++)
45 for(int j=0;j<limit;j++)
46 read(ma.m[i][j]);
47 matrix ans=fast_martix_pow(ma,n);
48 int out=0;
49 for(int i=0;i<limit;i++)
50 out+=ans.m[i][i]%mod;
51 printf("%dn",out%mod);
52 }
53 return 0;
54 }
- Flash/Flex学习笔记(48):反向运动学(下)
- 如何给sublime text3安装汉化包?so easy 哦
- Flash/Flex学习笔记(47):反向运动学(上)
- 汉诺塔问题算法介绍
- Flash/Flex学习笔记(36):自己动手实现一个滑块控件(JimmySilder)
- 使用Zabbix服务端本地邮箱账号发送报警邮件及指定报警邮件操作记录
- fckeditor上传问题的解决
- Flash/Flex学习笔记(46):正向运动学
- 异步Socket处理的一些测试值
- .Net中DES加密的细节问题
- 分布式监控系统Zabbix--完整安装记录 -添加web页面监控
- .Net中使用带UI的OCX的方法
- 2017数据科学领域15大热门GitHub项目
- 简单分页的性能优化
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- jQuery限制复选框checkbox的选中次数
- jQuery点击切换增加和删除class类
- Vue使用props和emit父子组件通信
- 听说Mysql你很豪横?-------------呕心沥血深入解析mysql备份与恢复!!!
- 排障集锦:九九八十一难之第十难!mysq备份恢复,Could not read entry at offset *: Error in log format or read error.
- Vue使用ref父子组件通信
- 听说Mysql你很豪横?-------------MySQL5.7主从同步
- 听说Mysql你很豪横?-------------MySQL5.7主从复制!读写分离!
- Vue兄弟组件传值
- Vue设置浏览器的标题title和图标icon
- VantUI封装自定义Tabbar路由跳转
- 听说Mysql你很豪横?-------------分分钟带你玩转SQL高级查询语句(常用查询,正则表达式,运算符)
- jQuery实现点击添加样式同胞移除样式
- 微信小程序生命周期
- 听说Mysql你很豪横?-------------分分钟带你玩转SQL高级查询语句(库函数,存储过程)