3-bipcon/3-bipcon.tex

   1 %%%%%%%%%%%%
   2 % Zápisek tøetího semináre z grafových agoritmù - ze dne 20.3.2006
   3 % Zapsal Jiøí Peinlich - peinlich@seznam.cz a Michal Kùrka - michal.kurka@gmail.com
   4 %%%%%%%%%%%
   5
   6 \input ../sgr.tex
   7
   8 \prednaska{3}{Bipartitní párování a globální k-souvislost}{zapsali Jiøí Peinlich, Michal Kùrka}
   9
  10 \h{Maximální párování v $k$-regulárním bipartitním grafu}
  11
  12 \s{Operace Degree Split} provádí rozdìlení sudì-regulárního grafu na dva
  13 podgrafy s polovièní regularitou. Operaci Degree Split definujme na
  14 sudì-regulárním grafu takto: V grafu najdeme eulerovský tah. Sestrojíme dva
  15 grafy, které budou mít oba stejnou mno¾inu vrcholù jako graf pùvodní a ka¾dý
  16 bude mít polovinu jeho hran. Mno¾inu hran prvního grafu budou tvoøit sudé hrany
  17 z nalezeného eulerovského tahu, mno¾inu hran druhého grafu pak hrany liché. Tuto
  18 operaci lze jistì provést v lineárním èase (v lineárním èase najdeme eulerovský
  19 tah i rozklad na sudé a liché hrany).
  20
  21 Dále budeme pracovat s $2^d$-regulárními grafy. Operací Degree Split tedy získáme dva grafy, které budou $2^{d-1}$-regulární.
  22
  23 Provedeme-li operaci Degree Split $\log k = d$ krát (polovinu hran v¾dy
  24 zahodíme), získáme 1-regulární podgraf a tedy párování pro zadaný graf (ve
  25 skuteènosti mù¾eme graf rozkládat na párování a vyrobit si 1-faktorizaci
  26 grafu). Slo¾itost nalezení párování bude tedy pro $2^d$ regulární grafy $\O(2^d n
  27 d)$.
  28
  29 Pokud zadaný graf nebude $2^d$-regulární budeme muset pøidat hrany tak, aby nový
  30 graf tuto vlastnost mìl. Operaci Degree Split pak budeme provádìt tak, abychom
  31 se k párování blí¾ili.
  32
  33 Místo toho, abychom do grafu hrany jen pøidávali, budeme v pøípadech, kdy je to
  34 mo¾né, pouze zvìt¹ovat násobnost hrany. U ka¾dé hrany si tedy budeme pamatovat
  35 její násobnost.
  36
  37 \s{Degree Split s~násobnostmi:} Pro sudì-regulární grafy s násobnostmi zavedeme operaci Degree Split takto:
  38 Graf $G=(V,E)$ rozdìlíme na dva grafy $G_1$ a $G_2$, bude platit $V(G_1) = V(G_2) = V$. Hrany nyní pøidìlíme následujícím zpùsobem:
  39 \algo
  40 \:Pokud $e\in E$ v $G$ má sudou násobnost (znaèíme $n(e)$), umístíme ji do $E_1$ i do $E_2$ s násobností ${n(e) \over 2}$, v opaèném pøípadì pøidáme do obou grafù hranu s násobností $\lfloor {n(e) \over 2} \rfloor$.
  41 \:Graf se zbylých hran má v¹echny stupnì sudé a je bez multiplicit. Provedeme na nìj pùvodní operaci Degree Split a rozdìlené mno¾iny hran pøidáme do $G_1$ a $G_2$.
  42 \endalgo
  43
  44 Celý proces lze stihnout v èase $\O(m)$ ($m$ je poèet hran $G$), nebo» v první èásti u ka¾dé hrany pouze zjistíme, zda má sudou násobnost, pøidáme nové hrany v konstantním èase (upravíme násobnosti), a v druhé èásti se provede Split, který má té¾ lineární slo¾itost. Operace Degree Split má tedy slo¾itost $\O(m)$ i v grafu s násobnostmi.
  45
  46 Mìjme nyní $k$-regulární bipartitní graf. Zvolme $t$ tak aby $2^t\geq kn$.
  47 Zvolme dále
  48 $\alpha := \lfloor {2^t \over k} \rfloor$ a
  49 $\beta := 2^t \bmod k$.
  50 Do grafu pøidáme hrany a upravíme násobnosti hran tak, aby byl $2^t$ regulární. Dále pøidáme triviální párování ($i$-tý vrchol vlevo se spojí s $i$-tým vrcholem vpravo) s násobností $\beta$. Tuto mno¾inu hran oznaème $F$. Platí $\beta < k \Rightarrow \vert F \vert < 2^t$.
  51
  52 Takto získáme $2^t$-regulární graf. Na tento graf budeme aplikovat operaci Degree Split a vybereme si v¾dy tu polovinu, kde bude ménì hran z $F$. Tímto zpùsobem budem graf dìlit dokud budou stupnì vrcholù vìt¹í ne¾ jedna. Tedy $t$-krát. Poslední takto získaný graf bude 1-regulární (párování). V ka¾dém kroku se zbavíme alespoò poloviny hran z $F$. Provedeme to ceklem $t$-krát a tedy výsledné párování bude perfektní párování zadaného grafu.
  53
  54 Slo¾itost algoritmu je $\O(kn \log n)$, proto¾e inicializace algoritmu se dá provést v lineárním èase, provede se $\log (kn)$ iterací po $\O(m)$.
  55
  56 \h{Algoritmy na hledání globální k-souvislosti}
  57 \s{Hranová k-souvislost}
  58 Problém zji¹tìní stupnì souvislosti grafu lze pøevést na problém hledání minimálního øezu v~grafu. Algoritmus pro hledání minimálního øezu lze zkonstruovat napøíklad tak, ¾e pro ka¾dé dva vrcholy $s$, $t$ zjistíme minimální $st$-øez. Algoritmus mù¾eme zrychlit, pokud si uvìdomíme, ¾e jeden z vrcholù $s$ nebo $t$ lze zvolit pevnì. Pokud pracujeme s orientovanými grafy, musíme projít jak $s \rightarrow t$ øezy tak i $t \rightarrow s$ øezy. Algoritmus bude mít slo¾itost $\O(n^{{5/3}}m)$.
  59
  60 U vrcholové $k$-souvislosti to ov¹em tak snadno nepùjde. Pokud by toti¾ fixovaný vrchol byl souèástí nìjakého minimálního separátoru, algoritmus mù¾e selhat. Pøesto ale nemusíme procházet v¹echny dvojice vrcholù. Staèí si pamatovat, kolik vrcholù $s$ jsme u¾ pro v¹echny $t$ zkontrolovali a nejmen¹í zatím nalezený separátor. Kdy¾ bude poèet vrcholù vìt¹í ne¾ nejmen¹í separátor, tak u¾ jsme jistì na¹li jeden z minimálních øezù. Slo¾itost takového algoritmu pak bude $\O(\kappa (G) n^{3/2} m)$, kde $\kappa(G)$ je stupeò souvislosti $G$, který hledáme.
  61
  62 Bez tokù se dá ov¹em najít hranová k-souvislost efektivnìji (i v ohodnocených grafech). Uka¾me si, jak na to.
  63
  64 \h{Algoritmus pro nalezení minimálního globálního øezu v neorientovaných grafech (Namagochi, Ibaraki)}
  65
  66 \noindent Zavedeme znaèení:
  67
  68 \itemize\ibull
  69 \:{$r(u,v)$ buï kapacita minimálního $uv$-øezu}
  70 \:{$d(v)$ buï celková kapacita hran vedoucích z $v$}
  71 \:{$d(U)$ buï kapacita hran vedoucích z mno¾iny $U \subseteq V$, speciálnì $d(U)=d(\overline{U})$}
  72 \:{$d(P,Q)$ buï kapacita hran vedoucích mezi mno¾inami $P,Q$}
  73 \endlist
  74
  75 \s{Definice}
  76 {\it Legálním uspoøádáním vrcholù} (LU) budeme nazývat lineární uspoøádání vrcholù $v_1, v_2,... ,v_n$ takové, ¾e platí
  77 $d(\{v_1 \ldots v_{i-1}\},v_i) \geq d(\{v_1 \ldots v_{i-1}\},v_j)$ pro $1 \leq i<j\leq n$.
  78
  79 \s{Lemma} Je-li $v_1 ... v_n$ LU na $G$, pak $r(v_{n-1},v_n)=d(v_n)$.
  80
  81 \s{Dùkaz} Buï $C$ nìjaký øez oddìlující $v_{n-1}$ a $v_n$. Utvoøme posloupnost vrcholù $u_i$ takto:
  82
  83 \algo
  84 \:$u_0 := v_1$
  85 \:$u_i := v_j$ tak, ¾e $j>i$, $v_i$ a $v_j$ jsou oddìleny øezem $C$ a $j$ je minimální takové.
  86 \endalgo
  87
  88 Ka¾dé $u_{i-1}$ je tedy buï $u_i$, pokud jsou $v_i$ a $v_{i-1}$ na stejné stranì øezu, nebo $u_{i-1}$ je $v_i$ pokud jsou $v_i$ a $v_{i-1}$ na opaèné stranì øezu.
  89 Dostáváme tedy, ¾e $d(\{v_1...v_{i-1}\},u_i)\leq d(\{v_1...v_{i-1}\},u_{i-1})$, proto¾e buïto  $u_i=u_{i-1}$ a pak je nerovnost splnìna jako rovnost, nebo je $u_i=v_j$, $j>i$ a nerovnost plyne z LU vrcholù $v_i$.
  90
  91 Chceme ukázat, ¾e velikost libovolného øezu je alespoò taková jako velikost øezu kolem vrcholu $v_n$.
  92 Platí, ¾e $ \vert C \vert \geq \sum_{i=1}^{n-1} d(v_i,u_i)$. Uká¾eme, ¾e pravá strana je alespoò $d(v_n)$.
  93
  94 $$\sum_{i=1}^{n-1} d(v_i,u_i) = \sum_{i=1}^{n-1} d(\{v_1..v_i\},u_i) - d(\{v_1 \ldots v_{i-1}\},u_i) \geq \sum_{i=1}^{n-1} d(\{v_1 \ldots v_i\},u_i) - d(\{v_1 .. v_{i-1}\},u_{i-1}) = $$
  95
  96 $$ = d(\{v_1 \ldots v_{n-1}\},u_{n-1}) - d(\{v_1 \ldots v_0\},u_0) = d(\{v_1 \ldots v_{n-1}\},v_n) - 0.$$
  97 \qed
  98
  99 Dokázali jsme, ¾e libovolný øez separující $v_{n-1}$ a $v_n$ je vìt¹í ne¾ jednoduchý øez skládající se jen z hran kolem~$v_n$. Kdy¾ tedy sestavíme nìjakou LU posloupnost vrcholù, budeme mít k dispozici jednoduchý minimální øez $v_{n-1}$ a~$v_n$. Následnì vytvoøíme graf $G'$, v nìm¾ $v_{n-1}$ a $v_n$ skontrahujeme. Rekurzivnì najdeme minimální øez v $G'$. Hledaný minimální øez poté buïto oddìluje vrcholy $v_n$ a $v_{n-1}$ a potom je øez kolem vrcholu $v_n$ minimální, nebo vrcholy $v_n$ a $v_{n-1}$ neoddìluje, a v takovém pøípadì jej najdeme rekurzivnì. Hledaný øez je tedy men¹í z rekurzivnì nalezeného øezu a øezu kolem $v_n$.
 100
 101 Zbývá ukázat, jak konstruovat LU. Postaèí hladovì: Pamatujeme si $\forall v\neq v_1 \ldots v_{i-1}$ hodnotu $d(\{v_1 \ldots v_{i-1},v)$, oznaème ji $z_v$. V ka¾dém kroku vybereme vrchol $v$ s maximální hodnotou $z_v$ a prohlásíme ho za $v_i$ a pøepoèítáme~$z_v$.
 102
 103 Zde se hodí datová struktura, která doká¾e rychle hledat maxima a zvy¹ovat hodnoty prvkù,
 104 napøíklad Fibonacciho halda. Ta zvládne \<DeleteMax> v~èase $(\log n)$ a \<Increase> v~$\O(1)$
 105 amortizovanì. Celkem pak ná¹ algoritmus bude mít slo¾itost $\O(n(m+n\log n))$ pro obecné kapacity.
 106
 107 Pokud jsou kapacity malá celá èísla, mù¾eme vyu¾ít balíèkové struktury. Budeme
 108 si udr¾ovat obousmìrný seznam zatím pou¾itých hodnot $z_v$, ka¾dý prvek takového
 109 seznamu bude obsahovat v¹echny vrcholy se spoleènou hodnotou $z_v$. Kdy¾ budeme
 110 mít seznam seøazený, vybrání minimálního prvku znamená pouze podívat se na
 111 první prvek seznamu a z nìj odebrat jeden vrchol, pøípadnì celý prvek ze seznamu
 112 odstranit. Operace \<Increase> poté bude reprezentovat pouze pøesunutí vrcholu o
 113 malý poèet pøihrádek, pøípadnì zalo¾ení nové pøihrádky na správném místì.
 114 \<DeleteMax> i \<Increase> pak budou mít slo¾itost $\O(1)$ amortizovanì a celý algoritmus $\O(mn)$.
 115
 116 \bye