Codeforces 585E. Present for Vitalik the Philatelist(容斥)
好题!学习了好多
写法①:
先求出gcd不为1的集合的数量,显然我们可以从大到小枚举计算每种gcd的方案(其实也是容斥),或者可以直接枚举gcd然后容斥(比如最大值是6就用2^cnt[2]-1+3^cnt[3]-1-(6^cnt[6]-1),cnt[x]表示x的倍数的个数),用容斥计算的话可以发现系数是莫比乌斯函数的相反数,就可以线性筛了。下面会记录一种O(MAX*ln(MAX))的筛法...求cnt的话可以选择直接枚举倍数计算O(MAX*ln(MAX))或者分解质因数,因为1e7内最多有8个不同质因子,求出所有质因子再枚举子集给cnt加上贡献。
然后枚举每一个数计算它对答案的贡献,就是gcd不为1的集合的数量里去掉这个数的因数的贡献,分解质因数之后一个一个去掉即可。
这种写法就不写代码啦~(其实写了但是写炸了懒得调了
O(MAX*ln(MAX))莫比乌斯函数筛:
miu[]=;
for(int i=;i<=mx;i++)
for(int j=i<<;j<=mx;j+=i) miu[j]-=miu[i];
写法②:
这是一种更简单的写法,尝试把上面统计数量后计算每个数的贡献这两个过程合二为一。
考虑去掉一个数的因数对答案的贡献怎么更简单地做,可以发现一个数的所有因数的cnt里都有这个数的贡献,所以我们实际上可以在统计答案的时候直接计算所有数的贡献顺便去掉这个数的贡献,也就是以下的公式:
中间部分是gcd为x的时候的方案数,右边是gcd为x时对此式子有贡献的数的数量。
我们求莫比乌斯函数的时候用上上面的O(MAX*ln(MAX))版筛法就可以把代码长度做到极短了。
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#define MOD(x) ((x)>=mod?(x)-mod:(x))
using namespace std;
const int maxn=, maxm=1e7+, inf=1e9+, mod=1e9+;
int n, x, mx, ans;
int mi[maxn], cnt[maxm], miu[maxm];
void read(int &k)
{
int f=; k=; char c=getchar();
while(c<'' || c>'') c=='-' && (f=-), c=getchar();
while(c<='' && c>='') k=k*+c-'', c=getchar();
k*=f;
}
int main()
{
read(n); mi[]=; for(int i=;i<=n;i++) mi[i]=(1ll*mi[i-]<<)%mod;
for(int i=;i<=n;i++) read(x), cnt[x]++, mx=max(mx, x);
miu[]=;
for(int i=;i<=mx;i++)
{
for(int j=i<<;j<=mx;j+=i) miu[j]-=miu[i], cnt[i]+=cnt[j];
if(!cnt[i] || !miu[i]) continue;
int delta=1ll*miu[i]*(mi[cnt[i]]-)%mod*(cnt[i]-n)%mod;
delta=MOD(delta+mod); ans=MOD(ans+delta);
}
printf("%d\n", ans);
}
写法③:
是在写法②的基础上的,可以省去求莫比乌斯函数的过程。
可以倒着枚举数i,一般可以较为简单的计算i的倍数的总贡献,然后删去i的倍数(此时不包括i)的贡献(此时倍数的贡献已经计算好了)来得到i的贡献。注意!!!这里i的答案指的是如果计算i的话答案是多少,但实际上可能是需要减去的,它会被它的因数减去而得到它的因数的答案。
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#define MOD(x) (x>=mod?(x)-mod:(x))
using namespace std;
const int maxn=, inf=1e9, maxm=1e7+, mod=1e9+;
int n, x, mx, anss;
int ans[maxm], mi[maxn], cnt[maxm];
inline void read(int &k)
{
int f=; k=; char c=getchar();
while(c<'' || c>'') c=='-'&&(f=-), c=getchar();
while(c<='' && c>='') k=k*+c-'', c=getchar();
k*=f;
}
int main()
{
read(n);
for(int i=;i<=n;i++) read(x), cnt[x]++, mx=max(mx, x);
mi[]=; for(int i=;i<=n;i++) mi[i]=(1ll*mi[i-]<<)%mod;
for(int i=mx;i>=;i--)
{
int sum=cnt[i];
for(int j=i<<;j<=mx;j+=i) sum+=cnt[j], ans[i]=MOD(ans[i]-ans[j]+mod);
ans[i]+=1ll*(mi[sum]-)*(n-sum)%mod; ans[i]=MOD(ans[i]);
anss=MOD(anss+ans[i]);
}
printf("%d\n", anss);
}
Codeforces 585E. Present for Vitalik the Philatelist(容斥)的更多相关文章
- Codeforces 585E - Present for Vitalik the Philatelist(简单莫反+狄利克雷前缀和)
Codeforces 题目传送门 & 洛谷题目传送门 一道不算太难的 D1E 罢--虽然我不会做/kk u1s1 似乎这场 Div1 挺水的?F 就是个 AC 自动机板子还被评到了 3k2-- ...
- 【CF 585E】 E. Present for Vitalik the Philatelist
E. Present for Vitalik the Philatelist time limit per test 5 seconds memory limit per test 256 megab ...
- 【CodeForces】585 E. Present for Vitalik the Philatelist
[题目]E. Present for Vitalik the Philatelist [题意]给定n个数字,定义一种合法方案为选择一个数字Aa,选择另外一些数字Abi,令g=gcd(Ab1...Abx ...
- CF585E. Present for Vitalik the Philatelist [容斥原理 !]
CF585E. Present for Vitalik the Philatelist 题意:\(n \le 5*10^5\) 数列 \(2 \le a_i \le 10^7\),对于每个数\(a\) ...
- 「CF585E」 Present for Vitalik the Philatelist
「CF585E」 Present for Vitalik the Philatelist 传送门 我们可以考虑枚举 \(S'=S\cup\{x\}\),那么显然有 \(\gcd\{S'\}=1\). ...
- CF 585 E Present for Vitalik the Philatelist
CF 585 E Present for Vitalik the Philatelist 我们假设 $ f(x) $ 表示与 $ x $ 互质的数的个数,$ s(x) $ 为 gcd 为 $ x $ ...
- CF585E:Present for Vitalik the Philatelist
n<=500000个2<=Ai<=1e7的数,求这样选数的方案数:先从其中挑出一个gcd不为1的集合,然后再选一个不属于该集合,且与该集合内任意一个数互质的数. 好的统计题. 其实就 ...
- Codeforces Round #345 (Div. 1) A - Watchmen 容斥
C. Watchmen 题目连接: http://www.codeforces.com/contest/651/problem/C Description Watchmen are in a dang ...
- CodeForces 559C Gerald and Gia (格路+容斥+DP)
CodeForces 559C Gerald and Gia 大致题意:有一个 \(N\times M\) 的网格,其中有些格子是黑色的,现在需要求出从左上角到右下角不经过黑色格子的方案数(模 \(1 ...
随机推荐
- PHP序列化serialize()和反序列化unserialize()
所谓的序列化,就是把保存在内存中的各种对象状态或属性保存起来,在需要时可以还原出来. serialize() 可处理除了 resource 之外的任何类型返回字符串,此字符串包含了表示 value 的 ...
- Linux建立互信关系(ssh公钥登录)
Linux有多种登录方式,比如telnet.ssh.支持ssh登录方式:口令登录和公钥登录 ssh登录方式:ssh [-l login_name] [-p port] [user@]hostname ...
- Python 日志记录与程序流追踪(基础篇)
日志记录(Logging) More than print: 每次用 terminal debug 时都要手动在各种可能出现 bug 的地方 print 相关信息来确认 bug 的位置: 每次完成 d ...
- windows下对python的pip更新到最新版本
1->打开windows的命令窗口. 2->进入到pip.exe所在的文件夹下,我安装的python在G:\python3.6文件夹下,pip.exe则在G:\python3.6\Scri ...
- 如何更改Arcmap里经纬度小数点后面的位数?
customize>arcmap option>data view >round coordinate to 改成想要显示的小数位数
- Linux sync命令的作用分析
Sync命令 在用reboot命令启动unix系统后,系统提示出错信息,部分应用程序不能正常工作.经仔细检查系统文件,并和初始的正确备份进行比较,发现某些文件确实被破坏了,翻来覆去找不到文件遭破坏 ...
- 单源最短路——Bellman-Ford算法
1.Dijkstra的局限性 Dijkstra算法是处理单源最短路径的有效算法,但它局限于边的权值非负的情况,若图中出现权值为负的边,Dijkstra算法就会失效,求出的最短路径就可能是错的. 列如以 ...
- Java 数组转字符
public static String toString(int[] arr){ String temp = ""; for(int i = 0;i<arr.length; ...
- [pascal入门]数组
一.本节目标 本节我们将要讲述数组.本节目标: 一维数组 二维数组 字符数组 二.一维数组 我们通过一个案例来简单的理解数组.班主任要计算班级里面50个同学数学成绩的平均成绩,道理上讲这是一个比较简单 ...
- WCF RestFull提交数据超出限额解决方法
最近在使用wcf restfull时出现了超大数据提交出错的问题. 服务端会返回错误:服务器处理请求时遇到错误.有关构造有效服务请求的内容,请参阅服务帮助页.异常消息为“反序列化对象 属于类型 Yes ...