Изменения

|definition='''Каталог''' (англ. ''catalog''){{---}} упорядоченный массив из элементов, на которых введено [[Отношение порядка|отношение порядка]]. В данной статье предполагается, что массив упорядочен по неубыванию.

| <center><tex> k </tex> [[Целочисленный двоичный поиск|бинарных поисков ]] </center> || <center><tex> \ThetaO(n) </tex></center> || <center><tex> O(k \log n) </tex></center>

* 1) # Мы можем проводить ссылки из каталога номер <tex> i </tex> в <tex> (i + 1) </tex>-ый й каталог таким образом, что разница между элементами, соединенными ссылками минимальна, что, очевидно, в некоторых случаях уменьшит время поискаэлемента в следующем каталоге, так как область поиска сожмется. <br>* 2) # Мы можем для оптимизации пункта 1 создать модифицированные каталоги <tex> M_i </tex>, где <tex> i </tex>-ый й каталог будет представлять каталог <tex> C_i </tex> слитый с <tex> M_{i + 1} </tex>

Введем определения: {{Определение|definition='''Подставной элемент''' {{---}} Будем называть подставным элементом такой элемент каталога <tex> M_i </tex>, который пришел из каталога <tex> M_{i + 1} </tex>. А также Сами каталоги <tex> M_i </tex> будем называть '''модифицированными каталогами'''.}}[[Файл:FCT_pic2.jpg|600px800px|left|thumb|Построение модифицированных каталогов]]'''Первый этап построения''':

'''Второй этап построения''':* В каждом ''модифицированном каталоге'' для каждого элемента заведем две ссылки. Для ''неподставных элементов'' это будут ссылки на максимальный''подставной элемент'' меньше текущего и на минимальный любого типа больше текущего. Если ''элемент подставной'', то ссылки будут на минимальный ''неподставной элемент" больше текущего и на максимальный ''неподставной элемент'' меньше текущего. Назовем их ссылками влево и вправо.

<br> Для наглядности заведем таблицу, где в <tex>i</tex>-ой й строке <tex> j </tex>-ая я ячейка будет окрашена в зеленый цвет, если она присутствует в каталоге <tex> C_i </tex>. Тогда результатом построения будет таблица, которая представлена на рисунке.Для упрощения рисунка ссылки вправо из неподставных элементов не были отображены, их следует воспринимать как следующий справа от рассматриваемого элемент в ряду таблицы любого цвета.

'''struct node ''' Node: element '''T''' key node '''Node''' left, right, down '''bool ''' is_alien

'''for''' i = k - 1 '''todownto''' 1 '''int''' pointer_in_C = 1 <font color=green>// указатель на самый левый элемент каталога C[i], который еще не рассмотрели </font> '''int''' pointer_in_next_M = 1 <font color=green>// указатель на самый левый элемент каталога M[i + 1], который еще не рассмотрели </font> '''int''' pointer_in_M = 1 <font color=green>// указатель на самый левый элемент каталога M[i], в который будем добавлять элемент </font> node '''Node''' last_non_alien = NULL <tex> \varnothing </tex> <font color=green>// указатель на последний ''неподставной элемент'' для текущей позиции </font> node '''Node''' last_alien = NULL <tex> \varnothing </tex> <font color=green>// указатель на последний ''подставной элемент'' для текущей позиции</font> '''while'''('''true''') '''if''' (pointer_in_next_M > M[i + 1].size && '''and''' pointer_in_C > C[i].size)

'''if''' (pointer_in_next_M > M[i + 1].size || '''or''' M[i + 1][pointer_in_next_M] <tex> \geqslant </tex>= C[i][pointer_in_C])

M[i][pointer_in_M].key left= last_alien

M[i][pointer_in_M].is_alien = '''true'''

lst_alien last_alien = M[i][pointer_in_M]

last_non_alien = (!'''if''' '''not''' M[i][M[i].size].is_alien ? last_non_alien = M[i][M[i].size] : NULL) '''else''' last_non_alien = <tex> \varnothing </tex> <font color=green>// теперь это last_non_alien указатель на первый справа ''неподставной элемент'' для текущей позиции</font> '''for''' j = M[i].size - 1 '''todownto''' 1 '''if''' (M[i][j].is_alien)

Из построения понятно, что мы тратим <tex> O(n_k) </tex> на построение последнего каталога, <tex> O(n_{k-1} + n_k / 2) </tex> на построение предпоследнего и т.д. Пусть <tex> p = 2^{\lfloor log n \rfloor- 1} </tex>. Тогда получаем оценку <tex> O(n_k (1 + 1/2 + 1/4 + ... \dots + 1/p) + n_{k - 1} (1 + 1/2 + 1/4 + ... \dots + 1/(p/2)) + ... \dots + n_1 ) </tex> <tex> = O(2 n_k + 2 n_{k -1} + ... \dots n_1) = O(n) </tex> памяти. По алгоритму понятно, что такая же оценка верна и для времени на предподсчет.

*В первом каталоге ответ на запрос найдем с помощью [[Целочисленный двоичный поиск|бинарного поиска ]] по <tex> C_1 M_1 </tex>. Пусть ответом для этого каталога будет ячейка <tex> cell </tex>, тогда если <tex> cell </tex> {{---}} подставная вершина, то перейдем по ссылке влево.

**Для того, чтобы перейти в новый ''модифицированный каталог'' мы перейдем из <tex> cell </tex> по ссылке влево, чтобы попасть в ''подставную вершину'', а потом из нее перейдем по ссылке вниз, чтобы попасть в следующий каталог. **Если теперь <tex> cell </tex> {{---}} ''неподставная вершина'', то нам достаточно рассмотреть двух ее соседей справа в <tex> C_i </tex>, так как <tex> cell.key <= \leqslant x </tex>, а каждая вторая ячейка из <tex> M_i </tex> попадает в <tex> M_{i - 1} </tex>, т.е. мы бы встретили ее ранее и перешли мы вниз по ней, но этого не случилось.**Обновив <tex> cell </tex> максимальным из подходящих значений нужно проверить, является ли она ''подставным элементом'', если да, то перейдем по ссылке влево, попав в ответ для текущего каталога, иначе это и будет ответ.

node '''Node''' cell = binary_search(M[1], x) '''if (''' cell.is_alien)

ans[1] = cell.key; <font color=green>// ans[i] - ответ на текущий запрос для каталога С[i] </font>

cell = cell.left.down '''if''' (cell.right <= tex> \leqslant </tex> x) <font color=green>// Попытка сдвинуться к большему элементу </font>

'''if''' (cell.right <= tex> \leqslant </tex> x) <font color=green>// Попытка сдвинуться к большему элементу </font> cell = cell.right <font color=green>// Замечание: по построению, если мы стоим в ''неподставном элементе'', то при сдвиге вправо мы можем оказаться в элементе любого типа</font> '''if''' (cell.is_alien) <font color=green>// Для этого есть проверка </font>

Как можно видеть, для того, чтобы найти ответ для первого каталога необходимо сделать один [[Целочисленный двоичный поиск|бинарный поиск]], что требует <tex> O(\log n) </tex> времени, после чего необходимо <tex> O(k) </tex> времени, чтобы ответить на запрос для всех остальных каталогов. Суммарное время работы <tex> O(\log n + k) </tex>.<br>  ====Примеры ответа на запрос===={| border="0" cellpadding="5" cellspacing="0" center|[[Файл:FCT_pic3.jpg|600px|center|thumb|Ответ на запрос <tex> x = 9 </tex>]]|[[Файл:FCT_pic4.jpg|600px|center|thumb|Ответ на запрос <tex> x = 6 </tex>]]|}Рассмотрим, как будет происходить ответ на запрос для <tex> x = 9 </tex> (картинка справа) и для <tex> x = 6 </tex> (картинка слева). Каталоги взяты из примера для [[Техника частичного каскадирования#Построение | построения]]. Оставлены только ссылки, по которым осуществляется переход, а элементы пронумерованы в порядке обхода.<br> ===Дополнительно===Данная техника может использоваться для ускорения некоторых алгоритмов, где требуется ответить на запрос на отрезке <tex> [L, R]</tex>, где <tex> L, R \in R^n, n \in \mathbb N </tex>. Однако иногда наблюдается замедление, о чем можно почитать <ref>[http://codeforces.com/blog/entry/21892?locale=en тутFractional cascading is in fact slow? ifsmirnov's blog]</ref>. 

*[[Пересечение прямоугольника с множеством непересекающихся отрезков (segment tree)]]*[[Перечисление точек в произвольном прямоугольнике за n * log ^(d - 1) n (range tree)]]== Примечания ==<references />