首先,考虑容斥,我们所要的答案是并集至少有\( 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;
}