bzoj 2839: 集合计数

最新推荐文章于 2019-03-05 15:47:00 发布

fyc_kabuto

最新推荐文章于 2019-03-05 15:47:00 发布

阅读量370

点赞数 1

分类专栏：容斥原理排列组合

本文链接：https://blog.csdn.net/qq_36808030/article/details/78986531

版权

容斥原理同时被 2 个专栏收录

27 篇文章 0 订阅

订阅专栏

排列组合

21 篇文章 0 订阅

订阅专栏

###题意：
一个有N个元素的集合有 $2^N$ 个不同子集（包含空集），现在要在这 $2^N$ 个集合中取出若干集合（至少一个），使得它们的交集的元素个数为K，求取法的方案数，答案模1000000007
###题解：
好题。
一开始觉得应该很简单，然而很快就证明了我很天真。
设 $f[k]=C_n^k*\sum_{i=1}^{2^{n-k}}C_{2^{n-i}}^i=C_n^k*(2^{2^{n-k}}-1)$
这个大概就是表示规定选k个，在剩下 $2^{n-k}$ 个集合中任意选的方案数。
然而这样算不但有交集大于k的情况，重复的也一大堆。
于是要容斥，考虑容斥系数。
f[k+1]显然是 $C_{k+1}^k$
f[k+2]的第一项是 $C_{k+2}^k$ 然而多减了f[k+1]的情况，所以要再容斥一下，加上 $C_{k+2}^{k+1}C_{k+1}^{k}$ 。
根据 $C_N^MC_M^S=C_N^S*C_{N-S}^{M-S}$
最后可以化简得： $g[k]=\sum_{i=k}^n(-1)^{n-i}C_{i}^kf[i]$

update

显然 $f[i]=\sum_{j=i}^nC_j^ig[j]$
直接二项式反演即可
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
const LL mod=1000000007;
LL n,K,f[1000010],inv[1000010],t=1,k[1000010];
void pre()
{
	f[0]=f[1]=inv[0]=inv[1]=k[0]=1;
	for(LL i=2;i<=1000000;i++) f[i]=(f[i-1]*i)%mod,inv[i]=(mod-mod/i)*inv[mod%i]%mod;
	for(LL i=1;i<=1000000;i++) (inv[i]*=inv[i-1])%=mod,k[i]=k[i-1]*2%(mod-1);
}
LL pow(LL a,LL b)
{
	LL ans=1;
	while(b)
	{
		if(b&1) (ans*=a)%=mod;
		(a*=a)%=mod;b>>=1;
	}
	return ans;
}
LL C(LL n,LL m)
{
	if(m>n) return 0;
	return f[n]*inv[n-m]%mod*inv[m]%mod;
}
int main()
{
	pre();
	scanf("%lld %lld",&n,&K);
	LL ans=0;
	for(LL i=K;i<=n;i++,t=-t)
		ans+=C(n,i)*(pow(2,k[n-i])-1)%mod*C(i,K)%mod*t,ans%=mod;
	(ans+=mod)%=mod;
	printf("%lld",ans);
}