bzoj1042: [HAOI2008]硬币购物 容斥原理

神一般的容斥原理，主要是方法难想

题意：给出四种面值的硬币c1,c2,c3,c4。n个询问。每次询问用d1、d2、d3、d4个相应的硬币能够拼出多少种总和为s？

思路：（1）首先，用完全背包求出f[i]表示四种硬币的数量无限制拼出i的方案数。

（2）接着我们来理解x=f[s]-f[s-(d1+1)*c1]的含义：x表示c1硬币的数量不超过d1个而其他三种硬币的数量不限制拼成s的方案数。我们举着例子来说明，假设现在有两种硬币，面值分别为1和2，那么我们求出f:f[0]=1,f[1]=1,f[2]=2,f[3]=2,f[4]=3,f[5]=3,f[6]=4。其中f[3]的两种分别为3=1+1+1=1+2,f[6]的四种为:6=1+1+1+1+1+1=1+1+1+1+2=1+1+2+2=2+2+2。加入我们现在求第一种硬币最多使用两个，第二种硬币无限制的方案数，按照我们说的x=f[6]-f[6--(2+1)*1]=f[6]-f[3]=2。也就是6=1+1+2+2=2+2+2两种。我们发现我们删除了1+1+1+1+1+1和1+1+1+1+2两种，为什么能够通过减去f[3]删掉这两种？我们来看f[3]，3=1+1+1=1+2，我们发现6中被删掉的两种正是通过这个f[3]增加3个1得到的。

#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <queue>

using namespace std;
typedef long long sint;
sint f[110000];
int c[10],n;
int s[10],sum;
sint f2(int k)
{
    if(k>=0) return f[k];
    return 0;
}
int main()
{
    for(int j=1;j<=4;j++) scanf("%d",&c[j]);
    scanf("%d",&n);
    f[0]=1;
    for(int i=1;i<=4;i++)
    {
        for(int j=1;j<=100100;j++)
        {
            if(j>=c[i])
            {
                f[j]+=f[j-c[i]];
            }
        }
    }
    sint ans;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=4;j++) scanf("%d",&s[j]);
        scanf("%d",&sum);
        ans=f[sum];
        ans-=f2(sum-(s[1]+1)*c[1]);
        ans-=f2(sum-(s[2]+1)*c[2]);
        ans-=f2(sum-(s[3]+1)*c[3]);
        ans-=f2(sum-(s[4]+1)*c[4]);
        ans+=f2(sum-(s[1]+1)*c[1]-(s[2]+1)*c[2]);
        ans+=f2(sum-(s[1]+1)*c[1]-(s[3]+1)*c[3]);
        ans+=f2(sum-(s[1]+1)*c[1]-(s[4]+1)*c[4]);
        ans+=f2(sum-(s[2]+1)*c[2]-(s[3]+1)*c[3]);
        ans+=f2(sum-(s[2]+1)*c[2]-(s[4]+1)*c[4]);
        ans+=f2(sum-(s[3]+1)*c[3]-(s[4]+1)*c[4]);
        ans-=f2(sum-(s[1]+1)*c[1]-(s[2]+1)*c[2]-(s[3]+1)*c[3]);
        ans-=f2(sum-(s[1]+1)*c[1]-(s[2]+1)*c[2]-(s[4]+1)*c[4]);
        ans-=f2(sum-(s[1]+1)*c[1]-(s[3]+1)*c[3]-(s[4]+1)*c[4]);
        ans-=f2(sum-(s[2]+1)*c[2]-(s[3]+1)*c[3]-(s[4]+1)*c[4]);
        ans+=f2(sum-(s[1]+1)*c[1]-(s[2]+1)*c[2]-(s[3]+1)*c[3]-(s[4]+1)*c[4]);
        printf("%lld\n",ans);
    }
    return 0;
}