RMQ问题详解 --Sparse-Table算法实现和线段树实现

RMQ:范围最小值问题
给出有n个元素的数组a1,a2,…,an.设计一个数据结构支持多次查询Query（L,R):计算min(aL,aL+1,…,aR).

Sparse-Table算法：
令d[i][j]表示下标从i开始长度为2^j的一段元素中的最小值。所以我们可以得到递推公式：
d[i][j]=min(d[i][j-1],d[i+2^j-1][j-1]).可以看出d[i][j]可以有他前面的一个数和左下角的一个数求出来。所以我们再求d[i][j]时先求出前面所以得列然后求行。
求d[i][j]:

void init(const vector<int>& A) {
    int n = A.size();
    for(int i = 0; i < n; i++) d[i][0] = A[i];
    for(int j = 1; (1<<j) <= n; j++)
      for(int i = 0; i + (1<<j) - 1 < n; i++)
        d[i][j] = max(d[i][j-1], d[i + (1<<(j-1))][j-1]);
  }

那我们如何求（L,R）的最小值？令k是满足2^k<=R-L+1的最大整数。以L开头和以R结尾的两个长度为2^k的区间刚好覆盖了查询区间[L,R].这两个区间的最小值就是整个区间的最小值。

int query(int L, int R) {
    int k = 0;
    while((1<<(k+1)) <= R-L+1) k++; // 如果2^(k+1)<=R-L+1，那么k还可以加1
    return max(d[L][k], d[R-(1<<k)+1][k]);
  }

例题：（UVA 11235 Frequent Values）
题目链接：http://acm.hust.edu.cn/vjudge/problem/23846

题意：
给你一个非降序的整数数组，你的任务是对于一系列查询（i,j）,回答：ai,ai+1,…,aj中出现最大的值出现的次数？
分析：
由于是非降序的，所以所有相同的值都聚集在一起。每一组相同的值我们看成是一段。我们用count[i]表示第i段数的个数。如果查询刚好包括几个连续的段，那么问题就是求区间最大值问题了。但是查询的区间可能将某一段的后半部分，中间几段和最后一段的前一部分。那么我们就要将最前面一部分和最后一部分的值要进行更正。我们用num[p],left[p],right[p]分别表示位置p所在段的编号和左右端点的位置。所以查询(L,R)的结果就变成了三部分的最大值：从L到L段结束处元素的个数（right[L]-L+1),从R段的开始到R处元素的个数（R-left[R]+1),中间num[L]+1段到num[R]-1段的count的最大值。
特殊情况：当L，R在同一段时答案为R-L+1;

AC代码：

#include<cstdio>
#include<algorithm>
#include<vector>
using namespace std;

const int maxn = 100000 + 5;
const int maxlog = 20;

// 区间最*大*值
struct RMQ {
  int d[maxn][maxlog];
  void init(const vector<int>& A) {
    int n = A.size();
    for(int i = 0; i < n; i++) d[i][0] = A[i];
    for(int j = 1; (1<<j) <= n; j++)
      for(int i = 0; i + (1<<j) - 1 < n; i++)
        d[i][j] = max(d[i][j-1], d[i + (1<<(j-1))][j-1]);
  }

  int query(int L, int R) {
    int k = 0;
    while((1<<(k+1)) <= R-L+1) k++; // 如果2^(k+1)<=R-L+1，那么k还可以加1
    return max(d[L][k], d[R-(1<<k)+1][k]);
  }
};

int a[maxn], num[maxn], left[maxn], right[maxn];
RMQ rmq;
int main() {
  int n, q;
  while(scanf("%d%d", &n, &q) == 2) {
    for(int i = 0; i < n; i++) scanf("%d", &a[i]);
    a[n] = a[n-1] + 1; // 哨兵
    int start = -1;
    vector<int> count;
    for(int i = 0; i <= n; i++) {
      if(i == 0 || a[i] > a[i-1]) { // 新段开始
        if(i > 0) {
          count.push_back(i - start);
          for(int j = start; j < i; j++) {
            num[j] = count.size() - 1; left[j] = start; right[j] = i-1;
          }
        }
        start = i;
      }
    }
    rmq.init(count);
    while(q--) {
      int L, R, ans;
      scanf("%d%d", &L, &R); L--; R--;
      if(num[L] == num[R]) ans = R-L+1;
      else {
        ans = max(R-left[R]+1, right[L]-L+1);
        if(num[L]+1 < num[R]) ans = max(ans, rmq.query(num[L]+1, num[R]-1));
      }
      printf("%d\n", ans);
    }
  }
  return 0;
}