Построение компонент рёберной двусвязности

Основные понятия

• Реберная двусвязность
• Компонента реберной двусвязности

Построение компонент реберной двусвязности будет осуществляться с помощью обхода в глубину.

Двупроходный алгоритм

Первый способ найти искомые компоненты - сначала определить критерий перехода в новую компоненту реберной двусвязности, а затем покрасить вершины графа в нужные цвета.

Первый проход определяет для каждой вершины [math]v[/math] две величины: [math]enter(v)[/math] - время входа поиска в глубину в вершину, [math]return(v)[/math] - минимальное из времен входа вершин, достижимых из [math]v[/math] по дереву [math]dfs[/math] и не более, чем одному обратному ребру. [math]return(v)[/math] находится как [math]min(enter(v), return(u), enter(w))[/math] для всех [math]u[/math] - сыновей [math]v[/math] в дереве [math]dfs[/math], [math]w[/math] - соседей [math]v[/math] по обратным ребрам. Важно, что ребро к родителю дерева [math]dfs[/math] не является обратным ребром обхода.

Псевдокод первого прохода:

 void dfs(v, родитель):
   увеличиваем текущее время
   enter(v) := текущее время    
   return(v) := enter(v)
   для всех вершин u, смежных v:
     если enter(u) равен нулю (вершина не посещена):
       dfs(u, v)
       return(v) := min(return(v), return(u))
     иначе если u не родитель:
       return(v) := min(return(v), enter(u))
 ...
 обнуляем массив enter
 текущее время := 0
 для всех вершин v графа:
   если enter(v) = 0:
     dfs(v, null)

Определим критерий перехода к новой компоненте. Воспользуемся ранее доказанной леммой.

Основываясь на этом, определим алгоритм окраски вершин графа. Перешли по мосту, следовательно началась новая компонента.

Псевдокод второго прохода:

 void paint(v, цвет):
   colors(v) := цвет
   для всех вершин u, смежных v:
     если colors(u) равен нулю (вершина не покрашена):
       если return(u) = enter(u):
         увеличиваем максимальный цвет
         paint(u, максимальный цвет)
       иначе:
         paint(u, цвет)
 ...
 обнуляем массив colors
 максимальный цвет := 0
 для всех вершин v графа:
   если colors(v) = 0:
     увеличиваем максимальный цвет
     paint(v, максимальный цвет)

Вершины каждой из компонент реберной двусвязности окажутся окрашенными в свой цвет.

Время работы алгоритма будет время работы двух запусков dfs, то есть 2 * [math] O(|V| + |E|)[/math], что есть [math] O(|V| + |E|)[/math].

Однопроходный алгоритм

Можно также искать компоненты реберной двусвязности путем конкатенации циклов. Воспользуемся тем, что реберная двусвязность является отношением эквивалентности на вершинах графа; тогда, если у двух циклов существует хоть одна общая вершина, все вершины, располагающиеся на этих циклах, принадлежат одной компоненте. Более того, две вершины [math]u[/math] и [math]v[/math] лежат в одной компоненте реберной двусвязности тогда и только тогда, когда существует последовательность простых циклов [math]c_1 c_2 ... c_n[/math], причем [math]u \in c_1[/math], [math]v \in c_n[/math], и [math]c_i[/math] имеет с [math]c_{i + 1}[/math] хотя бы одну общую вершину для всех [math]i \in {1 ... n - 1}[/math]. Действительно, если зафиксировать один путь от [math]u[/math] до [math]v[/math], а затем искать точки пересечения второго, не имеющего одинаковых ребер с первым, пути с ним, то получится последовательность циклов, точками сочленения между которыми будут как раз точки пересечения путей. И наоборот, последовательность простых циклов легко превратить в два реберно непересекающихся пути.

Совокупность компонент реберной двусвязности будем хранить как систему непересекающихся множеств вершин.

Псевдокод:

 int dfs(v, родитель): (возвращает 0, если у v и ее потомков нет обратных ребер, и представителя множества, содержащего цикл с v, в обратном случае)
   seen(v) = true
   value = 0, result = 0
   для всех вершин u, смежных v:
     если не seen(u):
       value = dfs(u, v)
       если value > 0:
         color(v) = value
         result = value
       иначе если u не родитель:
         color(v) = color(u)
         result = color(v)
   return result
 ...
 обнуляем массив seen
 нумеруем вершины графа натуральными числами от 1 до мощности множества вершин графа
 для всех вершин v графа:
   color(v) = номер вершины (номер цвета соответствует номеру вершины-представителя в множестве)
 для всех вершин v графа:
   если не seen(v):
     dfs(v, null)

Осталось лишь сопоставить всем вершинам отдельно взятой компоненты единственного представителя.

Псевдокод:

 int relax(v): (возвращает нового представителя)
   если color(v) не равен номеру v:
     color(v) = relax(color(v))
   return color(v)
 ...
 для всех вершин v графа:
   relax(v)

Теперь две вершины имеют одинаковый цвет тогда и только тогда, когда они принадлежат одной компоненте реберной двусвязности.

Время работы dfs [math] O(|V| + |E|)[/math]. А время восстановления [math] O(|V|)[/math], так как снм. Итоговое время работы алгоритма [math] O(|V| + |E|)[/math].

См. также

• Oбхода в глубину
• Использование обхода в глубину для поиска точек сочленения
• Построение компонент вершинной двусвязности
• Использование обхода в глубину для поиска мостов
• Визуализация построение компонент реберной двусзяности

Литература

Седжвик Роберт. Фундаментальные алгоритмы на C++. Часть 5: Алгоритмы на графах: Пер. с англ./Роберт Седжвик. — СПб.: ООО «ДиаСофтЮП», 2002. — С. 123-128

В.А.Кузнецов, А.М.Караваев. "Оптимизация на графах" - Петрозаводск, Издательство ПетрГУ 2007