2021-07-17

多重背包问题

有 N 种物品和一个容量是 V 的背包。

第 i 种物品最多有 si 件，每件体积是 vi，价值是 wi。

求解将哪些物品装入背包，可使物品体积总和不超过背包容量，且价值总和最大。
输出最大价值。

输入格式
第一行两个整数，N，V，用空格隔开，分别表示物品种数和背包容积。

接下来有 N 行，每行三个整数 vi,wi,si，用空格隔开，分别表示第 i 种物品的体积、价值和数量。

输出格式
输出一个整数，表示最大价值。

数据范围
0<N≤1000
0<V≤2000
0<vi,wi,si≤2000

与完全背包问题的区别无非多了一个每类物品的最大上限，从这一点入手可以很简单的得到

朴素解法如下：

1.#include<bits/stdc++.h>

using namespace std;
const int N=1005;
int f[N][N],w[N],v[N],s[N];
int main()
{
    int n,V;
    cin>>n>>V;
    for(int i=1;i<=n;++i)
        cin>>v[i]>>w[i]>>s[i];
    for(int i=1;i<=n;++i)
        for(int j=v[i];j<=V;++j)
        {
            for(int k=0;k<=s[i]&&k*v[i]<=j;++k)
                f[i][j]=max(f[i][j],f[i-1][j-k*v[i]]+k*w[i]);
        }
    cout<<f[n][V];
}
2.#include<bits/stdc++.h>

using namespace std;
const int N=1005;
int f[N][N],w[N],v[N],s[N];
int main()
{
    int n,V;
    cin>>n>>V;
    for(int i=1;i<=n;++i)
        cin>>v[i]>>w[i]>>s[i];
    for(int i=1;i<=n;++i)
        for(int j=v[i];j<=V;++j)
        {
            f[i][j]=f[i-1][j];
            for(int k=1;k<=s[i]&&k*v[i]<=j;++k)
                f[i][j]=max(f[i][j],f[i-1][j-k*v[i]]+k*w[i]);
        }
    cout<<f[n][V];
}

两个解法几乎一样，但循环的思路有所区别，第一个是第二个的一般化把，k=0就是f[i-1][j];因此2是可行的
但这两个方法的时间复杂度是大O的（nv^2）计算次数也超过了

那么来试一下能不能用前面的优化方法优化这一个
观察f[i][j-v[i]]=max(f[i-1][j-v[i],f[i-1][j-v[i]*2]+w[i]+f[i-1][j-v[i]3]+2w[i]+…+f[i-1][j-s[i]*v[i]]+(n-1)w[i]+s[i+1]v[i]+nw[i])
而f[i][j]=max(f[i-1][j],f[i-1][j-v[i]]+w[i],f[i-1][j-v[i]2]+2w[i]+f[i-1][j-v[i]3]+3w[i]+…+f[i-1][j-nv[i]]+(n)*w[i];
区别在于这里f[i][j-v[i]]比完全背包中的多了一项，且 没有办法确定这一项是不是最终决定了f[i][j-v[i]]的属性
完全背包中后面有几项由这个j来独立决定而这里s【i】也参与了决定，完全背包中必有f[i][j]比f[i][j-v[i]]多一项f[i-1][j]也即是他们
起始不同但结尾同，就是j不够了
而在这里如果j相对于v[i]足够大会使s[i]成为决定性因素，此时不再有上述性质，反而这里f[i][j],f[i][j-v[i]]都是有s【i】项

所以前面的优化行不通

下面采用常用的二进制优化方法将s次循环改变为logs次
数据2000时间复杂度降低100倍大概是1e7可行

所谓二进制优化就是，起始只需要算出logs个状态就可以表示所有s个状态了，从二进制可以表示所有整数可以简单看出
那么试着把s打包乘几组，然后组合起来一定能够达到不漏！！
这就变成了01背包问题
把N开成1000*log2000=11000多点

ac:
#include<bits/stdc++.h>

using namespace std;
const int N=15000;

int w[N],v[N],f[N],cnt=1;

int main()
{
    int n,V;
    cin>>n>>V;
    for(int i=1;i<=n;++i)
    {
        int tv,tw,ts;
        cin>>tv>>tw>>ts;
        for(int k=1;2*k<=ts;k*=2){
            v[cnt]=k*tv;
            w[cnt++]=k*tw;
            ts -=k ;
        }
        if(ts){
            v[cnt]=ts*tv,w[cnt++]=ts*tw;
        }
    }
    for(int i=1;i<=cnt;++i)
        for(int j=V;j>=v[i];--j)
            f[j]=max(f[j],f[j-v[i]]+w[i]);
    cout<<f[V];
}