动态规划之矩阵连乘问题求解

Input：输入两行标准输入，第一行表示矩阵的个数，第二行是这些矩阵的维数。
Output：输出两行标准输出，第一行为这些矩阵的最小数乘次数，第二行为这些矩阵加括号后的表达式。

输入样例1：

6
30 35 15 5 10 20 25

输出样例1：

15125
((A1(A2A3))((A4A5)A6))

具体使用的数据结构：
输入n表示矩阵的个数
int *p = new int[n+1]用来表示矩阵的维数。例第一个矩阵的维数是(p[0],p[1])，第i个矩阵的维数是(p[i-1],p[i])

二维数组m[i][j]用来表示第i个矩阵到第j个矩阵所需要的最小数乘次数
二维数组s[i][j]=k用来表示第i个矩阵到第j个矩阵在第k的矩阵后断开。

子问题的递归关系：

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
using namespace std;
int m[101][101];//m[i][j]表示矩阵i到j需要的最小数乘次数
int s[101][101];//s[i][j]表示矩阵i到j从哪里断开
void MatrixChain(int n,int *p)
{
    for(int i=1;i<=n;i++){
        m[i][i] = 0;//第一次遍历，对角线元素全部置0
    }
    for(int r=1;r<n;r++)
        for(int i=1;i<=n-r;i++){
            int j=i+r;
            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 temp = m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];
                if(temp<=m[i][j])
                    {m[i][j]=temp;s[i][j]=k;}
            }
        }
}
//递归打印出最优解表达式
void PrintAnswer(int i,int j)
{
    if(i==j) cout<<"A"<<i;
    else{
        cout<<"(";
        PrintAnswer(i,s[i][j]);
        PrintAnswer(s[i][j]+1,j);
        cout<<")";
    }
}
int main()
{
    int n;//表示有n个矩阵相乘
    cin>>n;
    int *p = new int[n+1];//数组下标从0开始,n个矩阵有n+1个维度参数
    for(int i=0;i<=n;i++){
        cin>>p[i];//依次输入矩阵维度
    }
    MatrixChain(n,p);
    cout<<m[1][n]<<endl;
    PrintAnswer(1,n);
    return 0;
}

运行结果：

总结：