Range Minimum Query (RMQ)

解法一：Sparse Table (ST) algorithm
对2^k 的长度的子数组进行动态规划。我们将使用数组M[0, N-1][0, logN]进行保存，其中M[i][j]是以i 开始，长度为 2^j的子数组的最小值的索引。

一旦我们预处理了这些值，让我们看看怎样使用它们去计算RMQ_A(i, j)。思路是选择两个能够完全覆盖区间[i..j]的块并且找到它们之间的最小值。设k = [log(j - i + 1)].。为了计算RMQ_A(i, j) 我们可以使用下面的公式

                                  

这个算法的总的时间复杂度为 <O(N logN), O(1)>

 

 

const int MAXN = 100; int abs(int a, int b) { return a >= b ? a - b : b - a; } void STinit(int A[], int M[][MAXN], int N) { int i, j; for (i = 0; i < N; i++) M[i][0] = i; for (j = 1; (1 << j) <= N; j++) for (i = 0; i + (1 << j) - 1 < N; i++) if (A[M[i][j-1]] < A[M[i + (1 << (j - 1))][j-1]]) M[i][j] = M[i][j-1]; else M[i][j] = M[i + (1 << (j - 1))][j-1]; } int STquery(int A[], int M[][MAXN], int i, int j) { int len = abs(i, j) + 1; int l = 0; while ((1 << l) <= len) l++; l--; if (A[M[i][l]] < A[M[j - (1 << l) + 1][l]]) return A[M[i][l]]; else return A[M[j - (1 << l) + 1][l]]; }

解法二：segment trees

为了解决RMQ问题我们也可以使用线段树。线段树是一个类似堆的数据结构，可以在基于区间数组上用对数时间进行更新和查询操作。我们用下面递归方式来定义线段树的[i, j]区间：
注意具有N个区间元素的线段树的高度为[logN] + 1。下面是区间[0,9]的线段树：

 

线段树和堆具有相同的结构，因此我们定义x是一个非叶结点，那么左孩子结点为2*x,而右孩子结点为2*x+1。
使用线段树解决RMQ问题，我们应该使用数组 M[1, 2 * 2^{[logN] + 1}]，这里M[i]保存结点i区间最小值的位置。初始时M的所有元素为-1。树应当用下面的函数进行初始化(b和e是当前区间的范围)：#include <cstdio> void init(int node, int A[], int M[], int b, int e) { if (b == e) M[node] = b; else { init(2 * node, A, M, b, b + (e - b) / 2); init(2 * node + 1, A, M, b + (e - b) / 2 + 1, e); if (A[M[2 * node]] < A[M[2 * node + 1]]) M[node] = M[2 * node]; else M[node] = M[2 * node + 1]; } } int query(int node, int A[], int M[], int b, int e, int i, int j) { if (i > e || j < b) return -1; if (i <= b && e <= j) return M[node]; int r1 = query(2 * node, A, M, b, b + (e - b) / 2, i, j); int r2 = query(2 * node + 1, A, M, b + (e - b) / 2 + 1, e, i, j); if (r1 == -1) return M[node] = r2; if (r2 == -1) return M[node] = r1; if (A[r1] < A[r2]) return M[node] = r1; return M[node] = r2; } int main() { int A[] = {10, 9, 8, 7, 6, 5, 4, 3, 2, 1}; int LEN_A = sizeof(A) / sizeof(int), LOG_LEN_A = 0; while (1 << LOG_LEN_A <= LEN_A) LOG_LEN_A++; LOG_LEN_A--; int *M = new int[2 * (1 << (LOG_LEN_A + 1)) + 1]; // 线段树的高度为LOG_LEN_A + 1, // 所以树的节点最多有2 * (1 << (LOG_LEN_A + 1)) 个 init(1, A, M, 0, LEN_A - 1); printf("%d/n", A[query(1, A, M, 0, LEN_A - 1, 2, LEN_A - 1)]); printf("%d/n", A[query(1, A, M, 0, LEN_A - 1, 3, LEN_A - 1)]); printf("%d/n", A[query(1, A, M, 0, LEN_A - 1, 4, LEN_A - 1)]); printf("%d/n", A[query(1, A, M, 0, LEN_A - 1, 5, LEN_A - 1)]); return 0; }
调用init和query都应该使用node = 1，b = 0，e = N - 1来调用。
可以很容易的看出任何查询都可以在O(log N)内完成。注意当我们碰到完整的in/out区间时我们停止了，因此数中的路径最多分裂一次。用线段树我们获得了<O(N), O(logN)>的算法。线段树非常强大，不仅仅是因为它能够用在RMQ上，还因为它是一个非常灵活的数据结构，它能够解决动态版本的RMQ问题和大量的区间搜索问题。