题目大意
定义两个结点数相同的图
G1
与图
G2
的异或为一个新的图
G
,其中如果
现在给定
s
个结点数相同的图
题目分析
考虑容斥:我们用
Bn
的时间来枚举
n
个点的子集划分方案,然后限制不同子集的点一定不可以连边,同一子集任意。
这样我们就可以列出若干条
既然是容斥,我们就要想办法通过容斥系数来减去不合法的方案。
考虑如果一幅图实际上被划分成了
n
个连通块,它会被我计算多少次
其中 S 是第二类斯特林数。
考虑对于一个
首先显然 n=1 的时候这条式子等于 1 ,那么
∑m=1n(−1)m−1(m−1)!S(n,m)=∑m=1n−1(−1)m−1(m−1)![S(n−1,m−1)+mS(n−1,m)]+(−1)n(n−1)!=∑m=1n−1(−1)m−1(m−1)!S(n−1,m−1)+∑m=1n−1(−1)m−1m!S(n−1,m)+(−1)n(n−1)!=(−1)∑m=1n−2(−1)m−1m!S(n−1,m)+∑m=1n−1(−1)m−1m!S(n−1,m)+(−1)n−1(n−1)!(这一条式子运用了S(n,0)=0(n>0))=(−1)n−1−1(n−1)!S(n−1,n−1)+(−1)n−1(n−1)!=0
听说这个玩意儿叫做斯特林反演。
这就说明了最后我只会算入一个连通块的情况!
时间复杂度 O(Bnn2s) 。
代码实现
#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;
}