HDOJ1421 搬寝室 DP

搬寝室

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个物品选其中j对，所有对的重量之差平方和等于其疲劳值，求最小疲劳值

动态规划，先对所有重量进行排序，dp[i][j]二维数组表示前i个选j对用的最小疲劳度

用状态转移方程 dp[i][j]=min(dp[i-2][j-1]+(a[i]-a[i-1])*(a[i]*a[i-1]),dp[i-1][j]) (大神那儿学到的，收获很大）。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
#define N 2000
#define M 40000
int a[N],dp[N][N];
int main()
{
	int n,i,j,k,t,s;
	while(~scanf("%d %d",&n,&k))
	{
		if(2*k>n) continue;
		
		int  MIN=M;
		memset(a,0,sizeof(a));
		for(i=1;i<=n;i++)
		scanf("%d",&a[i]);
		sort(a+1,a+1+n);
		
		memset(dp,M,sizeof(dp));
		dp[0][0]=0;
		dp[1][0]=0;
		for(i=2;i<=n;i++)
		 {
		 	dp[i][0]=0;
		 	for(j=1;j<=k;j++)
		 	{
		 		if(j*2>i) break;
		 		dp[i][j]=min(dp[i-2][j-1]+(a[i]-a[i-1])*(a[i]-a[i-1]),dp[i-1][j]);
			 }
		 }
	printf("%d\n",dp[n][k]);
	}
	return 0;
 }