###题意:
一个有N个元素的集合有
2
N
2^N
2N个不同子集(包含空集),现在要在这
2
N
2^N
2N个集合中取出若干集合(至少一个),使得它们的交集的元素个数为K,求取法的方案数,答案模1000000007
###题解:
好题。
一开始觉得应该很简单,然而很快就证明了我很天真。
设
f
[
k
]
=
C
n
k
∗
∑
i
=
1
2
n
−
k
C
2
n
−
i
i
=
C
n
k
∗
(
2
2
n
−
k
−
1
)
f[k]=C_n^k*\sum_{i=1}^{2^{n-k}}C_{2^{n-i}}^i=C_n^k*(2^{2^{n-k}}-1)
f[k]=Cnk∗∑i=12n−kC2n−ii=Cnk∗(22n−k−1)
这个大概就是表示规定选k个,在剩下
2
n
−
k
2^{n-k}
2n−k个集合中任意选的方案数。
然而这样算不但有交集大于k的情况,重复的也一大堆。
于是要容斥,考虑容斥系数。
f[k+1]显然是
−
C
k
+
1
k
-C_{k+1}^k
−Ck+1k
f[k+2]的第一项是
−
C
k
+
2
k
-C_{k+2}^k
−Ck+2k然而多减了f[k+1]的情况,所以要再容斥一下,加上
C
k
+
2
k
+
1
C
k
+
1
k
C_{k+2}^{k+1}C_{k+1}^{k}
Ck+2k+1Ck+1k。
根据
C
N
M
C
M
S
=
C
N
S
∗
C
N
−
S
M
−
S
C_N^MC_M^S=C_N^S*C_{N-S}^{M-S}
CNMCMS=CNS∗CN−SM−S
最后可以化简得:
g
[
k
]
=
∑
i
=
k
n
(
−
1
)
n
−
i
C
i
k
f
[
i
]
g[k]=\sum_{i=k}^n(-1)^{n-i}C_{i}^kf[i]
g[k]=i=k∑n(−1)n−iCikf[i]
update
显然
f
[
i
]
=
∑
j
=
i
n
C
j
i
g
[
j
]
f[i]=\sum_{j=i}^nC_j^ig[j]
f[i]=∑j=inCjig[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);
}