【RMQ_ST算法】POJ 3264 Balanced Lineup

RMQ(Range Minimum/Maximum Query):区间最值查询

(最近才知道原来RMQ是一大类问题，就是区间最值查询的问题，可以用线段树解决)

另外就是用下面写的ST(Sparse Table[稀疏表])算法【ST算法用了倍增的思想】：

ST是一种在线算法，预处理部分O(nlogn)，查询部分O(1)

预处理部分：

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]);

预处理部分的复杂度：O(nlogn)

释意：

1. max_[ i ][ j ]，从区间的第 i 的位置开始，长度为2 ^ j 的区间上的最大值。
2. 我们要求区间的最大值，可以将区间从中间分成两个子区间，然后取两个子区间的最大值中的较大者。区间长度一分为二，后边第二个下标变为 j-1 是很好理解的。
3. max_[ i ] [ j-1 ]表示的区间是[ i , i+2 ^ (j-1)-1 ]
4. 而max_[ i + (1 << j-1) ][ j-1 ])表示的区间是[ i+2 ^ (j-1) , i+2 ^ j-1]
5. (可能有些人跟我一样不理解这个区间，想个例子其实超级好理解：比如从区间的1位置开始，长度为8，分成两个区间其实是[ 1 , 4 ]和[ 5 , 8 ]，和上面的式子联系一下就明白了)

查询部分：

int query_max(int l,int r)
{
    int k=log2(r-l+1);//向下取整
    return max(max_[l][k],max_[r-(1<<k)+1][k]);
}

查询部分复杂度：O(1)
释意：

1. 我们取 k 为区间长度以2为底的指数(向下取整)
2. 如果区间长度刚好是2的次幂，max_[ l ][ k ]和max_[ r-(1<<k)+1 ][ k ]两者表示的是一样的区间，并且已经预处理过了，我们是不是可以直接返回已经预处理好的值？
3. 如果不是呢，我们就向下取整log2( r-l+1 )，然后将区间分成有重合部分的长度都为 2 ^ k 的两个区间：[ l , l+2 ^ k-1 ] [ r-2^k+1 , r ]这时候我们手动分的这两个区间是已经预处理好的，是不是可以直接返回其中的较大值即可？

区间最小值与上所述同理

例题+代码：

POJ 3264 Balanced Lineup

• 题意：区间最小值
• 思路：模板

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <cmath>
#include <cstring>
#include <string>
#include <vector>
#include <set>
#include <stack>
#include <list>
#include <map>
#define INF 0x3f3f3f3f
#define P(x) x>0 ? x : 0
#define MID l + r >> 1
#define lsn rt << 1
#define rsn rt << 1 | 1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr

using namespace std;
typedef long long ll;
typedef vector<int>:: iterator VITer;
const int maxN=5e4+5;
int max_[maxN][20],min_[maxN][20];
int N,Q;

void RMQ_pre()
{
    for(int j=1;(1<<j)<=N;j++)
    {
        for(int i=1;i + (1 << j-1)<=N;i++)
        {
            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 query_min(int l,int r)
{
    int k=log2(r-l+1);//向下取整
    return min(min_[l][k],min_[r-(1<<k)+1][k]);
}

int query_max(int l,int r)
{
    int k=log2(r-l+1);
    return max(max_[l][k],max_[r-(1<<k)+1][k]);
}

int main()
{
    while(~scanf("%d%d", &N, &Q))
    {
        for(int i=1;i<=N;i++)
        {
            scanf("%d",&max_[i][0]);
            min_[i][0]=max_[i][0];
        }
        RMQ_pre();
        while(Q--)
        {
            int a,b;
            scanf("%d%d", &a, &b);
            printf("%d\n",query_max(a,b)-query_min(a,b));
        }
    }
    return 0;
}