Алгоритм двух китайцев — различия между версиями

Материал из Викиконспекты
Перейти к: навигация, поиск
(Реализация)
(Реализация)
Строка 80: Строка 80:
 
  Ребро - структура из трех чисел - {откуда ребро, куда ребро, вес ребра}.
 
  Ребро - структура из трех чисел - {откуда ребро, куда ребро, вес ребра}.
 
  root - текущий корень.
 
  root - текущий корень.
 +
особенность реал
 
   
 
   
 
  проверяем, можно ли дойти из <tex>v</tex> до остальных вершин. Если можно - запускаем <tex>findMST</tex>
 
  проверяем, можно ли дойти из <tex>v</tex> до остальных вершин. Если можно - запускаем <tex>findMST</tex>
Строка 86: Строка 87:
 
     int minEdge[n]; // создаем массив минимумов, входящих в каждую компоненту, инициализируем бесконечностью.
 
     int minEdge[n]; // создаем массив минимумов, входящих в каждую компоненту, инициализируем бесконечностью.
 
     for each <tex>e \in E</tex>
 
     for each <tex>e \in E</tex>
minEdge[e.to] = min(e.w, minEdge[e.to])
+
        minEdge[e.to] = min(e.w, minEdge[e.to])
 
     for each <tex>v \in V, v != root</tex>
 
     for each <tex>v \in V, v != root</tex>
 
         res += minEdge[v] //веса минимальных ребер точно будут в результате
 
         res += minEdge[v] //веса минимальных ребер точно будут в результате
 
     edge zeroEdges[] //создаем массив нулевых ребер
 
     edge zeroEdges[] //создаем массив нулевых ребер
 
     for each <tex>e \in E</tex>
 
     for each <tex>e \in E</tex>
if (текущее ребро равно по весу минимальному, входящему в эту вершину)
+
        if (текущее ребро равно по весу минимальному, входящему в эту вершину)
 
             zeroEdges.pushback(e);
 
             zeroEdges.pushback(e);
     вцвфцвфц
+
     dfs(root, zeroEdges) // проверяем, можно ли дойти до всех вершин по нулевым ребрам
// нашли все нулевые дуги
+
    if (можно дойти)
graph G1 = graph(G.n, G.root, zeroEdges);
+
        return res
zeroEdges.clear();
+
    int newComponents[n]; // будущие компоненты связности
if (isConnect(G1)){
+
    newComponents = Сondensation(zeroEdges)  
return res;
+
    edge newEdges[] //создаем массив ребер в новом графе с вершинами в сконденсированными компонентами
}
+
vector<pair<pair<int, int>, long long>> Edges = vector<pair<pair<int, int>, long long>>();
vector<int> comp1 = Components(G1);
 
//теперь у нас есть новые компоненты
 
G1.r.clear();
 
vector<pair<pair<int, int>, long long>> Edges = vector<pair<pair<int, int>, long long>>();
 
 
for (int i = 0; i < G.r.size(); i++){
 
for (int i = 0; i < G.r.size(); i++){
 
if (comp1[G.r[i].first.first] != comp1[G.r[i].first.second])
 
if (comp1[G.r[i].first.first] != comp1[G.r[i].first.second])

Версия 17:57, 11 декабря 2012

Алгоритм двух китайцев — алгоритм построения минимального остовного дерева во взвешенном ориентированном графе с корнем в заданной вершине. Был разработан математиками Чу Йонджином и Лю Цзенхонгом.

Постановка задачи

Дан взвешенный ориентированный граф [math]G(V, E)[/math] и начальная вершина [math]v[/math]. Требуется построить корневое остовное дерево в [math]G[/math] с корнем в вершине [math]v[/math], сумма весов всех ребер которого минимальна.

Алгоритм

Описание

Если хотя бы одна вершина графа [math]G[/math] недостижима из [math]v[/math], то требуемое дерево построить нельзя.

  1. Для каждой вершины [math]u \ne v[/math] графа [math]G[/math] произведём следующую операцию: найдём ребро минимального веса, входящее в [math]u[/math], и вычтем вес этого ребра из весов всех рёбер, входящих в [math]u[/math]. [math]m(u) = \min \limits_{tu \in E}w(tu), w'(tu) = w(tu) - m(u)[/math].
  2. Строим граф [math]K = (V,K_0)[/math], где [math]K_0[/math] — множество рёбер нулевого веса графа [math]G[/math] c весовой функцией [math]w'[/math]. Если в этом графе найдётся остовное дерево с корнем в [math]v[/math], то оно и будет искомым.
  3. Если такого дерева нет, то построим граф [math]C[/math] — конденсацию графа [math]K[/math]. Пусть [math]y[/math] и [math]z[/math] — две вершины графа [math]C[/math], отвечающие компонентам сильной связности [math]Y[/math] и [math]Z[/math] графа [math]K[/math] соответственно. Положим вес ребра между вершинами [math]y[/math] и [math]z[/math] равным минимальному среди весов рёбер графа [math]G[/math] с весовой функцией [math]w'[/math], идущих из [math]Y[/math] в [math]Z[/math].
  4. Продолжим с пункта 2, используя граф [math]C[/math] вместо [math]G[/math].
  5. В [math]C[/math] построено MST [math]T[/math]. Построим теперь MST [math]T'[/math] в [math]G[/math] с весовой функцией [math]w'[/math]. Добавим к [math]T'[/math] все вершины компоненты сильной связности графа [math]K[/math], которой принадлежит [math]v[/math] (по путям нулевого веса из [math]v[/math]). Пусть в [math]T[/math] есть ребро [math]yz[/math], где [math]y[/math] отвечает компоненте сильной связности [math]Y[/math], а [math]z[/math] — компоненте сильной связности [math]Z[/math] графа [math]K[/math]. Между [math]Y[/math] и [math]Z[/math] в графе [math]G[/math] с весовой функцией [math]w'[/math] есть ребро [math]y'z'[/math], вес которого равен весу ребра [math]yz[/math]. Добавим это ребро к дереву [math]T'[/math]. Добавим к [math]T'[/math] все вершины компоненты [math]Z[/math] по путям нулевого веса из [math]z'[/math]. Сделаем так для каждого ребра дерева [math]T[/math].
  6. Полученное дерево [math]T'[/math] — MST в графе [math]G[/math].



Пример

Исходный граф. КитайГраф1.png
Произведем спуск до нулевых ребер (Фаза 1, 2). КитайГраф2.png
По нулевым ребрам нельзя дойти до всех вершин из [math]v[/math], поэтому строим конденсацию и добавляем наименьшие ребра между компонентами (Фаза 3).

Найдем [math]MST[/math] для данного графа.

КитайГраф3.png
Произведем спуск до нулевых ребер (Фаза 1, 2). КитайГраф4.png
По нулевым ребрам нельзя дойти до всех вершин из [math]v[/math], поэтому строим конденсацию и добавляем наименьшие ребра между компонентами (Фаза 3).

Найдем [math]MST[/math] для данного графа.

КитайГраф5.png
Произведем спуск до нулевых ребер (Фаза 1, 2). По полученным нулевым ребрам можно дойти из корня до всех вершин. Тогда запускаем [math]dfs[/math] из корня и возвращаем ребра. КитайГраф6.png
Находим корень в каждой из компонент, из него запускаем [math]dfs[/math] по нулевым ребрам, возвращаем результат. КитайГраф7.png
Находим корень в каждой из компонент, из него запускаем [math]dfs[/math] по нулевым ребрам. Полученое дерево и есть [math]MST[/math]. КитайГраф8.png

Корректность



Замечания:

  • После перевзвешивания в каждую вершину, кроме [math]v[/math], входит по крайней мере одно ребро нулевого веса.
  • Пусть [math]T[/math] — искомое дерево в [math]G[/math] с весовой функцией [math]w[/math]. [math]w'(T) = w(T) - \sum \limits_{u \in V \setminus v}m(u)[/math], т.е. [math]T[/math] - MST в [math]G[/math] с весовой функцией [math]w[/math] тогда и только тогда, когда [math]T[/math] — MST в [math]G[/math] с весовой функцией [math]w'[/math].
Лемма:
Кратчайшее дерево путей [math]T'[/math] в графе [math]G[/math] можно получить, найдя кратчайшее дерево путей [math]T[/math] в графе [math]C[/math], а затем заменив в нем каждую компоненту сильной связности деревом, построенным из дуг нулевой длинны.
Доказательство:
[math]\triangleright[/math]
Зафиксируем любое дерево путей и покажем, что в графе [math]G[/math] найдется дерево не большей длины, имеющее такую структуру, как сказано в лемме. Для такой структуры дерева необходимо и достаточно, чтобы в каждое из подмножеств входило только по одному ребру. Меньше быть не может, иначе получится отдельная компонента связности. Если же в какое-то подмножество входит больше чем одно ребро, то все ребра кроме одного можно заменить ребрами нулевой длины, лежащими внутри подмножества, что разве лишь уменьшит длину дерева и не нарушит связности. Повторяя это преобразование нужное число раз мы добьемся искомой структуры дерева.
[math]\triangleleft[/math]

Из сделанных замечаний и леммы следует, что дерево [math]T'[/math] — MST в [math]G[/math].

Реализация



обозначения:
Граф хранится в виде множества ребер + индекс корня.
Множество ребер - список смежности.
Ребро - структура из трех чисел - {откуда ребро, куда ребро, вес ребра}.
root - текущий корень.
особенность реал

проверяем, можно ли дойти из [math]v[/math] до остальных вершин. Если можно - запускаем [math]findMST[/math]
findMST:
   int result = 0;
   int minEdge[n]; // создаем массив минимумов, входящих в каждую компоненту, инициализируем бесконечностью.
   for each [math]e \in E[/math]
       minEdge[e.to] = min(e.w, minEdge[e.to])
   for each [math]v \in V, v != root[/math]
       res += minEdge[v] //веса минимальных ребер точно будут в результате
   edge zeroEdges[] //создаем массив нулевых ребер
   for each [math]e \in E[/math]
       if (текущее ребро равно по весу минимальному, входящему в эту вершину)
           zeroEdges.pushback(e);
   dfs(root, zeroEdges) // проверяем, можно ли дойти до всех вершин по нулевым ребрам
   if (можно дойти)
       return res
   int newComponents[n]; // будущие компоненты связности
   newComponents = Сondensation(zeroEdges) 
   edge newEdges[] //создаем массив ребер в новом графе с вершинами в сконденсированными компонентами

vector<pair<pair<int, int>, long long>> Edges = vector<pair<pair<int, int>, long long>>(); for (int i = 0; i < G.r.size(); i++){ if (comp1[G.r[i].first.first] != comp1[G.r[i].first.second]) Edges.push_back(make_pair(make_pair(comp1[G.r[i].first.first], comp1[G.r[i].first.second]), G.r[i].second - minEdge[G.r[i].first.second])); } int compon(0); for (int i = 0; i < G.n; i++){ compon = max(compon, comp1[i]); } graph G2 = graph(compon + 1, comp1[G.root], Edges); Edges.clear(); res += findMST(G2); return res;

Сложность

Всего будет построено не более [math]V[/math] конденсаций. Конденсацию можно построить за [math]O(E)[/math]. Значит, алгоритм можно реализовать за [math]O(VE)[/math].

Источники

  • Романовский И. В. Дискретный анализ, 3-е изд., перераб. и доп. - СПб.:Невский Диалект; БХВ-Петербург, 2003. - 320 с.: ил. - ISBN 5-7940-0114-3
  • http://is.ifmo.ru