[POJ3264]Balanced Lineup(RMQ, ST算法)

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

典型RMQ，这道题被我鞭尸了三遍也是醉了…这回用新学的st算法。

st算法本身是一个区间dp，利用的性质就是相邻两个区间的最值的最值一定是这两个区间合并后的最值，这条性质决定了这个dp子问题的重叠。可以利用这个性质预处理出这张表，只不过步长是２的幂次。

查询的时候也是如此，但是未必会精准地选中两个区间，不要紧，因为两个区间重叠的部分也会被自动算在求最值的内部。这个时候如果算的是区间和的话，要减去这一部分。（区间和的话直接用前缀和不就好了嘛）

/*
━━━━━┒ギリギリ♂ eye！
┓┏┓┏┓┃キリキリ♂ mind！
┛┗┛┗┛┃＼○／
┓┏┓┏┓┃ /
┛┗┛┗┛┃ノ)
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┃┃┃┃┃┃
┻┻┻┻┻┻
*/
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
using namespace std;
#define fr first
#define sc second
#define cl clear
#define BUG puts("here!!!")
#define W(a) while(a--)
#define pb(a) push_back(a)
#define Rint(a) scanf("%d", &(a))
#define Rll(a) scanf("%lld", &a)
#define Rs(a) scanf("%s", a)
#define Cin(a) cin >> a
#define FRead() freopen("in", "r", stdin)
#define FWrite() freopen("out", "w", stdout)
#define Rep(i, len) for(int i = 0; i < (len); i++)
#define For(i, a, len) for(int i = (a); i < (len); i++)
#define Cls(a) memset((a), 0, sizeof(a))
#define Clr(a, x) memset((a), (x), sizeof(a))
#define Full(a) memset((a), 0x7f, sizeof(a))
#define lrt rt << 1
#define rrt rt << 1 | 1
#define pi 3.14159265359
#define RT return
#define lowbit(x) x & (-x)
#define onenum(x) __builtin_popcount(x)
typedef long long LL;
typedef long double LD;
typedef unsigned long long Uint;
typedef pair<int, int> pii;
typedef pair<LL, LL> pLL;
typedef pair<string, LL> psi;
typedef map<string, LL> msi;
typedef vector<LL> vi;
typedef vector<LL> vl;
typedef vector<vl> vvl;
typedef vector<bool> vb;

const int maxn =50050;
int n, q;
int h[maxn];
int dp[maxn][30][2];

void st() {
    for(int i = 1; i <= n; i++) dp[i][0][0] = dp[i][0][1] = h[i];
    for(int j = 1; (1 << j) <= n; j++) {
        for(int i = 1; i + (1 << j) - 1 <= n; i++) {
            dp[i][j][0] = min(dp[i][j-1][0], dp[i+(1<<(j-1))][j-1][0]);
            dp[i][j][1] = max(dp[i][j-1][1], dp[i+(1<<(j-1))][j-1][1]);
        }
    }
}

int query(int l, int r) {
    int j = 0;
    while((1 << (j + 1)) <= r - l + 1) j++;
    return max(dp[l][j][1], dp[r-(1<<j)+1][j][1]) - min(dp[l][j][0], dp[r-(1<<j)+1][j][0]);
}

int main() {
//    FRead();
    while(~scanf("%d%d",&n,&q)) {
        For(i, 1, n+1) Rint(h[i]);
        st();
        int l, r;
        W(q) {
            scanf("%d%d",&l,&r);
            printf("%d
",query(l,r));
        }
    }
    RT 0;
}