170
правок
Изменения
Нет описания правки
==Постановка задачи==
Даны матроиды <tex>M_1 = \langle S, I_1 \mathcal{I}_1 \rangle</tex> и <tex>M_2 = \langle S, I_2 \mathcal{I}_2 \rangle</tex>. Необходимо найти максимальное по мощности независимое множество в пересечении <tex>M_1</tex> и <tex>M_2</tex>.
==Алгоритм решения==
Пусть множество <tex>J \in (I_1 \mathcal{I}_1 \cap I_2\mathcal{I}_2)</tex>.
<br>Определим [[Граф замен для двух матроидов|граф замен]] <tex>D_{M_1, M_2}(J) = \langle S, A(J) \rangle</tex>, где
<tex>A(J) = \{(y, z) | y \in J, z \in S\setminus J, J - y + z \in I_1 \mathcal{I}_1 \} </tex> <tex>\cup \{ (z', y') | z' \in S \setminus J, y' \in J, J - y' + z' \in I_2 \mathcal{I}_2 \}</tex>.
Пусть <tex>X_1 = \{ z \in S \setminus J | J + z \in I_1 \mathcal{I}_1 \}</tex>, <tex>X_2 = \{ z \in S \setminus J | J + z \in I_2 \mathcal{I}_2 \}</tex>, <tex>P</tex> {{---}} кратчайший путь из <tex>X_1</tex> в <tex>X_2</tex> в графе <tex>D_{M_1, M_2}(J)</tex>. <tex>P</tex> может и не существовать.
{{Лемма
|statement =
{{Утверждение
|statement =
<tex>r_1 (U) \le leqslant |J \cap U|</tex>
|proof =
От противного. Пусть <tex>r_1 (U) > |J \cap U|</tex>, тогда <tex>\exists z \in U \setminus (J \cap U) : (J \cap U) + z \in I_1\mathcal{I}_1</tex> при том, что <tex>J + z \notin I_1\mathcal{I}_1</tex>. В противном случае <tex>J + z \in I_1\mathcal{I}_1, z \in X_1</tex>, то есть <tex>X_1 \cap U \ne \emptyset</tex>, что противоречит отсутствию пути из <tex>X_1</tex> в <tex>X_2</tex>. Так как <tex>(J \cap U) + z \in I_1\mathcal{I}_1</tex>, а <tex>J + z \notin I_1\mathcal{I}_1</tex>, <tex>\exists y \in J \setminus U : J - y + z \in I_1\mathcal{I}_1</tex>. Однако, тогда <tex>(y, z) \in A(J)</tex>, что противоречит факту <tex>\delta^- (U) = \emptyset</tex>.
}}
{{Утверждение
|statement =
<tex>r_2 (S \setminus U) \le leqslant |J \cap (S \setminus U)|</tex>
|proof =
От противного. Пусть <tex>\exists z \in (S \setminus U) \setminus J : J \cap (S \setminus U) + z \in I_2\mathcal{I}_2</tex>. Аналогично доказательству предыдущего утверждения <tex>\exists y \in J \setminus (S \setminus U) : J - y + z \in I_2\mathcal{I}_2</tex>. Однако <tex>J \setminus (S \setminus U) = J \cap U</tex>, то есть <tex>(z, y)</tex> {{---}} дуга в <tex>D_{M_1, M_2}(J)</tex>, поэтому <tex>z \in U</tex> (т.к. <tex>y \in U</tex>). Противоречие.
}}
Так как <tex>|J| = |J \cap U| + |J \setminus U| \ge geqslant r_1 (U) + r_2 (S \setminus U), |J| = r_1 (U) + r_2 (S \setminus U)</tex>. Следовательно, <tex>J</tex> {{---}} максимальное по мощности независимое множество в пересечении <tex>M_1</tex> и <tex>M_2</tex>.
}}
{{Лемма
|statement =
<tex>J' = J \bigtriangleup V(P) \in I_1 \mathcal{I}_1 \cap I_2\mathcal{I}_2</tex>
|proof =
[[Файл:Graph_DY.png|400px|thumb|right|Граф замен]]
Пусть <tex>P = z_0, y_1, z_1, ..., y_t, z_t; G = \{ z_1, ..., z_t \} \cup (J \setminus \{ y_1, ..., y_t \})</tex>. Тогда <tex>G \subseteq S, |G| = |J|</tex> и дуги из <tex>\{ y_1, ..., y_t \}</tex> в
<tex>\{ z_1, ..., z_t \}</tex> составляют единственное полное паросочетание в <tex>J \bigtriangleup G</tex>. То есть, согласно [[Лемма о единственном паросочетании в графе замен | лемме о единственном паросочетании в подграфе замен]], <tex>G \in I_1\mathcal{I}_1</tex>.К тому же, <tex>\forall i \ge geqslant 1\ z_i \notin X_1</tex>, иначе <tex>P</tex> {{---}} не кратчайший путь из <tex>X_1</tex> в <tex>X_2</tex>. Это означает, что <tex>z_i + J \notin I_1\mathcal{I}_1</tex>, то есть <tex>r_1 (J \cup G) = r_1 (J) = r_1 (G) = |G| = |J|</tex>. Так как <tex>J + z_0 \in I_1\mathcal{I}_1, G + z_0 \in I_1\mathcal{I}_1</tex> (т.е. <tex>J' = \{ z_0, z_1, ..., z_t \} \cup (J \setminus \{ y_1, ..., y_t \}) \in I_1\mathcal{I}_1</tex>.Симметрично <tex>G = \{ z_0, ..., z_{t - 1} \} \cup (J \setminus \{ y_1, ..., y_t \}), J' \in I_2\mathcal{I}_2</tex> и, следовательно, <tex>J' \in (I_1 \mathcal{I}_1 \cap I_2\mathcal{I}_2)</tex>.
}}
'''while''' '''not''' isMaximal
построить [[Граф замен для двух матроидов|граф замен]] <tex>D_{M_1, M_2}(J)</tex>
<tex>X_1 \leftarrow \{ z \in S \setminus J | J + z \in I_1 \mathcal{I}_1 \}</tex> <tex>X_2 \leftarrow \{ z \in S \setminus J | J + z \in I_2 \mathcal{I}_2 \}</tex>
<tex>P</tex> <tex>\leftarrow</tex> кратчайший путь из <tex>X_1</tex> в <tex>X_2</tex>
'''if''' <tex>P \ne \emptyset</tex>
|about=
Эдмондса - Лоулера
|statement= Пусть <tex>M_1=\langle X, I_1\mathcal{I}_1\rangle</tex>, <tex>M_2=\langle X, I_2\mathcal{I}_2\rangle</tex> {{---}} матроиды. Тогда <br><tex>\max\limits_{I \in I_1 \mathcal{I}_1 \cap I_2 \mathcal{I}_2 } |I| = \min\limits_{A \subseteq X} \left(r_1(A) + r_2(X \setminus A)\right)</tex>.
Где <tex>r_1</tex> и <tex>r_2</tex> {{---}} ранговые функции в первом и втором матроиде соответственно.
|proof=
[[Файл:El_graph.png|thumb|140px|right|Завершение алгоритма]]
<div>
Докажем неравенство <tex>\max\limits_{I \in I_1 \mathcal{I}_1 \cap I_2 \mathcal{I}_2 } |I| \le leqslant \min\limits_{A \subseteq X} \left(r_1(A) + r_2(X \setminus A)\right)</tex>
Выберем произвольные <tex>I \in I_1 \mathcal{I}_1 \cap I_2\mathcal{I}_2</tex>, <tex>A \subseteq X</tex>, тогда
<tex>|I| = |I \cap A| + |I \cap (X \setminus A)|</tex>
<tex>|I| = r_1(I \cap A) + r_2(I \cap (X \setminus A))</tex>
Но <tex>r_1(I \cap A) \le leqslant r_1(A)</tex> и <tex>r_2(I \cap (X \setminus A)) \le leqslant r_2(X \setminus A)</tex>, значит
<tex>|I| \le leqslant r_1(A) + r_2(X \setminus A)</tex>
В силу произвольности <tex>I</tex> и <tex>A</tex> получаем
<tex>\max\limits_{I \in I_1 \mathcal{I}_1 \cap I_2 \mathcal{I}_2 } |I| \le leqslant \min\limits_{A \subseteq X} \left(r_1(A) + r_2(X \setminus A)\right)</tex>
Конструктивно построим <tex>\forall M_1, M_2</tex> такие <tex>I \in I_1 \mathcal{I}_1 \cap I_2\mathcal{I}_2</tex> и <tex>A \subseteq X</tex>, что <tex>|I| = r_1(A) + r_2(X \setminus A)</tex>.
Обозначим <tex>S = \left\{x|I \cup \{x\} \in I_1\mathcal{I}_1\right\}</tex>, <tex>T = \left\{x|I \cup \{x\} \in I_2\mathcal{I}_2\right\}</tex>. Если <tex>S \cap T \ne \varnothing</tex>, добавим их пересечение в <tex>I</tex>.
Построим [[Граф замен для двух матроидов|граф замен]] <tex>G_I</tex>. Добавим вершину <tex>z</tex>, не влияющую на независимость в первом матроиде {{---}} из неё будут вести рёбра во все вершины множества <tex>S</tex>. Пусть <tex>p</tex> {{---}} кратчайший путь из <tex>S</tex> в <tex>T</tex>, <tex>p_1</tex> {{---}} путь <tex>p</tex> с добавленным в начало ребром из <tex>z</tex>. По [[Лемма о единственном паросочетании в графе замен|лемме о единственном паросочетании]] и [[Лемма о единственном паросочетании в подграфе замен, индуцированном кратчайшим путем|лемме о единственном паросочетании, индуцированном кратчайшем путём]] <tex>I \bigtriangleup p_1 \in I_2\mathcal{I}_2</tex>. Теперь добавим вершину <tex>u</tex>, не влияющую на независимость во втором матроиде {{---}} в неё будут вести рёбра из всех вершин множества <tex>T</tex>. Тогда <tex>p_2</tex> (путь <tex>p</tex> с добавленным ребром в <tex>u</tex>) — кратчайший путь из <tex>S</tex> в <tex>u</tex>. Аналогично, <tex>I \bigtriangleup p_2 \in I_1\mathcal{I}_1</tex>. Отсюда следует, что <tex>I \bigtriangleup p \in I_1 \mathcal{I}_1 \cap I_2\mathcal{I}_2</tex>, причём <tex>|I \bigtriangleup p| = |I| + 1</tex>.</div>
Будем таким образом увеличивать <tex>I</tex>, пока существует путь <tex>p</tex>. Рассмотрим момент, когда такого пути не нашлось.
Докажем, что <tex>r_1(A) = |I \cap A|</tex> от противного.
Пусть <tex>r_1(A) > |I \cap A|</tex>, тогда существует <tex>w \in A \setminus (I \cap A)</tex>, такое, что <tex>(I \cap A) \cup \{w\} \in I_1\mathcal{I}_1</tex>. Если <tex>I \cup \{w\} \in I_1\mathcal{I}_1</tex>, то <tex>w \in S</tex> и из <tex>S</tex> есть путь в <tex>A</tex>. Значит, <tex>I \cup \{w\} \notin I_1\mathcal{I}_1</tex>. Отсюда следует, что существует <tex>y \in I \setminus A</tex>, такое что <tex>I \setminus \{y\} \cup \{w\} \in I_1\mathcal{I}_1</tex>. Но тогда ребро <tex>yw</tex> имеется в графе, то есть из <tex>y</tex> существует путь в <tex>T</tex>, что противоречит условию <tex>y \in I \setminus A</tex>.
Следовательно, <tex>r_1(A) = |I \cap A|</tex>. Аналогично, <tex>r_2(\overline A) = |I \cap \overline A|</tex>. Отсюда <tex>r_1(A) + r_2(\overline A) = |I|</tex>, то есть при найденных <tex>I</tex> и <tex>A</tex> достигается равенство.
}}
== См. также==
* [[Пересечение_матроидов,_определение,_примеры]]
* [[Алгоритм_построения_базы_в_объединении_матроидов]]
== Источник Источники информации ==
''Chandra Chekuri'' — [http://www.cs.illinois.edu/class/sp10/cs598csc/Lectures/Lecture17.pdf '''Combinatorial Optimization''']
[[Категория:Алгоритмы и структуры данных]]
[[Категория:Матроиды]]