[Codeforces 1097F]Alex and a TV Show（bitset + 莫比乌斯反演）

最新推荐文章于 2021-07-24 22:00:47 发布

xyz32768

最新推荐文章于 2021-07-24 22:00:47 发布

阅读量219

点赞数

分类专栏： CF、AT 文章标签： bitset 莫比乌斯反演

本文链接：https://blog.csdn.net/xyz32768/article/details/88363511

版权

CF、AT 专栏收录该内容

18 篇文章 0 订阅

订阅专栏

Address

Solution

发现我们要查询的是某集合内某数的出现次数模 $2$ 的结果
可以对于每个可重集维护一个 bitset ，第 $i$ 个集合的 bitset 第 $j$ 位表示第 $i$ 个集合中 $j$ 的出现次数的奇偶性
但我们发现这样做无法解决第 $3$ 种操作
我们不妨把维护的东西换一下
$A [i] [j]$ 表示第 $i$ 个集合内 $j$ 的倍数的个数的奇偶性
对于 $1$ 操作，直接将所有 $v$ 的约数 $k$ 令 $A [v] [k] = 1$ ，其他的 $k$ 把 $A [v] [k]$ 改成 $0$ 即可
对于 $2$ 操作显然还是 $A[x]\leftarrow A[y]\bigoplus A[z]$ ，其中 $\bigoplus$ 为异或运算
对于 $3$ 操作，显然 $x|\gcd(a,b)$ 当且仅当 $x ∣ a$ 且 $x ∣ b$ 。于是 $A[x]\leftarrow A[y]\& A[z]$ ， $\&$ 为与运算
而对于 $4$ ，我们如何在知道 $A [x]$ 的情况下知道第 $x$ 个集合中 $v$ 的个数的奇偶性呢？
对于 $A [x]$ ，进行莫比乌斯反演，得
$ans=\sum_{v|k}A[x][k]\mu(\frac kv)$
由于我们只关心奇偶性，所以我们可以把 $\mu$ 的值看成只有 $0$ 和 $1$ 两种
再开 bitset $\mu[v][k]$ 表示：如果 $v ∣ k$ 则 $\mu[v][k]=|\mu(\frac kv)|$ ，否则 $\mu[v][k]=0$
这样询问结果就是 $A[x]\& mu[v]$ 中 $1$ 的个数的奇偶性
复杂度 $O(\frac {(n+q+v)v}{64})$

Code

#include <bits/stdc++.h>

inline int read()
{
	int res = 0; bool bo = 0; char c;
	while (((c = getchar()) < '0' || c > '9') && c != '-');
	if (c == '-') bo = 1; else res = c - 48;
	while ((c = getchar()) >= '0' && c <= '9')
		res = (res << 3) + (res << 1) + (c - 48);
	return bo ? ~res + 1 : res;
}

const int N = 1e5 + 5, M = 7005;

int n, q, miu[M], tot, pri[M];
bool mark[M];

std::bitset<M> times[N], fy[M];

void sieve()
{
	mark[0] = mark[miu[1] = 1] = 1;
	for (int i = 2; i <= 7000; i++)
	{
		if (!mark[i]) miu[pri[++tot] = i] = 1;
		for (int j = 1; j <= tot; j++)
		{
			if (i * pri[j] > 7000) break;
			mark[i * pri[j]] = 1;
			if (i % pri[j] == 0) break;
			else miu[i * pri[j]] = miu[i];
		}
	}
	for (int i = 1; i <= 7000; i++)
		for (int j = i; j <= 7000; j += i)
			fy[i][j] = miu[j / i];
}

int main()
{
	n = read(); q = read();
	sieve();
	int op, x, y, z;
	while (q--)
	{
		op = read(); x = read(); y = read();
		if (op == 1)
		{
			int S = sqrt(y);
			times[x] = 0;
			for (int i = 1; i <= S; i++)
			{
				if (y % i) continue;
				times[x][i] = times[x][y / i] = 1;
			}
		}
		else if (op == 2) z = read(), times[x] = times[y] ^ times[z];
		else if (op == 3) z = read(), times[x] = times[y] & times[z];
		else printf("%d", (times[x] & fy[y]).count() & 1);
	}
	return puts(""), 0;
}