第38題
(RMQ 區(qū)間最值問題)給定序列a0,?,an-1,和m次詢問,每次詢問給定l,r,求max{al,?,ar}。
為了解決該問題,有一個算法叫the Method of Four Russians,其時間復雜度為O(n+m),步驟如下:
1) 建立Cartesian(笛卡爾)樹,將問題轉(zhuǎn)化為樹上的LCA(最近公共祖先)問題。
2) 對于LCA問題,可以考慮其Euler序(即按照DFS過程,經(jīng)過所有點,環(huán)游回根的序列),即求Euler序列上兩點間一個新的RMQ問題。
3) 注意新的問題為±1 RMQ,即相鄰兩點的深度差一定為1。
下面解決這個±1 RMQ問題,“序列”指Euler序列:
1) 設t為Euler序列長度。取。將序列每b個分為一大塊,使用ST表(倍增表)處理大塊間的RMQ問題,復雜度。
2)(重點)對于一個塊內(nèi)的RMQ問題,也需要O(1)的算法。由于差分數(shù)組2b?1種,可以預處理出所有情況下的最值位置,預處理復雜度O(b2b),不超過O(n)。
3) 最終,對于一個查詢,可以轉(zhuǎn)化為中間整的大塊的RMQ問題,以及兩端塊內(nèi)的RMQ問題。
試補全程序。
#include <iostream>
#include <cmath>
using namespace std;
const int MAXN = 100000, MAXT = MAXN << 1;
const int MAXL = 18, MAXB = 9, MAXC = MAXT / MAXB;
struct node {
int val;
int dep, dfn, end;
node *son[2]; // son[0], son[1] 分別表示左右兒子
} T[MAXN];
int n, t, b, c, Log2[MAXC + 1];
int Pos[(1 << (MAXB - 1)) + 5], Dif[MAXC + 1];
node *root, *A[MAXT], *Min[MAXL][MAXC];
void build() { // 建立 Cartesian 樹
static node *S[MAXN + 1];
int top = 0;
for (int i = 0; i < n; i++) {
node *p = &T[i];
while (top && S[top]->val < p->val)
①;
if (top)
②;
S[++top] = p;
}
root = S[1];
}
void DFS(node *p) { // 構(gòu)建 Euler 序列
A[p->dfn = t++] = p;
for (int i = 0; i < 2; i++)
if (p->son[i]) {
p->son[i]->dep = p->dep + 1;
DFS(p->son[i]);
A[t++] = p;
}
p->end = t - 1;
}
node *min(node *x, node *y) {
return ③ ? x : y;
}
void ST_init() {
b = (int)(ceil(log2(t) / 2));
c = t / b;
Log2[1] = 0;
for (int i = 2; i <= c; i++)
Log2[i] = Log2[i >> 1] + 1;
for (int i = 0; i < c; i++) {
Min[0][i] = A[i * b];
for (int j = 1; j < b; j++)
Min[0][i] = min(Min[0][i], A[i * b + j]);
}
for (int i = 1, l = 2; l <= c; i++, l <<= 1)
for (int j = 0; j + l <= c; j++)
Min[i][j] = min(Min[i - 1][j], Min[i - 1][j + (l >> 1)]);
}
void small_init() { // 塊內(nèi)預處理
for (int i = 0; i <= c; i++)
for (int j = 1; j < b && i * b + j < t; j++)
if (④)
Dif[i] |= 1 << (j - 1);
for (int S = 0; S < (1 << (b - 1)); S++) {
int mx = 0, v = 0;
for (int i = 1; i < b; i++) {
⑤;
if (v < mx) {
mx = v;
Pos[S] = i;
}
}
}
}
node *ST_query(int l, int r) {
int g = Log2[r - l + 1];
return min(Min[g][l], Min[g][r - (1 << g) + 1]);
}
node *small_query(int l, int r) { // 塊內(nèi)查詢
int p = l / b;
int S = ⑥;
return A[l + Pos[S]];
}
node *query(int l, int r) {
if (l > r)
return query(r, l);
int pl = l / b, pr = r / b;
if (pl == pr) {
return small_query(l, r);
} else {
node *s = min(small_query(l, pl * b + b - 1), small_query(pr * b, r));
if (pl + 1 <= pr - 1)
s = min(s, ST_query(pl + 1, pr - 1));
return s;
}
}
int main() {
int m;
cin >> n >> m;
for (int i = 0; i < n; i++)
cin >> T[i].val;
build();
DFS(root);
ST_init();
small_init();
while (m--) {
int l, r;
cin >> l >> r;
cout << query(T[l].dfn, T[r].dfn)->val << endl;
}
return 0;
}
① 處應填( )