洛谷 P2261 [CQOI2007]余数求和

前两天肝了整整两天的FFT，博客暂咕
今天打算先肝个莫队，于是复习一下整除分块

题面：
给出$n,k(1\le n,k\le 10^9)$，求
$$\sum _{i=1}^{n}k\mathrm{m}\mathrm{o}\mathrm{d}i$$
推一波公式：
$$\begin{array}{rl}& \sum _{i=1}^{n}k\mathrm{m}\mathrm{o}\mathrm{d}i\\ =& \sum _{i=1}^{n}k-i\ast \lfloor \frac{k}{i}\rfloor \\ =& nk-\sum _{i=1}^{n}i\ast \lfloor \frac{k}{i}\rfloor \end{array}$$
推到这里，观察一下$\lfloor \frac{k}{i}\rfloor$，发现：到了后面$i$大的时候$\lfloor \frac{k}{i}\rfloor$可能会保持不变！
重点来看$\lfloor \frac{k}{i}\rfloor$会在什么区间内保持不变：
设$\lfloor \frac{k}{i}\rfloor =\lfloor \frac{k}{j}\rfloor =q(i<j)$
则有$qi+x=qj+y=k(0\le x<i,0\le y<j)$
最极端的情况便是$y=0$
则$qi+x=qj=k$
肉眼可见$j=\frac{k}{q}$
考虑到可能$q\nmid j$
所以$j=\lfloor \frac{k}{q}\rfloor$
对于每一个$x$，它所在的$\lfloor \frac{k}{i}\rfloor$区间上界就是$\frac{k}{\lfloor \frac{k}{x}\rfloor }$
这样就可以通过等差数列对每一块的$i\ast \lfloor \frac{k}{i}\rfloor$进行计算了
每块计算公式：
$\frac{(l+r)\ast (r-l+1)}{2}\ast (k/l)$
肉眼可见初始块一定是$[1,1]$
递推式：$r=\lfloor \frac{k}{\lfloor \frac{k}{l}\rfloor }\rfloor ,l=r+1$

有一个细节：
如果$l>k$，那么$\lfloor \frac{k}{i}\rfloor =0$，会导致除以零
由此我们要加一个特判

#include<cstdio>
using namespace std;
inline int min(const int a,const int b)
{return a<b?a:b;}
int main()
{
	int n,k;
	scanf("%d%d",&n,&k);
	long long ans=0;
	for(int l=1,r=1;l<=n;l=r+1,r=k>l?min(n,k/(k/l)):n)
		ans+=1ll*(l+r)*(r-l+1)/2*(k/l);
	printf("%lld",1ll*n*k-ans);
	return 0;
}

这就是整除分块的一种