生活随笔
收集整理的這篇文章主要介紹了
BZOJ.1190.[HNOI2007]梦幻岛宝珠(分层背包DP)
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
題目鏈接
把重量表示為\(a\times2^b\)的形式,然后按\(b\)排序。
從高到低枚舉每一位,\(f[i]\)表示當前位容量為\(i\)時的最大價值(容量即\(a\times2^{bit}\))。對于同一位,直接\(01\)背包就行了。
如何轉移到下一位?\(f[i]\)轉移到\(f[i\times2+\text{W在這一位是否為1}]\)。注意到每一位的容量不會超過\(n\times a_{max}=1000\),所以再對\(1000\)取\(\min\)即可。
//840kb 64ms
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
typedef long long LL;
const int N=1005;struct Node
{int a,b,val;bool operator <(const Node &x)const{return b>x.b;}
}A[N];inline int read()
{int now=0,f=1;register char c=gc();for(;!isdigit(c);c=='-'&&(f=-1),c=gc());for(;isdigit(c);now=now*10+c-'0',c=gc());return now*f;
}
void Solve(int n,int m)
{static int f[N],g[N];int mx=0;for(int i=1,a,b; i<=n; ++i){a=read(), A[i].val=read();for(b=0; b<30&&!(a&1); ++b,a>>=1);mx=std::max(mx,a), A[i].a=a, A[i].b=b;}std::sort(A+1,A+1+n);const int lim=n*mx;memset(f,-0x3f,sizeof f);f[0]=0, A[n+1].b=-1;for(int i=30,now=1; ~i; --i){memset(g,-0x3f,sizeof g);for(int j=0,p; j<=lim; ++j)p=std::min((j<<1)+(m>>i&1),lim), g[p]=std::max(g[p],f[j]);memcpy(f,g,sizeof f);for(; A[now].b==i; ++now)for(int ai=A[now].a,vi=A[now].val,j=0; j<=lim-ai; ++j)f[j]=std::max(f[j],f[j+ai]+vi);//f[i]:剩余容量為i }int ans=0;for(int i=0; i<=lim; ++i) ans=std::max(ans,f[i]);printf("%d\n",ans);
}int main()
{int n,m;while(n=read(),m=read(),n!=-1&&m!=-1) Solve(n,m);return 0;
}
轉載于:https://www.cnblogs.com/SovietPower/p/9805221.html
總結
以上是生活随笔為你收集整理的BZOJ.1190.[HNOI2007]梦幻岛宝珠(分层背包DP)的全部內容,希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。