HDU1421 搬寝室,经典dp

1.题目叙述:

搬寝室

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


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
 

Author
xhd
 

Source
 

Recommend
lcy   |   We have carefully selected several similar problems for you:   1058  1159  1069  2571  1978 

3.思路分析:

首先,要怎么搬呢?即每一对要怎么取?如果有abcd四个数,且a<b<c<d,应该是取ab,cd好呢还是ac,bd好?抑或是bc,ad好呢?答案是第一种,因为:

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

即每对物品都应是重量最为接近的物品,也就是说对n件物品排序后,每对物品都应该是连续的。

定义数组w[i]为搬第i对物品所消耗的疲劳值;数组dp[n][k]来表示在n件物品中搬k对的最佳状态,而达到这一状态的决策可能为:

  1. 第n件物品不搬,即在前n - 1件物品中搬k对,那么疲劳值仍为dp[n - 1][k];
  2. 第n件物品要搬,那么根据上面所证,第n - 1件物品也要同时搬,即在前n - 2件物品中搬k - 1对物品,再搬最后一对物品,那么疲劳值为dp[n - 2][k - 1] + w[n - 1],n - 1是因为对数必然比总物品数少1。

为使疲劳值最小,因此最佳策略为取两种决策中的最小值,即应使:

dp[n][k] = min(dp[n - 1][k], dp[n - 2][k - 1] + w[n - 1])

4.AC代码:

#include 
   
   
    
    
#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
         #include 
        
          #include 
         
           #include 
          
            #include 
           
             #include 
             #include 
             
               #include 
              
                #define INF 0x3f3f3f3f #define maxn 2222 #define eps 1e-6 #define pi acos(-1.0) #define e 2.718281828459 #define mod (int)1e9 + 7; using namespace std; typedef long long ll; int w[maxn], dp[maxn][maxn]; int main() { #ifndef ONLINE_JUDGE freopen("in.txt", "r", stdin); freopen("out.txt", "w", stdout); long _begin_time = clock(); #endif int n, k; while (scanf("%d%d", &n, &k) != EOF) { for (int i = 1; i <= n; i++) scanf("%d", &w[i]); sort(w + 1, w + n + 1); for (int i = 1; i < n; i++) { w[i] = w[i + 1] - w[i]; w[i] *= w[i]; } memset(dp, 0, sizeof(dp)); for (int j = 1; j <= k; j++) for (int i = 2 * j; i <= n; i++) { if (i == j * 2) dp[i][j] = dp[i - 2][j - 1] + w[i - 1]; else dp[i][j] = min(dp[i - 1][j], dp[i - 2][j - 1] + w[i - 1]); } printf("%d\n", dp[n][k]); } #ifndef ONLINE_JUDGE long _end_time = clock(); printf("time = %ld ms\n", _end_time - _begin_time); #endif return 0; } 
               
              
            
           
          
         
       
      
      
     
     
    
    
   
   
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值