hdu 1421 dp

搬寝室

Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 18191 Accepted Submission(s): 6170

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

 1 /*
 2 dp(i,j)表示：i个数内,选出j对数的最优策略
 3 dp(i,j)=max(dp(i-1,j),dp(i-2,j-1)+(a[i]-a[i-1])^2)
 4 */
 5 #include <iostream>
 6 #include <cstdio>
 7 #include <cstring>
 8 using namespace std;
 9 
10 typedef __int64 LL;
11 const LL INF=1e15;
12 const int maxn=2005;
13 int a[maxn];
14 LL dp[maxn][maxn];
15 inline LL min(LL a,LL b){return a<b?a:b;}
16 
17 void quicksort(int l,int r)
18 {
19     if(l<r)
20     {
21         int t=a[l],i=l,j=r;
22         while(i!=j)
23         {
24             while(a[j]>=t && i<j) j--;
25             while(a[i]<=t && i<j) i++;
26             if(i<j) swap(a[i],a[j]);
27         }
28         a[l]=a[i];a[i]=t;
29         quicksort(l,i-1);
30         quicksort(i+1,r);
31     }
32 }
33 
34 int main()
35 {
36     int n,k,i,j;
37     while(~scanf("%d%d",&n,&k))
38     {
39         for(i=1;i<=n;i++) scanf("%d",a+i);
40         quicksort(1,n);
41         for(i=0;i<=n;i++)
42         for(j=0;j<=k;j++)
43             dp[i][j]=INF;
44         for(i=0;i<=n;i++) dp[i][0]=0;
45         for(i=2;i<=n;i++)
46         {
47             for(j=1;j<=k && j*2<=i;j++)
48                 dp[i][j]=min(dp[i-1][j],dp[i-2][j-1]+(a[i]-a[i-1])*(a[i]-a[i-1]));
49         }
50         printf("%I64d\n",dp[n][k]);
51     }
52     return 0;
53 }