【TC SRM 670】Gxor（斯特林反演，容斥原理）

题意
定义两张 n 个点带标号无向图的 xor：结果也是一张 n 个点
的带标号无向图：只包含恰好在一张图中出现的边。
给出 m 张图，问有多少方案选一个子集，xor 起来的图恰好
有 k 个连通块。

$O(Bell(n))$枚举$n$个点划分为几个块，求块间的边全部不存在的方案数，有$k$个块的方案数为$f(k)$。

那么考虑一个恰有$p$个联通块的方案$x$，对于$f(k)$，会计算$x$的次数是$S_2(p,k)$，也就是把$p$个联通块放入$k$个块中每个块不为空的方案数。

因为$x^n={\sum }_{i=0}^x{S}_2(n,i)x^{\underline{i}}$
这个等式是因为$x^n$的组合意义是把$n$个不同的球放到$x$个不同的盒子中。
我们可以枚举其中$i$个盒子非空即$\left(\genfrac{}{}{0ex}{}{x}{i}\right)=\frac{x^{\underline{i}}}{i!}$种方案。
然后把所有的球放入这$i$个不同的(所以要乘$i!$)盒子使得这$i$个盒子没有空的，也即$S_2(n,i)i!$
那么$x^n={\sum }_{i=0}^n{S}_2(n,i)x^{\underline{i}}$

又因为$x^{\overline{n}}={\sum }_{i=0}^n{S}_1(n,i)x^i$
证明，${\prod }_{i=0}^{n-1}(x+i)$中每个数选$x$，就代表第一类斯特林数中新建了一个圆排列，选$i$就代表插在了之前的$i$个数中的一个的前面。
所以$(-x{)}^{\overline{n}}={\sum }_{i=0}^n{S}_1(n,i)(-x{)}^i$
$(-x{)}^{\overline{n}}(-1{)}^n=x^{\underline{n}}={\sum }_{i=0}^n{S}_1(n,i)(-1{)}^{n-i}{x}^i$

所以$x^n={\sum }_{i=0}^n{S}_2(n,i){\sum }_{j=0}^i{S}_1(i,j)(-1{)}^{i-j}{x}^j$
$x^n={\sum }_{j=0}^n{x}^j{\sum }_{i=j}^n{S}_2(n,i)S_1(i,j)(-1{)}^{i-j}$
由此可得对于固定的$j$
${\sum }_{i=j}^n{S}_2(n,i)S_1(i,j)(-1{)}^{i-j}=[n=j]$

回到容斥，$f(k)$会计算恰有$p$个联通块的方案$S_2(k,p)$次，那么我们计算${\sum }_{i=1}^{n}f(i)S_1(i,1)(-1{)}^{i-1}$，根据斯特林反演，每个恰有$p$个联通块的方案的计算权和为${\sum }_{i=1}^p{S}_2(p,i)S_1(i,1)(-1{)}^{i-1}=[p=1]$
所以计算出$f(k)$即可。

$f(k)$的计算，对于一定不会出现的边可以列出方程，高斯消元求出自由元个数$r$，$f(k)=2^r$。

$\mathcal{A}CCode$

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

class Gxor{
	LL fac[10],a[55],b[55],ans;
	int T,n,c[10],lg[1<<10],dz[100] , cnt[100];
	void dfs(int u,int d){
		if(u == 0){
			
			int m = 0;
			for(int u=0,t=0;u<n;u++)
				for(int v=u+1;v<n;v++,t++)
					if(c[u] ^ c[v])
						a[m++] = b[t];
			a[m] = 0;
			int r = 0;
			
			for(int i=0;i<T;i++){
				if((a[r] >> i & 1) == 0)
					for(int j=r+1;j<m;j++)	
						if(a[j] >> i & 1){
							swap(a[r],a[j]);
							break;
						}
				for(int j=r+1;j<m;j++)
					if(a[j] >> i & 1)
						swap(a[r],a[j]) , a[j] ^= a[r];
				if(a[r]>>i&1) r++;
			}
			ans += fac[d] * (1ll << T >> r);
			return ;
		}
		for(int t=u;t;t=(t-1)&u) if(t & (u & -u)){
			for(int s=t;s;s&=(s-1))
				c[lg[s&-s]] = lg[u&-u];
			dfs(u-t,d+1);
		}
	}
	public:
	LL countsubs(vector<string>S){
		fac[1] = fac[0] = 1;
		for(int i=0;i<10;i++) lg[1<<i] = i; 
		for(int i=2;i<10;i++) fac[i] = fac[i-1] * (i-1);
		for(int i=2;i<10;i+=2) fac[i] = - fac[i];
		for(int i=1;i<=9;i++) dz[i * (i-1) / 2] = i;
		T = S.size() , n = dz[S[0].size()];
		for(int i=0;i<T;i++)
			for(int u=0,t=0;u<n;u++)
				for(int v=u+1;v<n;v++,t++)
					if(S[i][t] == '1')
						b[t] |= 1ll << i;
		dfs((1<<n)-1,0);
		return ans;
	}
};