【CodeForces】908 E. New Year and Entity Enumeration

【题目】E. New Year and Entity Enumeration

【题意】给定集合T包含n个m长二进制数，要求包含集合T且满足以下条件的集合S数：长度<=m，非和与的结果都在集合中。（详细的题意见原题）

【算法】数学（贝尔数）

【题解】这道题确实不太能理解这种做法，所以就简单写写了。

先不考虑S须包含集合T。

对于一个方案，按位考虑，所有含位 i 的数字and起来得到含该位的最小数字，记为f[i]。

对于f[x]≠f[y]，有f[x]&f[y]=0，证明：!(f[x]&f[y])&f[x]这个数字含有x位且<f[x]。

那么不同位的f[x]，要么相等要么不等且无交集，那么方案数对应1~m的集合划分数（贝尔数）。

贝尔数：B(n)=ΣC(n-1,k)*B(k)，k=0~n-1。求解复杂度O(m^2)。

最后考虑T，不同位如果竖着看的二进制数不同那么其f值一定不能相同，所以分成若干部分各自求解后再相乘即是答案。

#include<cstdio>
#include<map>
#define ll long long
using namespace std;
const int maxn=1010,MOD=1e9+7;
int c[maxn][maxn],f[maxn],n,m;
ll b[maxn];
map<ll,int>mp;
int main(){
    scanf("%d%d",&m,&n);
    for(int i=0;i<n;i++){
        for(int j=1;j<=m;j++){
            int x;
            scanf("%1d",&x);
            b[j]+=(1ll*x)<<i;
        }
    }
    for(int i=1;i<=m;i++)mp[b[i]]++;
    for(int i=0;i<=m;i++){
        c[i][0]=1;
        for(int j=1;j<=i;j++)c[i][j]=(c[i-1][j-1]+c[i-1][j])%MOD;
    }
    f[0]=1;
    for(int i=1;i<=m;i++){
        for(int j=0;j<i;j++)f[i]=(f[i]+1ll*c[i-1][j]*f[j]%MOD)%MOD;
    }
    int ans=1;
    for(map<ll,int>::iterator it=mp.begin();it!=mp.end();it++)ans=1ll*ans*f[it->second]%MOD;
    printf("%d",ans);
    return 0;
}