树形dp（洛谷1040 加分二叉树noip2003提高组第三题）

/*一旦根的位置在一个中序遍历之中确定，它的左边的序列就是它的左子树，右边就是右子树。那么整棵树的最优值就是左子树的最优值乘上右子树的最优值。
状态转移方程：f[i,j]表示中序遍历中i~j这棵子树能得到的最大加分 f[i,j]=max{f[i,k-1]*f[k+1,j]+f[k,k]}（i<=k<=j）
动态规划的顺序是自下而上，类似于区间动态规划。 
边界处理：将所有空子树的值都赋为1。《摘自网络》*/

#include<bits/stdc++.h>
#define maxn 200
using namespace std;
int a[maxn],b[maxn][maxn],dp[maxn][maxn];
int n;

void dfs(int from,int to)//递归打印树 
{
   if(from<=to)   
   {
      cout<<b[from][to]<<" ";
      dfs(from,b[from][to]-1);
      dfs(b[from][to]+1,to);
   }
}

int main()
{
	cin>>n;
	for(int i=0;i<=n+1;i++)
      for(int j=0;j<=n+1;j++)//初始化，注意全部 
        dp[i][j]=1;
	for(int i=1;i<=n;i++)
	{
	  cin>>a[i];
	  dp[i][i]=a[i];
	  b[i][i]=i;
    }
    
    for(int len=2;len<=n;len++)//这里枚举长度才无后效性 
	for(int i=1;i<=n-len+1;i++)
	  {
	      int j=i+len-1;
		  for(int k=i;k<=j;k++)
	      {
			   if(dp[i][j]<dp[i][k-1]*dp[k+1][j]+dp[k][k])
			   {
			     dp[i][j]=dp[i][k-1]*dp[k+1][j]+dp[k][k];
			     b[i][j]=k;
			   }
		  }
	  }
	cout<<dp[1][n]<<endl;
	dfs(1,n);
   return 0;
}