【“分块”算法知识点】
● 分块是用线段树的分区思想改良的暴力法。代码比线段树简单。效率比普通暴力法高。分块适合求解 m=n=10^5 规模的问题,或 m*sqrt(n)≈10^7 的问题。其中,n 为元素个数,m 为操作次数。
● “分块”算法的基本要素
(1)块的大小用 block 表示。通常,令 block=sqrt(n)。其中,n 为元素个数。
(2)块的数量用 cnt 表示。计算块的数量的代码如下:
int block=sqrt(n);
int cnt=n/block;
if(n % block) cnt++;
(3)定义 pos[i] 为第 i 个元素所在的块。
若下标从 1 开始,则有 pos[i]=(i-1)/block+1。其中,block=sqrt(n)。
若下标从 0 开始,则有 pos[i]=i/block。其中,block=sqrt(n)。
(4)块的左边界 le[] 及右边界 ri[]。
若用 le[i] 和 ri[i] 分别表示块 i 的第一个和最后一个元素的位置。
若下标从 1 开始,则有:
le[1]=1, ri[1]=block;
le[2]=block+1, ri[2]=2*block;
……
le[i]=(i-1)*block+1, ri[i]=i*block;
……
若下标从 0 开始,则有:
le[0]=0, ri[0]=block-1;
le[1]=block, ri[1]=2*block-1;
……
le[i]=i*block, ri[i]=(i+1)*block-1;
……
● 综上,“分块”算法 build() 函数的构建细节如下。
(1)下标从 0 开始,build() 函数的构建如下。
void build(int n) {
int block=sqrt(n);
int cnt=n/block;
if(n%block) cnt++;
for(int i=0; i<cnt; i++) {
le[i]=i*block;
ri[i]=(i+1)*block-1;
}
ri[cnt-1]=n-1;
for(int i=0; i<n; i++) pos[i]=i/block;
}
(2)下标从 1 开始,build() 函数的构建如下。
void build(int n) {
int block=sqrt(n);
int cnt=n/block;
if(n%block) cnt++;
for(int i=1; i<=cnt; i++) {
le[i]=(i-1)*block+1;
ri[i]=i*block;
}
ri[cnt]=n;
for(int i=1; i<=n; i++) pos[i]=(i-1)/block+1;
}
(3)上文中的两段 build() 函数,是下标从 0 及从 1 开始时,“分块”算法 build() 函数的基础代码。不同的问题,会在 build() 函数代码后添加一些不同的内容。详见:
https://blog.csdn.net/hnjzsyjyj/article/details/138903837
https://blog.csdn.net/hnjzsyjyj/article/details/138926594
https://blog.csdn.net/hnjzsyjyj/article/details/138863063
特别注意,针对不同的问题进行分析时,下图将具有极高的应用价值。其对整块及碎块的处理一目了然。
【参考文献】
https://blog.csdn.net/hnjzsyjyj/article/details/138903837
https://blog.csdn.net/hnjzsyjyj/article/details/138926594
https://blog.csdn.net/hnjzsyjyj/article/details/138863063