【DP计划】11.6——[BZOJ]硬币购物（背包+容斥）MEDIUM

Description

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

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

每次的方法数
Sample Input
1 2 5 10 2

3 2 3 1 10

1000 2 2 2 900
Sample Output
4

27

---

这道题的思路很清奇。

首先，最原始的多重背包肯定会T。但由于硬币总共只有4种，所以可以先把这道题当做完全背包来做，然后进行容斥。

具体来说，我们可以枚举一个二进制状态，上面的1表示用这种硬币，0表示不选。然后进行容斥。最先的状态是0，所以答案加上$f[tot-0]$，也就是完全背包做出来的总方案数。然后接下来的状态时1，所以减去$f[tot-c[1]\ast (d[1]+1)]$，为什么要$+1$，可以这么理解：我们需要减去取多于$d[i]$个的方案，所以最小的多于$d[i]$个的状态量就是$tot-c[i]\ast (d[i]+1)$，以此类推，容斥下去即可。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll read(){
    char c;ll x;while(c=getchar(),c<'0'||c>'9');x=c-'0';;
    while(c=getchar(),c>='0'&&c<='9') x=x*10+c-'0';return x;
}
ll ans,c[5],T,tot,d[5],f[100005],ct[20];
ll calc(ll x){
    ll res=0;
    for(ll p=1;x;x/=2,p++){
        if(x&1) res+=c[p]*(d[p]+1);
    }
    return tot-res<0?0:f[tot-res];
}
int main()
{
    c[1]=read();c[2]=read();c[3]=read();c[4]=read();T=read();
    for(ll i=1;i<(1<<4);i++) ct[i]=ct[i>>1]+(i&1);
    while(T--){
        memset(f,0,sizeof(f));f[0]=1;ans=0;
        d[1]=read();d[2]=read();d[3]=read();d[4]=read();tot=read();
        for(ll i=1;i<5;i++)
         for(ll j=c[i];j<=tot;j++)f[j]+=f[j-c[i]];
        for(ll i=0;i<(1<<4);i++){
            if(ct[i]&1) ans-=calc(i);
            else ans+=calc(i);
        }
        printf("%lld\n",ans);
    }
    return 0;
}