标签：F1 F2 limits int Camera Security ABC220H subseteq oplus

https://atcoder.jp/contests/abc220/tasks/abc220_h

题解

考虑折半搜索，将 \(n\) 个点分为大小为 \(\dfrac{n}{2}\) 的两个集合 \(S, T\)

设 \(F1[s]\ (s\subseteq S)\) 表示如果选了 \(s\) 中的点安装摄像头，那么被监视的边的数量的奇偶性

设 \(F2[t]\ (t\subseteq T)\) 表示如果选了 \(t\) 中的点安装摄像头，那么被监视的两端都在 \(T\) 中的边的数量的奇偶性

设 \(G[s]\ (s\subseteq S)\) 表示所有 (和 \(\{S\) \ \(s\}\) 中的点之间有奇数条边的) 的 \(T\) 中的点的集合

以上都可以在 \(O(n*2^{\frac{n}{2}})\) 的时间内计算

那么一对 \(s,t\) 满足题目条件当且仅当：

\[F1[s] \oplus F2[t] \oplus \ (\operatorname{popcount}(G[s]\& T)\&1) = 0 \]

考虑如何计数

枚举 \(G[s]\& T\)，考虑计算

\[H[P][0/1] = \sum\limits_{G[s]\& T = P} [F1[s]\oplus F2[t] = 0/1] \]

\[Ans=\sum\limits_{P\subseteq T} H[P][\operatorname{popcount}(P)] \]

发现上面第一个式子长得像个集合交卷积

枚举 \(v1=0/1, v2=0/1\)，计算 \(C1[P]=\sum\limits_{G[s]=P} [F1[s]=v1], C2[P]=[F2[P]=v2]\)，

那么 \(H[P][v1\oplus v2] = \sum\limits_{s\& t = P} C1[s]*C2[t]\)

使用 \(FWT\) 在 \(O(n*2^{\frac{n}{2}})\) 的时间内计算

代码

#include <bits/stdc++.h>
#define N 50
#define M (1<<20)+5
#define pb push_back
#define lb(x) (x&-x)
using namespace std;
typedef long long ll;

int n, m, G[M], cnt[M]; 
bool F1[M], F2[M];
ll C1[M], C2[M];
vector<int> E[N];
inline int chk(int s, int x) { return (s>>x)&1; }

void FWT(ll *F, int _n, int tp) {
	for (int i = 1; i < (1<<_n); i <<= 1) {
		for (int j = 0; j < (1<<_n); j += i+i) for (int k = 0; k < i; k++) {
			F[j+k] += tp*F[i+j+k];  
		}
	}
}

int main() {
	scanf("%d %d", &n, &m);
	for (int i = 1, u, v; i <= m; i++) {
		scanf("%d %d", &u, &v); 
		E[u].pb(v); E[v].pb(u);
	}
	int t1 = n/2, t2 = n-t1;
	for (int s = 1; s < (1<<t1); s++) {
		int t = s^lb(s), x = log2(lb(s))+1;
		F1[s] = F1[t]; G[s] = G[t];
		for (auto y : E[x]) {
			if ((y<=t1&&!chk(s,y-1))||y>t1) F1[s]^=1;
			if (y > t1) G[s] ^= (1<<(y-t1-1));
		}
	}
	for (int s = 1; s < (1<<t2); s++) {
		int t = s^lb(s), x = log2(lb(s))+t1+1;
		F2[s] = F2[t]; cnt[s] = cnt[t]+1;
		for (auto y : E[x]) {
			if (y > t1 && !chk(s,y-t1-1)) F2[s]^=1;
		}
	}
	ll ans = 0; int mx = (1<<t1)-1;
	for (int v1 = 0; v1 <= 1; v1++) for (int v2 = 0; v2 <= 1; v2++) {
		memset(C1, 0, sizeof(C1)); memset(C2, 0, sizeof(C2));
		for (int s = 0; s < (1<<t1); s++) if (F1[s] == v1) ++C1[G[mx^s]];
		for (int s = 0; s < (1<<t2); s++) if (F2[s] == v2) ++C2[s];
		FWT(C1, t2, 1); FWT(C2, t2, 1);
		for (int s = 0; s < (1<<t2); s++) C1[s] = C1[s]*C2[s];
		FWT(C1, t2, -1);
		for (int s = 0; s < (1<<t2); s++) {
			if (!((cnt[s]&1)^v1^v2)) ans += C1[s];
		}
	}
	printf("%lld\n", ans);
	return 0;
}

标签：F1,F2,limits,int,Camera,Security,ABC220H,subseteq,oplus
来源： https://www.cnblogs.com/ak-dream/p/AK_Dream122.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。