bzoj 2839: 集合计数【容斥原理+组合数学】

首先,考虑容斥,我们所要的答案是并集至少有\( k \)个数的方案数减去并集至少有\( k+1 \)个数的方案数加上并集至少有\( k \)个数的方案数……
在n个数中选i个的方案数是\( C_{n}^{i} \),n种集合的组合方案数为\( 2^n \)
并集至少有i个元素的方案数即为选\( i \)个元素的方案数\( C_{n}^{i} \),乘上剩下\( n-i \)个元素任意组合的方案数\( 2{2{n-i}-1} \)
然后乘上容斥系数\( (-1)^{i-k} \),再乘上在并集的\( i \)个元素中选择\( k \)个元素的方案数\( C_{i}^{k} \)
答案即为:\( ans=\sum_{i=k}{i<=n}(-1){i-k}*C_{n}{i}*C_{i}{k}*2{2{n-i}-1} \),ans可能为负数,记得最后\( ans=(ans\%mod+mod)\%mod \)

#include<iostream>
#include<cstdio>
using namespace std;
const long long mod=1e9+7,N=1000005;
long long n,k,inv[N],fac[N],ans;
long long ksm(long long a,long long b)
{
	long long r=1ll;
	while(b)
	{
		if(b&1)
			r=r*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return r;
}
long long C(long long n,long long m)
{
	return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
int main()
{
	scanf("%lld%lld",&n,&k);
	fac[0]=1;
	for(long long i=1;i<=n;i++)
		fac[i]=fac[i-1]*i%mod;//,cout<<fac[i]<<" ";
	inv[n]=ksm(fac[n],mod-2);
	for(long long i=n-1;i>=0;i--)
		inv[i]=inv[i+1]*(i+1)%mod;//,cout<<inv[i]<<endl;;
	for(long long i=n,b=2;i>=k;i--,b=b*b%mod)
		ans=(ans+((((i-k)&1)?-1:1)*C(n,i)%mod*C(i,k)%mod*(b+mod-1)%mod))%mod;
	printf("%lld",(ans%mod+mod)%mod);
	return 0;
}
posted @ 2018-01-06 10:08  lokiii  阅读(228)  评论(0编辑  收藏  举报