hdu1521-dp

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 11 3 

 

Sample Output

 4 

 

首先排序，使得两个数之间差的平方最小

若f[i-1][j]表示一共i-1件物品，要搬2*j件，若加一件物品a[i]，则有两种情况：

（1）不需要搬这件物品，则f[i][j]=f[i-1][j];

（2）需要搬这件物品，由于是排序好的，则一定会搬最后两件物品，最小疲劳相当于除掉这两件物品的最小疲劳再加上这两件物品的疲劳，即f[i][j]=f[i-2][j-1]+(a[i]-a[i-1])*(a[i]-a[i-1]);

状态方程：f[i][j]=min(f[i-1][j],f[i-2][j-1]+(a[i]-a[i-1])*(a[i]-a[i-1]);

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
 using namespace std;
int dp[2010][1005];
int a[2010];
 int mi(int x,int y)
 {
     if(x<y)
        return x;
     else
        return y;
 }
 int main()
 {
     int n,k;
     while(~scanf("%d%d",&n,&k))
     {
         memset(dp,0,sizeof(dp));
         for(int i=1;i<=n;i++)
            scanf("%d",&a[i]);
         sort(a+1,a+n+1);
         for(int i=2;i<=n;i++)
            for(int j=1;j*2<=i;j++)
         {
             if(i==2*j)
                dp[i][j]=dp[i-2][j-1]+(a[i]-a[i-1])*(a[i]-a[i-1]);
             else
                dp[i][j]=mi(dp[i-1][j],dp[i-2][j-1]+(a[i]-a[i-1])*(a[i]-a[i-1]));
         }
         printf("%d\n",dp[n][k]);
     }
     return 0;
 }