HDU1421 搬寝室，经典dp

1.题目叙述：

搬寝室

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

Total Submission(s): 27076    Accepted Submission(s): 9346


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

 

Source

ACM暑期集训队练习赛（二）

 

Recommend

lcy   |   We have carefully selected several similar problems for you:  1058 1159 1069 2571 1978 

3.思路分析：

首先，要怎么搬呢？即每一对要怎么取？如果有abcd四个数，且a<b<c<d，应该是取ab，cd好呢还是ac，bd好？抑或是bc，ad好呢？答案是第一种，因为：

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

即每对物品都应是重量最为接近的物品，也就是说对n件物品排序后，每对物品都应该是连续的。

定义数组w[i]为搬第i对物品所消耗的疲劳值；数组dp
[k]来表示在n件物品中搬k对的最佳状态，而达到这一状态的决策可能为：

第n件物品不搬，即在前n - 1件物品中搬k对，那么疲劳值仍为dp[n - 1][k]；
第n件物品要搬，那么根据上面所证，第n - 1件物品也要同时搬，即在前n - 2件物品中搬k - 1对物品，再搬最后一对物品，那么疲劳值为dp[n - 2][k - 1] + w[n - 1]，n - 1是因为对数必然比总物品数少1。
为使疲劳值最小，因此最佳策略为取两种决策中的最小值，即应使：

dp[n][k] = min(dp[n - 1][k], dp[n - 2][k - 1] + w[n - 1])

4.AC代码：

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define INF 0x3f3f3f3f
#define maxn 2222
#define eps 1e-6
#define pi acos(-1.0)
#define e 2.718281828459
#define mod (int)1e9 + 7;
using namespace std;
typedef long long ll;
int w[maxn], dp[maxn][maxn];
int main()
{
#ifndef ONLINE_JUDGE
freopen("in.txt", "r", stdin);
freopen("out.txt", "w", stdout);
long _begin_time = clock();
#endif
int n, k;
while (scanf("%d%d", &n, &k) != EOF)
{
for (int i = 1; i <= n; i++)
scanf("%d", &w[i]);
sort(w + 1, w + n + 1);
for (int i = 1; i < n; i++)
{
w[i] = w[i + 1] - w[i];
w[i] *= w[i];
}
memset(dp, 0, sizeof(dp));
for (int j = 1; j <= k; j++)
for (int i = 2 * j; i <= n; i++)
{
if (i == j * 2)
dp[i][j] = dp[i - 2][j - 1] + w[i - 1];
else
dp[i][j] = min(dp[i - 1][j], dp[i - 2][j - 1] + w[i - 1]);
}
printf("%d\n", dp
[k]);
}
#ifndef ONLINE_JUDGE
long _end_time = clock();
printf("time = %ld ms\n", _end_time - _begin_time);
#endif
return 0;
}