[HAOI2008]硬币购物

人生第一次开博客——一道动规题献上

（谢谢大家！谢谢大家对我的支持）
.小游戏——>>>2048—->>>flappy bird

题目：[HAOI2008]硬币购物http://www.lydsy.com/JudgeOnline/problem.php?id=1042

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

主线思路：dp预处理+容斥原理

设F[i]为不考虑每种硬币的数量限制的情况下，得到面值i的方案数。
为避免方案重复，要以k为阶段递推，边界条件为F[0]=1，这样预处理的时间复杂度就是O(S)。。。。可是没卵用啊Orz！！！（还是要靠我们的神牛沈大大）QWQ。。。。

接下来是一个神奇的方法（高能预警）：

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

当第1种硬币超过限制时，只要要用到D[1]+1枚硬币，剩余的硬币可以任意分配，所以方案数为 F[ S – (D[1]+1)C[1] ]，当且仅当(S – (D[1]+1)C[1])>=0，否则方案数为0。其余情况类似，每次询问只用问16次，所以询问的时间复杂度为O(1)。
（懵逼了吧！！！给你个眼神，体会去吧··················）

贴代码咯：

#include<cstdio>
#include<cstring>
using namespace std;
typedef long long LL;
LL c[5],d[5],f[110000],tot,ans;
void dfs(LL x,LL k,LL sum)
{
    if(sum<0)return;
    if(x==5)
    {
        if(k%2==1)ans-=f[sum];
        else ans+=f[sum];
        return;
    }
    dfs(x+1,k+1,sum-(d[x]+1)*c[x]);
    dfs(x+1,k,sum);
}
int main()
{
    for(int i=1;i<=4;i++)scanf("%d",&c[i]);
    scanf("%d",&tot);
    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]];
    while(tot--)
    {
        for(int i=1;i<=4;i++)scanf("%d",&d[i]);
        LL s;scanf("%d",&s);
        ans=0;
        dfs(1,0,s);
        printf("%lld\n",ans);
    }
    return 0;
}

欢迎评论！！！（谢谢大家！谢谢大家对我的支持~~~）