题意:
有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
*/