BUPT OJ96 矩阵幂

题目描述

​给你一个n*n的矩阵，

，

求其矩阵的k次幂，即Pk

输入格式

第一行，一个整数T（0<T<=10），表示要求矩阵的个数。

接下来有T组数据，每组数据格式如下：

第一行：两个数据n（2<=n<=10）、k（1<=k<=5），两个数字之间用一个空格隔开，其中n表示状况空间的总数，k表示待求的转移概率矩阵的步数。接下来有n行n列个正整数，其中，第i行第j列表示pij，（0<=pij<=10）。另外，数据保证最后结果不会超过10^8。

输出格式

输出为T组数据。

每组数据为已知矩阵的k次幂，格式为：

n行n列个正整数，每行数之间用空格隔开，注意，每行最后一个数后面不应该有多余的空格。

输入样例

32 29 89 33 34 8 49 3 03 5 75 24 0 3 0 10 0 5 8 58 9 8 5 39 6 1 7 87 2 5 7 3

输出样例

153 96108 811216 1248 7081089 927 5041161 1151 73947 29 41 22 16147 103 73 116 94162 108 153 168 126163 67 112 158 122152 93 93 111 97

矩阵幂的运算, 鉴于这题数据小的可怜, 就直接暴力了

下次会补一个高效的代码

注意格式输出

/*USER_ID: test#birdstormPROBLEM: 96SUBMISSION_TIME: 2014-03-01 23:10:01*/#include <stdio.h>#include <stdlib.h>#define For(i,m,n) for(i=m;i<n;i++) int a[11][11], b[11][11], c[11][11], n, k; void mul(){    int h, i, j, temp;    while(--k){        For(h,0,n) For(i,0,n){            temp=0;            For(j,0,n) temp+=b[h][j]*a[j][i];            c[h][i]=temp;        }        For(h,0,n) For(i,0,n) b[h][i]=c[h][i];    }} main(){    int i, j, t;    scanf("%d",&t);    while(t--){        scanf("%d%d",&n,&k);        For(i,0,n) For(j,0,n) scanf("%d",&a[i][j]),c[i][j]=b[i][j]=a[i][j];        mul();        For(i,0,n) For(j,0,n) printf("%d%c",c[i][j],j==n-1?'\n':' ');    }    return 0;}