九度OJ 1452 搬寝室 -- 动态规划

题目地址：http://ac.jobdu.com/problem.php?pid=1452

题目描述：

搬寝室是很累的,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件物品后的最佳状态是怎样的(也就是最低的疲劳度),请告诉他吧。

输入：

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

输出：

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

样例输入：

2 1
1 3

样例输出：

4

//在选定的最优方案中，每对物品都是重量相邻的一对物品
//dp[i][j]表示从前i件物品取j对物品的最优方案
//状态方程：dp[i][j]=min(dp[i-1][j],dp[i-2][j-1]+(goods[i]-goods[i-1])^2）
#include <stdio.h>
#include <stdlib.h>

#define MAXN 2001

int compare (const void * p, const void * q){
	return *(int *)p - *(int *)q;
}

int main(void){
	int n, k, m, i, j;
	int goods[MAXN];
	int dp[3][MAXN];

	while (scanf ("%d%d", &n, &k) != EOF){
		for (i=1; i<=n; ++i)
			scanf ("%d", &goods[i]);
		qsort (goods+1, n, sizeof(int), compare);
		for (i=0; i<3; ++i)
			for (j=0; j<MAXN; ++j)
				dp[i][j] = 0;
		for (i=2; i<=n; ++i){
			for (j=1; j<=k && j*2<=i; ++j){
				m = (goods[i] - goods[i-1]) * (goods[i] - goods[i-1]);
				if (j * 2 == i)
					dp[i%3][j] = dp[(i-2)%3][j-1] + m;
				else
					dp[i%3][j] = ((dp[(i-2)%3][j-1] + m) < dp[(i-1)%3][j]) ? (dp[(i-2)%3][j-1] + m) : dp[(i-1)%3][j];
			}
		}
		printf ("%d\n", dp[n%3][k]);
	}

	return 0;
}

HDOJ上相同的题目：http://acm.hdu.edu.cn/showproblem.php?pid=1421

HDOJ上相似的题目：http://acm.hdu.edu.cn/showproblem.php?pid=1500，CODE代码片

参考资料：http://zh.wikipedia.org/wiki/%E5%8A%A8%E6%80%81%E8%A7%84%E5%88%92

http://hawstein.com/posts/dp-novice-to-advanced.html