ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

Codeforces 1097 Alex and a TV Show

2019-01-17 09:00:55  阅读:264  来源: 互联网

标签:1097 Alex TV sum mu int maxn bitset 7005


传送门
除了操作 \(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;
}

标签:1097,Alex,TV,sum,mu,int,maxn,bitset,7005
来源: https://www.cnblogs.com/cjoieryl/p/10280400.html

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

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有