codeforces 547C 容斥原理

原创 2015年07月08日 17:44:34
/*
题意:
第一行有两个数n,q表示由n个数,q次询问
第二行有n个数a1,a2.....
接下来q行表示询问,每行一个数i表示ai在集合中是否存在(集合最初是空的)如果存在把
这个去掉以后问剩下的数两两互质的数有多少个。若果不存在把这个数加入集合,问两两互
质的有多少个。
注意:如果判断这个数在集合中存在的标志是这个数是否存在且下标是否一致,如果数存在但
下标不一致也应该把这个数加入到集合里头


思路:声明一个数x如果它不存在y>1且y的平方不是x的因数它就是好数
x=p1的a1次方*p2的a2次方*.....f(x)=a1+a2+a3+...换句话说a1,a2,...他们的值不是0就是1
然后声明d(x)为一个所有的数都有一个因数为x的集合。如果前面的集合中两两互质的已经算出
来了为ans,那么在加一个数,该是多少呢。ans=ans+d(1)-d(p1)-d(p2)-...d(pn)+d(p1p2)+
d(p1p3)+...-d(p1p2p3)-d(p1p2p4)....
d(1)的值应该为前面出现数的个数
d(p1)的值应该为前面出现有因数为P1的数的个数




上面的步骤就是容斥
*/
#include<bits/stdc++.h>
using namespace std;
const int maxn = 500005;
int a[200005],mul[80],cnt[maxn],has[200005];
vector<int> p[maxn];
long long ans;
void add(int val,int sign)
{
    for(int mask=0;mask<(1<<p[val].size());mask++)
    {
        if(!mask)
            mul[mask]=1;
        else
            mul[mask]=mul[mask&(mask-1)]*p[val][__builtin_ctz(mask)];
        if(sign<0)
            cnt[mul[mask]]--;
        if(__builtin_popcount(mask)&1)
            ans-=cnt[mul[mask]]*sign;
        else
            ans+=cnt[mul[mask]]*sign;
        if(sign>0)
            cnt[mul[mask]]++;
    }
}
int main()
{
    ios::sync_with_stdio(false);
    for(int i=2;i<maxn;i++)
        if(p[i].empty())
            for(int j=i;j<maxn;j+=i)
                p[j].push_back(i);
    int n,q;
    cin >> n >> q;
    for(int i=0;i<n;i++)
    {
        cin >> a[i];
        has[i]=1;
    }
    while(q--)
    {
        int x;
        cin >> x;
        x--;
        add(a[x],has[x]);
        has[x]*=-1;
        cout << ans << "\n";
    }
    return 0;
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

CodeForces 547C. Mike and Foam 莫比乌斯反演

用莫比乌斯反演定理,维护一个动态集合含有的互质的数对的个数

codeforces #305 547C C. Mike and Foam(莫比乌斯反演)

题目链接: 点击打开链接 题目大意: 给出一列数,最开始集合 为空,然后q次操作,每次给出一个x,如果第x个数存在,那么删去,不存在添加,问操作完互质的数有多少对 题目分析: 问互质的数的对...

*Codeforces Round #305 (Div. 1) C. Mike and Foam(容斥原理)

题目: http://codeforces.com/problemset/problem/547/C 题意: n种啤酒,每种啤酒有相应的毫升数,q个询问。一开始的酒架是空的,输入啤酒标号x,...

Codeforces Round #305 (Div. 1) C(容斥原理)

C. Mike and Foam time limit per test 2 seconds memory limit per test 256 megabytes ...

547B. Mike and Feet(codeforces Round #305)

B. Mike and Feet time limit per test 1 second memory limit per test 256 megabytes i...
  • caduca
  • caduca
  • 2015-06-19 15:32
  • 4303

codeforces #547D Mike and Fish 欧拉回路

题目大意:给定平面上的n个点,要求将每个点染成红色/蓝色,使得每行/每列的红色点数和蓝色点数之差≤1\leq 1将每一个横坐标/纵坐标看做一个点,每个点看做一条连接两个坐标的边 现在我们要将每条边染...

codeforces 547B【单调栈】

题意: 有一个长度为n的序列,序列有长度为1...n的连续子序列, 一个连续子序列里面最小的值称作这个子序列的子序列的strength, 要求出每种长度的连续子序列的最大的strength。 ...

547A . Mike and Frog(codeforces Round 305)

A. Mike and Frog time limit per test 1 second memory limit per test 256 megabytes i...
  • caduca
  • caduca
  • 2015-06-18 10:19
  • 4014

CodeForces 547B - Mike and Feet(单调栈)

题意:给定一个长度为 n (1 ≤ n ≤ 2 × 105) 的序列,求某一固定长度的连续子序列(长度为1 ~ n)中最小元素的最大值。 对于某一元素,找出它为最小值的最大连续长度区间。...

CodeForces 547E Mike and Friends [后缀树组+主席树]

给出N个串,M个询问,对于每个询问L,R,K,给出串L到串R的所有串中串K作为子串出现了几次。
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)