莫比乌斯反演
yuanxinyu402
这个作者很懒,什么都没留下…
展开
-
HYSBZ/BZOJ 2154 Crash的数字表格 - 莫比乌斯反演
题目描述由于不想码推导过程,粘一个by Liu Junhao最开始推到最后的式子,却傻傻的没有意识到可以分块加速,还以为自己推错了。。=_= 表示以后看到向下取整的东东要注意了。 注意:Mod要写就写完,不要懒,谁知道会不会爆??。#include<cstdio> #include<algorithm> using namespace std; #define MAXN 10000000 #def原创 2016-01-27 21:18:55 · 468 阅读 · 0 评论 -
SPOJ VLATTICE Visible Lattice Points - 莫比乌斯反演
题目描述一开始把坐标想成(1,1,1)开始,却在算斜率的时候忘记-(1,1,1),至此陷入歧途。最后忍无可忍看了题解。。。看的这个,就粘这个了想懂了还是很清楚的。#include<cstdio> #define MAXN 1000000 #define MAXP 80000 int mu[MAXN+10],cntpr,prime[MAXP+10],sum[MAXN+10],a; bool ispri原创 2016-01-27 16:03:31 · 288 阅读 · 0 评论 -
SPOJ - PGCD Primes in GCD Table - 莫比乌斯反演
题目描述1’.枚举prime#include<cstdio> #include<algorithm> using namespace std; #define MAXN 10000000 #define MAXP 800000int prime[MAXP+10],cntpr,mu[MAXN+10],a,b,sum[MAXN+10]; bool isprime[MAXN+10];void CalMob原创 2016-01-27 14:32:15 · 285 阅读 · 0 评论 -
BZOJ 2301: [HAOI2011]Problem b - 莫比乌斯反演
题目描述经过这道题的洗礼后,现在这个题就是差不多的了。 注意容斥原理,向上向下取整的问题。=_=#include<cstdio> #include<cstring> #include<algorithm> using namespace std; #define MAXN 50000 int mu[MAXN+10],prime[MAXN+10],cntpr,sum[MAXN+10],a,b,c,d原创 2016-01-27 12:22:28 · 277 阅读 · 0 评论 -
BZOJ(本校) 2525 公约数 - 莫比乌斯反演
题目描述 给出一个数N,求1 <= x, y <= N,且gcd(x, y)为素数的数对x, y的数量。 输入 一行一个数N。 输出 一个数表示答案 样例输入 4 样例输出 4 提示 【数据范围】 N ≤ 10000000 来源 HZOI gcd(x,y)为质数,1<=x,y<=n <=> gcd(x,y)=1 1<=x,y<=n/prime[i] prime[i]需原创 2016-01-27 10:22:59 · 510 阅读 · 0 评论 -
BZOJ 2440 完全平方数-莫比乌斯函数
题目描述 1.n的不大确定+f(n)([1,n]中不是平方数的倍数的数的个数)递增 -> 二分 2.容斥原理:f(x)=1的倍数-一个质数的倍数+2个质数乘积的倍数-… ; f(x)=sigma(mu[i]*floor(x/(i*i))) ,1<=i<=floor(sqrt(x)) #include<cstdio> #include<cstring> #include<cmath原创 2016-01-26 19:56:41 · 323 阅读 · 0 评论