矩阵链乘法

A1,A2....An，n个矩阵，矩阵Ai的维数为p[i-1]*p[i];

数据矩阵链乘法加全部括号，使得计算乘法的次数最少

对矩阵Ai，Ai+1,其维数分别为p[i-1]*p[i],p[i]*p[i+1]，Ai*Ai+1要执行乘法运算的次数为p[i-1]*p[i]*p[i+1];

m[i][j]为矩阵Ai*Ai+1*...*Aj最少的执行乘法运算的次数

假设我们在Ak后面加括号(i<=k<j),那么m[i][j]=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];

观察递推式发现m[i][j]仅取决于长度比它小的链，因此， 我们按长度由小到大求m[i][j]即可

#include <iostream>
#include<cstdio>
#define N 1005
#define INF 1<<29
using namespace std;
int p[N],m[N][N],s[N][N];
void Matrix_chain_order(int n){
    int i,j,k,l,q;
    for(i=1;i<=n;i++)
    m[i][i]=0;
    for(l=2;l<=n;l++){
        for(i=1;i<=n-l+1;i++){
            j=i+l-1;
            m[i][j]=INF;
            for(k=i;k<j;k++){
                q=m[i][k]+m[k][j+1]+p[i-1]*p[k]*p[j-1];
                if(q<m[i][j]){
                    m[i][j]=q;
                    s[i][j]=k;
                }
            }
        }
    }
}
void print_optimal_parens(int i,int j){
    if(i==j)
    printf("A%d",i);
    else {
        printf("(");
        print_optimal_parens(i,s[i][j]);
        print_optimal_parens(s[i][j]+1,j);
        printf(")");
    }
}
int main()
{
    int n,i;
    while(~scanf("%d",&n)){
        for(i=1;i<=n;i++){
            scanf("%d%d",&p[i-1],&p[i]);
        }
        Matrix_chain_order(n);
        print_optimal_parens(1,n);
    }
    return 0;
}