ST表·RMQ问题

ST表实际是一种动态规划思想求解RMQ问题的方法。离线预处理O(nlogn)， 在线查询O(1).

以最大值为例：

定义一个二维数组f[i][j]，表示从i开始向后2^j个数的最大值.处理完毕之后我们要求[L, R]内的最大值就可以比较f[L][k]以及f[R- 2^K+1][k]的最大值，也就是从L开始到L+2^K的最大值和R之前2^K区间的最大值比较，虽然是有重复区域，但是只要求最大值所以并没有影响；

定义：f[i][j]表示i到i+2^j-1这段区间的最大值。

预处理：f[i][0]=a[i]。即i到i区间的最大值就是a[i]。

状态转移：将f[i][j]平均分成两段，一段为f[i][j-1]，另一段为f[i+2^(j-1)][j-1]。

两段的长度均为2^j-1。f[i][j]的最大值即这两段的最大值中的最大值。

得到f[i][j]=max(f[i][j-1],f[i+2^(j-1)][j-1])。

模板：

以POJ3264为例：

#include <map>
#include <set>
#include <stack>
#include <cmath>
#include <queue>
#include <bitset>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <fstream>
#include <cstdlib>
#include <sstream>
#include <cstring>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;
//#define   maxd          1010
#define lson l , m , rt << 1
#define rson m + 1 , r , rt << 1 | 1
#define   mc(x, y)     memcpy(x, y, sizeof(x))
#define   ms(x,y)      memset(x,y,sizeof(x))
#define   rep(i,n)      for(int i=0;i<(n);i++)
#define   repf(i,a,b)   for(int i=(a);i<=(b);i++)
#define   PI           pair<int,int>
//#define   mapp            make_pair
#define   FI            first
#define   SE            second
#define   IT            iterator
#define   PB            push_back
#define   Times         10
typedef   long long     ll;
typedef   unsigned long long ull;
typedef   long double   ld;
typedef   pair<int,int > pce;
//#define N 100
const double eps = 1e-10;
const double  pi = acos(-1.0);
const  ll    mod = 1e9+7;
const  int   inf = 0x3f3f3f3f;
//const  ll    INF = (ll)1e18+300;
const int   maxd = 50500 + 10;

int n;
int Q;
int ac[maxd];
int max_[maxd][20];
int min_[maxd][20];

void RMQ() {
    for (int i = 1; i <= n; i++) {
        max_[i][0] = min_[i][0] = ac[i];
    }
    for (int j = 1; j <= 20; j++) {
    `    /*外层循环要从j开始，因为ST表我们是通过j-1来更新j的，如果是i为外层循环，第一次我们计算出了：
            1, 2, 4, 8,...这些位置，而下一次更新的时候，我们要的是max([1,4], [5, 8])但是我们没有计算出这两个区间的最大值，因为一次层的时候
            没有计算3, 5, 6, 7,...的值，因此必须j为第一层循环。

        */
        for (int i = 1; i <= n; i++) {
            if(i + (1<<j) - 1 <= n) {
                max_[i][j] = max(max_[i][j-1], max_[i + (1<<(j-1))][j - 1]);
                min_[i][j] = min(min_[i][j-1], min_[i + (1<<(j-1))][j - 1]);
            }
        }
    }
}

int main() {
    while(~scanf("%d%d", &n, &Q)) {
        for (int i = 1; i <= n; i++) {
            scanf("%d", &ac[i]);
        }
        RMQ();
        for (int i = 1; i <= Q; i++) {
            int a, b;
            scanf("%d%d", &a, &b);
            int k = (int)(log(b - a + 1)/log(2.0));
            /*查询操作就像图片中所示，取超过（r - l）/2的长度，分区间比较，而我们就可以取log2并向上取整.
            */
            int max_ans = max(max_[a][k], max_[b - (1<<k) + 1][k]);
            int min_ans = min(min_[a][k], min_[b - (1<<k) + 1][k]);
            cout << max_ans - min_ans << endl;
        }
    }
}