【Algothrim】 动态规划实例 (金明的预算方案 )

Description

金明今天很开心，家里购置的新房就要领钥匙了，新房里有一间金明自己专用的很宽敞的房间。更让他高兴的是，妈妈昨天对他说：“你的房间需要购买哪些物品，怎么布置，你说了算，只要不超过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]。（其中*为乘号）请你帮助金明设计一个满足要求的购物单。

Input 

输入文件的第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是所属主件的编号）
Output

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

Sample Input
1000 5
800 2 0
400 5 1
300 5 1
400 3 0
500 2 0

Analysis

有依赖的背包问题
简化的问题
这种背包问题的物品间存在某种“依赖”的关系。也就是说，i依赖于j，表示若选物品i，则必须选物品j。为了简化起见，我们先设没有某个物品既依赖于别的物品，又被别的物品所依赖；另外，没有某件物品同时依赖多件物品。
算法
这个问题由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的算法解决问题了。

这个题多了一些限制条件，所以每个物品的决策不止两个了。但由于题目里有这句话“每个主件可以有0个、1个或2个附件”，

所以可以这样做：

如果第i号物品为附件，就直接54了它。

如果第i号物品为主件，那么就有这么几种决策：不放入它；仅放入它；仅放入它和它的第一个附件；仅放入它和它的第二个附件；放入它和它的所有附件。

#define _CRT_SECURE_NO_WARNINGS#include<stdio.h>#include<string.h>int  max(int x, int y){if (x > y)return x;elsereturn y;}int N, M;int V[65], P[65], Q[65];int F[32005];int main(){memset(V,0,sizeof(V));memset(P, 0, sizeof(P));memset(Q, 0, sizeof(Q));memset(F, 0, sizeof(F));scanf("%d %d", &N, &M);int i, j,k;int k1, k2;int t1, t2;int l = 0;for (i = 1; i <= M; i++){scanf("%d %d %d", &V[i],&P[i],&Q[i]);}for (i = 1; i <= M; i++){k1 = 0;k2 = 0;t1 = 0;t2 = 0;if (Q[i] == 0){for (k = i + 1; k <= M; k++){if (Q[k] == i){t1 = k;k1 = 1;break;}}for (k = t1 + 1; k <= M;k++){if (Q[k] == i){t2 = k;k2 = 1;break;}}for (j = N; j >= V[i]; j--){F[j] = max(F[j - V[i]] + V[i] * P[i], F[j]);//放主if (j-V[i]-V[t1]>=0 && k1==1)//放付1F[j] = max(F[j - V[i]-V[t1]] + V[i] * P[i]+V[t1]*P[t1], F[j]);if (j - V[i] - V[t2] >= 0 && k2 == 1)//放付2F[j] = max(F[j - V[i] - V[t2]] + V[i] * P[i] + V[t2] * P[t2], F[j]);if (j - V[i] - V[t1] - V[t2] >= 0 && k1 == 1 && k2 == 1)//都放F[j] = max(F[j - V[i] - V[t2]-V[t1]] + V[i] * P[i] + V[t1] * P[t1] + V[t2] * P[t2], F[j]);}}}printf("%d\n", F[N]);return 0;}