bzoj1042: [HAOI2008]硬币购物

题意不用说了吧。一看大概就知道是背包。。
但是用背包做了半天还是没做出来
经过一位神犇的指导知道了这题要用到容斥原理
先两个for搞出4种硬币没有限制的情况下的方案数
我们要求的答案 Ans 就是 f[S] - (硬币1超限制的方案数) - (硬币2超限制的方案数) - (硬币3超限制的方案数) - (硬币4超限制的方案数) + (硬币1,2超限制的方案数) + (硬币1,3超限制的方案数) + (硬币1,4超限制的方案数) + …. - (硬币1,2,3超限制的方案数) - … + (硬币1,2,3,4超限制的方案数) 。　　

#include<cstdio>
#include<cstring>
using namespace std;
typedef long long LL;
LL c[5],d[5],f[110000];
int main() {
    for(int i=1;i<=4;i++)
        scanf("%d",&c[i]);
    memset(f,0,sizeof(f));f[0]=1;
    for(int i=1;i<=4;i++)
        for(int j=c[i];j<=100000;j++)
            f[j]+=f[j-c[i]];
    int n,s;
    scanf("%d",&n);
    while(n--) {
        for(int i=1;i<=4;i++)scanf("%d",&d[i]);
        scanf("%d",&s);LL ans=f[s];
        for(int i=1;i<=4;i++)
            if(s>=(d[i]+1)*c[i])
                ans-=f[s-(d[i]+1)*c[i]];
        for(int i=1;i<=3;i++)
            for(int j=i+1;j<=4;j++)
                    if(s>=((d[i]+1)*c[i]+(d[j]+1)*c[j]))
                        ans+=f[s-((d[i]+1)*c[i]+(d[j]+1)*c[j])];
        for(int i=1;i<=2;i++)
            for(int j=i+1;j<=3;j++)
                    for(int k=j+1;k<=4;k++)
                            if(s>=((d[i]+1)*c[i]+(d[j]+1)*c[j]+(d[k]+1)*c[k]))
                                ans-=f[s-((d[i]+1)*c[i]+(d[j]+1)*c[j]+(d[k]+1)*c[k])];
        if(s>=((d[1]+1)*c[1]+(d[2]+1)*c[2]+(d[3]+1)*c[3]+(d[4]+1)*c[4]))
            ans+=f[s-((d[1]+1)*c[1]+(d[2]+1)*c[2]+(d[3]+1)*c[3]+(d[4]+1)*c[4])];
        printf("%lld\n",ans);
    }
    return 0;
}