POJ 3264 Balanced Lineup 最基本的RMQ

其实这个题早就过了，只是当时为了偷懒拿线段树做的。当时肤浅的认为RMQ问题线段树都能胜任，现在才发现毕竟土洋。

不过线段树做多了，RMQ就显得很容易理解了。

其实RMQ很多算法，不过最通用的应属ST算法，总体思想就是一个很简单的dp。

以求区间最大值为例，数组dp[ i ][ j ]记录的是区间[ i , i + (1<<j)-1]内的最大值。

也就是说是从第 i 个元素开始，连续的 (1<<j)-1 个数内的最大值。

不难发现，dp[ i ][ 0 ]  == num[ i ] （num为存储输入数据的数组）。

而连续的 (1<<j)-1 个元素又可以分成左右两段，对于[ i , i + (1<<j)-1] 来说  就是 [ i , i + (1<<(j-1))-1] 和 [ i + (1<<(j-1)) , i + (1<<j)+1];

显然 区间[ i , i + (1<<j)-1]内的最值是 [ i , i + (1<<(j-1))-1] 和 [ i + (1<<(j-1)) , i + (1<<j)+1] 两个区间的最值中较大的一个。

即 dp[ i ][ j ] = max( dp[ i ][ j-1 ]  ,  dp[ i +(1<<(j-1))][j-1] )。

以上即为该算法的主要思路，而在编码实现过程中还是要注意到那个一万年不许变的问题——边界问题。

当 i +(1<<(j-1))  > n  (n为数据量)时，dp[ i ][ j ] =  dp[ i ][ j-1 ] ，因为此时没有必要分成两段了。

上面为初始化的思路，下面说一下查询的过程。

对于区间[ L , R ]内的最值可以视为[L,s1] , [s1+1,s2],[ s2+1,s3].......[sn+1,R]这一串连续的小区间的最值中较大的一个。

而 si 显然要满足 L <= si <= R,且每一段小区间的长度为 2^j (j为非负整数)。

即：

1. 每次都要找一个尽量大的满足上述要求的 j ，然后查询[L,L+(1<<j)-1]区间内的最值。

2. 然后更新L ，L += （1<<j）。

3.若L 〉 R，求值过程结束，否则返回第一步。

上面这些即为查询过程。

下面附上AC代码，效率略低唉。。

#include <iostream>
#include <algorithm>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <queue>
#include <stack>

#pragma comment(linker, "/STACK:1024000000");
#define LL long long int

using namespace std;

int dpmin[50010][20],dpmax[50010][20];
int num[50010];

int main()
{
    int i,j,m,n,l,r,Max,Min;
    while(scanf("%d %d",&n,&m) != EOF)
    {
        for(i = 1;i <= n; ++i)
        {
            scanf("%d",&num[i]);
        }

        for(i = 1;i <= n; ++i)
        {
            dpmin[i][0] = num[i];
            dpmax[i][0] = num[i];
        }

        for(j = 1;(1<<(j-1)) <= n; ++j)
        {
            for(i = 1;i <= n; ++i)
            {
                if(i+(1<<(j-1)) <= n)
                {
                    dpmax[i][j] = max(dpmax[i][j-1],dpmax[i+(1<<(j-1))][j-1]);
                    dpmin[i][j] = min(dpmin[i][j-1],dpmin[i+(1<<(j-1))][j-1]);
                }
                else
                {
                    dpmax[i][j] = dpmax[i][j-1];
                    dpmin[i][j] = dpmin[i][j-1];
                }
            }
        }

        while(m--)
        {
            scanf("%d %d",&l,&r);
            Max = num[l],Min = num[l];
            while(l <= r)
            {
                for(j = 0;l + (1<<j)-1 <= r; ++j)
                    ;
                --j;
                Max = max(Max,dpmax[l][j]);
                Min = min(Min,dpmin[l][j]);
                l += (1<<j);
            }
            printf("%d\n",Max-Min);
        }
    }
    return 0;
}