完全背包问题 解题报告

Posted ppprseter

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了完全背包问题 解题报告相关的知识,希望对你有一定的参考价值。

完全背包问题

(n)种物品,物品的体积分别为(V_1,V_2,dots,V_n),且每种物品的数量都可以看做是无限多的。现在有(m)次询问,每次询问给定一个容量为取的背包,请你回答是否存在一种物品选择方案,使得背包恰好能被完全装满(仅考虑体积,忽略长、宽、高等其他因素)。同时,要求所有选出的物品中,体积不小于(L)的物品总数量不能超过(C)件。

输入格式

第一行为两个正整数(n)(m),分别表示物品的种数以及询问的次数。
第二行为(n)个正整数(V_1,V_2,dots,V_n(V_i le 10000)),分别表示这(n)种物品的体积。
第三行为两个非负整数 (L(L le 20000 ))(C(Cle 30)),表示在选择方案中对大体积物品的数量限制。
接下来(m)行,每行一个正整数(W_i),表示这次询问中背包的容量。

输出格式

输出共(m)行,每行一个字符串,表示对应询问的答案。
对于每次询问,如果存在一种合法的方案,请输出(Yes),否则输出(No)

数据规模与约定

对于(10\%)的数据:(nle 8,W_ile 100)

对于(30\%)的数据:(W_ile 10000)

对于(60\%)的数据:(nle 30,mle 200)

对于另外(10\%)的数据:(n=2)

对于(100\%)的数据:(nle 50,mle 100000,W_ile 10^{18})


考虑暴力

(dp_{i,j,k})表前(i)个物品总体积为(k)选择了(j)件大物品的是否合法。

复杂度(O(nmax W_i c))

考虑修补这个想法。设(V_0)为最小的体积

  • 如果(V_0ge L),因为可选的物品总体积不会太大,所以考虑直接暴力,复杂度(O(nc^2sum V_i)),用( t{bitset})

优化一下就可以过了。

  • 如果(V_0 < L),因为(V_0)可以无限选,所以我们不妨求出在同余于(V_0)情况下的可能,这样就简化了状态。

(dp_{i,j,k})表示前(i)个物品选择了(j)件大物品选择的总体积同余于(V_0)下为(k)的最小总体积。

保证了最小总体积,我们就可以判断(W_i)是否大于等于(dp_{n,forall jle c,W_i mod V_0})来看是否合法

考虑转移

(dp_{i,j,k}=min(dp_{i,j,(v_0-v_i)mod v_0}+v_i,dp_{i-1,j,k}),v_i < L)

(dp_{i,j,k}=min(dp_{i,j-1,(v_0-v_i)mod v_0}+v_i,dp_{i-1,j,k}),v_i ge L)

发现第一个转移有环,可以建图解决。

源点(S)连接所有(0 sim v_0-1)的点,边权为(dp_{i,j,v_0}),然后其他点每个点(v)((v+v_i)mod v_0),边权为(v_i),然后跑一边最短路就可以了。

发现这个图比较特殊,可以通过寻找性质(O(n))解决。


Code:

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <bitset>
#define ll long long
int n,c,m,v[52],L;ll w;
int min(int x,int y){return x<y?x:y;}
namespace work1
{
std::bitset <1000001> dp[31];
void work()
{
    dp[0][0]=1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=c;j++)
            dp[j]|=dp[j-1]<<v[i];
    for(int i=1;i<=m;i++)
    {
        scanf("%lld",w);
        int ans=0;
        for(int j=0;j<=n;j++) ans|=dp[j][w];
        if(ans) puts("Yes");
        else puts("No");
    }
}
}
namespace work2
{
const int N=10010;
const int inf=0x3f3f3f3f;
int dp[31][N],used[32][N],mi,id;
#define dec(a,b) (((a-b)%v[1]+v[1])%v[1])
#define to ((now+v[p])%v[1])
void dfs0(int p,int k,int now)
{
    if(used[k][now]) return;
    used[k][now]=1;
    if(mi>dp[k][now]) mi=dp[k][now],id=now;
    dfs0(p,k,to);
}
void dfs1(int d,int p,int k,int now,int anc)
{
    if(now==anc) return;
    dp[k][now]=min(dp[k][now],d+v[p]);
    dfs1(dp[k][now],p,k,to,anc);
}
void topo(int p)
{
    memset(used,0,sizeof(used));
    for(int i=0;i<=c;i++)
        for(int j=0;j<v[1];j++)
            if(!used[i][j])
            {
                mi=inf;
                dfs0(p,i,j);
                dfs1(dp[i][id],p,i,(id+v[p])%v[1],id);
            }
}
void work()
{
    memset(dp,0x3f,sizeof(dp));
    dp[0][0]=0;
    for(int i=2;i<=n;i++)
    {
        if(v[i]<L) topo(i);
        else
        {
            for(int k=0;k<=c;k++)
                for(int j=0;j<v[1];j++)
                    dp[k][j]=min(k?dp[k-1][dec(j,v[i])]+v[i]:inf,dp[k][j]);
        }
    }
    for(int i=1;i<=m;i++)
    {
        scanf("%lld",&w);
        int ans=0;
        for(int j=0;j<=c;j++)
            if(dp[j][w%v[1]]!=inf)
                ans|=dp[j][w%v[1]]<=w;
        if(ans) puts("Yes");
        else puts("No");
    }
}
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",v+i);
    scanf("%d%d",&L,&c);
    std::sort(v+1,v+1+n);
    if(v[1]>=L) work1::work();
    else work2::work();
    return 0;
}

2018.10.30

以上是关于完全背包问题 解题报告的主要内容,如果未能解决你的问题,请参考以下文章

解题报告:hdu2191汶川地震 - 多重背包模板

动态规划之01背包详解解题报告

济南集训20191001解题报告

19.05.09 解题报告

题解报告:NYOJ #311完全背包(恰好装满)

动态规划-第二节:动态规划之背包类型问题