HDU 1421 搬寝室

搬寝室
Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u

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个物品，每个物品有重量，从中选出m对，使得这m对物品重量差的平方和最小。疲劳度：m对物品重量差的平方和。

分析

先对n种物品的重量排序令dp[j][i]表示前j个物品中选i对的最小疲劳度。则dp[j][i]中的i对物品可能含有第j个物品(这种情况下,第j种物品一定是和第j-1个物品配对)，那么dp[j][i] = dp[j-2][i-1]+(a[j-1]-a[j])*(a[j-1]-a[j])，dp[j][i]的i对也可能不含有第j个物品，此时有  dp[j][i]=dp[j-1][i]。

所以，状态转移方程为  dp[j][i]=min{dp[j-2][i-1]+(a[j-1]-a[j])*(a[j-1]-a[j]),dp[j-1][i]}

AC代码如下

#include <cstdio>
#include <cstring>
#include <algorithm>
#define maxn 2000+10
#define min(a,b) a<b?a:b
using namespace std;
int dp[maxn][maxn],a[maxn];
int main()
{
int n,m;
while(~scanf("%d%d",&n,&m))
{
for(int i = 1 ; i <= n ; i++)
scanf("%d",&a[i]);
sort(a + 1, a + n + 1);
memset(dp,0,sizeof(dp));
for(int i = 1 ; i <= m; i++)
{
for(int j = 2 * i ; j <= n ; j++)
{
if(j == 2 * i) dp[j][i] = dp[j-2][i-1] + (a[j-1] - a[j])*(a[j-1] - a[j]);
else
{
dp[j][i] = min(dp[j-1][i],dp[j-2][i-1] + (a[j-1] - a[j])*(a[j-1] - a[j]));
}
}
}
printf("%d\n",dp
[m]);
}
return 0;
}