加分二叉树(区间DP，记忆化搜索）

设一个n个节点的二叉树tree的中序遍历为（1,2,3,…,n），其中数字1,2,3,…,n为节点编号。

每个节点都有一个分数（均为正整数），记第i个节点的分数为didi，tree及它的每个子树都有一个加分，任一棵子树subtree（也包含tree本身）的加分计算方法如下：

subtree的左子树的加分 × subtree的右子树的加分 ＋ subtree的根的分数

若某个子树为空，规定其加分为1。叶子的加分就是叶节点本身的分数，不考虑它的空子树。

试求一棵符合中序遍历为（1,2,3,…,n）且加分最高的二叉树tree。

要求输出：

（1）tree的最高加分

（2）tree的前序遍历

输入格式
第1行：一个整数n，为节点个数。 

第2行：n个用空格隔开的整数，为每个节点的分数（0<分数<100）。

输出格式
第1行：一个整数，为最高加分（结果不会超过int范围）。     

第2行：n个用空格隔开的整数，为该树的前序遍历。如果存在多种方案，则输出字典序最小的方案。

数据范围
n<30
输入样例：
5
5 7 1 2 10
输出样例：
145
3 1 2 4 5

区间DP,二叉树的遍历) O(n3)
状态表示：f[i][j] 表示中序遍历是 w[i ~ j] 的所有二叉树的得分的最大值。
状态计算：f[i][j] = max(f[i][k - 1] * f[k + 1][j] + w[k])，
即将f[i][j]表示的二叉树集合按根节点分类，则根节点在 k 时的最大
得分即为 f[i][k - 1] * f[k + 1][j] + w[k]，则f[i][j]即为遍历
 k 所取到的最大值。

在计算每个状态的过程中，记录每个区间的最大值所对应的根节点编号。

那么最后就可以通过DFS求出最大加分二叉树的前序遍历了。

时间复杂度
状态总数是 O(n的2次方)，计算每个状态需要 O(n) 的计算量，
因此总时间复杂度是 O(n的3次方)。

#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

const int N = 30;

int n;
int w[N];
int f[N][N], g[N][N];

void dfs(int l, int r)
{
    if (l > r) return;
    int k = g[l][r];
    cout << k << ' ';
    dfs(l, k - 1);//递归处理左半边
    dfs(k + 1, r);//递归处理右半边
}

int main()
{
    cin >> n;

    for (int i = 1; i <= n; i ++ ) cin >> w[i];

    for (int len = 1; len <= n; len ++ )
        for (int l = 1; l + len - 1 <= n; l ++ )
        {
            int r = l + len - 1;
            if (len == 1) f[l][r] = w[l], g[l][r] = l;
            else
            {
                for (int k = l; k <= r; k ++ )
                {
                    int left = k == l ? 1 : f[l][k - 1];
                    int right = k == r ? 1 : f[k + 1][r];
                    int score = left * right + w[k];
                    if (score > f[l][r])
                    {
                        f[l][r] = score;
                        g[l][r] = k;
                    }
                }
            }
        }

    cout << f[1][n] << endl;

    dfs(1, n);

    return 0;
}

//记忆化搜索
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

const int N = 30;

int n;
int w[N];
int f[N][N], g[N][N];

void dfs2(int l, int r)
{
    if (l > r) return;
    int k = g[l][r];
    cout << k << ' ';
    dfs2(l, k - 1);//递归处理左半边
    dfs2(k + 1, r);//递归处理右半边
}
int dfs1(int l,int r){
	int& v=f[l][r];
	if(v!=-1)return v; 
	if(l>r)return 1;
	if(l==r)return v=w[l];
	
	v=-1e9;
	for(int k=l;k<=r;k++){
		int t=dfs1(l,k-1)*dfs1(k+1,r)+w[k];
		if(t>v){
			v=t;
			g[l][r]=k;
		}
	}
	return v;
} 
int main()
{
    cin >> n;

    for (int i = 1; i <= n; i ++ ) cin >> w[i];

        
    memset(f,-1,sizeof f);
    
    for (int i = 1; i <= n; i ++ ){
        g[i][i]=i;//根就是本身
        f[i][i]=w[i];
    }
    dfs1(1, n);
    
    cout << f[1][n] << endl;

    dfs2(1, n);
    
    return 0;
}

//记忆化搜索,简洁版
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

const int N = 30;

int n;
int w[N];
int f[N][N], g[N][N];

void dfs2(int l, int r)
{
    if (l > r) return;
    int k = g[l][r];
    cout << k << ' ';
    dfs2(l, k - 1);//递归处理左半边
    dfs2(k + 1, r);//递归处理右半边
}
int dfs1(int l,int r){
	int& v=f[l][r];
	if(v)return v; 
	if(l>r)return 1;
    if(l==r){//叶节点
        g[l][r]=r;
        return w[l];
    }
	
	for(int k=l;k<=r;k++){
		int t=dfs1(l,k-1)*dfs1(k+1,r)+w[k];
		if(t>v){
			v=t;
			g[l][r]=k;
		}
	}
	return v;
} 
int main()
{
    cin >> n;

    for (int i = 1; i <= n; i ++ ) cin >> w[i];

   
    dfs1(1, n);
    
    cout << f[1][n] << endl;

    dfs2(1, n);
    
    return 0;
}