hduoj 1421搬寝室【动态规划】

搬寝室

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 11 3

 

Sample Output

4

 

【分析】

首先证明：从4个数中 a b c d  依次递增；选取相邻的两个数一定是最小得

及:(a-b)^2+(c-d)^2<(a-c)^2+(b-d)^2

    &&(a-b)^2+(c-d)^2<(a-d)^2+(b-c)^2;

然后：令dp[i][j]表示前i个物品中选j对的最小疲劳度。

dp[i][j]可能含有第i个物品(这种情况下,第i种物品一定是和第i-1个物品配对)，则dp[i][j]=dp[i-2][j-1]+(val[i]-val[i-1])*(val[i]-val[i-1]);
dp[i][j]的j对也可能不含有第i个物品，此时有dp[i][j]=dp[i-1][j];

AC代码：

#include<cstdio>#include<cstring>#include<algorithm>#include<iostream>using namespace std;const int maxn = 20005;const int MAX_INT = 99999999;int a[maxn];int vis[maxn]; int dp[maxn][maxn/2];//dp[i][j]含有i或i-1个物品，配成j对时的最小疲劳度 int main(){int n, k;while(scanf("%d %d", &n, &k) != EOF){int i, j;a[0] = dp[0][0] = 0;for(i = 0; i < n; i++)cin>>a[i];sort(a, a+n);for(i = 0; i <= n; i++)for(j = 0; j <= k; j++)dp[i][j] = MAX_INT;//初始化配对数(j>0)时为较大数 for(i = 0; i <= n; i++)dp[i][0] = 0;//配对为0时,疲劳度为0 for(i = 2; i <= n; i++)for(j = 1; j*2 <= i; j++)dp[i][j] = min(dp[i-2][j-1]+(a[i-1]-a[i-2])*(a[i-1]-a[i-2]), dp[i-1][j]);printf("%d\n", dp[n][k]);}return 0;}