注意事项:
本题是"区间DP—环形石子合并"的扩展题,处理“环”的思路完全一样,就不多详细讲了,可以去那篇文章看,下面讲一下这题的dp。
题目:
在 Mars 星球上,每个 Mars 人都随身佩带着一串能量项链,在项链上有 N 颗能量珠。
能量珠是一颗有头标记与尾标记的珠子,这些标记对应着某个正整数。
并且,对于相邻的两颗珠子,前一颗珠子的尾标记一定等于后一颗珠子的头标记。
因为只有这样,通过吸盘(吸盘是 Mars 人吸收能量的一种器官)的作用,这两颗珠子才能聚合成一颗珠子,同时释放出可以被吸盘吸收的能量。
如果前一颗能量珠的头标记为 m,尾标记为 r,后一颗能量珠的头标记为 r,尾标记为 n,则聚合后释放的能量为 m×r×n(Mars 单位),新产生的珠子的头标记为 m,尾标记为 n。
需要时,Mars 人就用吸盘夹住相邻的两颗珠子,通过聚合得到能量,直到项链上只剩下一颗珠子为止。
显然,不同的聚合顺序得到的总能量是不同的,请你设计一个聚合顺序,使一串项链释放出的总能量最大。
例如:设 N=4,4 颗珠子的头标记与尾标记依次为 (2,3)(3,5)(5,10)(10,2)。
我们用记号 ⊕ 表示两颗珠子的聚合操作,(j⊕k) 表示第 j,k 两颗珠子聚合后所释放的能量。则第 4、1 两颗珠子聚合后释放的能量为:(4⊕1)=10×2×3=60。
这一串项链可以得到最优值的一个聚合顺序所释放的总能量为 ((4⊕1)⊕2)⊕3)=10×2×3+10×3×5+10×5×10=710。
输入格式
输入的第一行是一个正整数 N,表示项链上珠子的个数。
第二行是 N 个用空格隔开的正整数,所有的数均不超过 1000,第 i 个数为第 i 颗珠子的头标记,当 i<N 时,第 i 颗珠子的尾标记应该等于第 i+1 颗珠子的头标记,第 N 颗珠子的尾标记应该等于第 1 颗珠子的头标记。
至于珠子的顺序,你可以这样确定:将项链放到桌面上,不要出现交叉,随意指定第一颗珠子,然后按顺时针方向确定其他珠子的顺序。
输出格式
输出只有一行,是一个正整数 E,为一个最优聚合顺序所释放的总能量。
数据范围
4≤N≤100,
1 ≤ E ≤ 2.1×1e9
输入:
4
2 3 5 10
输出:
710
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 210;
int n, m, w[N]; //w[i]为第i个元素的值
int f[N][N];
int main() {
//读入,并将链复制接到尾部,例如123,变为123123
cin >> n;
for (int i = 1; i<=n; i++) {
cin >> m;
w[i] = w[i+n] = m;
}
//经典的区间dp写法(枚举区间长度->枚举左端点->枚举分界线)
for (int len = 3; len <= n+1; len++) { //从3个元素开始因为最少需要两个矩阵才能释放能量,例如 123=(1,2)(2,3)
for (int l = 1; l+len-1 <= n*2; l++) { //处理"环"开2*n
int r = l+len-1;
for (int k = l+1; k < r; k++){ //同理, 一个矩阵是需要2个元素才能确定的, 所以分界线从l+1开始
f[l][r] = max(f[l][r], f[l][k] + f[k][r] + w[l]*w[k]*w[r]);
}
}
}
//枚举每种缺口的情况(通过链复制接到尾部解决“环”),取最大值即为答案
int res = 0;
for (int i = 1; i<=n; i++) res = max(res, f[i][i+n]);
cout << res;
return 0;
}
思路:
经典y式dp法
1.状态表示
f[L][R]
: 对于L到R这个区间中的所有珠子(包含L和R)的所有合并方法,
属性为Max(释放的能量值)
2.状态计算
很明显如果直接按题里给的元组那样存储太麻烦了,可以直接将数据化简:
比如题目中给的例子: (2,3)(3,5)(5,10)(10,2)
实际上有效数字只有5个:2,3,5,10,2
那只存储这5个数字即可,需要用的时候直接2和3,3和5,5和10…即可。
然后处理“环”的方法这里就不说了,注意事项里提过了。
最后就是状态计算,按照分界线K
:
1.L到K的所有珠子释放最大的能量值为f[L][K]
2.K到R的所有珠子释放最大的能量值为f[K][R]
(注意这里是K而不是K+1,因为数组存的是有效数字而不是元组)
3.合并f[L][K]和f[K][R]两个珠子释放的能量值为w[L]*w[K]*w[R]
如果有所帮助请给个免费的赞吧~有人看才是支撑我写下去的动力!
声明:
算法思路来源为y总,详细请见https://www.acwing.com/
本文仅用作学习记录和交流