考慮有序選擇各子集，最後除以m!即可。設f[i]為選i個子集的合法方案數。

　　對f[i]考慮容斥，先只滿足所有元素出現次數為偶數。確定前i-1個子集後第i個子集是確定的，那麽方案數為A(2ⁿ-1,i-1)。

　　顯然不能為空集，於是去掉前i-1個已經滿足限制的方案，也即f[i-1]。

　　然後去掉第i個子集和之前重復的情況。顯然如果有重復，將這兩個去掉後仍然是合法的。那麽方案數為f[i-2]*(i-1)*(2ⁿ-1-(i-2))。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include
 
<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=0,f=1;char c=getchar();
    while (c<‘0‘||c>‘9‘) {if (c==‘-‘) f=-1;c=getchar();}
    while (c>=‘0‘&&c<=‘9‘) x=(x<<1)+(x<<3)+(c^48),c=getchar();
    return x*f;
}
 
#define P 100000007
#define N 1000010
int n,m,f[N],inv[N],p,A[N];
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj2339.in","r",stdin);
    freopen("bzoj2339.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),m=read();
    p=1;for (int i=1;i<=n;i++) p=(p<<1
 
)%P;p--;
    inv[1]=1;
    for (int i=2;i<=m;i++) inv[i]=P-1ll*(P/i)*inv[P%i]%P;
    for (int i=2;i<=m;i++) inv[i]=1ll*inv[i]*inv[i-1]%P;
    A[0]=1;for (int i=1;i<=m;i++) A[i]=1ll*A[i-1]*(p-i+1+P)%P;
    f[0]=1;f[1]=0;
    for (int i=2;i<=m;i++) f[i]=((A[i-1]-f[i-1]+P)%P-1ll*f[i-2]*(i-1)%P*(p-i+2+P)%P+P)%P;
    cout<<1ll*f[m]*inv[m]%P;
    return 0;
}

BZOJ2339 HNOI2011卡農（動態規劃+組合數學）