基础练习：矩阵乘法 -深度解析

气场五十米

---

Tip：

矩阵乘法万能算法：

int A[100], B[100], C[100] //计算A*B=C
for(int i=0;i<n;i++)
	for(int j=0;j<n;j++)
		for(int k=0;k<n;k++)
		C[i][j]+=A[i][k]*B[k][j];

---

知道上面的就简单了

求A得N次幂 就是将A*1 累加N次

思路

#include 
using namespace std;
int a[50][50],zero[50][50],c[50][50];
int main()
{
	int n,m;
	cin>>n>>m;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++) cin>>a[i][j];
	for(int i=0;i<n;i++) zero[i][i]=1;
	while(m--)
	{
		memset(c,0,sizeof(c));
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				for(int k=0;k<n;k++) c[i][j]+=zero[i][k]*a[k][j];
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++) zero[i][j]=c[i][j];
	}
	for(int i=0;i<n;i++,cout<<endl)
		for(int j=0;j<n;j++)
		{
			if(j!=0) cout<<" ";
			cout<<zero[i][j];
		}
	return 0;
}