AcWing 320. 能量项链（区间dp）

题意：

给定 N 个能量石，每个能量石有一个二元属性 (wi,1，wi,2)

其中 wi,1 表示第 i 个能量石和第 i−1 个能量石融合产生的 能量 的其中一个参数

当然 wi,2 表示第 i 个能量石和第 i+1 个能量石融合产生的 能量 的其中一个参数

魔法石是顺序且环形摆放的，每次可以融合相邻两个魔法石

融合两个能量石 i,i+1 所产生的能量为（题目保证，相邻能量石的参数一致，首尾一致）

$Ei,i+1=wi,1\times (wi,2或wi+1,1)\times wi+1,1$

融合后 左侧魔法石的第一个参数 和 右侧魔法石的第二个参数 合并成为一颗新的 魔法石

具体如下图所示：

求最终把 所有石头融合成一个石头时，产生的最大能量值

思路：

本题可以把区间长度len作为搜索的阶段来进行记忆化搜索，因此我们也可以采用 区间DP 的方式来处理

这题和 环形石子合并 十分相似，但又不尽相同

在 环形石子合并 中，每个石头只有单一的参数，而本题有两个参数，且最后一步是 将长度len为 3的区间合并成长度为2的区间，这也意味着我们需要在细节上做出改变

经过观察我们发现，合并两个石头 (a,b)，(b,c) 的操作就像是 矩阵乘法 一样，合并完后就变成了 (a,c)
因此我们可以离散的来存储每个参数，具体如下所示：

这样 状态表示 就更新为：当前合并的石子堆的左端石头的左参数是 l，右端石头的右参数是 r 的方案

这样对应的 初始状态 本来应该是 1个二元属性的石头，现在就变成了 长度为 2 的区间

这样合并区间后，需要记录的 新石头的参数 也刚好是 区间的两端对应的参数，如下图所示：

而且这里我们的转移方程也要修改为 $fl,r=max(fl,k+fk,r+El,r)$
以往的 区间DP 我们是把 区间 [a,b] 拆分为 [a,k] 和 [k+1,b]

因为 同一个石子 只会被合并到 一个石子堆 里

注意，本题在合并魔法石时，分割点 k 要被分到 左侧石子堆的右端点 和 右侧石子堆的左端点 中

因此，参数 k 要作为 两个区间的共同端点 来使用，即 [a,k] 和 [k,b]
此外我们原来只需要合并 n 个石头，这样转换后就要合并 n+1 个石头了 （本题可以经过一个周期回到起点，因此len要循环枚举到n+1）

具体分析如下所示：

$fl,r$状态表示：

集合：所有 当前合并的石子堆的左端石头的左参数是 l，右端石头的右参数是 r 的方案

属性：最大费用

$fl,r$状态计算：

$fl,r=max(fl,k+fk,r+El,r)(l<k<r)$

初始状态： $fl,l+1=0(1\le l\le n)$

目标状态： $f1,n+1$

关于如何解决 环的问题 不做额外阐述，有需要的可以参考之前写的这篇 环形石子合并

代码：

#include<bits/stdc++.h>

using namespace std;

const int N = 110, M = N<<1, inf = 0x3f3f3f3f;
int n;
int dp[M][M], w[M];
int mul[M];

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

//长度len应该从3开始枚举，因为长度是2的时候表示只有一个矩阵，不能释放能量，置为0，
//而且由于本题可以经过一个周期回到起点，因此len要循环枚举到n+1
        for(int len = 3; len <= n+1; ++len)
        {
                for(int l = 1; l + len - 1 <= n << 1; ++l)
                {
                        int r = l + len - 1;
                        for(int k = l+1; k < r; ++k)
                        {
                                dp[l][r] = max(dp[l][r], dp[l][k] + dp[k][r] + w[l] * w[k] * w[r]);
                        }
                }
        }

        int MAX = -inf;
        for(int i=1; i<=n; ++i)
        {
                MAX = max(MAX, dp[i][i+n]);
        }
        cout<<MAX<<'\n';

        return 0;
}