#include<iostream>
#define N 100
using namespace std;
int p[N];
int m[N][N];
int s[N][N];
int n;
void MatrixChain()
{
for(int i=1;i<=n;i++)
{
m[i][i]=0;
}
for(int r=2;r<=n;r++)
{
for(int i=1;i<=n-r+1;i++)
{
int j=i+r-1;
m[i][j]=m[i+1][j]+p[i-1]*p[i]*p[j];
//为了后面比较大小,而付的初始值;
s[i][j]=i;
//断开点位置;
for(int k=i+1;k<j;k++)
{
int t=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];
if(t<m[i][j])
{
m[i][j]=t;
s[i][j]=k;
}
}
}
}
}
void Traceback(int i,int j)
{
if(i==j)
{
cout<<"A["<<i<<"]";
return;
}
cout<<"(";
Traceback(i,s[i][j]);
Traceback(s[i][j]+1,j);
cout<<")";
}
int main()
{
cout<<"请输入矩阵的个数n : "<<endl;
cin>>n;
int i,j;
cout<<"请依次输入每个矩阵的行数和最后一个矩阵的列数:"<<endl;
for(i=0;i<=n;i++)
{
cin>>p[i];
}
MatrixChain();
Traceback(1,n);
cout<<endl;
cout<<"最小计算量的值为:"<<m[1][n]<<endl;
return 0;
}
矩阵连乘-动态规划
最新推荐文章于 2020-03-31 16:27:03 发布