整除分块学习笔记

参考资料

https://zhuanlan.zhihu.com/p/513784838

https://oi-wiki.org/math/number-theory/sqrt-decomposition/

常见化整除的公式

$$\begin{gathered} amodb=a-\lfloor \frac{a}{b}\rfloor \ast b \\ \lceil \frac{a}{b}\rceil =\lfloor \frac{a-1}{b}\rfloor +1 \end{gathered}$$

数论分块结论

对常数num，使得式子
$$\lfloor \frac{num}{i}\rfloor =\lfloor \frac{num}{j}\rfloor$$
成立的最大满足$i<=j<=n$的 $j$ 值为 $ \lfloor \frac{num}{\lfloor{\frac{num}{i}}\rfloor}\rfloor$

即$\lfloor \frac{num}{i}\rfloor$所在块的右端点为
$$\lfloor \frac{num}{\lfloor \frac{num}{i}\rfloor }\rfloor$$

应用

对形如
$$\sum _{i=st}^{ed}f(i)\lfloor \frac{num}{i}\rfloor$$
的式子，$\lfloor \frac{num}{i}\rfloor$呈块状分布，若能线性求出$\sum f(i)$，那么复杂可以通过分块降到$O(\sqrt{n})$

代码

//计算区间[st,ed]式子值
int f(int l,int r)//若f(i) = i,求[l,r]内f(i)的和
{
    return r - l + 1;
}
int block(int st, int ed, int num) {
	int R = 0;//块的右边界
	int ans = 0;
	ed = min(ed, num);//大于num下整除都是0
	for (int i = st; i <= ed; i = R + 1) {
		R = min(ed,num / (num / i));  //该区间的最后一个数
		ans += f(i,R)*(num / i);//区间[i,L]的num/i 都是一个值
	}
	return ans;
}