矩阵连乘

#include<iostream>
using namespace std;

#define NUM 51

int p[NUM];   //矩阵连乘阶数
int m[NUM][NUM];  //每步的最优值
int s[NUM][NUM];    //断开处
void MatrixChain(int n)
{
    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;
    else
    {
        cout<<"(";
        TraceBack(i,s[i][j]);
        TraceBack(s[i][j]+1,j);
        cout<<")";
    }
}
int main()
{
    p[0]=50;
    p[1]=10;
    p[2]=40;
    p[3]=30;
    p[4]=5;
    p[5]=20;
    p[6]=15;
    MatrixChain(6);
    for(int i=1; i<7; i++)
    {
        for(int j=1; j<7; j++)
        {
            cout<<m[i][j]<<" ";
        }
        cout<<endl;
    }

    for(int i=1; i<7; i++)
    {
        for(int j=1; j<7; j++)
        {
            cout<<s[i][j]<<" ";
        }
        cout<<endl;
    }

    TraceBack(1,6);
    return 0;
}