洛谷P1450 [HAOI2008]硬币购物【背包+容斥】

最新推荐文章于 2022-02-08 15:21:50 发布

niiick

最新推荐文章于 2022-02-08 15:21:50 发布

阅读量227

点赞数

分类专栏：动态规划--背包

本文链接：https://blog.csdn.net/niiick/article/details/83662000

版权

动态规划--背包专栏收录该内容

6 篇文章 1 订阅

订阅专栏

本文探讨了一种特殊的背包问题——硬币购物问题，通过使用完全背包预处理并结合容斥原理来解决限制条件下的支付方案数量计算。文章详细解释了如何在给定硬币面额、携带数量和购买价值的情况下，利用容斥原理消除非法支付方案，最终得出有效的支付方式数目。

摘要由CSDN通过智能技术生成

Time Limit: 10 Sec Memory Limit: 162 MB

Description

硬币购物一共有4种硬币。面值分别为c1,c2,c3,c4。某人去商店买东西，去了tot次。每次带di枚ci硬币，买si的价值的东西。请问每次有多少种付款方法。

Input

第一行 c1,c2,c3,c4,tot 下面tot行 d1,d2,d3,d4,s,其中di,s<=100000,tot<=1000

Output

每次的方法数

题目分析

第一眼单调队列优化多重背包，但再一看数据范围emmmm

首先假如没有物品个数的限制，那么这就是一个完全背包，这样的情况可以直接预处理
现在考虑加入限制，当前 $dp[S_i]$ 中包含的超过 $d_i$ 限制的方案数为 $dp[S_i-c_i*(d_i+1)]$
相当于先强制往包中放入 $d_i+1$ 个 $c_i$ ，剩余 $S_i-c_i*(d_i+1)$ 的容积不管怎么放都已经是不合法的了

先令 $ans=dp[S_i]$ ，然后 $ans-=\sum_{i=1}^4dp[S_i-c_i(d_i+1)]$
但是这时候我们又发现 $c_i,c_j(i!=j)$ 同时不合法的方案数多减了，我们再加回来
这时叒发现 $c_i,c_j,c_k(i!=j!=k)$ 同时不合法的方案数多加了，我们再减掉
以及叕加回四个同时不合法的方案数

所以就是个容斥问题
一般容斥都可以转化成二进制简化运算
复杂度约为 $O(4*maxn+T*2^4)$

#include<iostream>
#include<vector>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
typedef long long lt;

int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

const int maxn=100010;
int T,n;
int ci[5],di[5];
lt dp[maxn],S;

int main()
{
    for(int i=1;i<=4;++i) ci[i]=read();
    
    dp[0]=1;
    for(int i=1;i<=4;++i)
    for(int j=ci[i];j<maxn;++j)
    dp[j]+=dp[j-ci[i]];
    
    T=read();
    while(T--)
    {
        for(int i=1;i<=4;++i) di[i]=read();
        S=read(); lt res=dp[S];
        for(int i=1;i<=(1<<4)-1;++i)
        {
            int tt=S,cnt=0;
            for(int j=1;j<=4;++j)
            if(i&(1<<j-1)) cnt++,tt-=ci[j]*(di[j]+1); 
            if(tt<0) continue;
            if(cnt&1) res-=dp[tt];
            else res+=dp[tt];
        }
        printf("%lld\n",res);
    }
    return 0;
}