UVA 11149 - Power of Matrix(矩阵乘法)
2015-09-18 12:50
405 查看
UVA 11149 - Power of Matrix
题目链接题意:给定一个n*n的矩阵A和k,求∑kiAi
思路:利用倍增去搞。∑kiAi=(1+Ak/2)∑k/2iAi,不断二分就可以
代码:
#include <cstdio> #include <cstring> const int N = 45; int n, k; 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]) % 10; } } } return ans; } mat operator + (mat c) { mat ans; for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) ans.v[i][j] = (v[i][j] + c.v[i][j]) % 10; return ans; } } A; 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; } mat solve(mat x, int k) { if (k == 1) return x; mat ans; for (int i = 0; i < n; i++) ans.v[i][i] = 1; if (k == 0) return ans; ans = (ans + pow_mod(x, k>>1))* solve(x, k>>1); if (k&1) ans = ans + pow_mod(x, k); return ans; } int main() { while (~scanf("%d%d", &n, &k) && n) { for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) { scanf("%d", &A.v[i][j]); A.v[i][j] %= 10; } A = solve(A, k); for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) printf("%d%c", A.v[i][j], (j == n - 1 ? '\n' : ' ')); printf("\n"); } return 0; }
相关文章推荐
- 大嘴鱼游戏
- 路径问题
- 图的那些事儿——Dijkstra和Floyd
- [LeetCode]题解(python):013-Roman to Integer
- SQL SERVER 2012数据库:开启防火墙导致外部无法连接数据库解决办法
- 转-MAC 下安装PHONEGAP开发环境
- HDOJ4832Chess【dp+组合数学】
- Cocos2d-x项目脚本的创建
- 使用位操作
- 快速获取Windows系统上的国家和地区信息
- POJ 1847 Tram
- Dapper.net ORM
- LR录制手机APP脚本的一种方法
- 中缀表达式求值
- 第3周项目2 建设“顺序表”算法库
- iOS学习重要知识点整理02-进程和线程的一个简单解释
- 创建使用动态链接库
- 计算机网络之概念解析
- 杭电acm1005
- Codeforces Round #320 (Div. 2)(A~E)