[CSP-S模拟测试]:gcd(莫比乌斯反演)
题目描述
有$n$个正整数$x_1\sim x_n$,初始时状态均为未选。有$m$个操作,每个操作给定一个编号$i$,将$x_i$的选取状态取反。每次操作后,你需要求出选取的数中有多少个互质的无序数对。
输入格式
第一行两个整数$n,m$。第二行$n$个整数$x_1\sim x_n$。接下来$m$行每行一个整数。
输出格式
$m$行,每行一个整数表示答案。
样例
样例输入:
4 5
1 2 3 4
1
2
3
4
1
样例输出:
0
1
3
5
2
数据范围与提示
对于$20\%$的数据,$n,m\leqslant 1,000$。
对于另外$30\%$的数据,$x_i\leqslant 100$。
对于$100\%$的数据,$n,m\leqslant 200,000$,$x_i\leqslant 500,000$,$1\leqslant i\leqslant n$。
题解
我们先来设三个量:
$\alpha.s(i)$表示为$i$的倍数的数的个数。
$\beta.g(i)$表示 $gcd$为$i$的倍数的数个数。
$\gamma.f(i)$表示$gcd$为$i$的数的个数。
$s(i)$很好就能求出,而$g(i)=\frac{s(i)\times (s(i)-1))}{2}$,但是我们需要的是$f(i)$,该怎么办呢?
显然,$g(i)=\sum \limits_{i|d}f(d)$,那有又什么用呢?
这里就需要用到一个神奇的东东了:第二类莫比乌斯反演(详见信息学奥赛之数学一本通P145中间)。
于是这个式子便变成了:$f(i)=\sum \limits_{i|d}\mu(\frac{d}{i})g(d)$。
现在我们需要考虑的就只有修改操作了,每次插入或删除一个数的时候只要暴力枚举其因数即可。
时间复杂度:$\Theta(m\sqrt{\max x_i})$。
期望得分:$100$分。
实际得分:$100$分。
代码时刻
- #include<bits/stdc++.h>
- using namespace std;
- int n,m;
- long long a[200001];
- long long s[500001],g[500001],f[500001];
- long long mu[500001],prime[500001],cnt;
- bool vis[200001],v[500001];
- long long ans,mx;
- void pre_work()
- {
- mu[1]=1;
- for(int i=2;i<=mx;i++)
- {
- if(!v[i])mu[prime[cnt++]=i]=-1;
- for(int j=0;j<cnt&&i*prime[j]<=mx;j++)
- {
- v[i*prime[j]]=1;
- if(i%prime[j])mu[i*prime[j]]=-mu[i];
- else{mu[i*prime[j]]=0;break;}
- }
- }
- }
- int main()
- {
- scanf("%d%d",&n,&m);
- for(int i=1;i<=n;i++)
- scanf("%lld",&a[i]),mx=max(mx,a[i]);
- pre_work();
- while(m--)
- {
- int x,flag;
- scanf("%d",&x);
- flag=a[x];
- if(vis[x])
- {
- for(int i=1;i*i<=flag;i++)
- if(!(flag%i))
- {
- s[i]--;
- ans-=mu[i]*g[i];
- g[i]=s[i]*(s[i]-1)/2;
- ans+=mu[i]*g[i];
- if(flag/i!=i)
- {
- s[flag/i]--;
- ans-=mu[flag/i]*g[flag/i];
- g[flag/i]=s[flag/i]*(s[flag/i]-1)/2;
- ans+=mu[flag/i]*g[flag/i];
- }
- }
- vis[x]=0;
- }
- else
- {
- for(int i=1;i*i<=flag;i++)
- if(!(flag%i))
- {
- s[i]++;
- ans-=mu[i]*g[i];
- g[i]=s[i]*(s[i]-1)/2;
- ans+=mu[i]*g[i];
- if(flag/i!=i)
- {
- s[flag/i]++;
- ans-=mu[flag/i]*g[flag/i];
- g[flag/i]=s[flag/i]*(s[flag/i]-1)/2;
- ans+=mu[flag/i]*g[flag/i];
- }
- }
- vis[x]=1;
- }
- printf("%lld\n",ans);
- }
- return 0;
- }
rp++
[CSP-S模拟测试]:gcd(莫比乌斯反演)的更多相关文章
- [BZOJ 2820] YY的gcd(莫比乌斯反演+数论分块)
[BZOJ 2820] YY的gcd(莫比乌斯反演+数论分块) 题面 给定N, M,求\(1\leq x\leq N, 1\leq y\leq M\)且gcd(x, y)为质数的(x, y)有多少对. ...
- HDU1695 GCD(莫比乌斯反演)
传送门 看了1个多小时,终于懂了一点了 题目大意:给n,m,k.求gcd(x,y) = k(1<=x<=n, 1<=y<=m)的个数 思路:令F(i)表示i|gcd(x,y)的 ...
- hdu 1695 GCD 莫比乌斯反演入门
GCD 题意:输入5个数a,b,c,d,k;(a = c = 1, 0 < b,d,k <= 100000);问有多少对a <= p <= b, c <= q <= ...
- 洛谷P2257 YY的GCD 莫比乌斯反演
原题链接 差不多算自己推出来的第一道题QwQ 题目大意 \(T\)组询问,每次问你\(1\leqslant x\leqslant N\),\(1\leqslant y\leqslant M\)中有多少 ...
- HYSBZ - 2818 Gcd (莫比乌斯反演)
莫比乌斯反演的入门题,设 \(F(x): gcd(i,j)\%x=0\) 的对数,\(f(x): gcd(i,j)=x\)的对数. 易知\[F(p) = \lfloor \frac{n}{p} \rf ...
- 【BZOJ2818】Gcd [莫比乌斯反演]
Gcd Time Limit: 10 Sec Memory Limit: 256 MB[Submit][Status][Discuss] Description 给定整数N,求1<=x,y&l ...
- Luogu P2257 YY的GCD 莫比乌斯反演
第一道莫比乌斯反演...$qwq$ 设$f(d)=\sum_{i=1}^n\sum_{j=1}^m[gcd(i,j)==d]$ $F(n)=\sum_{n|d}f(d)=\lfloor \frac{N ...
- BZOJ 2818 Gcd (莫比乌斯反演 或 欧拉函数)
2818: Gcd Time Limit: 10 Sec Memory Limit: 256 MB Submit: 2534 Solved: 1129 [Submit][Status][Discu ...
- BZOJ 2820 luogu 2257 yy的gcd (莫比乌斯反演)
题目大意:求$gcd(i,j)==k,i\in[1,n],j\in[1,m] ,k\in prime,n,m<=10^{7}$的有序数对个数,不超过10^{4}次询问 莫比乌斯反演入门题 为方便 ...
随机推荐
- poj2236Wireless Network
Description An earthquake takes place in Southeast Asia. The ACM (Asia Cooperated Medical team) have ...
- Chrome开启多线程下载
Chrome多线程下载也和标签页预览一样属于Google测试中的功能,可通过在地址栏输入chrome://flags/,然后在搜索框中输入Parallel downloading,选择enabled, ...
- word中迅速将表格一分为二 拆分表格快捷键ctrl+shift+enter 重复上一个命令快捷键f4
这里说的是将一个表格拆分为两个表格 选择要拆分的行,快捷键ctrl+shift+enter,就拆分为两个表格了,是不是很快! 在多个表格需要拆分的时候,做一次这样的操作,然后不停的移动.F4,就可以了 ...
- 跨域(SpringBoot)
概念: 在 HTML 中,<a>, <form>, <img>, <script>, <iframe>, <link>等标签以及 ...
- 移动端web整理 移动端问题总结,移动web遇到的那些坑
meta基础知识 H5页面窗口自动调整到设备宽度,并禁止用户缩放页面 忽略将页面中的数字识别为电话号码 忽略Android平台中对邮箱地址的识别 当网站添加到主屏幕快速启动方式,可隐藏地址栏,仅针对i ...
- 那些年很脑残的bugs
1.老师给了前端界面,我们用java写后台. 我改了表单form的action属性,让它跳到自己写的servlet上面去.自己在servlet里面对数据库一顿操作猛如虎,然后让servlet跳回原来页 ...
- HDU Stealing Harry Potter's Precious(状压BFS)
状压BFS 注意在用二维字符数组时,要把空格.换行处理好. #include<stdio.h> #include<algorithm> #include<string.h ...
- mybatis where 中in的使用
当我们使用mybatis时,在where中会用到 in 如: where name in ('Jana','Tom'); 我们可以在sql中直接写 name in ('Jana','Tom') 或者 ...
- redis限流器的设计
1.定义注解 import java.lang.annotation.Documented; import java.lang.annotation.ElementType; import java. ...
- protocols - 协议定义文件
描述 该文件为普通 ASCII 文件,它描述了 TCP/IP 子系统中各类 DARPA internet 协议族.你应该参考这个文件, 而不是使用 ARPA 的 include 文件中的号码, 更不必 ...