2013-3-9矩阵幂算法

矩阵乘法的C++实现

原创于 2013-03-09 17:18:19 发布 · 594 阅读

0 ·

CC 4.0 BY-SA版权

C++ 专栏收录该内容

31 篇文章

订阅专栏

本文介绍了一种使用C++编程语言实现矩阵乘法的方法，包括内存分配、初始化、计算过程以及结果输出。

#include <iostream>
using namespace std;
int** multiply(int **p1,int** p2,int n)
{
int **p;
int i,j;
p=new int*[n];
for(i=0;i<n;i++)
p[i]=new int[n];
for(i=0;i<n;i++)
{
  for (j=0;j<n;j++)
  {
   p[i][j]=0;
  }
}
for(i=0;i<n;i++)
{
  for (j=0;j<n;j++)
  {
   int h=i,l=j;
   for(int k=0;k<n;k++)
   {
    p[i][j]+=p1[i][k]*p2[k][j];
   }
  }
}
return p;
}
int main()
{
int T,n,k;
cin>>T;
while(T)
{
  cin>>n>>k;
  int i=0,j=0;
//  int **p=(int *)malloc(sizeof(int)*n)
  int **p;
  p=new int*[n];
  for(i=0;i<n;i++)
  p[i]=new int[n];
  for(i=0;i<n;i++)
  {
   for (j=0;j<n;j++)
   {
    cin>>p[i][j];
   }
  }
  int **result;
  result=new int*[n];
  for(i=0;i<n;i++)
  result[i]=new int[n];
  for(i=0;i<n;i++)
  {
   for (j=0;j<n;j++)
   {
    result[i][j]=p[i][j];
   }
  }
  int t=T;
  int **r=new int[t];
  while(k-1)
  {
    result=multiply(result,p,n);
    k--;
  }
  r[t-T]=result;
  T--;
}
for(i=0;i<n;i++)
{
  for (j=0;j<n;j++)
  {
   cout<<result[i][j]<<" ";
  }
  cout<<endl;
}
return 0;
}