整除分块属于数论的一部分。这里尽可能脱离数论的知识。看一道例题。
已知正整数n,求
∑ i = 0 n = ⌊ n i ⌋ \sum_{i=0}^n=\lfloor\frac{n}{i}\rfloor i=0∑n=⌊in⌋
以n=8为例,结果如下表
i | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
---|---|---|---|---|---|---|---|---|
⌊ n i ⌋ \lfloor\frac{n}{i}\rfloor ⌊in⌋ | 8 | 4 | 2 | 2 | 1 | 1 | 1 | 1 |
可以看到上面的表中相连着的同一数值的区间不少,由这些数值可以将这张表划分为多个区间块,这就是整除分块了。整除的性质使得从1到n的数组表可根据数值划分为不同的分块,且分块数远远小于n。利用这种性质,我们如果能推导出每个分块具体的左右端点位置在哪,这个问题就可以快速求解出来了。
假设我们已知某一个分块的左端点l,要求解出该分块的右端点r。设该分块的数值为k,对于该分块中的每个数i,有k = ⌊ n i ⌋ \lfloor\frac{n}{i}\rfloor ⌊in⌋= ⌊ n l ⌋ \lfloor\frac{n}{l}\rfloor ⌊ln⌋即i k ≤ n ,稍加思索一番我们可得使i k ≤ n 成立的最大的i 值即是我们所求的右端点r,因此我们可以得到公式
r = ⌊ n k ⌋ = ⌊ n ⌊ n l ⌋ ⌋ r=\lfloor\frac{n}{k}\rfloor=\lfloor\frac{n}{\lfloor\frac{n}{l}\rfloor}\rfloor r=⌊kn⌋=⌊⌊ln⌋n⌋
代码如下:
ans = 0;
for(int l = 1, r; l <= n; l = r + 1)
{
r = n / (n / l);
ans += n / l * (r - l + 1);
}
简单的整除分块是不是很简单,如果你以为你已经掌握了精髓,那可能就是大错特错了。下面请看新的题目。
已知正整数n,a,b,求 ∑ i = 0 n = ⌊ n a i + b ⌋ \sum_{i=0}^n=\lfloor\frac{n}{ai+b}\rfloor i=0∑n=⌊ai+bn⌋
两种解法:1.直接进行推导, k = ⌊ n a l + b ⌋ k=\lfloor\frac{n}{al+b}\rfloor k=⌊al+bn⌋, r = m a x ( i ) , ( a i + b ) k < = n r=max(i),(ai+b)k<=n r=max(i),(ai+b)k<=n得到 r = ⌊ n − b k a k ⌋ r=\lfloor\frac{n-bk}{ak}\rfloor r=⌊akn−bk⌋
2.整体代换
r ∗ = ⌊ n ⌊ n a l + b ⌋ ⌋ r*=\lfloor\frac{n}{\lfloor\frac{n}{al+b}\rfloor}\rfloor r∗=⌊⌊al+bn⌋n⌋ r ∗ = a l + b r*=al+b r∗=al+b