[FROM LUOGU][HAOI2008]硬币购物

[HAOI2008]硬币购物

传送门

SOL
最初可能会考虑一个完全背包，然而会T
考虑完全背包是不是只用做一次？
答案是肯定的
对于一个限制d[i]，我们用d[i]+1个i硬币就不合法
那么对于一个询问s，i硬币不合法的情况数就是f[s-(d[i]+1)*c[i]]
也就是说我们先用够总面值为s-(d[i]+1)*c[i]的硬币，剩下的再用c[i]补足，那么就一定不合法
对于其他几种硬币显然，然后再考虑几种限制同时不满足
剩下的就是容斥了
另外，可以用位运算优化子集枚举

代码：

#include<bits/stdc++.h>
using namespace std;
#define re register
inline char nc(){
    static char buf[10000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,10000,stdin),p1==p2)?EOF:*p1++;
}
inline int rd(){
	int re data=0;static char ch=0;ch=nc();
	while(!isdigit(ch))ch=nc();
	while(isdigit(ch))data=(data<<1)+(data<<3)+(ch^48),ch=nc();
	return data;
}
int c[4],d[4],tot,s,k[16]={1,-1,-1,1,-1,1,1,-1,-1,1,1,-1,1,-1,-1,1};long long f[100005],ans;
signed main(){
	for(int re i=0;i^4;++i)c[i]=rd();f[0]=1,tot=rd();
	for(int re i=0;i^4;++i)for(int re j=c[i];j^100001;++j)f[j]+=f[j-c[i]];
	while(tot--){
		for(int re i=0;i^4;++i)d[i]=rd()+1;s=rd(),ans=f[s];
		for(int re sum,i=1;i^16;++i){
			sum=0;
			for(int re tmp=i,j=0;tmp;tmp>>=1,++j)if(tmp&1)sum+=d[j]*c[j];//前面已经将d[j]处理成d[j]+1 
			if(s>=sum)ans+=k[i]*f[s-sum];
		}printf("%lld\n",ans);
	}
	return 0;
}