Решение RMQ с помощью разреженной таблицы

Разреженная таблица (англ. sparse table) позволяет решать задачу online static RMQ за [math]O(1)[/math] на запрос, с предподсчётом за [math]O(N \log N)[/math] и использованием [math]O(N \log N)[/math] памяти.

Постановка задачи RMQ

Дан массив [math]A[1..N][/math] действительных чисел. Поступают запросы вида [math](l, r)[/math], на каждый запрос требуется найти минимум в массиве [math]A[/math], начиная с позиции [math]l[/math] и заканчивая позицией [math]r[/math].

Разреженная таблица

Разреженная таблица — двумерная структура данных [math]ST[i, j][/math], для которой выполнено следующее: [math]ST[i][j]=\min\left(A[i], A[i+1], ..., A[i+2^{j}-1]\right),\quad j \in [0 .. \log N][/math]. Иначе говоря, в этой таблице хранятся минимумы на всех отрезках, длины которых равны степеням двойки. Объём, занимаемый таблицей, равен [math]O(N \log N)[/math], и заполненными являются только те элементы, для которых [math]i+2^j \le N [/math].

Простой метод построения таблицы заключён в следующем реккурентном соотношении:

[math] ST[i][j] = \left\{ \begin{array}{rcl} \min\left(ST[i][j-1], ST[i+2^{j-1}][j-1]\right), j \gt 0 \\ A[i], j = 0 \\ \end{array} \right. [/math] .

Идемпотентность

Такая простота достигается за счет идемпотентности операции минимум: [math]\min(a, a)=a[/math]. Это один из ключевых моментов этого метода, так как она позволяет нам корректно считать минимум в области пересечения отрезков.

Таким образом мы получаем целый класс задач, которые могут решаться разреженной таблицей: вместо минимума может быть любая идемпотентноя бинарная операция [math] f(x, x) = x, x \in X [/math], где [math]X[/math] — некое множество, над которым задан массив [math]A[/math].

Применение к задаче RMQ

Дан запрос [math](l, r)[/math]. По нему найдём [math]k = \max \{j| 2^j \le r - l + 1\}[/math], т.е. логарифм длины запрашиваемого отрезка. Заметим, что [math]\min(A[l], A[l+1], ..., A[r]) = \min\left(ST[l][k], ST[r-2^k+1][j]\right)[/math]. Заметим [math]k[/math] зависит лишь от длины отрезка. Предпосчитать эту величину за [math]O(N\log N)[/math] можно введением функции [math]fl[l] = k[/math], для которой верно [math]fl[1] = 0, fl[x] = fl[\lfloor \frac{x}{2}\rfloor] + 1[/math]. Теперь мы можем находить [math]k[/math] за [math]O(1)[/math]. Таким образом, ответ на запрос происходит за константное время.

См. также

• Сведение задачи LCA к задаче RMQ
• Алгоритм Фарака-Колтона и Бендера
• Сведение задачи RMQ к задаче LCA

Источники

• Bender, M.A., Farach-Colton, M. et al. — Lowest common ancestors in trees and directed acyclic graphs. — J. Algorithms 57(2) (2005) — с. 75–94.