9.1 RMQ问题
RMQ(Range Minimum/Maximum Query)问题,就是对于给定数组,在其下标范围[i, j]内给出的最小值(或最大值,下文都称最小值)的问题。
RMQ和LCA问题可以互相转化,一个对RMQ和LCA问题的总结如下表[i]:
算法 | 处理方式 | 复杂度 | 备注 |
转化算法 | |||
LCA =>±1RMQ | N/A | O(n) | 引理1,规模O(n) - O(2n-1) |
RMQ => LCA-CT | N/A | O(n) | 引理2,规模不变 |
朴素算法 | |||
LCA-Naive | online | O(n^2) - O(1) | 动态规划 |
RMQ-Naive | online | O(n^2) - O(1) | 直接求解 |
经典算法 | |||
LCA-Tarjan | offline | O(na(n)) | |
RMQ-ST | online | O(nlogn) - O(1) | |
改进算法 | |||
RMQ-ST-Block | online | O(nloglogn) - O(1) | 将RMQ-ST分段处理 |
±1RMQ-ST-Block | online | O(n) - O(1) | 控制RMQ-ST-Block中分段的段种数,得到O(n)算法。 |
快速算法 | |||
RMQ-Fast | online | O(n) - O(1) | RMQ => LCA-CT => ±1RMQ |
折衷算法 | |||
RMQ-IT | online | O(n) - O(logn) | 线段树直接处理 |
RMQ-CT-Tarjan | offline | O(na(n)) | RMQ => LCA-CT再用LCA-Tarjan解决。 |
我们可以将一般RMQ的值得应用的算法列表:
算法 | 处理方式 | 复杂度 | 备注 |
RMQ-IT | online | O(n)-O(logn) | 询问不多时,竞赛首选。 |
RMQ-CT-Tarjan | offline | O(na(n)) | 由于是离线的,而且还要转化成CT,某些地方不应用。 |
RMQ-ST | online | O(nlogn)-O(1) | 询问多时,竞赛首选。 |
RMQ-ST-Block | online | O(nloglogn)-O(1) | 时间要求特别严格时采用。 |
1.2 Sparse Table算法
RMQ的一个经典的算法是SparseTable(ST)算法。它可以在O(nlgn)的预处理时间后,提供O(1)时间的查询。
定义d[i, j]为下标在区间[i, i+2^j-1]中的数组S的元素的最小值。即:
d[i, j] = min{S[k] | i ≤ k ≤ i+2^j-1}
有:
d[i, j-1] = min{S[k] | i ≤ k ≤ i+2^(j-1)-1}
d[i+2^(j-1), j-1] = min{S[k] | i+2^(j-1) ≤ k ≤ i+2^(j-1)+2^(j-1)-1} =min{S[k] | i+2^(j-1) ≤ k ≤ i+2^j-1}
故此有:
d[i, j] = min{S[k] | i ≤ k ≤ i+2^j-1} = min{d[i, j-1], d[i+2^(j-1),j-1]}
可以观察到区间[i, i+2^j-1]被分为区间[i, i+2^(j-1)-1]和区间[i+2^(j-1), i+2^j-1]。
对于所有{(i, j) | 0 ≤ i ≤ n-1, 0 ≤i+2^j-1 ≤ n-1 }都求出d[i, j]。首先依据定义有d[i, 0] = S[i]。又由于有上述的递推式,所以在已求出d[i, j-1]和d[i+2^(j-1), j-1]的前提下,d[i, j]可以在常数时间内求出。又由于{(i, j) | 0 ≤ i ≤ n-1,0 ≤ i+2^j-1 ≤ n-1} = {(i, j) | 0 ≤ i ≤ n-1-(2^j-1), 0 ≤ 2^j-1 ≤ n-1},所以可以在[0, n-1]区间内遍历2^j-1,并对每个j,求出每个合要求的i对应的d[i, j]。如果采用这种方法,在求解d[i, j]之前,d[i, j-1]和d[i+2^(j-1), j-1]都已经求出,故此求d[i, j]可以在常数时间内完成。由于j有O(lgn)个选择,而对每个j,i有O(n)个选择,故此求出数组d需要的时间是O(n lgn)。
假设要在[i, j]区间内查询最小值,即RMQ(i, j),则令k = max{l | 2^l ≤ (i-j+1)},即k为log(i-j+1)取下限。由此2^k*2 ≥ i-j+2,即i+2^k-1 ≥ j-2^k+1,所以[i, j]可被子区间[i, i+2^k-1]和[j-2^k+1, j]覆盖。而两个区间内的最小值分别是d[i, k]和d[j-2^k+1, k]。所以RMQ(i, j) = min{d[i, k], d[j-2^k+1,k]}。所以任意的查询都可以在O(1)时间内完成。
1.2.1实例
PKU JudgeOnline, 3264, Balanced Lineup.
1.2.2问题描述
给出一组数据,求出一个坐标区间内的最大值和最小值的相差范围。
1.2.3输入
63
1
7
3
4
2
5
15
46
22
1.2.4输出
6
3
0
1.2.5分析
这是经典的RMQ问题,用线段树也可以做。
可以发现,可以将二维数组d的两维互换,这样可以提高访问速度。
1.2.6程序
#include<stdio.h> #include<string.h> #include<iostream> using namespace std; #define maxNum 50002 int Dmax[16][maxNum];//注意这里为了提高访问速度,将数组的两维互换 int Dmin[16][maxNum]; int cow[maxNum]; int n; void bulidST() { int k; int l; int i, j; for(i = 0;i < n; i++){ Dmin[0][i] = cow[i]; Dmax[0][i] = cow[i]; } for(j = 1;(k = (1<<j) - 1) < n; j++){ for(i =0; i < n - k; i++){ l = (k + 1)>>1; if(Dmin[j- 1][i] < Dmin[j - 1][i + l]){ Dmin[j][i] = Dmin[j - 1][i]; }else{ Dmin[j][i] = Dmin[j - 1][i +l]; } if(Dmax[j- 1][i] > Dmax[j - 1][i + l]){ Dmax[j][i] = Dmax[j - 1][i]; }else{ Dmax[j][i] = Dmax[j - 1][i +l]; } } } } int checkST(int i, intj) { int k; int d; int RMaxQ; int RMinQ; k = 0; d = 1; while(d<< 1 <= j - i + 1){ d <<= 1; k++; } if(Dmin[k][i]< Dmin[k][j - d + 1]){ RMinQ = Dmin[k][i]; }else{ RMinQ = Dmin[k][j - d + 1]; } if(Dmax[k][i]> Dmax[k][j - d + 1]){ RMaxQ = Dmax[k][i]; }else{ RMaxQ = Dmax[k][j - d + 1]; } // cout << RMaxQ<< " : " << RMinQ << endl; returnRMaxQ - RMinQ; } int main() { int q; int i; int l, r; scanf("%d%d",&n, &q); for(i = 0;i < n; i++){ scanf("%d",&cow[i]); } bulidST(); for(i = 0;i < q; i++){ scanf("%d%d",&l, &r); printf("%d\n",checkST(l - 1, r - 1)); } } 本文章欢迎转载,请保留原始博客链接http://blog.csdn.net/fsdev/article[i] http://adn.cn/blog/article.asp?id=20