矩阵乘法-1

问题描述

　　给定一个N阶矩阵A，输出A的M次幂（M是非负整数）

　　例如：

　　A =

　　1 2

　　3 4

　　A的2次幂

　　7 10

　　15 22

输入格式

　　第一行是一个正整数N、M（1<=N<=30,0<=M<=5），表示矩阵A的阶数和要求的幂数

　　接下来N行，每行N个绝对值不超过10的非负整数，描述矩阵A的值

输出格式

　　输出共N行，每行N个整数，表示A的M次幂所对应的矩阵。相邻的数之间用一个空格隔开

样例输入

2 2

1 2

3 4

样例输出

7 10

15 22

#include<stdio.h>#include<stdlib.h>int**a,**b;int n,m;void multi(){int i,j,k,temp,**c;c=(int**)malloc(sizeof(int*)*n);for(i=0;i<n;i++)c[i]=(int*)malloc(sizeof(int)*n);for(i=0;i<n;i++){for(j=0;j<n;j++){temp=0;for(k=0;k<n;k++){temp+=b[i][k]*a[k][j];}c[i][j]=temp;}}for(i=0;i<n;i++)for(j=0;j<n;j++)b[i][j]=c[i][j];}int main(){int i,j;scanf("%d %d",&n,&m);a=(int**)malloc(sizeof(int*)*n);b=(int**)malloc(sizeof(int*)*n);for(i=0;i<n;i++){a[i]=(int*)malloc(sizeof(int)*n);b[i]=(int*)malloc(sizeof(int)*n);for(j=0;j<n;j++){scanf("%d",a[i]+j);b[i][j]=a[i][j];}}for(i=0;i<m-1;i++)multi();for(i=0;i<n;i++){for(j=0;j<n;j++)printf("%d ",b[i][j]);putchar('\n');}return 0;}