数学板块学习之整除分块

整除分块

题什么的看here

学莫比乌斯反演总看见整除分块，所以！。。先来学一下整除分块
直接通过给出问题看
求
$$\sum _{i=1}^n[\frac{n}{i}],n\le 10^{12}$$

首先来说暴力肯定跑不了
然后发现对于$[\frac{n}{i}]$最多只有$2\sqrt{n}$个值
并且如果$[\frac{n}{i^{\prime }}]=[\frac{n}{i}]$则$i^{\prime }$的最大值为$[\frac{n}{[\frac{n}{i}]}]$
用法就是设$L=1$，令$R=[\frac{n}{[\frac{n}{L}]}]$。将$(R-L+1)\ast [\frac{n}{L}]$累加，再令$L=R+1$

代码如下

for(int l = 1,r;l <= n; l = r+1) {
	r = (n/(n/l));
	ans += (r-l+1)*(n/l);
}

小技巧：

分块$[\frac{n}{i}]$ $$r=[\frac{n}{[\frac{n}{l}]}]即r=(n/(n/l))$$
分块$[\frac{n}{i^2}]$ $$r=\left[\sqrt{[\frac{n}{\frac{n}{l^2}}]}\right],即r=sqrt(n/(n/l/l))$$
分块$\lfloor \sqrt{\lfloor \frac{n}{i}\rfloor }\rfloor$ $$r=\lfloor \frac{n}{(\lfloor \sqrt{\lfloor \frac{n}{i}\rfloor }\rfloor {)}^2}\rfloor 即r=n/(sqrt(n/l)\ast sqrt(n/l))$$
证明直接上图，用正在学的latex写数学公式好爽，但是好慢。。。字丑了点（何止一点）QAQ敲代码敲得都不会写字了