算法导论 动态规划 矩阵链乘法

#include<iostream>
using namespace std;
const int N=100;
int p
={0};
int m

={0};
int s

={0};
void MATRIX_CHAIN_ORDER(int a)
{
int n=a;
int i,j,k,q,l;
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
m[i][j]=0;
for(l=2;l<=n;l++)
{
for(i=1;i<=n-l+1;i++)
{
j=i+l-1;
m[i][j]=INT_MAX;
for(k=i;k<=j-1;k++)
{
q=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];
if(q<m[i][j])
{
m[i][j]=q;
s[i][j]=k;
}
}
}
}
}
int main()
{
int La=0,i,j;
cout<<"Enter La: ";
cin>>La;
for(i=0;i<=La;i++)
cin>>p[i];
MATRIX_CHAIN_ORDER(La);
for(i=1;i<=La;i++)
for(j=La;j>=i;j--)
{
if(j%6==0)
cout<<endl;
cout<<m[i][j]<<' ';
}
cout<<endl;
for(i=1;i<La;i++)
for(j=La;j>=i+1;j--)
{
if(j%6==0)
cout<<endl;
cout<<s[i][j]<<' ';
}
return 0;
}