题目地址:
https://www.acwing.com/problem/content/322/
在Mars星球上,每个Mars人都随身佩带着一串能量项链,在项链上有 N N N颗能量珠。能量珠是一颗有头标记与尾标记的珠子,这些标记对应着某个正整数。并且,对于相邻的两颗珠子,前一颗珠子的尾标记一定等于后一颗珠子的头标记。因为只有这样,通过吸盘(吸盘是Mars人吸收能量的一种器官)的作用,这两颗珠子才能聚合成一颗珠子,同时释放出可以被吸盘吸收的能量。如果前一颗能量珠的头标记为 m m m,尾标记为 r r r,后一颗能量珠的头标记为 r r r,尾标记为 n n n,则聚合后释放的能量为 m ∗ r ∗ n m*r*n m∗r∗n(Mars单位),新产生的珠子的头标记为 m m m,尾标记为 n n n。需要时,Mars人就用吸盘夹住相邻的两颗珠子,通过聚合得到能量,直到项链上只剩下一颗珠子为止。显然,不同的聚合顺序得到的总能量是不同的,请你设计一个聚合顺序,使一串项链释放出的总能量最大。例如:设 N = 4 N=4 N=4, 4 4 4颗珠子的头标记与尾标记依次为 ( 2 , 3 ) , ( 3 , 5 ) , ( 5 , 10 ) , ( 10 , 2 ) (2,3), (3,5), (5,10), (10,2) (2,3),(3,5),(5,10),(10,2)。我们用记号 ⊕ ⊕ ⊕表示两颗珠子的聚合操作, ( j ⊕ k ) (j⊕k) (j⊕k)表示第 j j j, k k k两颗珠子聚合后所释放的能量。则第 4 4 4、 1 1 1两颗珠子聚合后释放的能量为: ( 4 ⊕ 1 ) = 10 ∗ 2 ∗ 3 = 60 (4⊕1)=10*2*3=60 (4⊕1)=10∗2∗3=60。这一串项链可以得到最优值的一个聚合顺序所释放的总能量为 ( ( 4 ⊕ 1 ) ⊕ 2 ) ⊕ 3 ) = 10 ∗ 2 ∗ 3 + 10 ∗ 3 ∗ 5 + 10 ∗ 5 ∗ 10 = 710 ((4⊕1)⊕2)⊕3)= 10*2*3+10*3*5+10*5*10=710 ((4⊕1)⊕2)⊕3)=10∗2∗3+10∗3∗5+10∗5∗10=710。
输入格式:
输入的第一行是一个正整数
N
N
N,表示项链上珠子的个数。第二行是
N
N
N个用空格隔开的正整数,所有的数均不超过
1000
1000
1000,第
i
i
i个数为第
i
i
i颗珠子的头标记,当
i
<
N
i<N
i<N时,第
i
i
i颗珠子的尾标记应该等于第
i
+
1
i+1
i+1颗珠子的头标记,第
N
N
N颗珠子的尾标记应该等于第
1
1
1颗珠子的头标记。至于珠子的顺序,你可以这样确定:将项链放到桌面上,不要出现交叉,随意指定第一颗珠子,然后按顺时针方向确定其他珠子的顺序。
输出格式:
输出只有一行,是一个正整数
E
E
E,为一个最优聚合顺序所释放的总能量。
数据范围:
4
≤
N
≤
100
4≤N≤100
4≤N≤100
1
≤
E
≤
2.1
∗
1
0
9
1≤E≤2.1∗10^9
1≤E≤2.1∗109
思路是动态规划。可以采用处理环形问题的通用技巧,即复制一份接到后面。这里虽然输入是 N N N个数,但实际上我们要求的是长 N + 1 N+1 N+1的区间的最大能量,所以要将这 N N N个数重复一遍拼在后面,形成长 2 N 2N 2N的数组,这样就转化为了线性问题。设 f [ i ] [ j ] f[i][j] f[i][j]是合并区间 [ i , j ] [i,j] [i,j]所能产生的最大能量,可以按照最后一次合并的中间值来分类,其范围为 [ i + 1 : j − 1 ] [i+1:j-1] [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 ] \max_i f[i][i+N] maxif[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^3) O(n3),空间 O ( n 2 ) O(n^2) O(n2)。