动态规划计算矩阵连乘

#include<iostream> 
using namespace std; 

const int MAX = 100; 
int p[MAX+1],m[MAX][MAX],s[MAX][MAX]; 

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


//根据s[][]记录的各个子段的最优解，将其输出
void traceback(int i,int j)
{
	if(i==j)
	{
  		cout<<'M'<<i+1;
	  	return ; 
	}
	if(i<s[i][j])
	  	cout<<'(';
	traceback(i,s[i][j]);
	if(i<s[i][j])
	  	cout<<')';
	if(s[i][j]+1<j)
	  	cout<<'(';
	traceback(s[i][j]+1,j);
	if(s[i][j]+1<j)
	  	cout<<')';
}
int main()
{
	int n;
	cout<<"请输入矩阵的个数："<<endl;
	cin>>n;
	cout<<"输入各个矩阵的维数:"<<endl;
	for(int i=0;i<=n;i++)
   		cin>>p[i];
   	cout<<"各个矩阵的维数为："<<endl; 
	for(int i=1;i<=n;i++)
		cout<<"M"<<i<<": "<<p[i-1]<<'\t'<<p[i]<<endl;   
   	cout<<"输出结果如下："<<endl;
   	matrixChain(n);
	cout<<endl;
	traceback(0,n-1);
	cout<<endl;
	return 0;
}