标签：
时间限制：1 秒
内存限制：128 兆
特殊判题：否
提交：2864
解决：1112
给定一个n*n的矩阵，求该矩阵的k次幂，即P^k。
输入包含多组测试数据。 数据的第一行为一个整数T(0<T<=10)，表示要求矩阵的个数。 接下来有T组测试数据，每组数据格式如下： 第一行：两个整数n（2<=n<=10）、k（1<=k<=5），两个数字之间用一个空格隔开，含义如上所示。 接下来有n行，每行n个正整数，其中，第i行第j个整数表示矩阵中第i行第j列的矩阵元素Pij且（0<=Pij<=10）。另外，数据保证最后结果不会超过10^8。
对于每组测试数据，输出其结果。格式为： n行n列个整数，每行数之间用空格隔开，注意，每行最后一个数后面不应该有多余的空格。
3 2 2 9 8 9 3 3 3 4 8 4 9 3 0 3 5 7 5 2 4 0 3 0 1 0 0 5 8 5 8 9 8 5 3 9 6 1 7 8 7 2 5 7 3
153 96 108 81 1216 1248 708 1089 927 504 1161 1151 739 47 29 41 22 16 147 103 73 116 94 162 108 153 168 126 163 67 112 158 122 152 93 93 111 97
#include<stdio.h> #include<stdlib.h> typedef struct Matrix { int m[11][11]; int n; }Matrix; Matrix multiply(Matrix a,Matrix b) { int i,j,k; Matrix c; c.n=a.n; for(i=1;i<=a.n;i++) for(j=1;j<=b.n;j++) { c.m[i][j]=0; for(k=1;k<=a.n;k++) c.m[i][j]+=a.m[i][k]*b.m[k][j]; } return c; } int main() { int total; while(scanf("%d",&total)!=EOF) { while(total--) { int n,k; scanf("%d%d",&n,&k); int flag=0; Matrix matrix,matrixa; matrix.n=n; int i,j; for(i=1;i<=n;i++) for(j=1;j<=n;j++) scanf("%d",&matrix.m[i][j]); matrixa=matrix; if(k&1&&k!=1) flag=1; while(k!=1) { matrix=multiply(matrix,matrix); k/=2; } if(flag) matrix=multiply(matrix,matrixa); for(i=1;i<=matrix.n;i++) { for(j=1;j<=matrix.n;j++) { printf("%d",matrix.m[i][j]); if(j!=matrix.n) printf(" "); } printf("\n"); } } } return 0; }
九度OJ 1474 矩阵幂（高效算法）
原文地址：http://blog.csdn.net/ruzhuxiaogu/article/details/43950773
