Изменения

'''Рандомизированное бинарное дерево поиска ''' (англ. ''Randomized binary search tree, RBST'', RBST) {{---}} структура данных, реализующая бинарное дерево поиска.

Как известно, можно подобрать такую последовательность операций с [[Дерево поиска, наивная реализация|бинарным деревом поиска в наивной реализации]], что его глубина будет пропорциональна количеству ключей, а следовательно запрос будет операции будут выполняться за <tex>O(n)</tex>. Поэтому, если поддерживать инвариант "случайности" в дереве, то можно добиться того, что математическое ожидание глубины дерева будет небольшим.

# Если <tex>T</tex> пусто, то оно является '''рандомизированным бинарным деревом поиска'''.# Если <tex>T</tex> непусто (содержит <tex>n</tex> вершин, <tex>n > 0</tex>), то <tex>T</tex> {{---}} '''рандомизированное бинарное дерево поиска ''' тогда и только тогда, когда его левое и правое поддеревья (<tex>L</tex> и <tex>R</tex>) оба являются '''RBST''', а также выполняется соотношение<tex>P[size(L) = i] = \fracdfrac{1}n, i = 1..n</tex>.

Из определения следует, что каждый ключ в RBST размера <tex>n </tex> может оказаться корнем с вероятностью <tex>\dfrac{1}{n}</ntex>.

(Похожие идеи используются в [[Декартово дерево| декартовом дереве]], поэтому во многих русскоязычных ресурсах термин '''рандомизированное бинарное дерево поиска''' используется как синонимическое название декартового дерева и [[Декартово дерево по неявному ключу|декартового дерева по неявному ключу]]).

Рассмотрим рекурсивный алгоритм вставки ключа <tex>x</tex> в RBST, состоящее из <tex>n</tex> вершин. С вероятностью <tex dpi = "150">\fracdfrac{1}{n+1}</tex> вставим ключ в корень дерева(разделим дерево по данному ключу и подвесим получившиеся деревья к новому корню), используя процедуру ''insert_at_root''<tex>\mathrm{insertAtRoot}</tex>. С вероятностью <tex dpi = "150">1 - \fracdfrac{1}{n+1} = \fracdfrac{n}{n+1}</tex> вставим его в правое поддерево, если он больше корня, или в левое поддерево, если меньше. Ниже приведён псевдокод процедуры вставки ''<tex>\mathrm{insert''}</tex>, процедуры ''insert_at_root''<tex>\mathrm{insertAtRoot}</tex>, а также процедуры ''<tex>\mathrm{split(k)''}</tex>, разбивающей дерево на два поддерева, в одном из которых все ключи строго меньше <tex>k</tex>, а в другом больше, либо равны; приведена достаточно очевидная рекурсивная реализация. (через <tex>Node </tex> обозначен тип вершины дерева, дерево представляется как указатель на корень).

'''Node''' insert(t : '''insertNode''' (T, x: '''T'''): '''int''' r = '''random'''(0.<tex>\dots</tex> t.size(T)) '''if''' (r == n)t.size T t = insert_at_rootinsertAtRoot(Tt, x) '''if''' (x < roott.key) T t = insert(Tt.left, x)

T t = insert(Tt.right, x) t.size = 1 + t.size '''return''' Tt

Заметим, что если дерево пусто, то ''<tex>\mathrm{insert'' }</tex> с вероятностью 1 делает <tex>x</tex> корнем.

// вставляет ключ x в дерево T '''Node''' insertAtRoot(t : '''insert_at_rootNode''' (T, x: '''T''') : <font color="green">//вставляем в дерево t ключ x</ создать пустые L и Rfont> L <l, r> = RBST() R = RBST() split(Tt, x, L, R) // создать пустое T T = RBST() Tt.key = x Tt.left = Ll Tt.left right = Rr '''return''' Tt

// разделяет дерево T по '''<Node''', '''Node>''' split(t : '''Node''', x // результат: деревья L и R '''splitT''' (T, ): <font color="green"> // разделяет дерево t по x, L, R)результат {{---}} пара деревьев r и l</font> '''if''' (t.size(T) == 0) // создать пустые L и R L = RBST() R = RBST() '''elsereturn''' <''null'', ''null'' > '''else if''' (x < Tt.key) R <l, r> = T split (Tt.left, x, L, R) t.left = r t.size = 1 + t.left).size + t.right.size r = t '''return''' <l, r>

L <l, r> = T split (Tt.right, x, L) t.right = l t.size = 1 + t.left.size + t.right.size l = t '''return''' <l, R)r>

|statement= Пусть после операции ''<tex>\mathrm{split'' }</tex> от дерева <tex>T</tex> по ключу <tex>x</tex> были получены деревья <tex>T_{L}</tex> и <tex>T_{R}</tex>. Тогда если <tex>T</tex> было рандомизированным бинарным деревом поиска, содержащим множество ключей <tex>K</tex>, то деревья <tex>T_{L}</tex> и <tex>T_{R}</tex> {{---}} рандомизированные бинарные деревья поиска, содержащие соответственно множества ключей <tex>K_{L} = \{y \in K | \mid y < x\}</tex> и <tex>K_{R} = \{y \in K | \mid y > x\}</tex>.

Пусть <tex>n > 0</tex>, и лемма верна при всех меньших размерах дерева.. Пусть также <tex>y = T.key, L = T.left, R = T.right</tex>. Если <tex>x > y</tex>, то <tex>y</tex> {{---}} корень <tex>T_{L}</tex>, <tex>L</tex> {{---}} левое поддерево <tex>T_{L}</tex>, а ''<tex>\mathrm{split'' }</tex> рекурсивно вызовется от <tex>R</tex>, разделив его на <tex>R'</tex> {{---}} правое поддерево <tex>T_{L}</tex> {{---}}, и <tex>T_{R}</tex>, которые по предположению индукции будут рандомизированными бинарными деревьями поиска. Но <tex>L</tex> также является RBST, т.к. является поддеревом <tex>T</tex>.

Итак для того, чтобы доказать, что <tex>T_{L}</tex> {{---}} рандомизированное бинарное дерево поиска, осталось показать, что любая его вершина <tex>z</tex> с вероятностью <tex dpi = "150">\fracdfrac{1}{m}</tex> окажется в корне, где <tex>m</tex> {{---}} размер <tex>T_{L}</tex>. Действительно:

<tex dpi = "150">P[A | \mid y < x] = \fracdfrac{P[A \; \wedge \; y < x]}{P[y < x]} = \fracdfrac{1 / n}{m / n} = \fracdfrac{1}{m}</tex>

|statement= Если <tex>T</tex> {{---}} рандомизированное бинарное дерево поиска, содержащее множество ключей <tex>K</tex>, <tex>x \notin K</tex>, тогда процедура ''<tex>\mathrm{insert(x, T)'' }</tex> вернёт рандомизированное бинарное дерево поиска <tex>T</tex>, содержащее множество ключей <tex>K \cap cup x</tex>.

Применим индукцию по <tex>n</tex> {{---}} размеру дерева. Если <tex>n = 0</tex>, то теорема верна: после операции ''<tex>\mathrm{insert(x, T)'' }</tex> получим дерево с корнем <tex>x</tex> и двумя пустыми поддеревьями.

В первом случае правое и левое поддеревья <tex>x</tex> по лемме являются рандомизированными BST, а также вероятность того, что <tex>x</tex> окажется в корне, равна <tex dpi = "150">\fracdfrac{1}{n + 1}</tex>. Т.е. новое дерево {{---}} рандомизированное BST.

Во втором случае корень у дерева останется прежнем. Заметим, что для каждого <tex>y \in K</tex> вероятность быть корнем была <tex dpi = "150">\fracdfrac{1}{n}</tex>, а корнем он останется с вероятностью <tex dpi = "150">\fracdfrac{n}{n + 1}</tex>, тогда для каждого <tex>y \in K</tex> вероятность быть корнем равна <tex dpi = "150">\fracdfrac{1}{n} \cdot \fracdfrac{n}{n + 1} = \fracdfrac{1}{n + 1}</tex>. По предположению же индукции поддерево, в которое вставляется <tex>x</tex> становится рандомизированным бинарным деревом поиска; а т.к. другое поддерево корня было рандомизированным, то новое дерево {{---}} рандомизированное BST.

Алгоритм удаления использует операцию ''<tex>\mathrm{merge'' }</tex> {{---}} слияние двух деревьев, удовлетворяющих условию: все ключи в одном из деревьев меньше ключей во втором. Для того, чтобы удалить некоторый ключ <tex>x</tex> из RBST сначала найдём вершину с этим ключом в дереве, используя стандартный алгоритм поиска. Если вершина не найдена, то выходим из алгоритма; в противном случае сливаем правое и левое поддеревья <tex>x</tex> (заметим, что ключи в левом поддереве меньше ключей в правом), удаляем <tex>x</tex>, а корень образовавшегося дерева делаем новым сыном родителя <tex>x</tex>. Псевдокод процедур удаления и слияния приведён ниже.

// удаляет ключ x из дерева T '''Node''' remove(t : '''removeNode'''(T, x: '''T'''): <font color="green">// удаляет ключ x из дерева T</font> '''if''' (t.size(T) == 0) // выйти, вернув пустое дерево T t = RBST()''null'' '''return''' Tt <font color="green">// вернуть пустое поддерево</font> '''if''' (x < Tt.key) Tt.left = remove(Tt.left, x) '''else''' '''if''' (x > Tt.key) Tt.right = remove(Tt.right, x)

// создать пустое дерево Q Q = RBST() Q q = merge(Tt.left, Tt.right) T t = Qq '''return''' Tt

// сливает деревья L и R // результат'''Node''' merge(l : дерево T '''Node''' , r : '''mergeNode'''(L): <font color="green">// сливает деревья l и r, R)результат {{---}} дерево t</font> '''int''' m = Ll.size '''int''' n = Rr.size

'''if''' (total == 0) // вернуть пустое T T t = RBST()''null'' '''return''' Tt <font color="green">// вернуть пустое поддерево</font> '''int''' r = '''random'''(1..<tex>\dots</tex> total) '''if''' (r < m) l.right = merge(l.right, r) <font color="green">// с вероятностью m / (m + n)</font> Ll.right size = merge(L1 + l.left.size + l.right, R).size '''return''' Ll '''if''' r.left = merge(l, r .left) < m) font color="green">// с вероятностью m n / (m + n)</font> R r.left size = merge(L, R1 + r.left).size + r.right.size '''return''' Rr

|statement= Пусть <tex>L</tex> и <tex>R</tex> {{---}} рандомизированные бинарные деревья поиска, содержащие соответственно множества ключей <tex>K_{L}</tex> и <tex>K_{R}</tex>, причём <tex>K_{L} < K_{R}</tex> (то есть каждый элемент <tex>K_{L}</tex> меньше каждого элемента <tex>K_{R}</tex>). Тогда операция ''<tex>\mathrm{merge(L, R)'' }</tex> вернёт рандомизированное бинарное дерево поиска, содержащее множество ключей <tex>K</tex> = <tex>K_{L} \cap cup K_{R}</tex>.

Пусть <tex>m > 0</tex> и <tex>n > 0</tex>, пусть также <tex>L.key = a</tex> или <tex>L.key = b</tex>. Без потери общности делаем корнем <tex>a</tex>. После рекурсивного слияния правого поддерева <tex>L</tex> и <tex>R</tex> получим рандомизированное бинарное дерево поиска (которое является правым поддеревом нового дерева). Левое же поддерево нового дерева тоже рандомизированное. Также верно, что для любого <tex>x \in K_{L}</tex> вероятность быть корнем равна <tex dpi = "150">\fracdfrac{1}{m + n}</tex>: действительно, вероятность оказаться в корне в <tex>L</tex> до слияния равна <tex dpi = "150">\fracdfrac{1}{m}</tex>, вероятность того, что элемент останется корнем после слияния равна <tex dpi = "150">\fracdfrac{m}{m + n}</tex>; осталось применить правило умножения.

|statement= Если <tex>T</tex> {{---}} рандомизированное бинарное дерево поиска, содержащее множество ключей <tex>K</tex>, тогда процедура ''<tex>\mathrm{remove(x, T)'' }</tex> вернёт рандомизированное бинарное дерево поиска <tex>T'</tex>, содержащее множество ключей <tex>K \backslash \{x\}</tex>

Возможно два случая: если <tex>x</tex> {{---}} корень <tex>T</tex>, то по лемме, после удаления получим рандомизированное бинарное дерево поиска; если же <tex>x</tex> {{---}} не корень <tex>T</tex>, то <tex>x</tex> рекурсивно удаляется из поддерева исходного дерева, и по предположению индукции после удаления получаем рандомизированное BST. Осталось лишь показать, что для любого <tex>y \in T, y \neq x</tex> вероятность оказаться корнем после удаления равна <tex dpi = "150">\fracdfrac{1}{n - 1}</tex>.

событие <tex>B</tex> {{---}} <tex>x</tex> был корнем <tex>T</tex> (до операции ''<tex>\mathrm{remove''}</tex>);

событие <tex>C</tex> {{---}} <tex>y</tex> стал корнем <tex>T'</tex> после операции ''<tex>\mathrm{merge'' }</tex> (но до этого им не являлся);

событие <tex>D</tex> {{---}} <tex>y</tex> был корнем <tex>T</tex> (до операции ''<tex>\mathrm{remove''}</tex>);

<tex>P[A] = P[A|\mid B] \times cdot P[B] + P[A|\mid\neg B] \times cdot P[\neg B] = P[C] \times cdot 1/n + P[D|\mid\neg B] \times cdot (n - 1)/n = </tex><tex dpi = "150">\fracdfrac{1}{n - 1} \times cdot \fracdfrac{1}{n} + \fracdfrac{1}{n - 1} \times cdot \fracdfrac{n - 1}{n} = \fracdfrac{1}{n - 1}</tex>.

Очевидно, что время работы приведённых алгоритмов пропорционально глубине дерева. Но т.к. математическое ожидание глубины рандомизированного бинарного дерева поиска есть <tex>O (\log n)</tex>, где <tex>n</tex> {{---}} число вершин в дереве<ref>Thomas H. Cormen, Charles E. Leiserson, Ronald L. Rivest, Clifford Stei, "Introduction to Algorithms", Second Edition {{---}} Chapter 12.4</ref>, то математическое ожидание времени работы поиска, вставки и удаления {{---}} также <tex>O (\log n)</tex>.

==См. также==*[[Дерево поиска, наивная реализация2-3 дерево]]* [[B-дерево]]* [[Splay-дерево]]*[[Декартово АВЛ-дерево]]*[[Декартово Красно-черное дерево по неявному ключу]]

==СсылкиПримечания==*[http:<references //en.wikipedia.org/wiki/Treap Википедия {{---}} Дерамида и RBST]*[http://en.wikipedia.org/wiki/Random_binary_tree Википедия {{---}} случайное бинарное дерево]>

== Литература Источники информации ==* [http://en.wikipedia.org/wiki/Random_binary_tree Wikipedia {{---}} Random binary tree]* [http://en.wikipedia.org/wiki/Treap Wikipedia {{---}} Treap]