BZOJ4671:异或图

传送门
直接求连通的不好做,考虑容斥
g i g_i gi 表示至少有 i i i 个连通块的方案数, f i f_i fi 表示恰好有 i i i 个的
那么
g x = ∑ i = x n { x i } f i    ⟺    f x = ∑ i = x n ( − 1 ) i − x [ x i ] g i g_x=\sum_{i=x}^{n}\begin{Bmatrix}x \\ i\end{Bmatrix}f_i\iff f_x=\sum_{i=x}^{n}(-1)^{i-x}\begin{bmatrix}x \\ i\end{bmatrix}g_i gx=i=xn{xi}fifx=i=xn(1)ix[xi]gi
那么
f 1 = ∑ i = 1 n ( − 1 ) i − 1 ( i − 1 ) ! g i f_1=\sum_{i=1}^{n}(-1)^{i-1}(i-1)!g_i f1=i=1n(1)i1(i1)!gi
g g g
考虑枚举点的拆分,相当于是不同的集合之没有边,这部分直接用线性基求出方案

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

int n = 1, graph[65][15][15], m, id[15];
char ch[2333];
ll fac[15], ans, bc[65], v;

void Dfs(int x, int f) {
	register int i, j, k, tot, num;
	if (x > n) {
		memset(bc, 0, sizeof(bc)), num = 0;
		for (i = 1; i <= m; ++i) {
			for (v = tot = 0, j = 1; j <= n; ++j)
				for (k = j + 1; k <= n; ++k)
					if (id[j] != id[k]) v |= (ll)graph[i][j][k] << tot, ++tot;
			for (j = 0; j < tot; ++j)
				if (v >> j & 1) {
					if (!bc[j]) {
						bc[j] = v, ++num;
						break;
					}
					v ^= bc[j];
				}
		}
		ans += (ll)((f & 1) ? 1 : -1) * fac[f - 1] * (1ll << (m - num));
		return;
	}
	for (i = 1; i <= f + 1; ++i) id[x] = i, Dfs(x + 1, max(i, f));
}

int main() {
	register int i, j, k, len, cnt;
	for (scanf("%d", &m), i = 1; i <= m; ++i) {
		scanf(" %s", ch + 1), len = strlen(ch + 1);
		while (n * (n - 1) / 2 < len) ++n;
		for (cnt = 0, j = 1; j <= n; ++j)
			for (k = j + 1; k <= n; ++k) graph[i][j][k] = ch[++cnt] - '0';
	}
	for (fac[0] = 1, i = 1; i <= n; ++i) fac[i] = fac[i - 1] * i;
	Dfs(1, 0), printf("%lld\n", ans);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值