動態規劃 hdu 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 1
1 3



Sample Output
4

Solution

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
#define min(a,b) (a)<(b)?(a):(b)
#define sqr(a) (a)*(a)
#define MAX 2100
#define MAXK 1100
int n,k,i,j,a[MAX],f[MAX][MAXK];
int main()
{
	while (cin>>n>>k)
	{
		for (i=1;i<=n;i++)
			cin>>a[i];
		sort(a+1,a+n+1);
		f[0][0]=0;
		for (i=1;i<=n;i++)
		{
			f[i][0]=0;
			for (j=1;j<=k&&j*2<=i;j++)
				if (j*2==i)
					f[i][j]=f[i-2][j-1]+sqr(a[i]-a[i-1]);
				else
					f[i][j]=min(f[i-1][j],f[i-2][j-1]+sqr(a[i]-a[i-1]));
		}
		cout<<f[n][k]<<endl;
	}
	return 0;
}