codeforces 451E. Devu and Flowers 组合数+容斥

题意:n个不同盒子,每个盒子里有ai个小球,总共拿s个,问有多少种选法

如果盒子里求是无限个,直接插板法就可以求了

对于有了个数限制,可以容斥来做,需要去掉的是个数超限的,那么我先保证他个数超限(拿ai+1个),然后对剩下的小球(s-ai-1)按原来的插板法求放法就可以了


#include<bits/stdc++.h>
#define eps 1e-9
#define PI 3.141592653589793
#define bs 1000000007
#define bsize 256
#define MEM(a) memset(a,0,sizeof(a))
typedef long long ll;
using namespace std;
const int mo=1e9+7; 
ll qpow(ll a,ll k)
{
	ll ans=1;
	while(k)
	{
		if(k&1)
		ans=(ans*a)%mo;
		a=(a*a)%mo;
		k>>=1; 
	}
	return ans;
 } 
ll C(ll a,ll b)
{
	if(a<b)
	return 0;
	if(b>a-b)
	b=a-b;
	ll up=1,down=1;
	for(ll i=0;i<b;i++)
	{
		up=up*(a-i)%mo;
		down=down*(i+1)%mo;
	}
	return up*qpow(down,mo-2)%mo;
}
ll lucas(ll a,ll b)
{
	if(b==0)
	return 1;
	return C(a%mo,b%mo)*lucas(a/mo,b/mo)%mo; 
}
ll s,f[100];
int main()
{
	int n,i,j;
	cin>>n>>s;
	for(i=0;i<n;i++)
	cin>>f[i];
	ll ans=0;
	int up=(1<<n)-1;
	for(i=0;i<=up;i++)
	{
		int sta=1;
		ll now=s;
		for(j=0;j<n;j++)
		{
			if(i&(1<<j))
			{
				sta*=-1;
				now-=f[j]+1;
			}
		}
	//	cout<<lucas(now+n-1,n-1)<<" "<<now<<endl;
		if(now>=0)
		ans+=sta*lucas(now+n-1,n-1);
		ans=(ans+mo)%mo;
	}
	cout<<ans<<endl;

 }


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值