HDOJ 1421 搬寝室（基础DP方程）

搬寝室

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

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
  
  
  
  

   
   汉语题面，题意就不用解释了
  
  
  
  

   
   题解：
  
  
  
  

   
   看题之前要先明白一个小技巧：
  
  
  
  

   
   有4个数 a b c d  依次递增；
选取相邻的两个数的消耗一定是最小的
即:
  
  
  
  

   
   (a-b)^2+(c-d)^2<(a-c)^2+(b-d)^2；
  
  
  
  

   
   (a-b)^2+(c-d)^2<(a-d)^2+(b-c)^2;

  
  
  
  

   
   所以先给输入的数据升序排个序吧。。。
  
  
  
  

   
   知道了这些然后来看这道题，很明显的一道dp的基础题，dp的题难点就是动态转移方程的建立，现在来看这个题，这个题里要求答案要知道的就是从哪些里面取，取多少...所以可以开一个二维的dp[i][j]数组，表示从前j个数中选2*i个数的最小消耗是dp[i][j]，所以我们可以得到动态转移方程为：dp[i][j]=min(dp[i][j-1],dp[i-1][j-2]+cal(a[j-1],a[j]));这里的cal函数计算的是拿a[j-1]和a[j]两个物品的消耗。
  
  
  
  

   
   这里当dp[i][j]中i*2=j时是特殊情况可以先扫一遍先初始化这些就好了即：dp[i/2][i]=cal(a[i-1],a[i])+dp[(i-2)/2][i-2];
  
  
  
  

   
   这样代码就不难得出来了。
  
  
  
  

   
   题目链接：
   
   http://acm.hdu.edu.cn/showproblem.php?pid=1421
  
  
  
  

   
   代码：
  
  
  
  

   
   
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn=2020;
int dp[maxn/2][maxn];//dp[][]数组

int cal(int x,int y){//计算机拿x和y时的消耗
    return (y-x)*(y-x);
}

int main(){
    int n,k;
    while(scanf("%d%d",&n,&k)==2){
        int a[maxn];//用来存输入的物品
        for(int i=1;i<=n;i++) scanf("%d",&a[i]);
        sort(a+1,a+n+1);//升序排序
        memset(dp,0,sizeof(dp));//初始化dp[][]数组为零
        for(int i=1;i<=n;i++){
            if(i%2==0) dp[i/2][i]=cal(a[i-1],a[i])+dp[(i-2)/2][i-2];//初始化dp[i][j]中i*2=j时的dp[][]数组中的值
        }
        for(int i=1;i<=k;i++){
            for(int j=2;j<=n;j++){
                if(j==2*i) continue;
                dp[i][j]=min(dp[i][j-1],dp[i-1][j-2]+cal(a[j-1],a[j]));//利用前面推出的方程计算出dp[][]数组中的各个值
            }
        }
        printf("%d\n",dp[k][n]);//打印出从n个里面选2*k个的最小消耗
    }
    return 0;
}
大家要是还有什么不明白的欢迎提出来