矩阵乘法

Description

给定一个N阶矩阵A，输出A的M次幂（M是非负整数）
例如：
A =
1 2
3 4
A的2次幂
7 10
15 22

Input

第一行是一个正整数N、M（1<=N<=30, 0<=M<=5），表示矩阵A的阶数和要求的幂数
接下来N行，每行N个绝对值不超过10的非负整数，描述矩阵A的值

Output

输出共N行，每行N个整数，表示A的M次幂所对应的矩阵。每个数后输出一个空格。

Sample Input

2 2
1 2
3 4

Sample Output

7 10 
15 22 

#include<bits/stdc++.h>
using namespace std;
int main()
{
	int n,m;  
	int a[35][35],b[35][35],c[35][35];
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		int i,j,k,x;
        memset(a,0,sizeof(a));
        memset(b,0,sizeof(b));
        memset(c,0,sizeof(c));
		for(i=0;i<n;i++)
		{
			for(j=0;j<n;j++)
			scanf("%d",&a[i][j]);
		}
		memcpy(b,a,sizeof(a));   
		if(m==0)//0的矩阵的主对角线为1 其他为0 
		{
			for(i=0;i<n;i++)
			{
				for(j=0;j<n;j++)
				if(i==j)
				a[i][j]=1; 
				else
				a[i][j]=0;
			 } 
			 for(i=0;i<n;i++)
			{
				for(j=0;j<n;j++)
			  {
			  	printf("%d ",a[i][j]);	
			  }
			  printf("\n");
			}	
			 
		 } 
		 else
		 {
		 	for(x=1;x<m;x++)
			{
				memset(c,0,sizeof(c));
					for(i=0;i<n;i++)
				{
					for(j=0;j<n;j++)
					{
						for(k=0;k<n;k++)
						c[i][j]+=a[i][k]*b[k][j];
					}
				}
				memcpy(b,c,sizeof(c));  
			} 	
			for(i=0;i<n;i++)
			{
				for(j=0;j<n;j++)
			  {
			  	printf("%d ",b[i][j]);	
			  }
			  printf("\n");
			}	
		 }	    
	}
	return 0;
}

0次幂
1 0
0 1

注意输出格式 这道题是每个输出末尾都是空格