矩阵连乘问题

矩阵连乘问题

应用动态规划算法,将前面计算的重复值保存,一次来推进下一步的计算,直到得到m[1][n]的值,为最少的乘法次数

 

m[i][j]:表示计算A[i...j]所需的最少数乘次数

 

m[i][j]=min(i=<k<j){m[i][k]+m[k+1][j]+p[2*i]p[2*k+1]p[2*j+1]};

s[i][j]记录断开位置,也就是最少计算次数的中间位置。

 

 

 

 

#include<iostream>
using namespace std;

void MatrixChain(int *p,int n,int m[][6],int s[][6])
{
 int i,j;
 for(i=0;i<6;i++)
 {
  m[i][i]=0;
  s[i][i]=0;
 }

 int k=5;
 for(k=1;k<=5;k++)
 {
  for(i=0;i<n-k;i++)
  {
   int temp=0;
   m[i][i+k]=m[i+1][i+k]+p[2*i]*p[2*i+1]*p[2*(i+k)+1];
   s[i][i+k]=i;
   for(j=i+1;j<i+k;j++)
   {
    temp=m[i][j]+m[j+1][i+k]+p[2*i]*p[2*j+1]*p[2*(i+k)+1];
    if(m[i][i+k]>temp)
    {
     m[i][i+k]=temp;
     s[i][i+k]=j;
    }
   }
  }
 }
 
}

void Traceback(int i,int j,int s[][6])
{
 if(i==j) return;
 Traceback(i,s[i][j],s);
 Traceback(s[i][j]+1,j,s);
 cout<<"Multiply A"<<i<<","<<s[i][j];
 cout<<"and A"<<(s[i][j]+1)<<","<<j<<endl;
}

int main(void)
{
 int p[12]={30,35,35,15,15,5,5,10,10,20,20,25};
 int n=6,i,j;
 int m[6][6]={0};
 int s[6][6]={0};
 MatrixChain(p,n,m,s);
 for(i=0;i<6;i++)
 { 
  for(j=0;j<6;j++)
  {
   cout<<m[i][j]<<"\t";
  }
  cout<<endl;
 }
 for(i=0;i<6;i++)
 { 
  for(j=0;j<6;j++)
  {
   cout<<s[i][j]<<"\t";
  }
  cout<<endl;
 }

 Traceback(0,5,s);
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值