DP训练——背包DP

背包DP

POJ3093
题意
给定V(V<=30)个物品和背包容量D(D<=1000),求将背包尽量放满的方案数。(尽量放满指的是剩余物品中重量最小的物品都不能再放入背包)
题解
考虑到背包要放的尽量满,所以额外枚举一下没有放入背包的最小物品即可。
具体实现上,可以先将重量升序排序,然后枚举第k个物品最终没有放入背包。
则第1到第k-1的物品必然放入背包,同时对第k+1到第V的物品做01背包后,将满足条件的方案数统计下即可。
状态定义:d[j]表示去掉不放入背包的最小物品k后,背包装入总重量为j的方案数
目标:\sum_{j=max(maxD-v[k]+1,0)}^{maxD}d[j](maxD为前k-1个物品已放入背包后,背包的最大容量)
边界:d[0]=1
转移方程:d[j]=d[j]+d[j-v[i]],k+1\leq i\leq n
PS:本题存在如下坑点。
1 一个物品都放不进背包的情况被视为方案数为0,因此需要在开始位置特判。
2 所有物品全放入背包都可以的情况被视为方案数为1,同样需要特判。
3 当枚举到第k个物品放入背包后,背包恰好放满的情况下,当前的答案就是1,不再进行dp。
代码如下

/*

*/
#define method_1
#ifdef method_1
/*

*/
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iomanip>
#include<ctime>
#include<string>
#include<bitset>
#define D(x) cout<<#x<<" = "<<x<<"  "
#define E cout<<endl
using namespace std;
typedef long long ll;
typedef pair<int,int>pii;
const int maxv=30+5;
const int maxd=1000+5;
const int INF=0x3f3f3f3f;
int T,kase=0;
int V,D;
int v[maxv];
int sum[maxv];
ll d[maxd];
ll ans;
void init(){
    ans=0ll;
}
void dp(){
    for(int k=1;k<=V;k++){ //枚举剩余的最小物品k
        int maxD=D-sum[k-1]; //则前k-1个物品已经放入背包,背包的最大容量变为D-sum[k-1]
        if(maxD<0) break;
        if(maxD==0){
            ans++;
            continue;
        }
        memset(d,0,sizeof(d));
        d[0]=1;
        for(int i=k+1;i<=V;i++) for(int j=maxD;j>=v[i];j--) d[j]+=d[j-v[i]];
        for(int i=max(0,maxD-v[k]+1);i<=maxD;i++) ans+=d[i];
        /*
        for(int i=max(0,maxD-v[k]+1);i<=maxD;i++){
            D(k);D(maxD);D(i);D(d[i]);E;
        }
        */
    }
    if(sum[V]<=D) ans++;
}
int main() {
    ios::sync_with_stdio(false);
    //freopen("POJ3093.in","r",stdin);
    cin>>T;
    while(T--){
        cout<<++kase<<" ";
        cin>>V>>D;
        init();
        for(int i=1;i<=V;i++) cin>>v[i];
        sort(v+1,v+V+1);
        if(v[1]>D){
            cout<<0<<endl;
            continue;
        }
//        reverse(v+1,v+V+1);
        for(int i=1;i<=V;i++) sum[i]=sum[i-1]+v[i];
        dp();
        cout<<ans<<endl;
    }
    return 0;
}
#endif
#ifdef method_2
/*

*/

#endif
#ifdef method_3
/*

*/

#endif

BZOJ2748
题意
给你初始数的大小beginlevel和所能达到最大数的大小maxlevel(maxlevel<=1000)
再给你n(n<=50)个数让你每次加上这个数或减去这个数,但不能超出区间[0,maxlevel],求最后得到的最大结果。
题解
因为每一层都有两个转移方向,同时又存在一层内没有任何一个合法的状态转移的情况。
因此,01背包时不能压缩掉第一个维度。(因为需要确定第n层是否能转移到)
其他部分就是01背包的模板。
状态定义:d[i,j]表示处理完第i个物品,背包装入总重量为j的存在性
目标:max(j),d[n,j]=1
边界:d[0][beginlevel]=1
转移方程:d[i,j]=\left\{\begin{matrix}d[i,j]|d[i-1,j-v[i]],0\leq j-v[i] & & \\ d[i,j]|d[i-1,j+v[i]],j+v[i]\leq maxlevel & & \end{matrix}\right.
代码如下

/*

*/
#define method_1
#ifdef method_1
/*

*/
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iomanip>
#include<ctime>
#include<string>
#include<bitset>
#define D(x) cout<<#x<<" = "<<x<<"  "
#define E cout<<endl
using namespace std;
typedef long long ll;
typedef pair<int,int>pii;
const int maxn=50+5;
const int maxl=1000+5;
const int INF=0x3f3f3f3f;
int n,bl,ml,c[maxn];
int d[maxn][maxl];
void dp(){
    memset(d,0,sizeof(d));
    d[0][bl]=1;
    for(int i=1;i<=n;i++) for(int j=0;j<=ml;j++){
        if(j+c[i]<=ml) d[i][j]|=d[i-1][j+c[i]];
        if(j-c[i]>=0) d[i][j]|=d[i-1][j-c[i]];
    }
    for(int j=ml;j>=0;j--) if(d[n][j]){
        cout<<j;
        return;
    }
    cout<<-1;
}
int main() {
    ios::sync_with_stdio(false);
//  freopen("2748.in","r",stdin);
    cin>>n>>bl>>ml;
    for(int i=1;i<=n;i++) cin>>c[i];
    dp();
    return 0;
}
#endif
#ifdef method_2
/*

*/

#endif
#ifdef method_3
/*

*/

#endif

BZOJ1190
题意
n(n<=100)个物品,选取若干个,使得重量不超过W(W<=2^{32})的情况下,总价值最大。
数据保证每个物品的重量符合a*2^b(a<=10,b<=30)
题解
这是一道按位划分后,从低位向高位做分层背包的题目。
首先,考虑利用题目数据的条件,即每个物品的质量都可以表示为a*2^b的形式。
这意味着,如果按照b的值,对所有物品进行分组后(b=i的物品属于第i组),每组内所有a的值的和最大为1000(a<=10,n<=100)
我们将组号设为dp数组的第一维(范围[0,30]),组内所有a的值设为dp数组的第二维(范围[0,1000])。
于是,dp数组的空间便可以承受。
状态定义:d[i,k]表示在只使用第i组内的物品,背包已经占用的重量为k的最大价值
目标:所有d值,用于下一步计算
边界:d[0,0]=0
转移方程:d[i,k]=max(d[i,k-G[i][j]]+V[i,j])(j循环枚举第i组内的所有物品,向量G保存第i组内的物品的a值)
接下来考虑如何用上述dp数组,将所有组合并起来计算答案。
状态定义:d[i,j]表示在前i组中,背包已经占用的重量为j*2^i加上W二进制第i位以下的体积的最大价值(可近似理解为对上述状态求前缀和)
目标:d[len,1](len表示2^{len}<=W的最大的len,这就意味着需要考虑所有在转移方程中考虑W几位上本来就有的1,即对应转移方程中的(W>>i-1)\&1和递推w前缀和数组时的向上取整)
边界:d[0,0]=0
转移方程:d[i,j]=max(d[i,j-k]+d[i-1,min(w[i-1],k*2+((W>>i-1)\& 1))])(k循环枚举在i-1组内选择的物品体积,则j-k为在第i组内选择的物品体积)
其他细节详见代码注释。
代码如下

/*

*/
#define method_1
#ifdef method_1
/*

*/
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iomanip>
#include<ctime>
#include<string>
#include<bitset>
#define D(x) cout<<#x<<" = "<<x<<"  "
#define E cout<<endl
using namespace std;
typedef long long ll;
typedef pair<int,int>pii;
const int maxn=100+5;
const int maxa=10+5;
const int maxb=30+5;
const int maxm=1000+5;
const int INF=0x3f3f3f3f;
int n,W;
vector<int>G[maxb],V[maxb];
int w[maxb];
int d[maxb][maxm];
int len;
void init(){
    len=0;
    memset(d,0,sizeof(d));
    memset(w,0,sizeof(w));
    memset(G,0,sizeof(G));
    memset(V,0,sizeof(V));
}
int main() {
//  ios::sync_with_stdio(false);
    //freopen("BZOJ1190.in","r",stdin);
    while(scanf("%d%d",&n,&W)){
        if(n==-1&&W==-1) break;
        init();
        int weight,value;
        for(int i=1;i<=n;i++){
            scanf("%d%d",&weight,&value);
            int a,b=0;
            while(!(weight&1)) weight>>=1,b++;
            a=weight;
            len=max(len,b);
            G[b].push_back(a);
            w[b]+=a;
            V[b].push_back(value);
            //按照b的值分组,在每一组内记录对应的a作为重量,value作为价值
        }  
        //首先对每一组进行01背包
        //d[i,k]表示在第i组内,背包已经占用的重量为k的最大价值
        //这时len是组数
        for(int i=0;i<=len;i++){
            for(int j=0;j<G[i].size();j++){
                for(int k=w[i];k>=G[i][j];k--){
                    d[i][k]=max(d[i][k],d[i][k-G[i][j]]+V[i][j]);
                }
            }
        }
        while(W>>len) len++;
        len--;
        //这时len是表示2^len<=W的最大的len
//        D(len);E;
        //d[i,j]表示在第0~i组中,背包已经占用的重量为$j*2^_i$加上W二进制第i位以下的体积的最大价值
        for(int i=1;i<=len;i++){
            w[i]+=(w[i-1]+1)/2; //w[i]=w[i]+ceil(w[i-1]/2)
            for(int j=w[i];j>=0;j--){ //w[i]表示前i组背包的最大加权重量
                for(int k=0;k<=j;k++){
                    d[i][j]=max(d[i][j],d[i][j-k]+d[i-1][min(w[i-1],k*2+((W>>i-1)&1))]); //k*2是因为进了一位
                }
            }
        }
        //按照d的状态定义,答案就是d[len,1]
        printf("%d\n",d[len][1]);
    }
    return 0;
}
#endif
#ifdef method_2
/*

*/

#endif
#ifdef method_3
/*

*/

#endif
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 216,402评论 6 499
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,377评论 3 392
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 162,483评论 0 353
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,165评论 1 292
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,176评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,146评论 1 297
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,032评论 3 417
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,896评论 0 274
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,311评论 1 310
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,536评论 2 332
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,696评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,413评论 5 343
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,008评论 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,659评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,815评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,698评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,592评论 2 353

推荐阅读更多精彩内容