动态规划算法解矩阵连乘（动态规划）

【问题描述】使用动态规划算法解矩阵连乘问题，具体来说就是，依据其递归式自底向上的方式进行计算，在计算过程中，保存子问题答案，每个子问题只解决一次，在后面计算需要时只要简单查一下得到其结果，从而避免大量的重复计算，最终得到多项式时间的算法。

【输入形式】在屏幕上输入第1个矩阵的行数和第1个矩阵到第n个矩阵的列数，各数间都以一个空格分隔。

【输出形式】矩阵连乘A1...An的最少数乘次数和最优计算次序。

【样例1输入】

30 35 15 5 10 20 25

【样例1输出】

15125

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

【样例说明】

 输入：第1个矩阵的行数和第1个矩阵到第6个矩阵的列数，以一个空格分隔。

 输出：矩阵连乘A1...An的最少数乘次数为15125，最优计算次序为((A1(A2A3))((A4A5)A6))。

【评分标准】根据输入得到准确的输出。

代码实现：

#include<bits/stdc++.h>
using namespace std;
int m[100][100]={0};
int s[100][100]={0};
void matrixchain(int *p,int n)
{
    int q,j,k;
    for(int i=1;i<=n;i++)
    m[i][i]=0;
    for(int len=2;len<=n;len++)
    {
        for(int i=1;i<=n-len+1;i++)
        {
            j=i+len-1;
            m[i][j]=1000000;
            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;
                }
            }
        }
    }
}
void optimalparens(int i,int j)
{
    if(i==j)
    cout<<"A"<<i;
    else
    {
        cout<<"(";
        optimalparens(i,s[i][j]);
        optimalparens(s[i][j]+1,j);
        cout<<")";
    }
}
int main()
{
    int number,p[100],n=0;
    while(cin>>number)
    {
        p[n]=number;
        n++;
    }
    matrixchain(p,n-1);
    cout<<m[1][n-1]<<endl;

    optimalparens(1,n-1);
    return 0;
}