码迷,mamicode.com
首页 > 其他好文 > 详细

Space Elevator [POJ2392] [DP][优化]

时间:2018-06-23 18:53:13      阅读:217      评论:0      收藏:0      [点我收藏+]

标签:无限   namespace   使用   ace   har   pre   elevator   while   main   

题目大意

n件物品,第i件hi高,有ci件,最高的一件不能超过ai的高度。问最高能堆多高

输入:

第一行,一个n

接下来每一行,为hi,ai,ci

输出,最高堆多高

样例输入:

3
7 40 3
5 23 8
2 52 6

样例输出:

48 (从下到上:3个2号,3个1号,6个3号)

分析:

我们很容易想到要先放限制高度小的,那我们就先按限制高度从小到大排序。

然后我们可以发现,这个和多重背包很像,“物重”“价值”都为hi,数量为ci

设dp[i][j]为用前i件,花费高度j的盒子,最高能堆多高(显然dp[][j]=j/0)。

那么状态转移方程为dp[i][j]=max{dp[i-1][j-k*hi]+hi*k}其中限制条件是(dp[i-1][j-k*hi]!=0或j==k*hi )

但我们发现这个复杂度是不优秀的,我们需要降复杂度。

以下是优化方法{

我们知道,0/1背包和完全背包只需要改变循环方向就行了。

而完全背包从小到大循环就是保证了在更新了标号小的后,计算标号大的时候利用小的就可以使用更新后的值,从而实现“物品无限”。

而我们多重背包是有限的个数,如何在完全背包上动点手脚?

我们在dp[j]相对记录一个usd[j]表示已经用了这个物品几次了。

在保证条件dp[j-hi]!=0或j==hi 和dp[j-hi]+hi>dp[j]时

还需满足usd[j-hi]+1<=ci,再把usd[j]更新成usd[j-hi]+1,这样就能保证该物品使用次数小于等于ci。

最后输出结果小心坑,找到第一个dp[i]不等于0,输出;特判ans==0的情况

代码:

技术分享图片
 1 #include<set>
 2 #include<map>
 3 #include<queue>
 4 #include<stack>
 5 #include<cmath>
 6 #include<cstdio>
 7 #include<cstring>
 8 #include<iostream>
 9 #include<algorithm>
10 #define RG register int
11 #define rep(i,a,b)    for(RG i=a;i<=b;++i)
12 #define per(i,a,b)    for(RG i=a;i>=b;--i)
13 #define ll long long
14 #define inf (1<<29)
15 #define maxn 405
16 #define maxm 40005
17 using namespace std;
18 int n;
19 int dp[maxm],usd[maxm];
20 struct Dat{
21     int w,a,c;
22     inline int operator < (const Dat &tt)const{
23         return a<tt.a;
24     }
25 }p[maxn];
26 inline int read()
27 {
28     int x=0,f=1;char c=getchar();
29     while(c<0||c>9){if(c==-)f=-1;c=getchar();}
30     while(c>=0&&c<=9){x=x*10+c-0;c=getchar();}
31     return x*f;
32 }
33 int main()
34 {
35     n=read();
36     int mx=0;
37     rep(i,1,n)    p[i].w=read(),p[i].a=read(),p[i].c=read();
38     sort(p+1,p+1+n);
39     rep(i,1,n)
40     {
41         memset(usd,0,sizeof(usd));
42         rep(j,p[i].w,p[i].a)
43             if((j==p[i].w||dp[j-p[i].w]) && dp[j-p[i].w]+p[i].w>dp[j] && usd[j-p[i].w]+1<=p[i].c)
44                 dp[j]=dp[j-p[i].w]+p[i].w,
45                 usd[j]=usd[j-p[i].w]+1;
46     }
47     per(i,p[n].a,1)    if(dp[i]){cout<<dp[i];return 0;}
48     puts("0");
49     return 0;
50 }
View Code

 

Space Elevator [POJ2392] [DP][优化]

标签:无限   namespace   使用   ace   har   pre   elevator   while   main   

原文地址:https://www.cnblogs.com/ibilllee/p/9217535.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!