Description

期末考试结束了,班主任L老师要将成绩单分发到每位同学手中。L老师共有n份成绩单,按照编号从1到n的顺序叠

放在桌子上,其中编号为i的成绩单分数为w_i。成绩单是按照批次发放的。发放成绩单时,L老师会从当前的一叠

成绩单中抽取连续的一段,让这些同学来领取自己的成绩单。当这批同学领取完毕后,L老师再从剩余的成绩单中

抽取连续的一段,供下一批同学领取。经过若干批次的领取后,成绩单将被全部发放到同学手中。然而,分发成绩

单是一件令人头痛的事情,一方面要照顾同学们的心理情绪,不能让分数相差太远的同学在同一批领取成绩单;另

一方面要考虑时间成本,尽量减少领取成绩单的批次数。对于一个分发成绩单的方案,我们定义其代价为:

BZOJ4897: [Thu Summer Camp2016]成绩单【DP of DP】-LMLPHP

其中,k是方案中分发成绩单的批次数,对于第i批分发的成绩单,〖max〗_i是最高分数,〖min〗_i是最低分数。

a,b是给定的评估参数。现在,请你帮助L老师找到代价最小的分发成绩单的方案,并将这个最小的代价告诉L老师

。当然,分发成绩单的批次数k是由你决定的。

Input

第一行包含一个正整数n,表示成绩单的数量。

第二行包含两个非负整数a,b,表示给定的评估参数。

第三行包含n个正整数w_i,表示第i张成绩单上的分数。

Output

仅一个正整数,表示最小的代价是多少。

Sample Input

10

3 1

7 10 9 10 6 7 10 7 1 2

Sample Output

15

【样例数据说明】

第1批:第2至4份成绩单,落差值为1,剩余成绩单为76710712;

第2批:第4份成绩单,落差值为0,剩余成绩单为767712;

第3批:第1至4份成绩单,落差值为1,剩余成绩单为12;

第4批:剩余的2份成绩单,落差值为1。

总代价为4×3+(12+12)×1=15。

HINT

n<=50, a<=100, b<=10, w_i<=1000


思路

首先考虑dp

\(f_{l,r}\)表示把区间\([l,r]\)的所有成绩单删除的代价

然后就不会了,考虑整个区间,要么递归成两个子区间分别处理,或者直接考虑这个区间的最后一次删除

那么最后一次删除的如果是\([p,q]\)这个区间的数,记录\(g_{l,r,p,q}\)表示区间\([l,r]\)最后只剩下\([p,q]\)的最小代价

然后f的转移可以枚举\(p,q\)从g转移也可以划分成两个区间递归成f转移

g可以划分成两个区间,分别从\(g_{l,k,p,q}+g_{k+1,r,p,q},g_{l,k,p,q}+f_{k+1,r},f_{l,k}+g_{k+1,r,p,q}\)转移


#include<bits/stdc++.h>

using namespace std;

const int N = 51;
const int INF_of_int = 1e8; int f[N][N], g[N][N][N][N];
int a, b, n, m, w[N], pre[N]; int getg(int l, int r, int down, int up);
int getf(int l, int r); int getg(int l, int r, int down, int up) {
if (~g[l][r][down][up]) return g[l][r][down][up];
if (l == r) return g[l][r][down][up] = (down <= w[l] && w[l] <= up) ? 0 : INF_of_int;
int &cur = g[l][r][down][up];
cur = INF_of_int;
for (int k = l; k < r; k++) {
cur = min(cur, getg(l, k, down, up) + getg(k + 1, r, down, up));
cur = min(cur, getg(l, k, down, up) + getf(k + 1, r));
cur = min(cur, getf(l, k) + getg(k + 1, r, down, up));
}
return cur;
} int getf(int l, int r) {
if (~f[l][r]) return f[l][r];
if (l == r) return f[l][r] = a;
int &cur = f[l][r];
cur = INF_of_int;
for (int i = 1; i <= m; i++)
for (int j = i; j <= m; j++)
cur = min(cur, getg(l, r, i, j) + a + b * (pre[j] - pre[i]) * (pre[j] - pre[i]));
for (int k = l; k < r; k++)
cur = min(cur, getf(l, k) + getf(k + 1, r));
return cur;
} int main() {
#ifdef dream_maker
freopen("input.txt", "r", stdin);
#endif
memset(f, -1, sizeof(f));
memset(g, -1, sizeof(g));
scanf("%d %d %d", &n, &a, &b);
for (int i = 1; i <= n; i++) {
scanf("%d", &w[i]);
pre[i] = w[i];
}
sort(pre + 1, pre + n + 1);
m = unique(pre + 1, pre + n + 1) - pre - 1;
for (int i = 1; i <= n; i++)
w[i] = lower_bound(pre + 1, pre + m + 1, w[i]) - pre;
printf("%d", getf(1, n));
return 0;
}
05-18 08:10