Обсуждение участницы:Анна

Материал из Викиконспекты
Версия от 21:59, 19 ноября 2015; Анна (обсуждение | вклад) (Доказательство корректности гамма-алгоритма)
Перейти к: навигация, поиск

Алгоритм разделения АВЛ-дерева на два, где в первом дереве все ключи меньше заданного x, а во втором - больше

Пусть у нас есть дерево [math]T[/math]. Мы должны разбить его на два дерева [math]T_{1}[/math] и [math]T_{2}[/math] такие, что [math]T_{1} \leqslant x[/math] и [math]x \lt T_{2}[/math].

Предположим, что корень нашего дерева [math]\leqslant x[/math], в таком случае все левое поддерево вместе с корнем после разделения отойдет в дерево [math]T_{1}[/math]. Тогда рекурсивно спускаемся в правое поддерево и там проверяем это условие (так как часть правого поддерева тоже может содержать ключи [math]\leqslant x[/math]). Если же корень оказался [math]\gt x[/math], то мы спускаемся той же рекурсией, но только в левое поддерево и ищем там.

Пусть мы пришли в поддерево [math]S[/math], корень которого [math]\leqslant x[/math]. В таком случае этот корень со своим левым поддеревом должен отойти в дерево [math]T_{1}[/math]. Поэтому мы делаем следующее: запоминаем ссылку на правое поддерево [math]S[/math], удаляем корень, запоминая его значение (не меняя конфигурацию дерева, то есть просто делаем ссылки на него NULL'ами). Таким образом, мы отделяем сбалансированное АВЛ-дерево (бывшее левое поддерево [math]S[/math]). Делаем новую вершину со значением бывшего корня правым листом самой правой вершины [math]S[/math] и запускаем балансировку. Обозначим полученное дерево за [math]T'[/math]. Теперь нам нужно объединить его с уже построенным ранее [math]T_{1}[/math] (оно может быть пустым, если мы первый раз нашли такое дерево [math]S[/math]). Для этого мы ищем в дереве [math]T_{1}[/math] самое правое поддерево [math]P[/math] высоты, равной высоте [math]T'[/math] (спускаясь от корня всегда в правые поддеревья). Делаем новое дерево [math]K[/math], сливая [math]P[/math] и [math]T'[/math] (очевидно, все ключи в [math]T_{1}[/math] меньше ключей в [math]T'[/math], поэтому мы можем это сделать). Теперь в дереве [math]T_{1}[/math] у отца вершины, в которой мы остановились при поиске дерева [math]P[/math], правым поддеревом делаем дерево [math]K[/math] и запускаем балансировку. После нужно спуститься в правое поддерево бывшего дерева [math]S[/math] (по ссылке, которую мы ранее запомнили) и обработать его.

Если мы пришли в поддерево [math]Q[/math], корень которого [math]\gt x[/math], совершаем аналогичные действия: делаем NULL'ами ссылки на корень [math]Q[/math], запоминая ссылку на его левое поддерево. Делаем новую вершину со значением бывшего корня левым листом самой левой вершины [math]Q[/math] и запускаем балансировку. Объединяем полученное АВЛ-дерево с уже построенным ранее [math]T_{2}[/math] аналогичным первому случаю способом, только теперь мы ищем самое левое поддерево [math]T_{2}[/math].

Рассмотри пример (рис. 1). Цветом выделены поддеревья, которые после разделения должны отойти в дерево [math]T_{1}[/math]. [math]x = 76[/math].

Рис. 1. Разделение АВЛ-дерева на два.

Корень дерева [math]\leqslant x[/math], поэтому он со всем выделенным поддеревом должен отойти в дерево [math]T_{1}[/math]. По описанному выше алгоритму отделяем это поддерево с корнем и делаем из них сбалансированное АВЛ-дерево [math]T'[/math] (рис. 2). Так как это первая ситуация, в которой корень рассматриваемого поддерева был [math]\leqslant x[/math], [math]T'[/math] становится [math]T_{1}[/math]. Далее по сохраненной ссылке спускаемся в правое поддерево. Его корень [math]\gt x[/math]. Следовательно, строим из него и его правого поддерева [math]T_{2}[/math] и спускаемся в левое поддерево. Снова корень [math]\leqslant x[/math]. Строим новое [math]T'[/math] и объединяем его с уже существующим [math]T_{1}[/math] (рис. 3).

Рис. 2. Создание T'.
Рис. 3. Объединение T' и T1.

Далее действуем по алгоритму и в итоге получаем (рис. 4):

Рис. 4. АВЛ-деревья после разделения.

Данный алгоритм имеет сложность [math]O(\log^{2} n)[/math]. Рассмотрим решение, которое имеет сложность [math]O(\log{n})[/math].

Вернемся к примеру (рис. 1). Теперь рекурсивно спустимся вниз и оттуда будем строить деревья [math]T_{1}[/math] и [math]T_{2}[/math], передавая наверх корректные АВЛ-деревья. То есть для рис. 1 первым в дерево [math]T_{1}[/math] придет вершина [math]75[/math] с левым поддеревом (выделено светло-зеленым цветом), так как это корректное АВЛ-дерево, оно же и вернется из рекурсии. Далее мы попадем в вершину со значением [math]70[/math] и должны слить ее и ее левое поддерево (выделено светло-синим) с тем, что нам пришло. И сделать это нужно так, чтобы передать наверх корректное АВЛ-дерево. Будем действовать по такому алгоритму, пока не дойдем до вершины.

Пусть мы пришли в поддерево [math]S[/math] с корнем [math]\leqslant x[/math]. Тогда сольем его с уже построенным на тот момент [math]T_{1}[/math] ([math]T_{1}[/math] пришло снизу, а значит по условию рекурсии это корректное АВЛ-дерево, [math]S \leqslant T_{1}[/math] и [math]h(T_{1}) \leqslant h(S)[/math]). Но так как обычная процедура слияния сливает два АВЛ-дерева, а [math]S[/math] не является корректным АВЛ-деревом, мы немного ее изменим. Пусть мы в дереве [math]S[/math] нашли самое правое поддерево [math]K[/math], высота которого равна высоте [math]T_{1}[/math]. Тогда сделаем новое дерево [math]T'[/math], корнем которого будет вершина [math]S[/math] (без нее это дерево является сбалансированным), правым поддеревом — [math]T_{1}[/math], левым — [math]K[/math]. И подвесим [math]T'[/math] на то место, где мы остановились при поиске [math]K[/math]. Запустим балансировку. В случае, когда корень поддерева, в которое мы пришли, [math]\gt x[/math], все аналогично.

Разберем пример на рис. 1. Пусть мы рекурсивно спустились до узла [math]77[/math]. Ключ больше [math]x[/math], поэтому эта вершина станет деревом [math]T_{2}[/math] и передастся наверх. Теперь мы поднялись в узел [math]75[/math]. Он со своим левым поддеревом станет деревом [math]T_{1}[/math] и мы снова поднимемся наверх в узел [math]70[/math]. Он со своим левым поддеревом снова должен отойти в дерево [math]T_{1}[/math], и так как теперь дерево [math]T_{1}[/math] уже не пустое, то их надо слить. После слияния по описанному выше алгоритму получим (рис. 5)

Рис. 5.

После мы поднимемся в вершину с ключом [math]80[/math]. Она с правым поддеревом отойдет в дерево [math]T_{2}[/math] (рис. 6).

Рис. 6.

И на последней итерации мы поднимемся в корень дерева с ключом [math]50[/math], он с левым поддеревом отойдет в дерево [math]T_{1}[/math], после чего алгоритм завершится.

Пусть поддеревьев с ключами [math]\leqslant x[/math] оказалось больше, чем поддеревьев с ключами [math]\gt x[/math]. Докажем для них логарифмическую асимптотику. Дерево на последнем уровне имеет высоту [math]H_{k}[/math] (она может быть не равна [math]1[/math], если мы придём в [math]x[/math]). Его мы передаем наверх и вставляем в поддерево высотой [math]H_{k-1}[/math]. [math]H_{k} \leqslant H_{k-1}[/math], так как разница высот поддеревьев у любой вершины не больше [math]1[/math], и мы при переходе от [math]H_{k}[/math] к [math]H_{k-1}[/math] поднимаемся как минимум на одну вершину вверх. Слияние этих поддеревьев мы выполним за [math]H_{k-1} - H_{k}[/math], получим в итоге дерево высоты не большей, чем [math]H_{k-1}[/math]. Его мы передадим наверх, поэтому в следующий раз слияние будет выполнено за [math]H_{k-2} - H_{k - 1}[/math] и так далее. Таким образом мы получим [math](H - H_{1}) + (H_{1} - H_{2}) + (H_{2} - H_{3}) + \cdots + (H_{k - 1} - H_{k}) = H - H_{k} = O(\log{n})[/math].

Итоговая асимптотика алгоритма — [math]O(\log{n})[/math].

Гамма-алгоритм

Задача:
Определить, является ли граф планарным, и, если да, произвести его плоскую укладку.

Существует теорема Понтрягина-Куратовского, которая говорит, что граф планарен тогда и только тогда, когда он не содержит подграфов, гомеоморфных [math] K_{5} [/math] или [math] K_{3, 3} [/math]. Но этот критерий очень трудно проверить на практике, поэтому данная теорема представляет лишь теоретический интерес.

Чтобы проверить планарность графа и произвести его плоскую укладку, удобно пользоваться гамма-алгоритмом.

Входные данные

На вход алгоритму подаются графы со следующими свойствами:

  1. Граф связный,
  2. Граф содержит хотя бы один цикл,
  3. Граф не имеет мостов.

Если нарушено свойство [math]1[/math], то граф нужно укладывать отдельно по компонентам связности. Если нарушено свойство [math]2[/math], то граф — дерево и нарисовать его плоскую укладку тривиально.

Более подробно рассмотрим случай, когда в графе [math]G[/math] нарушено свойство [math]3[/math]. Сначала все мосты нужно убрать, далее произвести отдельную укладку всех компонент следующим образом: уложим одну компоненту связности, а следующую компоненту, связанную с первой в графе [math]G[/math] мостом, будем рисовать в той грани, в которой лежит вершина, принадлежащая мосту. Иначе может сложиться ситуация, когда концевая вершина моста будет находиться внутри плоского графа, а следующая компонента - снаружи. Таким образом мы сможем соединить мостом нужные вершины. Далее будем так поступать с каждой новой компонентой.

Описание алгоритма

Рассмотрим работу алгоритма, параллельно разбирая на примере каждый шаг. Пусть дан граф [math]G[/math] (рис. 1).

Рис. 1. Исходный граф.

[math]1.[/math] Первый этап - инициализация алгоритма.

В графе [math]G[/math] выбирается любой простой цикл и производится его укладка на плоскость. Пусть в примере это будет цикл [math]\{1, 2, 3, 4, 5, 6\}[/math]. После его укладки получаем две грани: [math]\Gamma_{1}[/math] и [math]\Gamma_{2}[/math] (рис. 2).

Рис. 2. Укладка цикла на плоскость.

Уже уложенную во время работы алгоритма часть будем обозначать [math]G_{plane}[/math]. В примере сейчас [math]G_{plane}[/math] — выбранный цикл [math]\{1, 2, 3, 4, 5, 6\}[/math].

[math]2.[/math] Второй этап - общий шаг.

Построим множество сегментов. Каждый сегмент [math]S[/math] относительно уже построенного [math]G_{plane}[/math] представляет собой одно из двух:

  1. ребро, оба конца которого принадлежат [math]G_{plane}[/math], но само оно не принадлежит [math]G_{plane}[/math],
  2. связную компоненту графа [math]G \backslash G_{plane}[/math], дополненную всеми ребрами графа [math]G[/math] такими, у которых один из концов принадлежит связной компоненте, а второй принадлежит графу [math]G_{plane}[/math].

Вершины, одновременно принадлежащие [math]G_{plane}[/math] и какому-либо сегменту, назовем контактными. На рис. 3 изображены сегменты для нашего примера. Контактные вершины обведены в квадрат.

Рис. 3. Выделение сегментов.

Пусть в каком-то сегменте нет ни одной контактной вершины. В таком случае граф до выделения [math]G_{plane}[/math] был несвязным, что противоречит условию. Пусть контактная вершина в сегменте только одна. Это значит, что в графе был мост, чего быть не может так же по условию. Таким образом, в каждом сегменте имеется не менее двух контактных вершин. Соответственно, в каждом сегменте есть цепь между любой парой контактных вершин.

Пусть грань [math]\Gamma[/math] вмещает сегмент [math]S[/math], если номера всех контактных вершин [math]S[/math] принадлежат этой грани, [math]S \subset \Gamma[/math]. Очевидно, таких граней может быть несколько. Множество таких граней обозначим [math]\Gamma(S)[/math], а их число — [math]|\Gamma(S)|[/math].

Итак, рассмотрим все сегменты [math]S_{i}[/math] и для каждого определим число [math]|\Gamma(S_{i})|[/math]. Если найдется такой номер [math]i[/math], что [math]|\Gamma(S_{i})| = 0[/math], то граф не планарен, алгоритм завершает работу. Иначе выбираем такой сегмент [math]S_{i}[/math], для которого число [math]|\Gamma(S_{i})|[/math] минимально. Если таких сегментов несколько, можно выбрать любой из них. Найдем в этом сегменте цепь между двумя контактными вершинами и уложим ее в любую грань из множества [math]|\Gamma(S_{i})|[/math], совместив контактные вершины сегмента с соответствующими вершинами грани. Выбранная грань разобьется на две. Выбранный сегмент после того, как из него взяли цепь, либо исчезнет, либо распадется на меньшие части, в которых будут новые контактные вершины, ведущие к вершинам обновленного [math]G_{plane}[/math].

В примере для любого [math]i[/math]: [math]S_{i} \subset \{\Gamma_{1}, \Gamma_{2}\}[/math], то есть [math]|\Gamma(S_{i})| = 2[/math]. Следовательно, можем выбрать любой [math]S_{i}[/math]. Пусть это будет сегмент [math]S_{1}[/math]. В нем есть цепь [math]\{1, 4\}[/math]. Вставим эту цепь в грань [math]\Gamma_{1}[/math], например, и этот сегмент исчезнет. После укладки цепи граф [math]G[/math] и сегменты будут выглядеть так (рис. 4):

Рис. 4. Граф и сегменты после второго этапа.


[math]3.[/math] Третий и последующие этапы аналогичны второму. Повторять вышеуказанные действия нужно либо до тех пор, пока не будет получена плоская укладка, то есть множество сегментов не останется пустым, либо пока не будет получено, что граф не планарен.

Разберем пример до конца. Повторим снова общий шаг. Теперь [math]|\Gamma(S_{1})| = |\Gamma(S_{3})| = 1[/math], [math]|\Gamma(S_{2})| = |\Gamma(S_{4})| = 2[/math]. Возьмем [math]S_{1}[/math]. В ней цепь [math]\{2, 5\}[/math], которую мы уложим в грань [math]\Gamma_{2}[/math], после чего этот сегмент исчезнет. Теперь картина будет следующая (рис. 5):

Рис. 5. Граф и сегменты после третьего этапа.

На следующем шаге [math]|\Gamma(S_{1})| = |\Gamma(S_{3})| = 2[/math], [math]|\Gamma(S_{2})| = 1[/math]. Выбираем сегмент [math]S_{2}[/math], содержащий цепь [math]\{3, 5\}[/math]. Уложим ее в грань [math]\Gamma_{2}[/math], после чего этот сегмент снова исчезнет (рис. 6).

Рис. 6. Граф и сегменты после четвертого этапа.

Теперь [math]|\Gamma(S_{1})| = 1[/math], а [math]|\Gamma(S_{3})| = 2[/math]. Уложим сначала цепь [math]\{2, 4\}[/math] из первого сегмента, он пропадет, потом уложим цепь [math]\{6, 7, 5\}[/math] из третьего. В результате граф будет полностью уложен на плоскость, множество сегментов останется пустым (рис. 7).

Рис. 7. Плоская укладка графа.

Таким образом, мы получили плоскую укладку исходного графа [math]G[/math].

Опишем алгоритм коротко и формально:

  1. Инициализация. Выбирается простой цикл в исходном графе и изображается на плоскости.
  2. Общий шаг. Этот шаг повторяется до тех пор, пока граф не будет уложен или пока не будет получено, что граф не планарен.
    • строится множество сегментов
    • для каждого сегмента вычисляется величина [math]|\Gamma(S)|[/math]. Если существует [math]i[/math]: [math]|\Gamma(S_{i})| = 0[/math], то граф не планарен, алгоритм завершает работу
    • выбирается сегмент с минимальным числом [math]|\Gamma(S_{i})|[/math]
    • в этом сегменте выбирается цепь между двумя контактными вершинами
    • эта цепь укладывается в любую грань, вмещающую данный сегмент
  3. Либо получена плоская укладка графа, либо граф оказался не планарен.

Доказательство корректности гамма-алгоритма

Перед доказательством корректности приведем ряд важных вспомогательных лемм.

Пусть два сегмента [math]S_{1}[/math] и [math]S_{2}[/math] называются конфликтующими относительно уже уложенной части, если:

  1. грань вмещает каждый из сегментов [math]S_{1}[/math] и [math]S_{2}[/math]
  2. в этих сегментах есть две цепи между контактными вершинами [math]L_{1}[/math] и [math]L_{2}[/math] соответственно такие, что их невозможно уложить в одну грань без пересечения

Например, на рис. 1 конфликтующими являются сегменты [math]S_{1}[/math] и [math]S_{2}[/math].

Рис. 1. Конфликтующие сегменты.
Лемма (1):
Конфликтующие сегменты [math]S_{1}[/math] и [math]S_{2}[/math] обладают следующим свойством: если [math]|\Gamma(S_{1})| \geq 2[/math] и [math]|\Gamma(S_{2})| \geq 2[/math], то [math]\Gamma(S_{1}) = \Gamma(S_{2}) = 2[/math].
Доказательство:
[math]\triangleright[/math]

Сначала докажем, что [math]\Gamma(S_{1}) = \Gamma(S_{2})[/math]. Предположим противное. Тогда по условию леммы найдутся три различные грани [math]\Gamma_{1}, \Gamma_{2}[/math] и [math]\Gamma_{3}[/math] такие, что [math]\Gamma_{1} \in \Gamma(S_{1}), \Gamma_{2} \in \Gamma(S_{2}), \Gamma_{3} \in Q = \Gamma(S_{1}) \cap \Gamma(S_{2}) \neq \emptyset[/math]. Тогда любые цепи [math]L_{1} \subset S_{1}[/math] и [math]L_{2} \subset S_{2}[/math] укладываются в [math]\Gamma_{1}[/math] и [math]\Gamma_{2}[/math] соответственно. Но это значит, что любая пара цепей [math]L_{1}[/math] и [math]L_{2}[/math] одновременно укладывается вне грани [math]\Gamma_{3}[/math]. Следовательно, они одновременно укладываются и внутри грани [math]\Gamma_{3}[/math], причем без пересечений. Но это противоречит тому, что [math]S_{1}[/math] и [math]S_{2}[/math] — конфликтующие сегменты. Таким образом, [math]\Gamma(S_{1}) = \Gamma(S_{2})[/math].

Теперь покажем, что [math]|Q| = 2[/math]. Доказательство снова поведем методом от противного. Пусть [math]|Q| \geq 3[/math]. Тогда снова существует три различные грани [math]\Gamma_{1}, \Gamma_{2}[/math] и [math]\Gamma_{3} \in Q[/math]. Аналогичными рассуждениями снова приходим к противоречию с тем, что [math]S_{1}[/math] и [math]S_{2}[/math] — конфликтующие сегменты.
[math]\triangleleft[/math]

Замечание. Пусть имеется множество сегментов таких, что они удовлетворяют следующей схеме: есть первый сегмент [math]S_{1}[/math], второй [math]S_{2}[/math], который конфликтует с [math]S_{1}[/math], третий [math]S_{3}[/math], который конфликтует с [math]S_{2}[/math], но не с [math]S_{1}[/math], и так далее, причем каждый из них укладывается в две грани. Тогда из леммы 1 следует, что эти грани являются общими для всех сегментов множества и можно укладывать их следующим образом: цепь [math]L_{1} \subset S_{1}[/math] укладывается в первую грань [math]\Gamma_{1}[/math], цепь [math]L_{2} \subset S_{2}[/math] — во вторую [math]\Gamma_{2}[/math], [math]L_{3} \subset S_{3}[/math] — снова в [math]\Gamma_{1}[/math] и так для всех элементов множества. Если цепочка сегментов замыкается в цикл четной длины, то проблем никаких нет. Если длина цикла нечетная, то два последних конфликтующих сегмента нужно будет уложить в одну грань без пересечений, что по определению невозможно. Следовательно, получить плоскую укладку нельзя.


Определение:
Частичной укладкой [math]G'[/math] планарного графа [math]G[/math] называется граф, который можно получить из какой-либо укладки графа [math]G[/math] на плоскости путем удаления некоторых ребер и вершин.


Таким образом во всякую частичную укладку планарного графа [math]G[/math] можно уложить оставшуюся часть, а именно недостающие вершины и ребра графа [math]G[/math].

Построим некоторый служебный граф [math]A(G')[/math] по следующей схеме: каждому сегменту в [math]G'[/math] сопоставим некоторую вершину в [math]A(G')[/math], причем две вершины будут смежны тогда и только тогда, когда соответствующие им сегменты являются конфликтующими.

Например, для графа на рис. 1 служебный граф будет иметь вид (рис. 2):

Рис. 2. Служебный граф. Вершины обозначены как соответствующие сегменты.
Лемма (2):
Если результатом некоторого шага работы гамма-алгоритма является частичная укладка [math]G'[/math] планарного графа [math]G[/math] такая, что [math]|\Gamma(S)| \geq 2[/math] для любого сегмента [math]S[/math] относительно [math]G'[/math], то [math]A(G')[/math]двудольный граф.
Доказательство:
[math]\triangleright[/math]
Докажем от противного. Пусть [math]A(G')[/math] — не двудольный. Тогда по теореме Кенинга в нем есть цикл нечетной длины. Этому циклу соответствует некоторая последовательность сегментов [math]S_{1}, S_{2}, \cdots S_{2m+1}, S_{1}[/math] относительно [math]G'[/math], в которой каждые соседние сегменты конфликтующие по определению. По лемме 1 [math]\Gamma(S_{i}) = \{\Gamma_{1}, \Gamma_{2}\}, i \in \{1 \cdots 2m+1\}[/math]. Так как [math]G'[/math] — частичная укладка графа, то все сегменты [math]S_{1}, S_{2}, \cdots S_{2m+1}[/math] могут быть уложены. А так как соседние сегменты этой последовательности конфликтующие, то они должны быть уложены в разные грани, что невозможно, так как число сегментов в последовательности нечетное. Получили противоречие. Следовательно, [math]A(G')[/math] — двудольный.
[math]\triangleleft[/math]
Теорема (О корректности гамма-алгоритма):
Гамма-алгоритм корректен, то есть если [math]G[/math] — планарный граф, то результатом каждого шага гамма-алгоритма является частичная укладка [math]G'[/math].
Доказательство:
[math]\triangleright[/math]

Докажем индукцией по числу шагов.

База индукции: полученный на этапе инициализации граф [math]G'_{0}[/math] является простым циклом, он будет присутствовать в любой укладке графа [math]G[/math]. Таким образом, [math]G'_{0}[/math] является частичной укладкой.

Шаг индукции: пусть граф [math]G'_{k-1}[/math], полученный на [math]k-1[/math]-м шаге работы алгоритма, является частичной укладкой. Докажем, что граф [math]G'_{k} = G'_{k-1} \cup L_{k}[/math], полученный на [math]k[/math]-м шаге присоединением цепи [math]L_{k}[/math], также является частичной укладкой.

Заметим, что на текущем шаге нет такого сегмента [math]S[/math] относительно [math]G'_{k-1}[/math], для которого бы выполнялось равенство [math]\Gamma(S) = \emptyset[/math], так как в противном случае существовала бы цепь этого сегмента, контактные вершины которой принадлежали бы разным граням и укладка которой была бы невозможна. Следовательно, нельзя было бы уложить [math]S[/math],что противоречит тому, что [math]G[/math] — планарный граф. Значит, мы можем рассматривать только следующие два случая:

[math]1.[/math] Существует сегмент [math]S[/math] для которого есть единственная вмещающая грань [math]\Gamma[/math], то есть [math]|\Gamma(S)| = 1[/math]. Так как только грани [math]\Gamma[/math] принадлежат все контактные вершины [math]S[/math], то укладка этого сегмента в эту грань неизбежна. Это значит, что помещая любую цепь [math]L \subset S[/math], снова получим частичную укладку графа.

[math]2.[/math] Для любого сегмента [math]S[/math] [math]|\Gamma(S)| \geq 2[/math].
[math]\triangleleft[/math]