【DP】zoj3164

这个题是背包的一个大综合题，大概出了有依赖关系的背包以外蕴含了所有背包类型，然后这里归一下类几种背包问题的解法。

01背包：暴力解法：f[i][j]代表考虑了前i个物品装了j的体积 转移就用f[i-1][j-v[i]]+w[i]；

               优化空间：1，f[j]代表考虑了j的体积 转移f[j-v[i]]+w[i]]，注意枚举体积时要倒序枚举，这样省空间的原理是f[j]                                          的状态更新好之后，在体积小于j的状态中并不会用到j来更新，所以不会出现重复装入。

                                 2，利用滚动数组的思想f[1][j]=f[0][j-v[i]]+w[i]，然后再把f[1]的值赋给f[0]，这样做的话即使体积是正                                        向枚举的也不会出现重复装入，然而这样做的话，细节要注意的是初始化的时候，f[1]和f[0]的                                          初始化的值要一样。

多重背包：暴力解法：枚举每次装入的数量，然后按照01背包来做；

                  优化时间：枚举数量的时候利用二进制拆分的思想，然后要注意的是拆分的方式，假设数量为k，假设k的二                                      进制位为1......，有n位，那么要拆分成2^(n-2)  (也就是1<<(n-2)) 也就是拆分后最大位应该是n-1                                      位，然后把2^(n-3),2^(n-4).....2^0逐次按照01背包放入，最后放入k-((1<<n-1)-1)

完全背包：更新的时候正向枚举体积即可

分组背包：不能一起加入的物品归为一个组，这类物品先进行一次背包，然后我的做法是额外开一个数组记录这组中                   加入某件物品后达到体积j的价值的最大值，也就是每次更新时用没有考虑过这一组的值进行更新，然后更                   新后数组单存起来，然后在这些单存起来的数组中找每个j体积对应的价值最大值。然后再把这个最大值数                   组作为考虑到此组后的背包的状态，这样就避免了组内选东西的重复，然而组内每一类物品无论是什么类                   型的背包都可以照常进行；

目测背包九讲里只讲了每个物品只有一个的分组背包，并没有讲泛化背包加分组背包的做法。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<vector>
using namespace std;
int kind[2010],w[2010],v[2010],k[2010],dp[3][2010],n,d,tail=0,g,ans;
vector < int > q[2010];
void limit_work(int x,int z)
{
     int i,j;
     for (i=1;i*2<k[x];i=i*2)
         for (j=d;j>=v[x]*i;j--)
             dp[z][j]=max(dp[z][j],dp[z][j-v[x]*i]+w[x]*i);
     if (k[x]-i>0) i=k[x]-(i-1);
     for (j=d;j>=v[x]*i;j--)
         dp[z][j]=max(dp[z][j],dp[z][j-v[x]*i]+w[x]*i);
     return ;
}
void unlimit_work(int x,int z)
{
    for (int i=v[x];i<=d;i++)
        dp[z][i]=max(dp[z][i],dp[z][i-v[x]]+w[x]);
}
int main()
{
    int x;
    char ch;
    while (scanf("%d %d",&n,&d)!=EOF)
    {
        memset(dp,-127,sizeof(dp));
        memset(kind,0,sizeof(kind));
        memset(k,0,sizeof(k));
        memset(w,0,sizeof(w));
        memset(v,0,sizeof(v));
        for (int i=1;i<=tail;i++)
            q[i].clear();
        tail=0;
        dp[0][0]=0;
        dp[2][0]=0;
        ans=0;
        for (int i=1;i<=n;i++)
            scanf("%d%d%d",&k[i],&w[i],&v[i]);
        scanf("%d",&g);
        for (tail=1;tail<=g;)
        {
            scanf("%d",&x);
            q[tail].push_back(x);
            kind[x]=tail;
            ch=getchar();
            if (ch!=' ') tail++;
        }
        tail--;
        for (int i=1;i<=n;i++)
            if (kind[i]==0)
            {
                if (k[i]==0) unlimit_work(i,0);
                else limit_work(i,0);
            }
        for (int i=1;i<=tail&&g>0;i++)
        {
            memset(dp[1],-127,sizeof(dp[1]));
            for (int j=0;j<q[i].size();j++)
            {
                for (int jj=0;jj<=d;jj++)
                   dp[2][jj]=dp[0][jj];
                if (k[q[i][j]]==0) 
                   unlimit_work(q[i][j],2);
                else 
                   limit_work(q[i][j],2);
                for (int j=0;j<=d;j++)
                   dp[1][j]=max(dp[1][j],dp[2][j]);
            }
            for (int j=0;j<=d;j++)
                dp[0][j]=dp[1][j];
        }
        if (dp[0][d]>=0) printf("%d\n",dp[0][d]);
        else printf("i'm sorry...\n");
    }
    return 0;
}

话说这是博客上第一道不是poj的题。。。。明天看final，让我先奶起来：tourist，你赢定了