BZOJ1042

最新推荐文章于 2019-09-13 22:33:37 发布

Le_ballon_rouge

最新推荐文章于 2019-09-13 22:33:37 发布

阅读量1k

点赞数 1

分类专栏：组合数学文章标签：容斥原理

本文链接：https://blog.csdn.net/Le_ballon_rouge/article/details/47758289

版权

组合数学专栏收录该内容

2 篇文章 0 订阅

订阅专栏

传送门：BZOJ1042

首先，计算出购买面值为i的物品的方案数f(i),这一步强制有序就可以了。
然后，每一次查询时ans=(f(s)-d1溢出方案-d2…+d1d2+d2d3+…-d1d2d3….+d1d2d3d4)，即容斥原理。
注意到d1溢出时，至少使用了d1+1个物品，于是剩下S-(d1+1)c1都可以随意分配，于是d1溢出的方案就是f(S-(d1+1)c1)

代码上的小细节见下。

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;

long long f[100005];
int n;
int c[5];
int W[10];
long long ans;
int tot;
int d[5],S;

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

void Dfs(int x,int S,int num)
{
    if(S<0)
        return;
    if(num==5){
        if(x&1)
            ans-=f[S];
        else
            ans+=f[S];
        return;
    }
    Dfs(x+1,S-(d[num]+1)*c[num],num+1);
    Dfs(x,S,num+1);
}

void Solve()
{
    for(int i=1;i<=tot;i++){
        scanf("%d%d%d%d%d",&d[1],&d[2],&d[3],&d[4],&S);
        ans=0;
        Dfs(0,S,1);
        cout<<ans<<endl;
    }
}

void Readdata()
{
    freopen("loli.in","r",stdin);
    scanf("%d%d%d%d%d",&c[1],&c[2],&c[3],&c[4],&tot);
}

void Close()
{
    fclose(stdin);
    fclose(stdout);
}

int main()
{
    Readdata();
    First();
    Solve();
    Close();
    return 0;
}