矩阵乘法模拟 北邮机试 矩阵幂

#include<bits/stdc++.h>
using namespace std;
int a[10][10];
int b[10][10];
int c[10][10];
int d[10][10];
int n,kk;
int main()
{
	//freopen("kkk.txt","r",stdin);
	cin>>n>>kk;
	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++) c[i][i]=1;
	while(kk>0)
	{
	    fill(d[0],d[0]+10*10,0);
			for(int i=0;i<n;i++)//控制左行 
			    for(int k=0;k<n;k++)//控制右列 
	                for(int j=0;j<n;j++)//控制走步 
	                    d[i][k]+=c[i][j]*a[j][k];
	        
	        for(int i=0;i<n;i++)
	            for(int j=0;j<n;j++)
	                c[i][j]=d[i][j];//c矩阵存储最后结果 
	    kk--;
    }
    for(int i=0;i<n;i++)
	{
	    for(int j=0;j<n;j++) cout<<c[i][j]<<" ";
	    cout<<endl;
	}
	return 0;
}