HDU-OJ-1421 搬寝室

搬寝室

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)

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
  
  

 

Author

xhd

 

————————————————————崭新的分割线———————————————————— 思路：分析一下，本题是要规划数对差的平方，使之总和最小。对于前 i 个数，考虑搬走一对物品、两对物品、三对物品…… 数组dp[i][j]表示前 i 个物品中搬走 j 对。 首先要有序。排个序即可。 情况一：2*j == i——这意味着全部都搬走。dp[i][j] = dp[i-2][j-1] + (a[i]-a[i-1])^2 情况二：搬走最后两对或者不搬走最后一个。dp[i][j] = min{dp[i-1][j], dp[i-2][j-1] + (a[i]-a[i-1])^2} 如果不搬走，疲劳值为0。dp[ ][ ]初始化为0。 代码如下：

/****************************************/
 #include <cstdio>
 #include <cstdlib>
 #include <cstring>
 #include <algorithm>
 #include <cmath>
 #include <stack>
 #include <climits>
 #include <queue>
 #include <vector>
 #include <map>
 #include <string>
 #include <iostream>
 using namespace std;
/****************************************/
#define S (a[i] - a[i-1])*(a[i] - a[i-1])
const int N = 2010;
int dp[N][N>>1], a[N];

int main()
{
	int n, k;
	while(~scanf("%d%d", &n, &k)) {
		for(int i = 1; i <= n; i++)
			scanf("%d", &a[i]);
		sort(a+1, a+n+1);
		for(int i = 0; i <= n; i++)
			for(int j = 0; j <= k; j++)
				dp[i][j] = 0;
		for(int i = 2; i <= n; i++)
			for(int j = 1; (j<<1) <= i; j++)
				if((j<<1) == i)
					dp[i][j] = dp[i-2][j-1] + S;
				else
					dp[i][j] = min(dp[i-1][j], dp[i-2][j-1] + S);//状态转移，有点类似背包的姿势。dp[i-1][j]意味着状态为不考虑第i个物体的j状态
		printf("%d\n", dp[n][k]);
	}
	return 0;
}