[动态规划]金明的预算方案

题目传送门:https://www.luogu.org/problemnew/show/P1064

---

题目描述

金明今天很开心，家里购置的新房就要领钥匙了，新房里有一间金明自己专用的很宽敞的房间。更让他高兴的是，妈妈昨天对他说：“你的房间需要购买哪些物品，怎么布置，你说了算，只要不超过N元钱就行”。今天一早，金明就开始做预算了，他把想买的物品分为两类：主件与附件，附件是从属于某个主件的，下表就是一些主件与附件的例子：

主件 附件

电脑 打印机，扫描仪

书柜 图书

书桌 台灯，文具

工作椅 无

如果要买归类为附件的物品，必须先买该附件所属的主件。每个主件可以有0个、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个非负整数

v p q （其中v表示该物品的价格（v<10000），p表示该物品的重要度（1~5），q表示该物品是主件还是附件。如果q=0，表示该物品为主件，如果q>0，表示该物品为附件，q是所属主件的编号）

输出格式：

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

---

题解

NOIP提高组的一道题，仔细想一想难度不算太大。
一道经典的背包问题，在01背包的基础上增加了一些难度。
首先，题目描述中说所有价格都是10的整数倍，所以在读入数据的时候我们就可以直接除以十，输出答案时再乘回来即可（防止超时超空间）。
状态设定：f[i]表示重量还剩i时所能取得的最大价值。
每个物品的价格就是重量，其重要度imp[i]乘以价格pri[i]就是价值。
对于每一件主件Mi会有不超过2件的附件，所以我们枚举每个附件是否购买即可。
那么就会有四个转移方程:

dp[j]=max(dp[j-w[mitem[i]]]+v[mitem[i]],dp[j]);//只取主件
dp[j]=max(dp[j-w[mitem[i]]-w[item[mitem[i]][1]]]+v[mitem[i]]+v[item[mitem[i]][1]],dp[j]);//取一号附件
dp[j]=max(dp[j-w[mitem[i]]-w[item[mitem[i]][2]]]+v[mitem[i]]+v[item[mitem[i]][2]],dp[j]);//取二号附件
dp[j]=max(dp[j-w[mitem[i]]-w[item[mitem[i]][1]]-w[item[mitem[i]][2]]]+v[mitem[i]]+v[item[mitem[i]][1]]+v[item[mitem[i]][2]],dp[j]);//两个附件都取

mitem[i]表示第i个主件的编号，item[mitem[i]][j]表示第i个主件的第j个附件的编号

注意判断剩余重量是否足够以及这个主件到底有没有1,2号附件即可。

---

代码:

#include<cstdio>
#include<iostream>
#include<ctime>
#include<string>
#include<algorithm>
#include<cmath>
#define maxn 100000
#define maxm 100000
#define inf 99999999
using namespace std;
int item[maxn][5];
int mitem[maxn];
int dp[maxn];
int w[maxn],v[maxn];
int n,m,mcnt,maxx;
int main()
{
    scanf("%d%d",&n,&m);
    n/=10;
    for(int i=1;i<=m;i++)
    {
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        if(c==0)
        {
            mitem[++mcnt]=i;
            w[i]=a/10;
            v[i]=a*b/10;
        }
        else
        {
            item[c][++item[c][0]]=i;
            w[i]=a/10;
            v[i]=a*b/10;
        }
    }
    for(int i=1;i<=mcnt;i++)
    {
        for(int j=n;j>=0;j--)
        {
            if(j-w[mitem[i]]>=0) 
                dp[j]=max(dp[j-w[mitem[i]]]+v[mitem[i]],dp[j]);
            if(item[mitem[i]][1]&&j-w[mitem[i]]-w[item[mitem[i]][1]]>=0) 
                dp[j]=max(dp[j-w[mitem[i]]-w[item[mitem[i]][1]]]+v[mitem[i]]+v[item[mitem[i]][1]],dp[j]);
            if(item[mitem[i]][2]&&j-w[mitem[i]]-w[item[mitem[i]][2]]>=0) 
                dp[j]=max(dp[j-w[mitem[i]]-w[item[mitem[i]][2]]]+v[mitem[i]]+v[item[mitem[i]][2]],dp[j]);
            if(item[mitem[i]][1]&&item[mitem[i]][2]&&j-w[mitem[i]]-w[item[mitem[i]][1]]-w[item[mitem[i]][2]]>=0) 
                dp[j]=max(dp[j-w[mitem[i]]-w[item[mitem[i]][1]]-w[item[mitem[i]][2]]]+v[mitem[i]]+v[item[mitem[i]][1]]+v[item[mitem[i]][2]],dp[j]);
        }
    }
    for(int i=0;i<=n;i++) maxx=max(maxx,dp[i]);
    printf("%d",maxx*10);
    return 0;
}