[BZOJ4671/JZOJ4591]异或图

题目大意

定义两个结点数相同的图$G_1$与图$G_2$的异或为一个新的图$G$，其中如果$(u,v)$在$G_1$与$G_2$中的出现次数之和为$1$，那么边$(u,v)$在$G$中，否则这条边不在$G$中。
现在给定$s$个结点数相同的图$G_{1...s}$，设$S=\{G_1,G_2,...,G_s\}$, 请问$S$有多少个子集的异或为一个连通图?

$2\le n\le 10,1\le s\le 60$

---

题目分析

考虑容斥：我们用$B_n$的时间来枚举$n$个点的子集划分方案，然后限制不同子集的点一定不可以连边，同一子集任意。
这样我们就可以列出若干条$s$个未知数的异或方程组，做一遍高斯消元求出自由元个数就可以知道相应的方案。
既然是容斥，我们就要想办法通过容斥系数来减去不合法的方案。
考虑如果一幅图实际上被划分成了$n$个连通块，它会被我计算多少次

$$\sum_{m=1}^nS(n,m)$$
其中$S$是第二类斯特林数。
考虑对于一个$m$个连通块的局面，我定义它的容斥系数是${(-1)}^{m-1}(m-1)!$，那么一个$n$个连通块的局面会被我计算多少次呢？
$$\sum_{m=1}^n(-1)^{m-1}(m-1)!S(n,m)$$
首先显然$n=1$的时候这条式子等于$1$，那么$n>1$时呢？

$$\begin{align} \sum_{m=1}^n(-1)^{m-1}(m-1)!S(n,m)&=\sum_{m=1}^{n-1}{(-1)}^{m-1}(m-1)!\left[S(n-1,m-1)+mS(n-1,m)\right]+{(-1)}^n(n-1)!\\ &=\sum_{m=1}^{n-1}{(-1)}^{m-1}(m-1)!S(n-1,m-1)+\sum_{m=1}^{n-1}{(-1)}^{m-1}m!S(n-1,m)+{(-1)}^n(n-1)!\\ &=(-1)\sum_{m=1}^{n-2}{(-1)}^{m-1}m!S(n-1,m)+\sum_{m=1}^{n-1}(-1)^{m-1}m!S(n-1,m)+(-1)^{n-1}(n-1)!\text{（这一条式子运用了$S(n,0)=0(n>0)$）}\\ &={(-1)}^{n-1-1}(n-1)!S(n-1,n-1)+{(-1)}^{n-1}(n-1)!\\ &=0 \end{align}$$
听说这个玩意儿叫做斯特林反演。
这就说明了最后我只会算入一个连通块的情况！
时间复杂度 $O(B_nn^2s)$。

---

代码实现

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>

using namespace std;

typedef long long LL;

const int MAXS=62;
const int N=12;

bool edge[MAXS][N][N];
int fact[N],bid[N];
char str[N*N];
int p[MAXS];
int n,s,cnt;
LL f[MAXS];
LL ans;

void ADD(LL vec)
{
    int cur=s-1;
    for (;cur>=0&&!((vec>>cur)&1);--cur);
    for (int i=0;i<cnt;++i)
        if (cur>p[i]) swap(f[i],vec),swap(cur,p[i]);
        else if (cur==p[i]) for (vec^=f[i];cur>=0&&!((vec>>cur)&1);--cur);
    if (cur>=0) f[cnt]=vec,p[cnt]=cur,++cnt;
}

LL calc()
{
    for (;cnt;--cnt) f[cnt-1]=p[cnt-1]=0;
    for (int i=1;i<n;++i)
        for (int j=i+1;j<=n;++j)
            if (bid[i]!=bid[j])
            {
                LL vec=0;
                for (int k=0;k<s;++k) vec|=(1ll<<k)*edge[k][i][j];
                ADD(vec);
            }
    return 1ll<<s-cnt;
}

void dfs(int x,int id)
{
    if (x==n+1)
    {
        ans+=(id&1?1:-1)*calc()*fact[id-1];
        return;
    }
    for (int i=1;i<=id;++i) bid[x]=i,dfs(x+1,id);
    bid[x]=++id,dfs(x+1,id);
}

int main()
{
    freopen("graph.in","r",stdin),freopen("graph.out","w",stdout);
    scanf("%d",&s),n=0;
    for (int i=0,l;i<s;++i)
    {
        scanf("%s",str),l=strlen(str);
        for (int j=1;!n;++j) if (j*(j-1)>>1==l) n=j;
        for (int j=1,x=0;j<n;++j)
            for (int k=j+1;k<=n;++k,++x)
                edge[i][j][k]=str[x]-'0';
    }
    fact[0]=1;
    for (int i=1;i<=n;++i) fact[i]=fact[i-1]*i;
    ans=0,dfs(1,0),printf("%lld\n",ans);
    fclose(stdin),fclose(stdout);
    return 0;
}