来自Peterwuyihong 的题单。

前置知识

前置芝士1 数论分块

UVA11526 H(n)

P2261 [CQOI2007]余数求和

P2260 [清华集训2012]模积和

  • 其中有一个式子需要注意一下:

    \[\sum_{i=1}^n\sum_{j=1}^m\left\lfloor\dfrac{n}{i}\right\rfloor\left\lfloor\dfrac{m}{i}\right\rfloor=\left(\sum_{i=1}^n\left\lfloor\dfrac{n}{i}\right\rfloor\right)\times\left(\sum_{j=1}^m\left\lfloor\dfrac{m}{i}\right\rfloor\right) \]

P3935 Calculating

  • 最重要的一部分就是计算 \(1\) 到 \(n\) 内每个数的约数个数之和,这样不是很好做,考虑枚举约束。

    \[\sum_{i=1}^{n}\sigma_0(i)=\sum_{d=1}^{n}\left\lfloor\dfrac{n}{d}\right\rfloor \]

前置芝士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

\[\begin{aligned} \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~\texttt{and}~x|j]\\ &=\sum_{x=1}^n\mu(x)\left\lfloor {n\over x}\right\rfloor\left\lfloor {m\over x}\right\rfloor \end{aligned} \]

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

形式2

\[\begin{aligned} \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}\varphi(x)\\ &=\sum_{x=1}^n\varphi(x)\sum_{i=1}^n\sum_{j=1}^m[x|i~\texttt{and}~x|j]\\ &=\sum_{x=1}^n\varphi(x)\left\lfloor {n\over x}\right\rfloor\left\lfloor {m\over x}\right\rfloor \end{aligned} \]

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

形式3

\[\begin{aligned} \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)\left\lfloor {n\over dk}\right\rfloor\left\lfloor {m\over dk}\right\rfloor\text{,然后令}dk=T\\ &=\sum_{T=1}^n\left\lfloor {n\over T}\right\rfloor\left\lfloor {m\over T}\right\rfloor\sum_{d|T}f(d)\mu\left({T\over d}\right) \end{aligned} \]

数论分块即可 \(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

杜教筛

\(\sum_{i=1}^nf(i),f\) 为积形函数

主要思路就是构造一个积形函数\(g\),使得它本身前缀和与

\[\sum_{i=1}^n\sum_{d|i}f(d)g\left(\frac{i}{d}\right) \]

好算一点

那么我们有

\[\sum_{i=1}^n\sum_{d|i}f(d)g\left(\frac{i}{d}\right)=\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\left(\frac{i}{d}\right)-\sum_{d=2}^ng(d)\sum_{k=1}^{\lfloor\frac{n}{d}\rfloor}f(k) \]

递归计算即可,当然可以预处理,具体见wiki写的

Min25筛

学会这个就能搞定一切积形函数前缀和了

题单放不下了我就放在这里

P1835 素数密度(请小题大做)

P5493 质数前缀统计(min25前置芝士)

P5325 【模板】Min_25筛(板子)

SP20173 DIVCNT2 - Counting Divisors (square)(1倍经验)

SP20174 DIVCNT3 - Counting Divisors (cube)(2倍经验)

SP34096 DIVCNTK - Counting Divisors (general)(3倍经验)

SP19975 APS2 - Amazing Prime Sequence (hard)(有技巧的min25,需要熟练了解其精髓)

SP22549 DIVFACT4 - Divisors of factorial (extreme)(用min25处理中间问题)

欧拉计划分区(可能要FQ

中文翻译

欧拉计划388

欧拉计划625

莫比乌斯反演 做题记录相关推荐

  1. 退役前的做题记录1.0

    退役前的做题记录1.0 租酥雨最近很懒qwq,具体表现在写题的时候不想发题解了. 但是想想这样也不太好,就决定发个一句话(半句话到几句话不等)题解上来. 2018-09.18-2018-09.28 [ ...

  2. 退役前的做题记录2.0

    退役前的做题记录2.0 最近在刷省选题......大致上是按照省份刷的. 不过上面的题目顺序是按照写题的顺序排列的,所以可能会有点乱哈. [BZOJ2823][AHOI2012]信号塔 最小圆覆盖,随 ...

  3. 2020.9月做题记录

    八月的做题记录因为是暑假所以鸽掉了. 离联赛真的不远了,要继续努力啊qwq- week -1 2020.08.30 2020.08.30 今天考试,修了20+次锅,修的我都没有心情做题了- 然后开始消 ...

  4. 退役前的做题记录5.0

    退役前的做题记录5.0 出于某种原因新开了一篇. [CodeChef]Querying on a Grid 对序列建立分治结构,每次处理\((l,mid,r)\)时,以\(mid\)为源点建立最短路树 ...

  5. 2020.7月做题记录

    转眼就到了2020的下半年了-前方仍是一片茫然. 长期计划 prufer 序列 2020.07.02-2020.07.04 Problem Finished P2624 [HNOI2008]明明的烦恼 ...

  6. 退役前的做题记录4.0

    退役前的做题记录4.0 最近主要在LOJ上写题 536. 「LibreOJ Round #6」花札 比较显然的二分图博弈模型,先手必胜当且仅当起始点一定在最大匹配中.连边可以对每种颜色以及数字建一个点 ...

  7. 概率期望题(期望 DP)做题记录

    概率期望题(期望 DP)做题记录 P3830 [SHOI2012]随机树 难点在于第二问:生成树的期望深度. 不 wei zhuo 捏,设 \(dp_{i,j}\) 表示已经有了 \(i\) 个叶子结 ...

  8. 数数题(计数类 DP)做题记录

    数数题(计数类 DP)做题记录 CF1657E Star MST 我们称张无向完全图是美丽的当且仅当:所有和 \(1\) 相连的边的边权之和等于这张完全图的最小生成树的边权之和. 完全图点数为 \(n ...

  9. CSDN 第六期编程竞赛做题记录

    CSDN 第六期编程竞赛做题记录 -- CSDN编程竞赛报名地址:https://edu.csdn.net/contest/detail/16 9.18周日闲来无视写一下 csdn 的编程题,每期编程 ...

最新文章

  1. python applymap_Python pandas.DataFrame.applymap函数方法的使用
  2. 2012自学CCNP路由与交换之四交换机初始化
  3. 爬虫为什么使用asyncio以及邮件系统为什么使用celery
  4. HD-SDI光端机是什么?其性能特点和技术参数有哪些?
  5. 【代码真相】之 开篇
  6. Fragment的保存
  7. linux命令本质,Linux 的命令机制
  8. 详细解读Spark的数据分析引擎:Spark SQL
  9. mate7 刷机 android 7,Mate7 四大版本完整稳定版刷机包大集合!
  10. Kubernetes详解(十八)——Pod就绪性探针实战
  11. 第十届全国大学生智能汽车竞赛获奖名单
  12. 前端css简易拾色器
  13. 路由与交换 实验一 路由器基本配置
  14. FFmpeg —— 14.示例程序(八):音频编码器(PCM编码为AAC)
  15. [MATLAB] 心形图
  16. 如何调试分布式系统:与微服务调试工具“Squash”创始人Idit Levine的对话
  17. Portapack应用开发教程 自定义开机图片
  18. 什么是AQS?AQS有什么作用?常见的基于AQS的组件有哪些?
  19. openbsd运行Linux应用程序,OpenBSD操作系统(1)——服务管理程序rcctl 类似linux下的service或syst...
  20. ChatGPT 辅助专利写作

热门文章

  1. vim 寄存器 操作_vim指令
  2. android 活动切换动画,android – 在使用ChangeImageTransform共享元素转换的两个活动之间动画化ImageView...
  3. vba 不等于_EXCEL中VBA基础语句(1)
  4. 哈希表(散列表)基础概念与经典题目(Leetcode题解-Python语言)之下——设计键
  5. Git入门教程(一)
  6. [Java网络编程]UDP通信程序练习
  7. [Java基础]自动装箱和拆箱
  8. poj2182 Lost Cows-暴力
  9. UVA - 839 Not so Mobile
  10. ./include/caffe/util/cudnn.hpp: error: too few arguments to function ‘SetPooling2dDescriptor