【ACWing】320. 能量项链

题目地址：

https://www.acwing.com/problem/content/322/

在Mars星球上，每个Mars人都随身佩带着一串能量项链，在项链上有$N$颗能量珠。能量珠是一颗有头标记与尾标记的珠子，这些标记对应着某个正整数。并且，对于相邻的两颗珠子，前一颗珠子的尾标记一定等于后一颗珠子的头标记。因为只有这样，通过吸盘（吸盘是Mars人吸收能量的一种器官）的作用，这两颗珠子才能聚合成一颗珠子，同时释放出可以被吸盘吸收的能量。如果前一颗能量珠的头标记为$m$，尾标记为$r$，后一颗能量珠的头标记为$r$，尾标记为$n$，则聚合后释放的能量为$m\ast r\ast n$（Mars单位），新产生的珠子的头标记为$m$，尾标记为$n$。需要时，Mars人就用吸盘夹住相邻的两颗珠子，通过聚合得到能量，直到项链上只剩下一颗珠子为止。显然，不同的聚合顺序得到的总能量是不同的，请你设计一个聚合顺序，使一串项链释放出的总能量最大。例如：设$N=4$，$4$颗珠子的头标记与尾标记依次为$(2，3),(3，5),(5，10),(10，2)$。我们用记号$\oplus$表示两颗珠子的聚合操作，$(j\oplus k)$表示第$j$，$k$两颗珠子聚合后所释放的能量。则第$4$、$1$两颗珠子聚合后释放的能量为：$(4\oplus 1)=10\ast 2\ast 3=60$。这一串项链可以得到最优值的一个聚合顺序所释放的总能量为$((4\oplus 1)\oplus 2)\oplus 3)=10\ast 2\ast 3+10\ast 3\ast 5+10\ast 5\ast 10=710$。

输入格式：
输入的第一行是一个正整数$N$，表示项链上珠子的个数。第二行是$N$个用空格隔开的正整数，所有的数均不超过$1000$，第$i$个数为第$i$颗珠子的头标记，当$i<N$时，第$i$颗珠子的尾标记应该等于第$i+1$颗珠子的头标记，第$N$颗珠子的尾标记应该等于第$1$颗珠子的头标记。至于珠子的顺序，你可以这样确定：将项链放到桌面上，不要出现交叉，随意指定第一颗珠子，然后按顺时针方向确定其他珠子的顺序。

输出格式：
输出只有一行，是一个正整数$E$，为一个最优聚合顺序所释放的总能量。

数据范围：
$4\le N\le 100$
$1\le E\le 2.1\ast 10^9$

思路是动态规划。可以采用处理环形问题的通用技巧，即复制一份接到后面。这里虽然输入是$N$个数，但实际上我们要求的是长$N+1$的区间的最大能量，所以要将这$N$个数重复一遍拼在后面，形成长$2N$的数组，这样就转化为了线性问题。设$f[i][j]$是合并区间$[i,j]$所能产生的最大能量，可以按照最后一次合并的中间值来分类，其范围为$[i+1:j-1]$，所以有： f [ i ] [ j ] = max ⁡ i + 1 ≤ k ≤ j − 1 { f [ i ] [ k ] + f [ k ] [ j ] + a [ i ] ∗ a [ k ] ∗ a [ j ] } f[i][j]=\max_{i+1\le k\le j-1}\{f[i][k]+f[k][j]+a[i]*a[k]*a[j]\} f[i][j]=i+1≤k≤j−1max​{f[i][k]+f[k][j]+a[i]∗a[k]∗a[j]}最后只需要返回${\max }_{i}f[i][i+N]$即可。代码如下：

#include <iostream>
using namespace std;

const int N = 210;
int n;
int a[N];
int f[N][N];

int main() {
    cin >> n;
    for (int i = 1; i <= 2 * n; i++)
        if (i <= n) cin >> a[i];
        else a[i] = a[i - n];

    int res = 0;
    for (int l = 3; l <= n + 1; l++)
        for (int i = 1; i + l - 1 <= 2 * n; i++) {
            int j = i + l - 1;
            for (int k = i + 1; k <= j - 1; k++)
                f[i][j] = max(f[i][j], f[i][k] + f[k][j] + a[i] * a[k] * a[j]);
            
            if (l == n + 1) res = max(res, f[i][j]);
        }

    cout << res << endl;

    return 0;
}

时间复杂度$O(n^3)$，空间$O(n^2)$。