【BZOJ2339】[HNOI2011]卡农
题解:虽然集合具有无序性,但是为了方便,我们先考虑有序的情况,最后将答案除以m!即可。
考虑DP。如果我们已经知道了前m-1个集合,那么第m个集合已经是确定的了。因为内层集合的n个元素可以随便出现,那么总数就是A(2^n-1,m-1)。但是可能存在不合法的情况。
1.在前m-1个集合中,n个数出现的次数已经都是偶数了,那么第m个集合为空,不合法,此时方案数为f[m-1]。
2.第m个集合与之前某个集合相同,那么我们不考虑这两个集合,剩下的方案数为f[i-2];该集合可能是第1...m-1个;该集合可能有2^n-1-(m-2)中情况。所以方案数为f[i-2]*(m-1)*(2^n-1-(m-2))
#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
typedef long long ll;
const ll mod=100000007;
ll n,m,n2,c,m1;
ll ine[1000010],f[1000010];
int main()
{
scanf("%lld%lld",&n,&m);
ll i;
ine[1]=1;
for(i=2;i<=m;i++) ine[i]=(mod-(mod/i)*ine[mod%i]%mod)%mod;
for(n2=i=1;i<=n;i++) n2=(n2<<1)%mod;
n2=(n2-1+mod)%mod,c=1,f[0]=1,f[1]=0,m1=1;
for(i=2;i<=m;i++)
{
c=c*(n2-i+2)%mod,m1=m1*ine[i]%mod;
f[i]=((c-f[i-1]-f[i-2]*(n2-i+2)%mod*(i-1))%mod+mod)%mod;
}
printf("%lld",f[m]*m1%mod);
return 0;
}