【bzoj2839】【2013福建】集合计数

一个有N个元素的集合有2^N个不同子集(包含空集),现在要在这2^N个集合中取出若干集合(至少一个),使得它们的交集的元素个数为K,求取法的方案数,答案模1000000007。

额很好的容斥题

额考虑选择所有的交集为C_{n}^{k}

这个很好理解但是剩下多少种集合:还有n-k个数

\sum_{i=0}^{n-k}C_{n-k}^{i}额稍有常识就知道这等价于2^{n-k}

然后从中选取集合这也是组合数

\sum_{i=0}^{2^{n-k}}C_{2^{n-k}}^{i}我们再次稍有常识一点

等价于2^{2^{n-k}}

额但是这个东西是什么:什么也不是他并不表示不少于k的集合选择数

对于长度为j的在这个式子中是重复出现了C_{j}^{k}

反向证明:在一个长为j的串出现了多少次?上式为一个很好的证明

然后容斥!

你看:设F_{i}=C_{n}^{i}*(2^{2^{n-i}}-1)

对于长度为i+1的,我们需要减去F_{i+1}C_{k+1}^{k}倍这好理解

因为对于这种情况实际上和上面的是一个原理

这里陷入一个递归

发现:对于F_{i+2}我们需要容斥掉C_{k+2}^{k}但是对于F_{i+1}C_{k+1}^{k}的被多减了所以乘上C_{k+2}^{k+1}

然后一顿操作(展开)加上C_{k+2}^{k}

所以答案为

\sum_{k}^{n}(-1)^{(n-k)}C_{i}^{k}

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
typedef int INT;
#define int long long
const int N=1000000+100;
int Quick_Pow(int x,int k,int mod){
	int ret=1;
	while(k){
		if(k%2){
			ret=ret*x%mod;
		}
		x=x*x%mod;
		k/=2;
	}
	return ret;
}
int fac[N+1];
int inv[N+1];
void Pre(){
	fac[0]=1;
	for(int i=1;i<N;++i){
		fac[i]=fac[i-1]*i%mod;
	}
	
	inv[N-1]=Quick_Pow(fac[N-1],mod-2,mod);
	for(int i=N-2;i>=0;--i){
		inv[i]=inv[i+1]*(i+1)%mod;
	}
}
int C(int n,int m){
	if(m>n)return 0;
	return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
int n,k;
int F(int x){
	return C(n,x)*((Quick_Pow(2,Quick_Pow(2,n-x,mod-1),mod)-1)%mod+mod)%mod;
}
int ans=0;
INT main(){
	cin>>n>>k;
	Pre();
	int cur=-1;
	for(int i=k;i<=n;++i){
		cur=cur*-1;
		ans=((ans+cur*F(i)*C(i,k)%mod)%mod+mod)%mod;
	}
	cout<<ans<<'\n';
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值