因为有了个数的限制,我们可以考虑容斥原理,首先对于整体求完全背包,dp值代表值为i的选择方式。
那么有一个超过限制的情况就是选择这个面值的多了一个,那么贡献就是c[i]*(d[i]+1),所以应该减去dp[s-(c[i]*(d[i]+1))]。
利用容斥,应该减去奇数个超限制的,加上偶数个超限制的。
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
int c[5],tot,d[5],s;
ll dp[100005],ans;
int main()
{
for(int i=1;i<=4;i++)scanf("%d",&c[i]);
dp[0]=1;scanf("%d",&tot);
for(int i=1;i<=4;i++)
{
for(int j=c[i];j<=100002;j++)
{
dp[j]+=dp[j-c[i]];
}
}
for(int i=1;i<=tot;i++)
{
for(int j=1;j<=4;j++)scanf("%d",&d[j]);ans=0;
scanf("%d",&s);
for(int j=0;j<=15;j++)
{
int tt=s;int cnt=0;
for(int k=1;k<=4;k++)
{
if((1<<(k-1))&j)
{
tt-=c[k]*(d[k]+1);cnt++;
}
}
if(tt<0)continue;
if(cnt&1)ans-=dp[tt];
else ans+=dp[tt];
}
printf("%lld\n",ans);
}
return 0;
}