2045: 矩阵乘法

2045: 矩阵乘法
时间限制: 1 Sec  内存限制: 128 MB
提交: 128  解决: 34
[提交][状态][讨论版]
题目描述
给定一个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 
提示
来源
基础练习

[提交][状态][讨论版]
中文  English  

心得：一开始不是超时就是不对，后来问了同学才知道可以用递归做，而且还要注意m是0的情况。

#include<bits/stdc++.h>
using namespace std;
int a[35][35],b[35][35],c[35][35];
int n;
void f(int m)
{
	if(m==1) return ;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			c[i][j]=0;
			for(int k=1;k<=n;k++)
			{
				c[i][j]+=b[i][k]*a[k][j];
			}
		}
	}
	memcpy(b,c,sizeof(c));
	f(m-1);
}
int main(void)
{
	int i,j,k,m;
	cin>>n>>m;
	for(i=1;i<=n;i++)
	for(j=1;j<=n;j++)
	cin>>a[i][j],b[i][j]=a[i][j];
	if(m==0)
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=n;j++)
			if(i==j) cout<<1<<" ";
			else cout<<0<<" ";
			cout<<endl;
		}
	else
	{
		f(m);
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=n;j++)
			{
				cout<<b[i][j]<<" ";
			}
			cout<<endl;
		}
	}
	return 0;
}