UVA 11551 - Experienced Endeavour(矩阵高速幂)
2016-03-30 08:51
447 查看
UVA 11551 - Experienced Endeavour
题目链接题意:给定一列数,每一个数相应一个变换。变换为原先数列一些位置相加起来的和,问r次变换后的序列是多少
思路:矩阵高速幂,要加的位置值为1。其余位置为0构造出矩阵,进行高速幂就可以
代码:
#include <cstdio> #include <cstring> const int N = 55; int t, n, r, a ; struct mat { int v ; mat() {memset(v, 0, sizeof(v));} mat operator * (mat c) { mat ans; for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { for (int k = 0; k < n; k++) ans.v[i][j] = (ans.v[i][j] + v[i][k] * c.v[k][j]) % 1000; } } return ans; } }; mat pow_mod(mat x, int k) { mat ans; for (int i = 0; i < n; i++) ans.v[i][i] = 1; while (k) { if (k&1) ans = ans * x; x = x * x; k >>= 1; } return ans; } int main() { scanf("%d", &t); while (t--) { scanf("%d%d", &n, &r); for (int i = 0; i < n; i++) scanf("%d", &a[i]); int x; mat Mat; for (int i = 0; i < n; i++) { scanf("%d", &x); int b; while (x--) { scanf("%d", &b); Mat.v[i][b] = 1; } } Mat = pow_mod(Mat, r); for (int i = 0; i < n; i++) { int ans = 0; for (int j = 0; j < n; j++) { ans = (ans + a[j] * Mat.v[i][j]) % 1000; } printf("%d%c", ans, (i == n - 1 ? '\n' : ' ')); } } return 0; }
相关文章推荐
- Socket连接与HTTP连接
- ie 7/8不支持trim方法
- Android源码编译环境搭建问题及解决方案小结
- 标准Web系统的架构分层
- 设计的目标
- android基础学习之通过Edittext为spinner添加内容
- 欢迎使用CSDN-markdown编辑器
- BaseServlet 继承 httpServlet
- 2016年3月26日--4月2日(共16小时,剩3036小时)
- algrothm_面向对象_封装
- mongodb_Mongo3.2.4框架及角色介绍
- 目标检测程序开发(一)——OpenCV2.4.11环境配置
- 分布式系统架构的基本原则和实践(转)
- 安全编程
- 基于jQuery实现Ajax验证用户名是否存在实例
- 深度优先搜索DFS(迷宫问题)
- 【BZOJ-1588】营业额统计 Splay
- CSS Hack解析
- 排序算法-堆排序
- algrothm_面向对象_递归