codeforces 547C 容斥原理
来源:互联网 发布:讯捷数据恢复 编辑:程序博客网 时间:2024/06/08 06:32
/*
题意:
第一行有两个数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的数的个数
上面的步骤就是容斥
*/
题意:
第一行有两个数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;}
1 0
- codeforces 547C 容斥原理
- Mike and Foam - CodeForces #547 C 容斥原理
- Codeforces Round #305 (Div. 1) C(容斥原理)
- 【CodeForces】651C - Watchmen(排序,容斥原理)
- cf 547C 容斥原理
- codeforces 630KIndivisibility(容斥原理)
- CodeForces-630 K. Indivisibility【容斥原理】
- 【CodeForces 630K】-Indivisibility 容斥原理
- 【CodeForces 630K】Indivisibility(容斥原理)
- Codeforces 547C Mike and Foam 容斥
- codeforces/#305 Div1/547/C Mike and Foam 【容斥】
- 【Codeforces】547C Mike and Foam 容斥
- *Codeforces Round #305 (Div. 1) C. Mike and Foam(容斥原理)
- Codeforces div2 #317C Lengthening Sticks(组合数学+ 容斥原理)
- 【CodeForces】559C - Gerald and Giant Chess(dp & 逆元 & 容斥原理)
- Codeforces #345 Div2 C Watchmen 容斥
- CF 547 C Mike and Foam(容斥原理)
- cf 341C (容斥原理)
- 黑马程序员-Java基础---集合
- canvas教程
- 七、Oracle分页
- HTTP状态码,超全的
- Codeforces 240E. Road Repairs 最小树形图+输出路径
- codeforces 547C 容斥原理
- 部分ubuntu系统安装时时区不对,导致系统总是在重启后时间会有相差
- 13 Fibonacci数
- php 不等待返回的实现方法(异步调用)
- Info.plist文件知多少(此处有表情)
- Nginx+Memcached+Tomcat集群配置实践
- 隐马尔科夫模型详解(一)
- 修改linux的文件时,如何快速找到要修改的内容
- Core Animation编程 Part1/2