HDU 1421-搬寝室

搬寝室

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


题目链接: 点击打开链接

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
  
  

 

分析：

我们需要在 n 件物品中选择出 2*k 件物品，但是为了让主人公在搬完这 2*k 件物品后尽可能的有一个好状态，所以我们要降低疲劳度，主人公每趟拿两件物品，疲劳度是和两手物品重量的差值的平方成正比的，所以我们要对 n 件物品的重量从小到大排序，因为两件重量最相邻的物品的差值最小，所以我们又可以再定义一个数组，用来保存两件相邻物品的重量差值，于是我们再开一个dp数组，dp[i][j] 代表在 i 件物品中选择 j 对物品所得的疲劳度，因为下标是从 0 开始的，所以dp[1][1]表示从两件物品中取走一对物品，对此，我们可以推到该动态转移方程为： dp[i][j] = min(dp[i - 1][j], dp[i - 2][j - 1] + b[i-1]);（i 从 2 到 n-1）（j 从 1 到 k），初始化将 dp[i][0] 都设为 0 ，其余的都设为 INF，让dp[1][1] = b[0],具体代码实现如下（在dp数组和b数组之间总是处理不好，最后好不容易出来的）： 

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

int const MAX = 2005;
int const INF = 0x3f3f3f3f;
int dp[MAX][MAX], a[MAX], b[MAX];///这三个数组分析中已经解释过了

int main()
{
    int n, k;
    while(scanf("%d %d", &n, &k) != EOF)
    {
        for(int i = 0; i < n; i++)
            scanf("%d", &a[i]);
        sort(a , a + n );
        for(int i = 0; i < n-1; i++)
            b[i] = (a[i + 1] - a[i]) * (a[i + 1] - a[i]);
        for(int i = 0; i <= n; i++)
        {
            for(int j = 1; j <= k; j++)
                dp[i][j] = INF;
            dp[i][0] = 0;///这里需注意，从i件物品中取0对物品的疲劳度是0
        }
        dp[1][1] = b[0];///因为下标是从0开始，所以这里的i为1，其实代表的是2件物品
        for(int i = 2; i < n; i++)
            for(int j = 1; j <= k; j++)
                    dp[i][j] = min(dp[i - 1][j], dp[i - 2][j - 1] + b[i-1]);
        printf("%d\n", dp[n-1][k]);
    }
}