2021/3/3 基础练习 矩阵乘法

基础练习 矩阵乘法

【练习时间】2020/3/3
【题目名称】基础练习 矩阵乘法
资源限制
时间限制：1.0s 内存限制：512.0MB
问题描述
　　给定一个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

【解题思路】

1. 特殊情况，矩阵的零次幂是对角矩阵；
2. 矩阵乘法规则是“左行右列”，详情请参见——>矩阵乘除法 。只需要对两个矩阵按行按列逐个计算即可，同时我们需要第三方矩阵帮助存储每次运算的结果，并复制给乘数矩阵再次参与运算。

【源代码】

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
int main(){
    int n,m;
	cin>>n>>m;
	int a[n][n],b[n][n],c[n][n];
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<n;j++)
		{
			cin>>a[i][j];
			b[i][j]=a[i][j];
		}
		
	} 
	if(m==0)//矩阵的零次幂是对角矩阵 
	{
		for(int i=0;i<n;i++)
		{
			for(int j=0;j<n;j++)
			{
				if(i==j)
				b[i][j]=1;
				else b[i][j]=0;
			}
		}
	}
	else{
		for(int x=0;x<m-1;x++)//相乘次数 
		{
			for(int i=0;i<n;i++)//控制第几行 
			{
				for(int k=0;k<n;k++)//控制第几列 
				{
					int sum=0;
					for(int j=0;j<n;j++)//某一列、行的所有数字 
					{
						sum+=a[j][k]*b[i][j];
					} 
					c[i][k]=sum;
				}
			}
			for(int i=0;i<n;i++)
			{
				for(int j=0;j<n;j++)
				b[i][j]=c[i][j];
			}
		}
		
	}
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<n;j++)
		{
			cout<<b[i][j]<<" "; 
		}
		cout<<endl; 
	}
    return 0;
}