BZOJ1042【HAOI2008】硬币购物（DP+容斥原理）

【HAOI2008】硬币购物（DP+容斥原理）

1042: [HAOI2008]硬币购物
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

多谢hzwer大神，他这题的博客(http://hzwer.com/5286.html)给我带来了很多收获，代码也是按照他思路打的，所以这篇其实是【半原创】。

这题的核心就是容斥原理，先dp预处理没有硬币使用数量限制的方法数，就是完全背包，这一部分的代码如下：

maxs=100000;
f[0]=1;
    for(int i=1;i<=4;i++)
        for(int j=c[i];j<=maxs;j++)
            f[j]+=f[j-c[i]];

现在我们已经得出了没有硬币使用数量限制的方法数，但是这个结果比正解要多许多不必要的东西，就是超过使用限制的方法数，那些是要减掉的，那么减掉什么呢？

根据容斥原理
答案为 ：
得到面值S的超过限制的方案数 – 第1种硬币超过限制的方案数 – 第2种硬币超过限制的方案数 – 第3种硬币超过限制的方案数 – 第4种硬币超过限制的方案数 + 第1,2种硬币同时超过限制的方案数 + 第1,3种硬币同时超过限制的方案数 + …… + 第1,2,3,4种硬币全部同时超过限制的方案数。

我相信许多人看了这一大坨文字就会跟我当初一样，一脸懵逼

其实很简单，我们有许多多余的方法数，在这很多的多余的方法数中有第1,2,3,4种硬币超过限制的方案数，我们要减去他，但如果减去了他们其实就会减去许多重复的东西，我们只是约定了第1种一定要超过限制，但并没有限制第2,3,4种的方案数超过限制，所以在第1种超过限制的方案数中，含有第2,3,4种超过限制的方案数，他们被我们减掉了两次，所以我们要加上他们，之后的同理。

这一部分的代码如下：

void dfs(int x,int k,int S)//k代表超过限制的硬币种数，S代表当前要完成的总钱数
{
    if(S<0) return;
    if(x==5)
    {
        if(k%2==1) ans-=f[S];
        else ans+=f[S];
        return;
    }
    dfs(x+1,k+1,S-(d[x]+1)*c[x]);
    dfs(x+1,k,S);
}

完整代码如下：

#include<cstdio>
#include<iostream>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxs=100000;
LL f[110000];
int tot,c[5],d[5];
LL ans,S;
void dfs(int x,int k,int S)
{
    if(S<0) return;
    if(x==5)
    {
        if(k%2==1) ans-=f[S];
        else ans+=f[S];
        return;
    }
    dfs(x+1,k+1,S-(d[x]+1)*c[x]);
    dfs(x+1,k,S);
}
int main()
{
    for(int i=1;i<=4;i++) scanf("%d",&c[i]);
    scanf("%d",&tot);
    f[0]=1;
    for(int i=1;i<=4;i++)
        for(int j=c[i];j<=maxs;j++)
            f[j]+=f[j-c[i]];
    for(int i=1;i<=tot;i++)
    {
        for(int j=1;j<=4;j++) scanf("%d",&d[j]);
        scanf("%d",&S);
        ans=0;dfs(1,0,S);
        printf("%lld\n",ans);
    }
}