hdu-1421-搬寝室 线性dp


搬寝室

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


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





设状态dp[i+1][j][k]表示对于第i个物品,选不选这个与前一个物品的组合, k=0表示不选i, k=1表示选i。
则有, dp[i+1][j][1]=min(dp[i][j][0], dp[i][j][1], dp[i][j-1][0])。


#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cmath>
#include <cstdlib>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>

using namespace std;

#define ll long long
#define rep(x, y, z) for(int x=y;x<z;x++)
const int maxn=2e3+7;
int n, k;
int a[maxn];
ll ans;
int dp[maxn][maxn/2][2];

void init(){
    memset(dp, -1, sizeof(dp));
    ans=0;
}
void solve(){
    sort(a, a+n);
    rep(i, 1, n){
        dp[i+1][1][1]= (a[i]-a[i-1])*(a[i]-a[i-1]);
        if(dp[i][1][1] != -1){
            if(dp[i][1][0] != -1) dp[i+1][1][0]=min(dp[i][1][1], dp[i][1][0]);
            else dp[i+1][1][0]=dp[i][1][1];
        }
        else dp[i+1][1][0]=dp[i][1][0];
        rep(j, 2, k+1){
            
            if(dp[i][j-1][0] != -1) dp[i+1][j][1]=dp[i][j-1][0]+(a[i]-a[i-1])*(a[i]-a[i-1]);
            if(dp[i][j][1] != -1){
                if(dp[i][j][0] != -1) dp[i+1][j][0]=min(dp[i][j][1], dp[i][j][0]);
                else dp[i+1][j][0]=dp[i][j][1];
            }
            else dp[i+1][j][0]=dp[i][j][0];
        }
    }
    if(dp[n][k][0] != -1) {
        if(dp[n][k][1]  != -1) printf("%d\n", min(dp[n][k][0], dp[n][k][1]));
        else printf("%d\n", dp[n][k][0]);
    }
    else printf("%d\n", dp[n][k][1]);
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    
    while(~scanf("%d%d", &n, &k)){
        init();
        rep(i, 0, n) scanf("%d", &a[i]);
        solve();
    }
    return 0;
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值