4897: [Thu Summer Camp2016]成绩单

4897: [Thu Summer Camp2016]成绩单

Time Limit: 40 Sec Memory Limit: 512 MB
Submit: 89 Solved: 51
[Submit][Status][Discuss]
Description

期末考试结束了，班主任L老师要将成绩单分发到每位同学手中。L老师共有n份成绩单，按照编号从1到n的顺序叠
放在桌子上，其中编号为i的成绩单分数为w_i。成绩单是按照批次发放的。发放成绩单时，L老师会从当前的一叠
成绩单中抽取连续的一段，让这些同学来领取自己的成绩单。当这批同学领取完毕后，L老师再从剩余的成绩单中
抽取连续的一段，供下一批同学领取。经过若干批次的领取后，成绩单将被全部发放到同学手中。然而，分发成绩
单是一件令人头痛的事情，一方面要照顾同学们的心理情绪，不能让分数相差太远的同学在同一批领取成绩单；另
一方面要考虑时间成本，尽量减少领取成绩单的批次数。对于一个分发成绩单的方案，我们定义其代价为：

其中，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+(1^2+0^2+1^2+1^2)×1=15。

HINT

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

Source

鸣谢Sengxian上传

[Submit][Status][Discuss]


式子里那个$k$可以当做开玩笑的，因为每一段的代价显然是互不影响的
定义状态$f[i][j][l][r]$表示对于区间$[i,j]$，取到剩余数字值域为区间$[l,r]$的最小代价
特别地，定义$f[i][j][0][0]$为取完的最小代价
对于一段区间$[i,j]$，最后和左端点一起删除的一定是某段从左起的序列
最后和右端点一起删除的一定是某段从右起的序列
那么就可以枚举一个分界点$k$来转移
于是$f[i][j][l][r] = min\{f[i][k][l][r] + f[k+1][j][l][r],f[i][k][0][0] + f[k+1][j][l][r],f[i][k][l][r] + f[k+1][j][0][0]\}$
显然，$f[i][j][0][0] = min\{f[i][j][l][r] + a + b * (r - l)^2\}$
暴力记搜就行了，复杂度远远不足$O(n^5)$

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int N = 53;
const int INF = 100000000;

int n,a,b,cur = 1,w[N],A[N],f[N][N][N][N];
bool vis[N][N][N][N];

#define min(a,b) ((a) < (b) ? (a) : (b))
#define max(a,b) ((a) > (b) ? (a) : (b))
#define update(x,y) ((x) = (x) > (y) ? (y) : (x))

inline int Dfs(int i,int j,int l,int r)
{
    if (vis[i][j][l][r]) return f[i][j][l][r];
    f[i][j][l][r] = INF; vis[i][j][l][r] = 1;
    if (!l && !r)
    {
        int Min = cur + 1,Max = 0;
        for (int x = i; x <= j; x++)
            Min = min(Min,w[x]),Max = max(Max,w[x]);
        f[i][j][l][r] = a + b * (A[Max] - A[Min]) * (A[Max] - A[Min]);
        for (int x = 1; x <= cur; x++)
            for (int y = x; y <= cur; y++)
            {
                int ret = Dfs(i,j,x,y); if (ret == INF) continue;
                update(f[i][j][l][r],ret + a + b * (A[y] - A[x]) * (A[y] - A[x]));
            }
        return f[i][j][l][r];
    }
    else
    {
        if (i == j)
        {
            if (l <= w[i] && w[i] <= r)
                f[i][j][l][r] = 0;
            return f[i][j][l][r];
        }
        for (int k = i; k < j; k++)
        {
            int L,R;
            L = Dfs(i,k,0,0);
            R = Dfs(k + 1,j,l,r);
            update(f[i][j][l][r],L + R);
            L = Dfs(i,k,l,r);
            R = Dfs(k + 1,j,0,0);
            update(f[i][j][l][r],L + R);
            L = Dfs(i,k,l,r);
            R = Dfs(k + 1,j,l,r);
            update(f[i][j][l][r],L + R);
        }
        return f[i][j][l][r];
    }
}

int main()
{
    #ifdef DMC
        freopen("DMC.txt","r",stdin);
    #endif

    cin >> n >> a >> b;
    for (int i = 1; i <= n; i++)
        scanf("%d",&w[i]),A[i] = w[i];
    sort(A + 1,A + n + 1);
    for (int i = 2; i <= n; i++)
        if (A[i] != A[i - 1]) A[++cur] = A[i];
    for (int i = 1; i <= n; i++)
        w[i] = lower_bound(A + 1,A + cur + 1,w[i]) - A;
    cout << Dfs(1,n,0,0) << endl;
    return 0;
}