【NOIP2006】 金明的预算方案

分类：非树形有依赖的背包问题

难度：提高+/省选-

---

题面

题目描述
金明今天很开心，家里购置的新房就要领钥匙了，新房里有一间金明自己专用的很宽敞的房间。更让他高兴的是，妈妈昨天对他说：“你的房间需要购买哪些物品，怎么布置，你说了算，只要不超过NN元钱就行”。今天一早，金明就开始做预算了，他把想买的物品分为两类：主件与附件，附件是从属于某个主件的，下表就是一些主件与附件的例子：
主件 附件
电脑 打印机，扫描仪
书柜 图书
书桌 台灯，文具
工作椅 无
如果要买归类为附件的物品，必须先买该附件所属的主件。每个主件可以有00个、1个或2个附件。附件不再有从属于自己的附件。金明想买的东西很多，肯定会超过妈妈限定的N元。于是，他把每件物品规定了一个重要度，分为5等：用整数1-5表示，第5等最重要。他还从因特网上查到了每件物品的价格（都是10元的整数倍）。他希望在不超过N元（可以等于N元）的前提下，使每件物品的价格与重要度的乘积的总和最大。
设第j件物品的价格为v[j]，重要度为w[j]，共选中了k件物品，编号依次为j1,j2,…,jk，则所求的总和为：v[j1] * w[j1] + v[j2] * w[j2] + …+ v[jk] * w[jk]。
请你帮助金明设计一个满足要求的购物单。

输入输出格式
输入格式：
第1行，为两个正整数，用一个空格隔开：
N m （其中N(<32000)表示总钱数，m(<60)为希望购买物品的个数。） 从第2行到第m+1行，第j行给出了编号为j-1的物品的基本数据，每行有3个非负整数vpq （其中v表示该物品的价格（v<10000），p表示该物品的重要度（1-5），q表示该物品是主件还是附件。如果q=0，表示该物品为主件，如果q>0，表示该物品为附件，q是所属主件的编号）

输出格式：
一个正整数，为不超过总钱数的物品的价格与重要度乘积的总和的最大值（<200000）。

输入输出样例
输入样例#1：
1000 5
800 2 0
400 5 1
300 5 1
400 3 0
500 2 0

输出样例#1：
2200

题解

关于分组背包问题和有依赖的背包问题详解参见背包九讲

分组背包问题伪代码：

for 所有的组k
for v=V…0
for 所有的i属于组k
f[v]=max{f[v],f[v-c[i]]+w[i]}

有依赖的背包问题：

这个问题由NOIP2006金明的预算方案一题扩展而来。遵从该题的提法，将不依赖于别的物品的物品称为“主件”，依赖于某主件的物品称为“附件”。由这个问题的简化条件可知所有的物品由若干主件和依赖于每个主件的一个附件集合组成。

按照背包问题的一般思路，仅考虑一个主件和它的附件集合。可是，可用的策略非常多，包括：一个也不选，仅选择主件，选择主件后再选择一个附件，选择 主件后再选择两个附件……无法用状态转移方程来表示如此多的策略。（事实上，设有n个附件，则策略有2^n+1个，为指数级。）

考虑到所有这些策略都是互斥的（也就是说，你只能选择一种策略），所以一个主件和它的附件集合实际上对应于P06中的一个物品组，每个选择了主件又选择了若干个附件的策略对应于这个物品组中的一个物品，其费用和价值都是这个策略中的物品的值的和。但仅仅是这一步转化并不能给出一个好的算法，因为物品组中的物品还是像原问题的策略一样多。

再考虑P06中的一句话： 可以对每组中的物品应用P02中“一个简单有效的优化”。 这提示我们，对于一个物品组中的物品，所有费用相同的物品只留一个价值最大的，不影响结果。所以，我们可以对主件i的“附件集合”先进行一次01背包，得到费用依次为0…V-c[i]所有这些值时相应的最大价值f’[0…V-c[i]]。那么这个主件及它的附件集合相当于V-c[i]+1个物品的物品 组，其中费用为c[i]+k的物品的价值为f’[k]+w[i]。也就是说原来指数级的策略中有很多策略都是冗余的，通过一次01背包后，将主件i转化为 V-c[i]+1个物品的物品组，就可以直接应用P06的算法解决问题了。

代码（C++）

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
struct cas {
    int v,p,q;
}a[60],pat[60][60];
int n,m,t[60],V[60][10],P[60][10],cnt[60],f[32000],ans;

int main() {
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++) {
        scanf("%d%d%d",&a[i].v,&a[i].p,&a[i].q);
        if(a[i].q) {
            t[a[i].q]++;
            pat[a[i].q][t[a[i].q]].v = a[i].v;
            pat[a[i].q][t[a[i].q]].p = a[i].p;
            pat[a[i].q][t[a[i].q]].q = a[i].q;
        }
    }
    for(int i=1;i<=m;i++) {
        if(t[i]) {
            memset(f,-1,sizeof(f));
            f[0] = 0;
            for(int j=1;j<=t[i];j++) {
                for(int k=n-a[i].v;k>=pat[i][j].v;k--) {
                    if(f[k]<f[k-pat[i][j].v]+pat[i][j].v*pat[i][j].p && f[k-pat[i][j].v]!=-1) {
                        f[k] = f[k - pat[i][j].v] + pat[i][j].v * pat[i][j].p;
                    }
                }
            }
            for(int j=0;j<=n-a[i].v;j++) {
                if(f[j]!=-1) {
                    cnt[i]++;
                    V[i][cnt[i]] = j+a[i].v;
                    P[i][cnt[i]] = f[j] + a[i].v * a[i].p;
                }
            }
        }
        if(!a[i].q) {
            cnt[i]++;
            V[i][cnt[i]] = a[i].v;
            P[i][cnt[i]] = a[i].v * a[i].p;
        }
    }
    memset(f,0,sizeof(f));
    for(int i=1;i<=m;i++) {
        for(int j=n;j>=0;j--) {
            for(int k=1;k<=cnt[i];k++) {
                if(j>=V[i][k]) {
                    f[j] = max(f[j],f[j-V[i][k]]+P[i][k]);
                }
            }
        }
    }
    for(int i=0;i<=n;i++) {
        ans = max(ans,f[i]);
    }
    printf("%d",ans);

    return 0;
}