矩阵连乘（C++实现）

矩阵连乘问题：n个二维矩阵A1、A2……An。请利用结合律改变计算次序，使得总的乘法计算量最小。

懒了，只放代码

#include <iostream>  
using namespace std;
#define MAX 50
int p[MAX+1];             //存储各个矩阵的列数以及第一个矩阵的行数（作为第0个矩阵的列数） 
int m[MAX][MAX];          //m[i][j]存储子问题的最优解  
int s[MAX][MAX];           //s[i][j]存储子问题的最佳分割点
int n;                    //矩阵个数
 
 void matrix(int n,int m[][n],int s[][n],int p[])
{  
   
    int i,j,k;  
    for(i=0;i<n;i++)  
        m[i][i]=0;                  //最小子问题仅含有一个矩阵 ，对角线全为0 
                               
    for(i=2;i<=n;i++)  
        for(j=0;j<n-i+1;j++)
        {                         
            m[j][j+i-1]=INT_MAX;   
            for(k=0;k<i-1;k++)
            {                                    //k代表分割点  
                if(m[j][j+i-1]>m[j][j+k]+m[j+k+1][j+i-1]+p[j]*p[j+k+1]*p[j+i])
                {  
                    m[j][j+i-1]=m[j][j+k]+m[j+k+1][j+i-1]+p[j]*p[j+k+1]*p[j+i];  
                    s[j][j+i-1]=k;                           //记录分割点  
                }  
            }  
        }  
}  
  
 void printmatrix(int leftindex,int rightindex)//递归打印输出  
{  
    if(leftindex==rightindex)  
        cout<<leftindex<<endl;  
    else{  
        cout<<"(" 
        printmatrix(leftindex,leftindex+s[leftindex][rightindex]);  
        printmatrix(leftindex+s[leftindex][rightindex]+1,rightindex);  
        cout<<")"<<endl;  
    }  
}  
int main()  
{  
    int i;
    cout<<"请输入矩阵相乘的矩阵个数:"<<endl;
    cin>>n;
    cout<<"请依次输入矩阵的行和烈（如A*B，A=20*30，B=30*40，即输入20 30 40)"<<endl;
    for(i=0;i<n+1;i++)
    {
        cin>>p[i])
    }
    matrix(n,m,s,p);  
    cout<<"矩阵连乘最小次数:"<<m[0][n-1]<<endl;  
    printmatrix(0,n-1);  
    cout<<endl;
    return 0;  
}