Poj 3264 Balanced Lineup (数据结构_RMQ)

题目链接：http://poj.org/problem?id=3264

题目大意：给定一个序列，n个数，给出m个询问，每次询问[li,ri]区间的最大值与最小值的差。

解题思路：一直想花点时间补下数据结构方面的知识，但因为状态dp和树形dp较难，花了大把时间。想放弃又害怕痛苦！到现在也没能搞定这两个家伙，今天早晨终于归零，开始学RMQ.

    这题应该是有两种做法，建两可线段树复杂度o(nlogn)，o(logn)查询,总复杂度O((n+m)logn，但这应该不会超时，没去敲。我用的是St算法，O(nlogn)初始化一个dp[MAX][MIN]数组，O（1）查询，St算法的优越之处就在于这个查询速度。

    那么它是怎么实现的呢？其实就是DP,运用倍增思想，把区间压缩。dp[i][j]表示从位置i开始（包括i自己）长度为1<<j的区间的最值，那么dp[i][0] = arr[i]，这是显然的，那么dp[i][j]该怎么转移呢？dp[i][j] = max（dp[i][j-1],dp[i+(1<<(j-1))][j-1]),把长度为1<<j的区间分成两段长度为1<<(j-1)的区间，取其中的最大值。长度长的从长度短的转移过来，所以在编码中要先写j的循环再写i循环。初始完成后就坐等高效查询。当需要查询[l,r]之间最值的时候，先计算出区间长度len，然后求k = log2len(即log（r-l+1)/log(2))，最后max(dp[l][k],dp[r-(1<<k)+1][k])即是答案。为什么可以用dp[l][k]和dp[r-(1<<k)+1][k]表示一个区间？他们表示的是两个至少一个交点的区间，如果大家不是很明白，可以想下查询【1，8】或【1，10】的情况，前者8个交点，后者6个交点。咳，口齿不清，表达得不清楚，一句话概括下--那是必然的。

测试数据：

6 3
1
7
3
4
2
5
1 5
4 6
2 2





代码：
#include <stdio.h>
#include <string.h>
#include <math.h>
#define MIN 20
#define MAX 51000
#define min(a,b) (a)<(b)?(a):(b)
#define max(a,b) (a)>(b)?(a):(b)


int n,m,arr[51000];

//模板
struct RMQ {

    int n,mmax[MAX][MIN],mmin[MAX][MIN];
    void Create(int flag);
    int Query(int x,int y,int flag);
}rmq;
void RMQ::Create(int flag) {
    //flag = 0建立最大rmq，flag = 1最小
    int i,j,k;

    if (flag == 0) {

        for (i = 1; i <= n; ++i)
            mmax[i][0] = arr[i];
        for (j = 1; (1<<j) <= n; ++j)
            for (i = 1; i+(1<<j)-1 <= n; ++i)
                mmax[i][j] = max(mmax[i][j-1],mmax[i+(1<<(j-1))][j-1]);
    }
    else {

        for (i = 1; i <= n; ++i)
            mmin[i][0] = arr[i];
        for (j = 1; (1<<j) <= n; ++j)
            for (i = 1; i+(1<<j)-1 <= n; ++i)
                mmin[i][j] = min(mmin[i][j-1],mmin[i+(1<<(j-1))][j-1]);
    }
}
int RMQ::Query(int l,int r,int flag){
    //flag = 0查询最大rmq，flag = 1最小
    int k = (int)(log(r-l+1.0)/log(2.0));
    if (flag == 0)
        return max(mmax[l][k],mmax[r-(1<<k)+1][k]);
    else
        return min(mmin[l][k],mmin[r-(1<<k)+1][k]);
}


int main()
{
    int i,j,k;


    scanf("%d%d",&n,&m);
    for (i = 1; i <= n; ++i)
        scanf("%d",&arr[i]);


    rmq.n = n;
    rmq.Create(0);
    rmq.Create(1);
    for (i = 1; i <= m; ++i) {

        scanf("%d%d",&j,&k);
        printf("%d\n",rmq.Query(j,k,0) - rmq.Query(j,k,1));
    }
}