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

# 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;
}