rq06金明的预算方案题解

古老的背包题了，动态规划的气息浓郁。

在著名的背包九讲里面有详细的介绍，我先搬运工一下。

7有依赖的背包问题

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

 7.2算法 这个问题由NOIP2006中“金明的预算方案”一题扩展而来。遵从该题的提法，将不依赖于别的物品的物品称为“主件”，依赖于某主件的物品称为“附件”。由这个问题的简化条件可知所有的物品由若干主件和依赖于每个主件的一个附件集合组成。 按照背包问题的一般思路，仅考虑一个主件和它的附件集合。可是，可用的策略非常多，包括：一个也不选，仅选择主件，选择主件后再选择一个附件，选择主件后再选择两个附件……无法用状态转移方程来表示如此多的策略。事实上，设有n个附件，则策略有2n+1个，为指数级。 考虑到所有这些策略都是互斥的（也就是说，你只能选择一种策略），所以一个主件和它的附件集合实际上对应于6中的一个物品组，每个选择了主件又选择了若干个附件的策略对应于这个物品组中的一个物品，其费用和价值都是这个策略中的物品的值的和。但仅仅是这一步转化并不能给出一个好的算法，因为物品组中的物品还是像原问题的策略一样多。 再考虑对每组内的物品应用2.3中的优化。我们可以想到，对于第k个物品组中的物品，所有费用相同的物品只留一个价值最大的，不影响结果。所以，可以对主件k的“附件集合”先进行一次01背包，得到费用依次为0...V−Ck所有这些值时相应的最大价值Fk[0...V−Ck]。那么，这个主件及它的附件集合相当于V−Ck+1个物品的物品组，其中费用为v的物品的价值为Fk[v−Ck]+Wk，v的取值范围是Ck≤v≤V。 也就是说，原来指数级的策略中，有很多策略都是冗余的，通过一次01背包后，将主件k及其附件转化为V−Ck+1个物品的物品组，就可以直接应用6的算法解决问题了。 

7.3较一般的问题 更一般的问题是：依赖关系以图论中“森林”3的形式给出。也就是说，主件的附件仍然可以具有自己的附件集合。限制只是每个物品最多只依赖于一个物品（只有一个主件）且不出现循环依赖。解决这个问题仍然可以用将每个主件及其附件集合转化为物品组的方式。唯一不同的是，由于附件可能还有附件，就不能将每个附件都看作一个一般的01背包中的物品了。若这个附件也有附件集合，则它必定要被先转化为物品组，然后用分组的背包问题解出主件及其附件集合所对应的附件组中各个费用的附件所对应的价值。 事实上，这是一种树形动态规划，其特点是，在用动态规划求每个父节点的属性之前，需要对它的各个儿子的属性进行一次动态规划式的求值。这已经触及到了“泛化物品”的思想。看完8后，你会发现这个“依赖关系树”每一个子树都等价于一件泛化物品，求某节点为根的子树对应的泛化物品相当于求其所有儿子的对应的泛化物品之和。 7.4小结 NOIP2006的那道背包问题我做得很失败，写了上百行的代码，却一分未得。后来我通过思考发现通过引入“物品组”和“依赖”的概念可以加深对这题的理解，还可以解决它的推广问题。用物品组的思想考虑那题中极其特殊的依赖关系：物品不能既作主件又作附件，每个主件最多有两个附件，可以发现一个主件和它的两个附件等价于一个由四个物品组成的物品组，这便揭示了问题的某种本质。 后来，我在《背包问题九讲》第一版中总结此事时说：“失败不是什么丢人的事情，从失败中全无收获才是。”之后的NOIP2007的比赛中，我得了满分。

上面是纯搬运的，如果是初接触动态规划，是很有必要把背包九讲看一遍。

这道题目简单的来说就是把附件考虑到主件里面来，循环的时候只考虑主件，而主件存在不拿，拿主件，拿主件加附件1，拿主件加附件2，拿主件加附件1和2一共五种情况。如果之前能把输入的数据处理好了，状态转移的部分几乎和01背包一样。

#include<stdio.h>
#include<stdlib.h>
struct hey
{   int num;
    int w[6];
    int v[6];
    int flag;
       }a[67];
int max(int a,int b,int c)
{
    if(a>=b&&a>=c)return(a);
    if(b>=a&&b>=c)return(b);
    if(c>=a&&c>=b)return(c);
    }
main()
{   
    int i,m,n,q[67],v[67],p[67];
    int j=1,M,N,k,f[67][3207]={0};
    scanf("%d %d",&N,&M);
    N=N/10;
    for(i=1;i<=M;i++)
          scanf("%d %d %d",&v[i],&p[i],&q[i]);
    for(i=1;i<=M;i++)
       {
        v[i]=v[i]/10;
        if(q[i]==0)
           {
             a[j].num=i;
             a[j].w[1]=v[i]*p[i];
             a[j].v[1]=v[i];   
             a[j].flag=1;
             j++;
            } 
        }
    m=j-1;
    for(i=1;i<=M;i++)
        if(q[i]!=0)
        {
           j=1;
           while(a[j].num!=q[i])j++;
           if(j>m)continue;
           a[j].flag++;
           a[j].w[a[j].flag]=a[j].w[1]+v[i]*p[i];
           a[j].v[a[j].flag]=a[j].v[1]+v[i];
           if(a[j].flag==3)
           {
              a[j].w[4]=a[j].w[2]+a[j].w[3]-a[j].w[1];
              a[j].v[4]=a[j].v[2]+a[j].v[3]-a[j].v[1];
              a[j].flag++;
           }
        }
    for(i=1;i<=m;i++)
       for(k=1;k<=a[i].flag;k++)
         for(j=N;j>=0;j--)
           if(j-a[i].v[k]>=0)
                f[i][j]=max(f[i-1][j-a[i].v[k]]+a[i].w[k],f[i-1][j],f[i][j]);
           else f[i][j]=max(0,f[i][j],f[i-1][j]);
     printf("%d\n",f[m][N]*10);    
               
      }