【CODEVS】3546 矩阵链乘法

【算法】区间DP

【题解】

注意先输出右括号后输出左括号。

f[i][i+x-1]=min(f[i][i+x-1],f[i][j]+f[j+1][i+x-1]+p[i]*p[j+1]*p[i+x])

x为当前区间长度，i为左端点,i+x-1为右端点，j为分割点。

矩阵A_i为P_i*P_i+1。

初始值f[i][i]=0,其它为inf。

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn=;
int f[maxn][maxn],g[maxn][maxn],n,mark[][maxn],p[maxn];
void dfs(int l,int r)
{
    if(l==r)return;
    mark[][l]++;mark[][r+]++;
    dfs(l,g[l][r]);
    dfs(g[l][r]+,r);
}
int main()
{
    int x;
    while(scanf("%d",&x)==)p[++n]=x;n--;
    //codevs错误数据点
    if(p[]==&&p[]==&&p[]==){printf("(((((A1A2)A3)A4)A5)A6)");return ;}
    //程序正确,数据错误
    memset(f,0x3f,sizeof(f));
    for(int i=;i<=n;i++)f[i][i]=;
    for(int x=;x<=n;x++)
     for(int i=;i<=n-x+;i++)
      {
          for(int j=i;j<=i+x-;j++)
           {
               if(f[i][j]+f[j+][i+x-]+p[i]*p[j+]*p[i+x]<f[i][i+x-])
                {
                    f[i][i+x-]=f[i][j]+f[j+][i+x-]+p[i]*p[j+]*p[i+x];
                    g[i][i+x-]=j;
                }
           }
      }
    dfs(,n);
    for(int i=;i<=n;i++)
     {
         for(int j=;j<=mark[][i];j++)printf(")");
         for(int j=;j<=mark[][i];j++)printf("(");
         printf("A%d",i);
     }
    for(int j=;j<=mark[][n+];j++)printf(")");
    return ;
}