codeforces 451E Devu and Flowers (容斥原理)

题意:

有n个花瓶,每个花瓶里面有数枝话,现在要从这些花瓶从取出s枝花,问有多少种方案。

题解:

直接计算这个组合数的哀家并不会。可以考虑用容斥原理,在每个花瓶中取花和将花放到每个花瓶里面的问题是等价的,那么就这样考虑:我们有s枝花,要放入n个空花瓶有多少方案。我们任意枚举花瓶的编号,开始不选择花瓶,那么要将这么花分配到各个花瓶的方案是C(s+n-1,n-1)(搁板法),那坑定有重复或者不可取的方案,因为有些花瓶的容量没那么大,那么第二次枚举一个花瓶同样计算,第三次两个花瓶。。。。

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
#define B(x) (1<<(x))
using namespace std;
typedef long long ll;
void cmax(int& a,int b){ if(b>a)a=b; }
void cmin(int& a,int b){ if(b<a)a=b; }
void cmax(ll& a,ll b){ if(b>a)a=b; }
void cmin(ll& a,ll b){ if(b<a)a=b; }
void add(int& a,int b,int mod){ a=(a+b)%mod; }
void add(ll& a,ll b,ll mod){ a=(a+b)%mod; }
const int oo=0x3f3f3f3f;
const ll MOD=1000000007;
int n;
ll s;
ll f[22];

ll quick_pow(ll a,ll k){
    ll ans=1;
    while(k){
        if(k&1) ans=(ans*a)%MOD;
        a=(a*a)%MOD;
        k>>=1;
    }
    return ans;
}

ll C(ll x,ll y){
    ll a=1,b=1;
    for(ll i=1;i<=y;i++){
        a=((x-i+1)%MOD*a)%MOD;
        b=(b*i)%MOD;
    }
    return a*quick_pow(b,MOD-2)%MOD;
}

ll C_mod(ll x,ll y){
    if(y==1)
        return 0;
    return C(x%MOD,y%MOD)*C_mod(x/MOD,y/MOD)%MOD;
}

ll solve(){
    int full=B(n)-1;
    ll ans=0;
    for(int i=0;i<=full;i++){
        ll sum=s,cnt=0;
        for(int j=0;j<n;j++)if(i&B(j)){
            cnt++;
            sum-=(f[j]+1);
        }
        if(sum<0)continue;
        if(cnt&1) add(ans,-C(sum+n-1,n-1),MOD);
        else add(ans,C(sum+n-1,n-1),MOD);
    }
    return (ans+MOD)%MOD;
}

int main()
{
    while(scanf("%d %I64d",&n,&s)!=EOF){
        for(int i=0;i<n;i++)scanf("%I64d",&f[i]);
        cout<<solve()<<endl;
    }
    return 0;
}
/**
2 3
1 3
2 4
2 2
3 5
1 3 2

2 270030023747
891135146290 437305641972
ans=30021858

20 4385085334307
273634411136 208521328637 450482376435 844118010709 197241285878 472126475472 2414038897 672334205413 809269727018 409013884362 739986692075 953956651947 462216461906 388007176838 245504550965 527140291750 632844435887 550532123833 757200390348 944901802640

*/




  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值