![](https://img-blog.csdnimg.cn/20201014180756922.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
Mobius反演
文章平均质量分 87
DZYO
Never stop
展开
-
BZOJ:2440: [中山市选2011]完全平方数(容斥原理)
传送门题意:给你一个数k,求从1开始第k个非完全平方数。定义一个数为完全平方数当且仅当这个数质因数分解后存在二次项。 如: 4=224=2^2为完全平方数,30=2∗3∗530=2*3*5为非完全平方数。题解: 1.从1开始第k个非完全平方数⇔\Leftrightarrowmin(n)min(n),比nn小的非完全平方数有k个。有了这个性质,可以二分。 2.nn小的非完全平方数个数=∑i=1n原创 2017-07-23 16:59:01 · 503 阅读 · 0 评论 -
BZOJ2154,BZOJ2693 (Mobius反演)
BZOJ2154:Crash的数字表格 BZOJ2693:jzptab BZOJ 2154 给n,mn,m,求∑i=1n∑j=1n[lcm(i,j)]\sum\limits_{i=1}^n\sum\limits_{j=1}^n[lcm(i,j)] 将原式适当变形得:∑i=1n∑j=1n[lcm(i,j)]====∑i=1n∑j=1ni⋅jgcd(i,j)∑d∑i=1⌊nd⌋∑i=1⌊md⌋d2⋅原创 2017-07-24 21:21:00 · 383 阅读 · 1 评论 -
BZOJ2820: YY的GCD(Mobius反演)
传送门题意 给xx,yy,求1≤i≤x,1≤j≤y1 \le i \le x,1\le j \le y且gcd(i,j)gcd(i,j)为质数的(i,j)(i,j)对数。题解 不妨设n<mn<m。 首先有重要等式: ⌊⌊ab⌋c⌋=⌊abc⌋\lfloor\frac{\lfloor \frac{a}{b}\rfloor}{c}\rfloor=\lfloor \frac{a}{bc}\rflo原创 2017-07-23 19:49:57 · 327 阅读 · 0 评论 -
BZOJ2005:能量采集(Mobius反演)
传送门题意: 求 ∑i=1n∑j=1mgcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^m gcd(i,j)题解: MobiusMobius反演。首先,n=∑d|nϕ(d)n=\sum_{d|n}\phi(d)然后phi(d)phi(d)可线性筛出。∑i=1n∑j=1mgcd(i,j)=∑i=1n∑j=1m∑d|i,d|jϕ(d)=∑dϕ(d)⌊nd⌋⌊md⌋\sum_{i=1原创 2017-09-15 10:26:29 · 266 阅读 · 0 评论 -
杜教筛
如果一个积性函数,能够快速求得其狄利克雷卷积的前缀和,则可以在O(n23)O(n^{\frac{2}{3}})的时间内快速筛出其前缀和,这种筛法叫做杜教筛。比如:∑d|nμ(d)=[u=1]⇒∑i=1n∑d|iμ(d)=1\sum\limits_{d|n}\mu(d)=[u=1] \Rightarrow \sum\limits_{i=1}^n\sum\limits_{d|i}\mu(d)=1那么如果原创 2017-09-17 16:28:42 · 1456 阅读 · 0 评论 -
51Nod1355:斐波那契的最小公倍数 (min-max容斥+Mobius反演)
传送门 题解: 对于fib数列有gcd(i,j)=fgcd(i,j)gcd(i,j)=fgcd(i,j)\gcd(i,j) = f_{\gcd(i,j)}(可用归纳法证明)。 那么对于gcd(f{T})gcd(f{T})\gcd(f_{\{T\}}) 显然等于fgcd{T}fgcd{T}f_{\gcd\{T\}}。 怎么求lcm? 直接min-max对指数容斥即可,易得: lcmf{T}...原创 2018-05-09 21:28:23 · 737 阅读 · 1 评论 -
Codeforces 809E: Surprise me!(Mobius反演)
传送门 题解: 对于每个iii,处理出:fi=∑a∑b[(vala,valb)==i]φ(vala)φ(valb)dis(a,b)f_i=\sum_a\sum_b[(val_a,val_b)==i]\varphi(val_a)\varphi(val_b)dis(a,b)fi=a∑b∑[(vala,valb)==i]φ(vala)φ(valb)dis(a,b) 那么ans=∑ifiiφ...原创 2018-11-13 22:02:47 · 333 阅读 · 0 评论