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
状态空间为 f ( i , j ) 表示从 i 件物品中取出 j 件的最小疲劳度
状态转移方程 f ( i , j ) = min ( f [ i - 1 ] [ j ] , f [ i - 2 ] [ j - 1 ] + a [ i ] )
f [ i - 1 ] [ j ] 表示第 i 件不取
f [ i - 2 ] [ j - 1 ] + a [ i ] 表示第 i 件取 (a [ i ] 表示的是第 i 件和第 i - 1 件的疲劳度)
2.
如果取 4 件 重量为 a , b , c , d 的物品( a < b < c < d)
那么 ab 和 cd 配对是疲劳度最小的选择
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#define Min(a,b) ((a)<=(b)?(a):(b))
#define ms(a) memset(a,0,sizeof(a))
using namespace std;
#define N 2005
int f[N][N];//状态空间
int w[N];//存各件物品的重量
int a[N];//表示的是第 i 件和第 i - 1 件的疲劳度
int main(){
int n,k;
while(cin>>n>>k){
ms(f);
ms(w);
ms(a);
for(int i=1;i<=n;i++){
cin>>w[i];
}
sort(w+1,w+n+1);//先排序
for(int i=2;i<=n;i++){
a[i]=(int)pow((w[i]-w[i-1]),2);
}
for(int i=1; i<=n; i++)
for(int j=i/2+1; j<=n; j++)
f[i][j] = INT_MAX;
for(int i=1;i<=n;i++) //初始条件
f[i][0]=0;
f[1][1]=INT_MAX;
f[2][2]=a[2];
for(int i=2;i<=n;i++){
for(int j=1;2*j<=i;j++){
f[i][j]=Min(f[i-1][j],f[i-2][j-1]+a[i]);
}
}
cout<<f[n][k]<<endl;
}
return 0;
}