Peter的莫比乌斯反演与各种筛法题单

前置芝士1——数论分块

UVA11526 H(n)

P2261 [CQOI2007]余数求和

P2260 [清华集训2012]模积和

P3935 Calculating

前置芝士2——线性筛

P2158 [SDOI2008]仪仗队

SP526 DIV - Divisors

P4626 一道水题 II

SP22461 SMALL - Smallest Number

CF1017F The Neutral Zone论如何选择筛法

以及你的算法优化技巧与数学能力

P6810 「MCOI-02」Convex Hull 凸包

P5495 Dirichlet 前缀和

P6788 「EZEC-3」四月樱花

莫比乌斯反演

就是把一个看起来只能暴力算的柿子化成一个可以一下子算出来或者数论分块可以算出来的东西

以下都默认 n\le m

形式1

\sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)=1]=\sum_{i=1}^n\sum_{j=1}^m\sum_{x|i}\sum_{x|j}\mu(x) =\sum_{x=1}^n\mu(x)\sum_{i=1}^n\sum_{j=1}^m[x|i][x|j] =\sum_{x=1}^n\mu(x)\lfloor {n\over x}\rfloor\lfloor {m\over x}\rfloor

数论分块即可 O(n) 预处理, O(\sqrt n) 单次询问

形式2

\sum_{i=1}^n\sum_{j=1}^m\gcd(i,j)=\sum_{i=1}^n\sum_{j=1}^m\sum_{x|i}\sum_{x|j}\phi(x) =\sum_{x=1}^n\phi(x)\sum_{i=1}^n\sum_{j=1}^m[x|i][x|j] =\sum_{x=1}^n\phi(x)\lfloor {n\over x}\rfloor\lfloor {m\over x}\rfloor

数论分块即可 O(n) 预处理, O(\sqrt n) 单次询问

形式3

\sum_{i=1}^n\sum_{j=1}^mf(\gcd(i,j))=\sum_{d=1}^nf(d)\sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)=d] =\sum_{d=1}^nf(d)\sum_{i=1}^{\lfloor{n\over d}\rfloor}\sum_{j=1}^{\lfloor{m\over d}\rfloor}[\gcd(i,j)=1] =\sum_{d=1}^nf(d)\sum_{k=1}^{\lfloor{n\over d}\rfloor}\mu(k)\lfloor {n\over dk}\rfloor\lfloor {m\over dk}\rfloor

然后是一个重要的思路,令dk=T

=\sum_{T=1}^n\lfloor {n\over T}\rfloor\lfloor {m\over T}\rfloor\sum_{d|T}f(d)\mu({T\over d})

数论分块即可 O(n) 预处理, O(\sqrt n)单次询问

莫比乌斯反演形式千千万,要多做才能做出感觉来。

莫比乌斯反演优化多次/单次询问

疯狂 LCM

LCMSUM - LCM Sum(2倍经验)

loj LCMSUM(3倍经验)

P2398 GCD SUM

P1390 公约数的和(2倍经验)

SP21615 NAJPWG - Playing with GCD

SP3871 GCDEX - GCD Extreme

AT5310 [ABC162E] Sum of gcd of Tuples (Hard)(难度绿???我大不服

莫反加整除分块

P3455 [POI2007]ZAP-Queries

P2522 [HAOI2011]Problem b

SP26017 GCDMAT - GCD OF MATRIX

SP26045 GCDMAT2 - GCD OF MATRIX (hard)(卡常,慎入!!

VLATTICE - Visible Lattice Points

杜教筛

主要思路就是构造一个积形函数$g$,使得它本身前缀和与 $$\sum_{i=1}^n\sum_{d|i}f(d)g(\frac{i}{d})$$ 好算一点 那么我们有 $$\sum_{i=1}^n\sum_{d|i}f(d)g(\frac{i}{d})=\sum_{d=1}^ng(d)\sum_{k=1}^{\lfloor\frac{n}{d}\rfloor}f(k)$$ $$g(1)\sum_{k=1}^nf(k)=\sum_{i=1}^n\sum_{d|i}f(d)g(\frac{i}{d})-\sum_{d=2}^ng(d)\sum_{k=1}^{\lfloor\frac{n}{d}\rfloor}f(k)$$ 递归计算即可,当然可以预处理,具体见wiki写的 ### [Min25筛](https://oi-wiki.net/math/min-25/) 学会这个就能搞定一切积形函数前缀和了 题单放不下了我就放在这里 [P1835 素数密度(请小题大做)](https://www.luogu.com.cn/problem/P1835) [P5493 质数前缀统计(min25前置芝士)](https://www.luogu.com.cn/problem/P5493) [P5325 【模板】Min_25筛(板子)](https://www.luogu.com.cn/problem/P5325) [SP20173 DIVCNT2 - Counting Divisors (square)(1倍经验)](https://www.luogu.com.cn/problem/SP20173) [SP20174 DIVCNT3 - Counting Divisors (cube)(2倍经验)](https://www.luogu.com.cn/problem/SP20174) [SP34096 DIVCNTK - Counting Divisors (general)(3倍经验)](https://www.luogu.com.cn/problem/SP34096) [SP19975 APS2 - Amazing Prime Sequence (hard)(有技巧的min25,需要熟练了解其精髓)](https://www.luogu.com.cn/problem/SP19975) [SP22549 DIVFACT4 - Divisors of factorial (extreme)(用min25处理中间问题)](https://www.luogu.com.cn/problem/SP22549) ### 欧拉计划分区(可能要翻墙 [中文翻译](https://pe-cn.github.io/problems/) [欧拉计划388](https://projecteuler.net/problem=388) [欧拉计划625](https://projecteuler.net/problem=625)

  1. P1829 - [集训队互测 2010] Crash的数字表格 / JZPTAB
  2. AT_agc038_c - [AGC038C] LCMs
  3. P3327 - [SDOI2015] 约数个数和
  4. UVA10214 - 树林里的树 Trees in a Wood.
  5. P1447 - [NOI2010] 能量采集
  6. P5176 - 公约数
  7. P4917 - 天守阁的地板
  8. CF1139D - Steps to One
  9. P3704 - [SDOI2017] 数字表格
  10. P7486 - 「Stoi2031」彩虹
  11. P2568 - GCD
  12. P2257 - YY的GCD
  13. SP4491 - PGCD - Primes in GCD Table
  14. P5221 - Product
  15. SP27942 - PROD1GCD - Product it again
  16. P6825 - 「EZEC-4」求和
  17. P4449 - 于神之怒加强版
  18. P2714 - 四元组统计
  19. P3911 - 最小公倍数之和
  20. P5518 - [MtOI2019] 幽灵乐团 / 莫比乌斯反演基础练习题
  21. P4466 - [国家集训队] 和与积
  22. P6156 - 简单题
  23. P6222 - 「P6156 简单题」加强版
  24. P4450 - 双亲数
  25. P7360 - 「JZOI-1」红包
  26. P4318 - 完全平方数
  27. SP4168 - SQFREE - Square-free integers
  28. CF235E - Number Challenge
  29. P4619 - [SDOI2018] 旧试题
  30. P4464 - [国家集训队] JZPKIL
  31. P4213 - 【模板】杜教筛
  32. P3172 - [CQOI2015] 选数
  33. P6055 - [RC-02] GCD
  34. SP26108 - TRENDGCD - Trending GCD
  35. SP19985 - GCDEX2 - GCD Extreme (hard)
  36. P6860 - 象棋与马
  37. P3768 - 简单的数学题
  38. P1587 - [NOI2016] 循环之美