Codeforces 1097 Alex and a TV Show

Cyhlnj

于 2019-01-17 08:51:06 发布

阅读量341

点赞数

分类专栏：反演与容斥 Bitset

本文链接：https://blog.csdn.net/oi_Konnyaku/article/details/86517620

版权

反演与容斥同时被 2 个专栏收录

37 篇文章 0 订阅

订阅专栏

Bitset

3 篇文章 0 订阅

订阅专栏

传送门
除了操作 $3$ 都可以 $b i t s e t$
现在要维护
$C_i=\sum_{gcd(j,k)=i}A_jB_k$
类比 $F W T$ ，只要求出 $A'_i=\sum_{i|d}A_d$
就可以直接按位相乘了
求答案就是莫比乌斯反演， $A_i=\sum_{i|d}\mu(\frac{d}{i})A'_i$
把每个数字的 $\mu$ 的 $b i t s e t$ 预处理出来，乘法就是 $a n d$
最后用 $c o u n t$ 统计答案

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

const int maxn(7005);

int mu[maxn], n, q, pr[maxn], tot;
bitset <maxn> bc[100005], bcm[7005], bcv[7005], ispr;

int main() {
    int i, j, op, l, r, v;
    mu[1] = 1;
    for (i = 2; i <= 7000; ++i) {
        if (!ispr[i]) pr[++tot] = i, mu[i] = -1;
        for (j = 1; j <= tot && pr[j] * i <= 7000; ++j) {
            ispr[pr[j] * i] = 1;
            if (i % pr[j]) mu[i * pr[j]] = -mu[i];
            else {
                mu[i * pr[j]] = 0;
                break;
            }
        }
    }
    for (i = 1; i <= 7000; ++i)
        for (j = i; j <= 7000; j += i) {
            bcv[j].set(i);
            if (mu[j / i]) bcm[i].set(j);
        }
    scanf("%d%d", &n, &q);
    for (i = 1; i <= q; ++i) {
        scanf("%d%d%d", &op, &l, &r);
        if (op == 1) bc[l] = bcv[r];
        else if (op == 2) scanf("%d", &v), bc[l] = bc[r] ^ bc[v];
        else if (op == 3) scanf("%d", &v), bc[l] = bc[r] & bc[v];
        else putchar(((bc[l] & bcm[r]).count() & 1) + '0');
    }
    return 0;
}