解题思路:
考虑容斥,枚举点的集合划分,强制两两集合间的点不连通,集合内的点任意连,若划分成m个集合,则最后至少有m个连通块。
而一个m的划分,在容斥时会被计算 ∑i=1mSim 次( S 为第二类斯特林数)。
所以列出容斥系数计算式:
打表找规律可得 fi=(−1)i−1(i−1)!
那么现在只用考虑若我们dfs出了一种划分,如何计算满足的方案数?
考虑把每张图看做一个二进制数,每一位代表一个属于不同集合的点对,要使所有这类点对不连通,那就是问该二进制数集合有多少子集异或和为0。
考虑线性基,那么先选出线性无关组后,加上任意剩余数组成的集合都有且仅有一个子集异或和为0,且互不相同。那么设线性无关组个数为cnt,方案数即为 2s−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;
}