知识点 - 整除分块

知识点 - 整除分块

1.约数个数

令$\tau (n)$为n的约数个数，求其前缀和。

解析

若n的质因数分解为$n=p_1^{e_1}\cdot p_2^{e_2}\cdots p_k^{e_k}$，则
$$\tau (n)=(e_1+1)\cdot (e_2+1)\cdots (e_k+1)$$
不过求这个前缀和时不需要此公式。

我们有
$$\sum _{i=1}^n\tau (i)=\sum _{i=1}^n\sum _{j=1}^n[j|i]=\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor$$
对于
$$\begin{gathered} y_1=\lfloor \frac{n}{x}\rfloor \\ {y}_2=\frac{n}{x} \end{gathered}$$
$n=12$时作图如下：

当$i\le \sqrt{n}$时，$\lfloor \frac{n}{i}\rfloor$显然只有$O(\sqrt{n})$个取值；当$i>\sqrt{n}$时，$\lfloor \frac{n}{i}\rfloor <\sqrt{n}$显然也只有$O(\sqrt{n})$个取值；
对于固定$\lfloor \frac{n}{i}\rfloor$,$i$的取值是一段连续的区间，这段区间是$[\lfloor \frac{n}{\lfloor \frac{n}{i}\rfloor +1}\rfloor +1,\lfloor \frac{n}{\lfloor \frac{n}{i}\rfloor }\rfloor ]$
几何上可以理解为先找到对应的红线，即求出i点的值$v=\lfloor \frac{n}{i}\rfloor$；然后求出红线的终点,即右区间$r=\lfloor \frac{n}{v}\rfloor$,以及它的起点$l=\lfloor \frac{n}{v+1}\rfloor$

将$\lfloor \frac{n}{\lfloor \frac{n}{i}\rfloor }\rfloor$作为分界线，所以有如下代码：

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

右区间的代数证明，?戳这里。

2.约数和

令$\sigma (i)$为n的约数和，求前缀和。

oeis A024916

约数和

解析

类似上面，有：
$$\sigma (n)=\frac{p_1^{e_1+1}-1}{p_1-1}\cdot \frac{p_2^{e_2+1}-1}{p_2-1}\cdots \frac{p_k^{e_k+1}-1}{p_k-1}$$
同上也没用。因为我们可以直接写出：
$$\sum _{i=1}^n\sigma (i)=\sum _{i=1}^n\sum _{j=1}^n[j|i]\cdot j=\sum _{i=1}^{n}i\cdot \sum _{j=1}^n[i|j]=\sum _{i=1}^{n}i\cdot \lfloor \frac{n}{i}\rfloor =\sum _{i=1}^n\frac{\lfloor \frac{n}{i}\rfloor \cdot (\lfloor \frac{n}{i}\rfloor +1)}{2}$$
最后的等式
$$\sum _{i=1}^{n}i\cdot \lfloor \frac{n}{i}\rfloor =\sum _{i=1}^n\frac{\lfloor \frac{n}{i}\rfloor \cdot (\lfloor \frac{n}{i}\rfloor +1)}{2}$$
可以看成对每个$\lfloor \frac{n}{i}\rfloor$固定的块，将 $\lfloor \frac{n}{i}\rfloor$提出来，然后对$i\in [\lfloor \frac{n}{\lfloor \frac{n}{i}\rfloor +1}\rfloor +1,\lfloor \frac{n}{\lfloor \frac{n}{i}\rfloor }\rfloor ]$等差数列求和。

int f(int n) {
	int ans = 0;
	for (int l = 1, r; l <= n; l = r + 1) {
		r = n / (n / l);
		ans += (r - l + 1)*(n / l)*(n / l + 1) / 2;
	}
	return ans;
}
signed main()
{
	int l, r;
	while (cin >> l>>r) {
		cout << f(r) - f(l - 1) << endl;
	}
	return 0;
}