HDU - problem Tr A 【快速幂 + 矩阵】
来源:互联网 发布:oppo手机root软件 编辑:程序博客网 时间:2024/06/05 07:43
Tr A
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 4488 Accepted Submission(s): 3377
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的内容。
每组数据的第一行有n(2 <= n <= 10)和k(2 <= k < 10^9)两个数据。接下来有n行,每行有n个数据,每个数据的范围是[0,9],表示方阵A的内容。
Output
对应每组数据,输出Tr(A^k)%9973。
Sample Input
22 21 00 13 999999991 2 34 5 67 8 9
Sample Output
22686
Author
xhd
Source
HDU 2007-1 Programming Contest
Recommend
linle | We have carefully selected several similar problems for you: 1588 2256 2604 2254 3117
#include <map>#include <set>#include <cstdio>#include <cstring>#include <algorithm>#include <queue>#include <iostream>#include <stack>#include <cmath>#include <string>#include <vector>#include <cstdlib>//#include <bits/stdc++.h>//#define LOACL#define space " "using namespace std;typedef long long LL;typedef __int64 Int;typedef pair<int, int> paii;const int INF = 0x3f3f3f3f;const double ESP = 1e-5;const double PI = acos(-1.0);const int MOD = 9973;const int MAXN = 15;int n, m;struct Matrix {LL m[MAXN][MAXN];int row, col;};Matrix ori, res, u;void init() {memset(res.m, 0, sizeof(res.m));ori.row = ori.col = n;}void scan_in() {for (int i = 1; i <= n; i++) {for (int j = 1; j <= n; j++) {scanf("%d", &ori.m[i][j]);}}}Matrix multi(Matrix x, Matrix y) {Matrix z;memset(z.m, 0, sizeof(z.m));z.row = x.row; z.col = y.col;for (int i = 1; i <= x.row; i++) {for (int k = 1; k <= x.col; k++) {for (int j = 1; j <= y.col; j++) {z.m[i][j] += x.m[i][k]*y.m[k][j]%MOD;}z.m[i][k] %= MOD;}}return z;}Matrix pow_mod(Matrix a, int x){Matrix b;b.col = a.col; b.row = a.row;memset(b.m, 0, sizeof(b.m));for (int i = 1; i <= n; i++) {b.m[i][i] = 1;} while(x){ if(x&1) b = multi(a,b); a = multi(a, a); x >>= 1; } return b;}int main() {int T;scanf("%d", &T);while (T--) {scanf("%d%d", &n, &m);init(); scan_in();res = pow_mod(ori, m);int ans = 0;for (int i = 1; i <= n; i++) {ans += res.m[i][i];}printf("%d\n", ans%MOD);}return 0;}
0 0
- HDU - problem Tr A 【快速幂 + 矩阵】
- hdu 1575Tr A 矩阵快速幂
- HDU 1575 Tr A(矩阵快速幂)
- HDU 1575 Tr A (矩阵快速幂)
- HDU 1575-Tr A(矩阵快速幂)
- HDU 1575 Tr A (矩阵快速幂)
- HDU - 1575 Tr A(矩阵快速幂)
- hdu 1575 Tr A(矩阵快速幂)
- Hdu 1575 Tr A【矩阵快速幂】
- HDU 1575 Tr A (矩阵快速幂)
- 【HDU 1575】Tr A (矩阵快速幂)
- HDU 1575 Tr A 矩阵快速幂
- 【HDU 1575 Tr A】+ 矩阵快速幂
- HDU-1575-Tr A【矩阵快速幂】
- hdu 1575 TR A 矩阵快速幂
- HDU 1575 - Tr A(矩阵快速幂)
- HDU 1575 Tr A 矩阵快速幂
- HDU.1575 Tr A ( 矩阵快速幂)
- 关于C语言的数据类型需要了解的基本概念
- VIM常用用法
- 手机控制的浇花神器
- Python神图
- C#之2.0
- HDU - problem Tr A 【快速幂 + 矩阵】
- 小白手动编译安装mysql过程
- 2016.11.1
- UVA 11136 Hoax or what(multiset)
- 管道(传统进程间通信方式)
- Mac 10.11安装MySQLdb-python时报错总结
- 关于结构体指针以及.和->区别的几个尝试
- 查询activity是否存在栈中
- android四大组件之service知识点最强总结