HDU 1241

搬寝室

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Submission(s): 9598    Accepted Submission(s): 3235

Problem Description

搬寝室是很累的,xhd深有体会.时间追述2006年7月9号,那天xhd迫于无奈要从27号楼搬到3号楼,因为10号要封楼了.看着寝室里的n件物品,xhd开始发呆,因为n是一个小于2000的整数,实在是太多了,于是xhd决定随便搬2*k件过去就行了.但还是会很累,因为2*k也不小是一个不大于n的整数.幸运的是xhd根据多年的搬东西的经验发现每搬一次的疲劳度是和左右手的物品的重量差的平方成正比(这里补充一句,xhd每次搬两件东西,左手一件右手一件).例如xhd左手拿重量为3的物品,右手拿重量为6的物品,则他搬完这次的疲劳度为(6-3)^2 = 9.现在可怜的xhd希望知道搬完这2*k件物品后的最佳状态是怎样的(也就是最低的疲劳度),请告诉他吧.

 

Input

每组输入数据有两行,第一行有两个数n,k(2<=2*k<=n<2000).第二行有n个整数分别表示n件物品的重量(重量是一个小于2^15的正整数).

 

Output

对应每组输入数据,输出数据只有一个表示他的最少的疲劳度,每个一行.

 

Sample Input

2 1 1 3

 

Sample Output

4

 

 

解题：

定义数组a[i]为搬第i对物品所消耗的疲劳值；数组f[n][k]来表示在n件物品中搬k对的最佳状态，而达到这一状态的决策可能为：

1. 第n件物品不搬，即在前n - 1件物品中搬k对，那么疲劳值仍为f[n - 1][k]；
2. 第n件物品要搬，那么根据上面所证，第n - 1件物品也要同时搬，即在前n - 2件物品中搬k - 1对物品，再搬最后一对物品，那么疲劳值为f[n - 2][k - 1] + a[n - 1]，n - 1是因为对数必然比总物品数少1。

为使疲劳值最小，因此最佳策略为取两种决策中的最小值，即应使：

f[n][k] = min(f[n - 1][k], f[n - 2][k - 1] + a[n - 1])

#include<stdio.h>
#include<stdlib.h> 

int cmp ( const void *a , const void *b )
{
    return *(int *)a - *(int *)b;
}

int n, k;
int i, j;
int a[2005];
int f[2005][1003];

int getMin(int a, int b)
{
    return (a<b?a:b);
}

    
int getCr(int n, int k)
{
    if(2*k > n)
    {
        return 1099999999;            /*因为是取最小的消耗值，所以设一个很大的数*/
    }
    if(k == 0)
    {
        return 0;
    }
    else
    {
        return f[n][k];
    } 
}

int main()
{        
    while(scanf("%d%d", &n, &k) != EOF)
    {
        for(i=1; i<=n; i++)
        {
            scanf("%d", &a[i]);
        }
        qsort(a, n+1, sizeof(int), cmp);    /*对差值的平方进行排序 */ 
        for(i=1; i<n; i++)                /*计算相邻两个数差的平方 */ 
        {
            a[i] = a[i+1] - a[i];
            a[i] *= a[i];
        }
        
        /* 计算最小消耗值*/ 
        for(i=1; i<=k; i++)    
        {
            for(j=2*i; j<=n; j++)
            {
                f[j][i] = getMin(getCr(j-1, i), getCr(j-2, i-1) + a[j-1]);  
            }
        }
        printf("%d\n", f[n][k]);    
    }
    return 0;
}