bzoj4671: 异或图【容斥原理+线性基】

解题思路：

考虑容斥，枚举点的集合划分，强制两两集合间的点不连通，集合内的点任意连，若划分成m个集合，则最后至少有m个连通块。

而一个m的划分，在容斥时会被计算$\sum\limits_{i=1}^mS_m^i$次($S$为第二类斯特林数)。

所以列出容斥系数计算式：$\sum\limits_{i=1}^mS_m^if_i=[m=1]$，

打表找规律可得$f_i=(-1)^{i-1}(i-1)!$

那么现在只用考虑若我们dfs出了一种划分，如何计算满足的方案数？

考虑把每张图看做一个二进制数，每一位代表一个属于不同集合的点对，要使所有这类点对不连通，那就是问该二进制数集合有多少子集异或和为0。

考虑线性基，那么先选出线性无关组后，加上任意剩余数组成的集合都有且仅有一个子集异或和为0，且互不相同。那么设线性无关组个数为cnt，方案数即为$2^{s-cnt}$。

也可以将该矩阵转置来算，线性无关组个数是一样的，但好像会快些。

#include<bits/stdc++.h>
#define ll unsigned long long
using namespace std;

int s,n;
int g[65][15][15],id[15];
ll ans,fac[11],base[65];
char ch[100];

void dfs(int x,int num)
{
    if(x==n+1)
    {
        int cnt=0;
        for(int k=1;k<=s;k++)
        {
            int l=0;ll t=0;
            for(int i=1;i<=n;i++)
                for(int j=i+1;j<=n;j++)
                    if(id[i]!=id[j])t|=(1ll<<l)*g[k][i][j],l++;
            for(int i=1;i<=cnt;i++)
                if((t^base[i])<t)t^=base[i];
            if(t)base[++cnt]=t;
        }
        ans+=1ll*((num&1)?1:-1)*fac[num-1]*(1ll<<(s-cnt));
        return;
    }
    for(int i=1;i<=num+1;i++)
    {
        id[x]=i;
        dfs(x+1,num+(i>num));
    }
}

int main()
{
    //freopen("lx.in","r",stdin);
    scanf("%d",&s);
    for(int k=1;k<=s;k++)
    {
        scanf("%s",ch);int l=strlen(ch);
        for(int i=1;!n;i++)if(i*(i-1)/2==l)n=i;
        l=0;
        for(int i=1;i<=n;i++)
            for(int j=i+1;j<=n;j++)
                g[k][i][j]=ch[l++]-'0';     
    }
    fac[0]=1;
    for(int i=1;i<=n;i++)fac[i]=fac[i-1]*i;
    dfs(1,0);
    cout<<ans;
    return 0;
}