【DP】硬币购物

【题意】
一共有4种硬币。面值分别为c1,c2,c3,c4。某人去商店买东西，去了tot次。每次带di枚面值为ci的硬币，买si的价值的东西。请问每次有多少种付款方法。
【思路】
设f[i]为不考虑硬币限制下，买价值为i的物品的方案数。那么很容易得出f[j]=f[j-c[1]]+f[j-c[2]]…的递推式。注意要推完一种硬币再推下一种硬币。

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

例如这样的推法就是错的。假设我们的四种硬币面值为2，3，100，1000，那么按照这个推法，f[5]=f[2]+f[3]，(2,3)与(3,2)被重复计数了。实际上的重复情况复杂的多，我们不可能用/2等方法解决它。因此，要推完一种硬币再推下一种硬币。

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

这样的方法就是对的。
对于每次询问，我们先假设答案为f[s]，但s>=(d[1]+1)*c[1]，即f[s]中包含了第一种硬币超过限制的情况。我们要减去这部分。
很明显，对于第一种硬币超过限制的情况，所用的第一种硬币数量>=d[1]+1。
容易发现，总和为s，第一种硬币至少d[1]+1的情况个数，等于总和为s-(d[1]+1)*c[1]，硬币个数没有限制的情况个数！
因为总和为s，第一种硬币至少d[1]+1的情况个数中，那d[1]+1个硬币是固定的，是不会变的，对于情况个数是不会有影响的。
所以，我们将原先为f[s]的答案减去f[s-(d[1]+1)*c[1]]即可得出正确答案。
这只是一种硬币超过限制的解决办法。对于两种硬币超过限制我们也同样这么减去，但重复减去了一部分：那些两种硬币同时超过限制的情况。
因此我们要加回两种硬币同时超过限制的情况。根据上面的式子，很容易就知道这样做就是

if(s>=(d[1]+1)*c[1]+(d[2]+1)*c[2])//说明f[s]包含了两种同时超过限制的情况
    ans+=f[s-(d[1]+1)*c[1]-(d[2]+1)*c[2]]

然后这样加回去又会多加了三种超过限制的，所以三种的也要减去。然后三种的又会多减了四种的，所以四种要加回。
【代码】

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 110000
#define LL long long
using namespace std;
LL read()
{
    LL ans=0,f=1;char x=getchar();
    while(x<'0'||x>'9')
    {
        if(x=='-')f=-1;
        x=getchar();
    }
    while(x>='0'&&x<='9')
    {
        ans=ans*10+x-'0';
        x=getchar();
    }
    return ans*f;
}
int m;
LL c[5],d[5],f[MAXN];
int main()
{
    for(int i=1;i<=4;i++)c[i]=read();

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

    scanf("%d",&m);
    while(m--)
    {
        for(int i=1;i<=4;i++)d[i]=read()+1;
        LL x=read(),ans = f[x];

        for(int i=1;i<=4;i++)
            if(x>=d[i]*c[i])
                ans-=f[x-d[i]*c[i]];

        for(int i=1;i<=3;i++)
            for(int j=i+1;j<=4;j++)
            {
                if(i==j)continue;
                if(x>=d[i]*c[i]+d[j]*c[j])
                    ans+=f[x-d[i]*c[i]-d[j]*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(i==j||j==k||i==k)continue;
                    if(x>=d[i]*c[i]+d[j]*c[j]+d[k]*c[k])
                        ans-=f[x-d[i]*c[i]-d[j]*c[j]-d[k]*c[k]];
                }

        if(x>=d[1]*c[1]+d[2]*c[2]+d[3]*c[3]+d[4]*c[4])
            ans+=f[x-d[1]*c[1]-d[2]*c[2]-d[3]*c[3]-d[4]*c[4]];
        printf("%lld\n",ans);

    }
    return 0;
}