Problem Description
搬寝室是很累的,xhd深有体会.时间追述2006年7月9号,那天xhd迫于无奈要从27号楼搬到3号楼,因为10号要封楼了.看着寝室里的n件物品,xhd开始发呆,因为n是一个小于2000的整数,实在是太多了,于是xhd决定随便搬2k件过去就行了.但还是会很累,因为2k也不小是一个不大于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
思路
1、先对a[N]从小到大排序;
2、动态方程:dp[i][j]=min(dp[i-1][j-2]+(a[j]-a[j-1])*(a[j]-a[j-1]),dp[i][j-1])。 dp[i][j]代表从i件物品里选取j对,当选取物品数为i(从2开始是因为物品件数最小为2)时,j从1到i/2遍历。
3、注意初始化问题, 最大值的选取问题。
代码
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std;
int a[2001] = {0};
int dp[2001][2001] = {0}; // dp[i][j]代表从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);
// 初始化dp状态值为最大, i=0&j=0时为1; 其他为最大值
int MAX_INT = pow(2, 31)-1;
for (int i=1; i<=n; i++)
for (int j=1; j<=n; j++)
dp[i][j] = MAX_INT;
for(int i=2; i<=n; i++)
{
for(int j=1; j*2<=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]);
// cout<<i<<"-"<<j<<":"<<dp[i][j]<<endl; // 调试用
}
}
printf("%d\n", dp[n][k]);
}
return 0;
}