Tr A
来源:互联网 发布:webvr源码 编辑:程序博客网 时间:2024/05/29 17:02
Time Limit : 1000/1000ms (Java/Other) Memory Limit : 32768/32768K (Java/Other)
Total Submission(s) : 48 Accepted Submission(s) : 28
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
#include<stdio.h>#include<algorithm>using namespace std;#define mod 9973int n, k;typedef struct node{int mat[15][15];}mat;mat A; // 定义 A 与 mat相同的结构体void In() // 输入矩阵 初始化{scanf("%d%d", &n, &k); // 注意要从 A.mat[0][0] 开始赋值for (int i = 0; i < n; i++)for (int j = 0; j < n; j++){scanf("%d", &A.mat[i][j]);}}mat mul(mat a,mat b) // 矩阵乘法{mat C;memset(C.mat, 0, sizeof(C.mat)); // 特别注意 等价于 C.mat[i][j] = 0;for (int i = 0; i<n; i++) { for (int j = 0; j<n; j++) { for (int k = 0; k < n; k++) {C.mat[i][j] =( C.mat[i][j]+a.mat[i][k] * b.mat[k][j] % mod)%mod; } } }return C;}mat Pow(mat A, int k) // 求 A^k {mat B;memset(B.mat, 0, sizeof(B.mat)); // 特别注意 等价于 B.mat[i][j] = 0;for (int i = 0; i<n; i++) // 单位矩阵初始化{B.mat[i][i] = 1;}while (k>0){if (k & 1) B = mul(B, A);A = mul(A, A);k >>= 1;}return B;}int main(){ int t;scanf("%d",&t);while (t--){int sum = 0;mat x; // 在定义一个 mat 结构体In(); // 输入矩阵x = Pow(A, k);for (int i = 0; i < n;i++){sum += x.mat[i][i];sum %= mod;} printf("%d\n", sum % mod);}return 0;}
0 0
- Tr A
- Tr A
- Tr A
- Tr A
- Tr A
- Tr A
- Tr A
- Tr A
- Tr A hdu 1575
- HDOJ 1575 - Tr A
- hdu 1575 Tr A
- hdu1575 Tr A
- Tr A(矩阵)
- 题目59:Tr A
- HDU 1575 Tr A
- Tr A hdu1575
- HDU 1575 Tr A
- hdu 1575 Tr A
- 关于yii的relations
- OC 基本数据类型与对象之间的转换方法
- input_dev结构体分析
- c++第七周作业(第五章)
- TextView文字+drawableLeft居中显示
- Tr A
- 牛腩新闻发布系统之感受
- windows8 安装IIS 和 添加网站
- 写的全球IP地址搜纯ASP(VBscript)
- uva 10594 - Data Flow(最小费用最大流)
- TTP224 4路电容式触摸开关(数字触摸传感器)模块
- 中文操作系统的VBscript中,使用的是unicode字符集
- iOS开发- UICollectionView详解+实例
- DRIVER_OBJECT & DEVICE_OBJECT