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
题意:计算搬n件物品中的k对的最小疲劳值(注意每搬一对,疲劳值在原来的基础上累加)
题解:动态规划
a[i]用来存储每一件物品的重量;
w[i]代表搬第i-1与第i件的疲劳值;
dp[i][j]代表前i件物品搬j对后的最小疲劳值。
接下来,对于第i件物品而言,寻找状态方程
1.第i件物品不搬,即在前i- 1件物品中搬j对,那么疲劳值为 dp[i - 1][j]
2.第i件物品搬,则一定是第i件物品和第i-1件物品一起搬(这样取差后平方值最小),即在前i - 2件物品中搬j - 1对物品,再搬最后一对物品,那么疲劳值为 dp[i - 2][j - 1] + w[i]
(w[i]代表搬第i-1与第i件的疲劳值)
dp[i][j] = min(dp[i - 1][j], dp[i - 2][j - 1] + w[i])
应该注意的是要考虑边界问题:
1.当2 * j > i时,即要搬的数量超过了物品总量,故此时令dp[i][j]为无穷大;
2.当j = 0时,即一对物品都没搬时,所需疲劳值应该是0,此时令dp[i][j] = 0。
#include<stdio.h>
#include<algorithm>
#define INF 0x3f3f3f3f
#define NUM 2010
using namespace std;
int dp[NUM][NUM/2];
int a[NUM];
int w[NUM];
int get(int i, int j)
{
if(j * 2 > i)
return INF;
if(j == 0)
return 0;
return dp[i][j]; //正常情况
}
int main()
{
int n, k;
while(scanf("%d%d", &n, &k) != EOF)
{
for(int i = 1; i <= n; ++i)
{
scanf("%d", &a[i]);
}
sort(a + 1, a + n + 1); //把a数组中1-n的元素排序,这样连续求出的疲劳值最小
for(int i = 2; i <=n; ++i)
{
w[i] =(a[i -1] - a[i])*(a[i - 1] - a[i]);
}
for(int i=2;i<=n;i++)
{
for(int j=1;2*j<=i;j++)
{
dp[i][j] = min(get(i - 1, j), get(i - 2, j - 1) + w[i]);
}
}
printf("%d\n", dp[n][k]);
}
return 0;
}