Остовные деревья Лекция 4

Задача «Минимальное остовное дерево» Дано: Граф G, веса c: E(G) R. Найти остовное дерево в G наименьшего веса или определить, что G несвязный.

Задача «Максимальный взвешенный лес» Дано: Граф G, веса c: E(G) R. Найти лес в G наибольшего веса.

Эквивалентные задачи Будем говорить, что задача P линейно сводится к задаче Q, если существуют функции f и g, вычислимые за линейное время, такие что f преобразует частную задачу x из P в частную задачу y из Q, и g преобразует решение f (x) в решение x. Eсли P линейно сводится к Q и Q линейно сводится к P, то обе задачи называются эквивалентными.

Эквивалентность Предложение 4.1 Задача «Минимальное остовное дерево» и задача «Максимальный взвешенный лес» эквивалентны.

Доказательство (G,c) исходный пример задачи «Максимальный взвешенный лес». Удалим все ребра отрицательного веса. Положим c'(e) = – c(e). Добавим минимальное множество ребер F, так чтобы полученный граф G' стал связным. (Веса можно взять любые.) Решим задачу «Минимальное остовное дерево» на примере (G',c'). Удалив из решения множество ребер F, получим решение исходной задачи. (G,c) исходный пример задачи «Минимальное остовное дерево». Положим c'(e) = K – c(e), где K = 1 + max e E(G) c(e). Решение задачи «Максимальный взвешенный лес» на примере (G',c') дает решение задачи «Минимальное остовное дерево» на исходном примере.

Условия оптимальности Теорема 4.2 Пусть (G,c) пример задачи «Минимальное остовное дерево», и пусть T остовное дерево в G. Тогда следующие условия эквивалентны: a) T оптимум. b)Для любого e = {x, y} E(G)\ E(T), все ребра из x-y-пути в T не дороже чем e. c)Для любого e E(T), e ребро наименьшей стоимости из (V(C)), где C связная компонента на T– e.

(c) (a) (с) Пусть T такое, что для любого e E(T), e ребро наименьшей стоимости из (V(C)), где C связная компонента на T– e. Пусть T* оптимальное решение, такое что E(T) E(T*) максимально возможное. Покажем, что T = T*. Пусть e = {x,y} E(T)\ E(T*). Пусть C связная компонента на T– e. T* + e содержит цикл D. Так как e E(D) δ(C), то существует f e, f E(D) δ(C). Тогда (T* + e) – f является остовным деревом. T* оптимум c(e) c(f) и (с) c(f) c(e). c(f) = c(e) и (T* + e) – f является оптимальным остовным деревом. Противоречие, так как в (T* + e) – f больше на одно общее ребро с T, чем в T*.

Алгоритм Краскала (1956) Input: Связный граф G, веса c: E(G) R. Output: Остовное дерево T наименьшего веса. 1) Сортируем ребра так, что c(e 1 ) c(e 2 ) …c(e m ). 2) Set T (V(G), ). 3) For i 1 to m do: If T+e i не содержит цикла then T T +e i.

Алгоритм Краскала (2) Теорема 4.3 Алгоритм Краскала находит оптимальное решение за O(mn).

Алгоритм Краскала (3) Теорема 4.4 Алгоритм Краскала можно реализовать за O(m log n).

Алгоритм Краскала (1956) Input: Связный граф G, веса c: E(G) R. Output: Остовное дерево T наименьшего веса. 1) Сортируем ребра так, что c(e 1 ) c(e 2 ) …c(e m ). 2) Set T (V(G), ). 3) For i 1 to m do: If T+e i не содержит цикла then T T +e i.

Как улучшить шаг 3 Основная цель шага 3 проверить не приведет ли добавление ребра e i = {v,w} к образованию цикла. Это эквивалентно проверки лежат ли v и w в одной связной компоненте. По ходу алгоритма будем строить дополнительный орлес B с V(B) = V(G), такой что связные компоненты B индуцированы на тех же вершинах, что и связные компоненты T.

Проверка Первоначально, B = (V(G), ø ) и h(v) = 0, для всех v V(G), где h(v) длина максимального пути из v в B. 3.1 Для ребра e i = {v,w} находим корни r v и r w ордеревьев в B, содержащих v и w. 3.2 Если r v = r w, то переходим к следующему ребру и идем на Если r v r w, то добавляем e i к T Если h(r v ) h(r w ), то добавляем дугу (r v, r w ) к B, иначе добавляем дугу (r w, r v ) к B.

Время работы шага 3 Время работы пропорционально длине r v -v-пути в B. Покажем, что любое ордерево в B с корнем r имеет по крайней мере 2 h(r) вершин. Когда B = (V(G), ø ) утверждение очевидно. Пусть алгоритм добавляет дугу (x,y) в B. Получаем новое ордерево с корнем в x и числом вершин 2 h(x) + 2 h(y). Если h(x) > h(y), то значение h(x) не меняется и утверждение справедливо. Если h(x) = h(y), то значение h(x) увеличивается на 1. Число вершин в новом ордереве 2 h(x) + 2 h(y) =2 h(x)+1. h(r) log n, и трудоемкость шага 3 mlog n.

Алгоритм Прима (1957) Input: Связный граф G, веса c: E(G) R. Output: Остовное дерево T минимального веса. 1) Выбрать v V(G). T ({v}, ). 2) While V(T) V(G) do: Выбрать ребро e G (V(T)) минимальной стоимости. T T +e.

Алгоритм Прима (2) Теорема 4.5 Алгоритм Прима находит решение за O(n 2 ) элементарных операций.

Как реализовать шаг 2 While V(T) V(G) do: Выбрать ребро e G (V(T)) минимальной стоимости. T T +e. Для каждой вершины v V(T) будем хранить самое дешевое ребро (кандидата) из v в V(T). Тогда выбор ребра минимальной стоимости и замена кандидатов может быть выполнена за O(n) элементарных операций.

Задача «Максимальный взвешенный ориентированный лес» Дано: Орграф G, веса c: E(G) R. Найти ориентированный лес в G наибольшего веса.

Задача «Минимальное остовное ориентированное дерево» Дано: Орграф G, веса c: E(G) R. Найти остовное ориентированное дерево в G наименьшего веса или определить, что оно не существует.

Задача «Минимальное остовное корневое ориентированное дерево» Дано: Орграф G, вершина r V(G), веса c: E(G) R. Найти остовное ориентированное дерево с корнем r в G наименьшего веса или определить, что оно не существует.

Эквивалентность трех задач Предложение 4.6. Задачи «Максимальный взвешенный ориентированный лес», «Минимальное остовное ориентированное дерево» и «Минимальное остовное корневое ориентированное дерево» эквивалентны. Упражнение 4.1 Доказать предложение 4.6.

Ориентированный лес Орграф называется ориентированным лесом, если соответствующий ему граф является лесом и каждая вершина v имеет не более одного входящего ребра.

Ориентированный лес и циклы Предложение 4.7. Пусть B орграф с для всех x V(B). Тогда B имеет цикл тогда и только тогда, когда соответствующий ему граф имеет цикл. Лемма 4.8. (Karp [1972]) Пусть B 0 подграф G максимального веса с для всех v V(B 0 ). Тогда существует оптимальный ориентированный лес B в G такой, что для каждого цикла C в B 0, |E(C)\ E(B)| = 1.

Доказательство леммы a1a1 b1b1 a2a2 b2b2 a3a3 b3b3 С B 0 E(C)\ E(B)={(a 1, b 1 ),…, (a k, b k )} Покажем, что в B есть b i -b i-1 -путь для всех i. Пусть B – оптимальный орлес в G имеющий макси -мально много ребер из B 0.

Покажем, что в B есть b i -b i-1 -путь для всех i. a i-1 b i-1 aiai bibi a i+1 b i+1 С B 0 e E(B) E(B):={(x,y) E(B)}\{e} U {(a i,b i )} P B B не орлес. [b i-1 a i ] B

Основная идея Найти B 0 ориентированный подграф G максимального веса, в котором в каждую вершину входит не больше одной дуги. Стянуть каждый цикл в B 0 в вершину. Перераспредилить веса в новом графе G 1, так чтобы любой оптимальный орлес в G 1 соответствовал оптимальному орлесу в G.

Алгоритм Эдмондса построения ориентированного леса максимального веса (1967) Input: орграф G, веса c: E(G) R +. Output: орлес максимального веса B of G. 1)Set i 0, G 0 G, и c 0 c. 2)Пусть B i подграф G максимального веса с для всех v B i. 3) If B i не содержит циклов then B B i и go to (5). 4)Построим (G i+1, c i+1 ) из (G i, c i ): do для каждого цикла C из B i. Стянем C к одной вершине v C в G i+1. For каждого ребра e = ( z, y) E( G i ) с z V(C), y V(C) do: Set c i+1 (e) c i (e) – c i ( (e,C)) + c i (e C ) и (e) e, где e ( z, v C ), (e,C)=(x,y) E(C), и e C самое дешевое ребро C. 5) If i = 0 then stop. 6) For каждого цикла C из B i-1 do: If есть ребро e ( z, v C ) E(B) then E(B) (E(B)\{e }) (e) ( E(C)\{ ( (e),C)}) else E(B) E(B) (E(C)\{e C }). Set V(B) V(G i-1 ), i i–1 и go to (5).

Шаг 4 z y x z vCvC С B i e α(e,C) For каждого ребра e = ( z, y) E(G i ) с z V(C), y V(C) do: Set c i+1 (e) c i (e) – c i ( (e,C)) + c i (e C ) и (e) e, где e ( z, v C ), (e,C)=(x,y) E(C), и e C самое дешевое ребро C. eCeC e

Шаг 6 z vCvC e E(B) z y x С B i α(e,C) eCeC e vCvC y x С B i α(e,C) eCeC

Алгоритм Эдмондса Теорема 4.9 Алгоритм Эдмондса находит оптимальное решение.

Доказательство Последовательно применяя шаг 4 алгоритма, мы получили последовательность (G i, c i ), i = 0,…, k. Ясно, что полученный на последнем применении шага 4 орлес B является оптимальным для (G k, c k ). Покажем, что на шаге 6 мы последовательно строим оптимальные решения для (G i, c i ), i = k–1,…,0. Мы хотим показать, что шаг 6 переводит ориентированный лес наибольшего веса B для G i в ориентированный лес наибольшего веса B* для G i–1.

Доказательство (2) Пусть B' i–1 произвольный орлес в G i–1, такой что |E(C)\ E(B' i–1 )| = 1 для любого C из B i–1. Пусть B' i получается из B' i–1 стягиванием циклов в B i–1. Тогда B' i орлес в G i.

Шаг 4 z y x z vCvC С B i e α(e,C) Set c i+1 (e) c i (e) – c i ( (e,C)) + c i (e C ) eCeC e

Индукция По индукции, B ориентированный лес наибольшего веса в G i. c i (B) c i ( B' i )

Шаг 6 z vCvC e E(B) z y x С B i α(e,C) eCeC e vCvC y x С B i α(e,C) eCeC

Упражнение 4.2 Пусть (V,T 1 ) и (V,T 2 ) два дерева на одном множестве вершин V. Доказать, что для любого ребра e T 1 существует ребро f T 2 такое, что (V,(T 1 \ {e}) U {f}) и (V,(T 2 \ {f}) U {e}) деревья.

Упражнение 4.3 Дан граф G с произвольными весами c: E(G) R. Найти остовный подграф в G наименьшего веса.

Упражнение 4.4 Дан граф G с произвольными весами c: E(G) R. Найти остовное дерево T в G, такое что вес максимального ребра в T наименьший (max{c(e)|e T} min).