整数分块——operator_

1.算法简介

这是一种很简单但好用且优雅的 $\sqrt{n}$ 分块。

考虑求 ${\sum }_{i=1}^n\lfloor \frac{n}{i}\rfloor$ ，很容易发现 $\lfloor \frac{n}{i}\rfloor$ 的值的种类不多，考虑按值域来求解。

事实上，$\lfloor \frac{n}{i}\rfloor$ 的值最多只有 $2\sqrt{n}$ 种(并不严谨，但反正是 $\sqrt{n}$ 级别的)，证明如下：

对于 $i>\sqrt{n}$ ，$\lfloor \frac{n}{i}\rfloor \le \frac{n}{i}<\sqrt{n}$ ，显然只有最多 $\sqrt{n}$ 个整数符合。

对于 $i\le \sqrt{n}$ ，每个 $i$ 的 $\lfloor \frac{n}{i}\rfloor$ 不一样也最多只有$\sqrt{n}$ 个。

然后就证毕了。

如果你想，你可以去 AC 一道CF绿题。

然后我们回到本来的式子上去，因为值域不连续，而区间是连续的，所以考虑枚举每一种取值相等的区间，具体来说就是递推左边界，并求出对应区间的右边界。显然地，在 $i\le \sqrt{n}$ 时左边界等于右边界。

求右边界的方法很简单：r=n/(n/l)(均为整除)。 其证明如下：

令 $a=\lfloor \frac{n}{l}\rfloor$

$\lfloor \frac{n}{r}\rfloor =a\iff a\le \frac{n}{r}\iff r\le \frac{n}{a}$

所以 $r_{max}=\lfloor \frac{n}{a}\rfloor$

代码简洁至极：

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

2.例题

本来不想放例题的，但洛谷上没标签，硬生生找了好久的题

例1.P3935 Calculating

原式 $={\sum }_{i=1}^{r}f(i)-{\sum }_{i=1}^{l-1}f(i)$

令 $g(x)={\sum }_{i=1}^{x}f(i)$ ，则原式 $=g(r)-g(l-1)$ (就是个前缀和罢了…)

考虑 $f(i)$ ，发现其实就是 $i$ 的因数个数，所以我们换个思路，不要对于每个数求因数个数，而是求每个数的贡献 ，所以 $g(x)={\sum }_{i=1}^x\lfloor \frac{x}{i}\rfloor$ ，然后就是一个板题了。

例2.P2261 [CQOI2007] 余数求和

原式 $={\sum }_{i=1}^{n}k-i\lfloor \frac{k}{i}\rfloor =nk-{\sum }_{i=1}^{n}i\lfloor \frac{k}{i}\rfloor$

还是按板子来，对于每一块相等的数作一个等差数列求和作为贡献就好了。

例3.[ARC068E] Snuke Line

容易发现，第 $i$ 个商品会对列车 $d$ 做出贡献当且仅当 $l_i\le k\ast d\le r_i$ ，等价于 $\frac{l_i}{d}\le k\le \frac{r_i}{d}$

因为 $k$ 是整数，所以有 $\lfloor \frac{l_i-1}{d}\rfloor <\lceil \frac{l_i}{d}\rceil \le k\le \lfloor \frac{r_i}{d}\rfloor$

即 $\lfloor \frac{l_i-1}{d}\rfloor <\lfloor \frac{r_i}{d}\rfloor$

这个东西发现满足条件的 $d$ 不超过 $\sqrt{n}$ 个连续区间，可以用整数分块+差分搞定。

例4.P2260 [清华集训2012] 模积和

先容斥掉 $i\ne j$ 的限制，然后按例2的方法做，得到一个比较恶心的式子：(记$g_{n,x}={\sum }_{i=1}^{n}i\lfloor \frac{x}{i}\rfloor ,minn=min(n,m))$

$$(n^2-g_{n,n})\ast (m^2-g_{m,m})-minn\cdot nm+n{g}_{minn,m}+m{g}_{minn,n}-\sum _{i=1}^{minn}{i}^2\lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{i}\rfloor$$

然后呢？就没有然后了，都是可以整数分块维护的。

特别注意一下最后那个东西，需要用到 ${\sum }_{i=1}^n{i}^2=\frac{n(n+1)(2n+1)}{6}$