X-Git-Url: http://mj.ucw.cz/gitweb/?a=blobdiff_plain;f=3-bipcon%2F3-bipcon.tex;h=362bd7b75a79cb904c7c8009f4040da281cb03ca;hb=b2484e9e3f714201f2c3543a2a265d015e6e5915;hp=81c3f0088983a3c55442cf49a8e6b4b7c60334ed;hpb=66de0031bc8248c38306aac0930ea56807f8d0d3;p=ga.git diff --git a/3-bipcon/3-bipcon.tex b/3-bipcon/3-bipcon.tex index 81c3f00..362bd7b 100644 --- a/3-bipcon/3-bipcon.tex +++ b/3-bipcon/3-bipcon.tex @@ -1,137 +1,150 @@ \input ../sgr.tex -\prednaska{3}{Bipartitní párování a globální k-souvislost}{zapsali Jiøí Peinlich, Michal Kùrka} - -\>V~minulé kapitole jsme se zabývali aplikacemi tokù na~hledání maximálního párování -a minimálního øezu. V~této si pøedvedeme dva algoritmy pro podobné problémy, -které se obejdou bez tokù. - -\h{Maximální párování v regulárním bipartitním grafu \cite{alon:matching}} - -Nejprve si nadefinujme operaci {\I Degree Split,} která dostane jako vstup libovolný -$2k$-regulární graf $G=(V,E)$ a rozdìlí ho na~podgrafy $G_1=(V,E_1)$ a $G_2=(V,E_2)$, které budou -oba $k$-regulární. Tuto operaci mù¾eme snadno provést v~lineárním èase tak, ¾e si graf -rozdìlíme na~komponenty, v~ka¾dé nalezneme eulerovský tah a jeho sudé hrany dáme do~$G_1$ -a liché do~$G_2$. - -To nám pomù¾e ke~snadnému algoritmu pro nalezení maximálního párování ve~$2^d$-regulárním -bipartitním grafu.\foot{V¹imnìte si, ¾e takové párování bude v¾dy perfektní (viz Hallova vìta).} -Staèí provést Degree Split na~dva $2^{d-1}$-regulární grafy, na~libovolný jeden z~nich -aplikovat dal¹í Degree Split atd., a¾ se dostaneme k~$1$-regulárnímu grafu, který -je perfektním párováním v~$G$. To v¹e jsme schopni stihnout v~lineárním èase, -jeliko¾ velikosti grafù, které splitujeme, exponenciálnì klesají. Také bychom -mohli rekurzivnì zpracovávat obì komponenty a tak se v~èase $\O(m\log n)$ dobrat -ke~kompletní 1-faktorizaci zadaného grafu.\foot{To je rozklad hran grafu na~disjunktní -perfektní párování a má ho ka¾dý regulární bipartitní graf.} - -Pokud zadaný graf nebude $2^d$-regulární, pomù¾eme si tím, ¾e ho novými hranami -doplníme na $2^d$-regulární a pak si pøi splitech budeme vybírat ten podgraf, -do~kterého padlo ménì nových hran, a uká¾eme, ¾e nakonec v¹echny zmizí. -Abychom graf pøíli¹ nezvìt¹ili, budeme sna¾it místo pøidávání úplnì nových -hran pouze zvy¹ovat násobnost hran existujících. Pro ka¾dou hranu $e$ si tedy -budeme pamatovat její násobnost $n(e)$. - -{\I Degree Split grafu s~násobnostmi} pak budeme provádìt následovnì: hranu~$e$ s~násobností $n(e)$ umístíme do~$G_1$ -i~do~$G_2$ s~násobností $\lfloor n(e)/2 \rfloor$ a pokud bylo $n(e)$ liché, pøidáme hranu do~pomocného grafu -$G^\prime$. V¹imnìte si, ¾e $G^\prime$ bude sudì-regulární graf bez násobností, tak¾e na~nìj mù¾eme aplikovat pùvodní -Degree Split a $G^\prime_i$ pøiøadit ke~$G_i$. To~v¹e zvládneme v~èase $\O(m)$. - -Mìjme nyní $k$-regulární bipartitní graf. Zvolme $t$ tak aby $2^t\geq kn$. -Zvolme dále parametry +\prednaska{3}{Bipartitní párování a globální k-souvislost}{} + +V~předeÅ¡lých kapitolách jsme se zabývali aplikacemi toků na~hledání maximálního párování +a minimálního $st$-řezu. Nyní si předvedeme dva algoritmy pro podobné problémy, +které se obejdou bez toků. + +\h{Maximální párování v regulárním bipartitním grafu \cite{alon:matching}} + +Nejprve si nadefinujme operaci {\I Å¡těpení grafu,} která zadaný graf $G=(V,E)$ +se vÅ¡emi vrcholy sudého stupně a sudým počtem hran v~každé komponentě souvislosti +rozloží na~disjunktní podgrafy $G_1=(V,E_1)$ a $G_2=(V,E_2)$, +v~nichž bude pro každý vrchol~$v$ platit ${\rm deg}_{G_1}(v) = {\rm +deg}_{G_2}(v) = {\rm deg}_G(v)/2$. Tuto operaci můžeme snadno +provést v~lineárním čase tak, že si graf rozdělíme na~komponenty, v~každé +nalezneme eulerovský tah a jeho hrany budeme přidávat střídavě do~$G_1$ a do~$G_2$. + +Å těpení nám pomůže ke~snadnému algoritmu pro nalezení maximálního párování ve~$2^t$-regulárním +bipartitním grafu.\foot{VÅ¡imněte si, že takové párování bude vždy perfektní (viz Hallova věta).} +Komponenty takového grafu mají určitě sudý počet hran, takže jej můžeme +rozÅ¡těpit na~dva $2^{t-1}$-regulární grafy. Libovolný jeden z~nich pak opět rozÅ¡těpíme +atd., až dostaneme $1$-regulární graf, který je perfektním párováním v~$G$. +To vÅ¡e jsme schopni stihnout v~lineárním čase, jelikož velikosti grafů, které +Å¡těpíme, exponenciálně klesají. Také bychom mohli rekurzivně zpracovávat obě +části a tak se v~čase $\O(m\log n)$ dobrat ke~kompletní 1-faktorizaci +zadaného grafu.\foot{To je rozklad hran grafu na~disjunktní perfektní párování +a má ho každý regulární bipartitní graf.} + +Pokud zadaný graf nebude $2^t$-regulární, pomůžeme si tím, že ho novými hranami +doplníme na $2^t$-regulární a pak si při Å¡těpeních budeme vybírat ten podgraf, +do~kterého padlo méně nových hran, a ukážeme, že nakonec vÅ¡echny zmizí. +Abychom graf příliÅ¡ nezvětÅ¡ili, budeme se snažit místo přidávání úplně nových +hran pouze zvyÅ¡ovat násobnost hran existujících. Pro každou hranu $e$ si tedy +budeme pamatovat její násobnost $n(e)$. + +{\I Å těpení grafu s~násobnostmi} pak budeme provádět následovně: hranu~$e$ +s~násobností $n(e)$ umístíme do~$G_1$ i~do~$G_2$ s~násobností $\lfloor n(e)/2 +\rfloor$ a pokud bylo $n(e)$ liché, přidáme hranu do~pomocného grafu +$G^\prime$. VÅ¡imněte si, že $G^\prime$ bude graf bez násobností, v~němž mají +vÅ¡echny vrcholy sudý stupeň, takže na~něj můžeme aplikovat původní Å¡těpící algoritmus +a $G^\prime_i$ přiřadit ke~$G_i$. To~vÅ¡e zvládneme v~čase $\O(m)$. + +Mějme nyní $k$-regulární bipartitní graf. Obě jeho partity jsou stejně velké, +označme si počet vrcholů v~každé z~nich~$n$. Zvolme $t$ tak aby $2^t\geq kn$. +Zvolme dále parametry $\alpha := \lfloor 2^t/k \rfloor$ a $\beta := 2^t \bmod k$. -Ka¾dé pùvodní hranì nastavíme násobnost~$\alpha$ a pøidáme triviální párování~$F$ -($i$-tý vrchol vlevo se spojí s~$i$-tým vrcholem vpravo) s~násobností~$\beta$. -V¹imnìte si, ¾e $\betai$, $v_i$ a $v_j$ jsou oddìleny øezem $C$ a $j$ je minimální takové. +\:$u_i := v_j$ tak, že $j>i$, $v_i$ a $v_j$ jsou odděleny řezem $C$ a $j$ je minimální takové. +[Tedy $v_j$ je nejbližší vrchol na~druhé straně řezu.] \endalgo -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. Dostáváme tedy, ¾e $d(\{v_1\ldots v_{i-1}\},u_i)\leq d(\{v_1\ldots -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$. - -Chceme ukázat, ¾e velikost libovolného øezu je alespoò taková, jako velikost øezu kolem vrcholu $v_n$. -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)$: +Každé $u_{i-1}$ je tedy buď rovno $u_i$, pokud jsou $v_i$ a $v_{i-1}$ na stejné straně řezu, nebo rovno $v_i$, pokud +jsou $v_i$ a $v_{i-1}$ na~stranách opačných. Z~toho dostáváme, že $d(\{v_1\ldots v_{i-1}\},u_i)\leq d(\{v_1\ldots +v_{i-1}\},u_{i-1})$, protože buďto $u_{i-1}=u_i$, a pak je nerovnost splněna jako rovnost, nebo je $u_{i-1}=v_i$ a +nerovnost plyne z~legálnosti uspořádání. +Chceme ukázat, že velikost naÅ¡eho řezu~$C$ je alespoň taková, jako velikost řezu kolem vrcholu $v_n$. +VÅ¡imneme si, že $\vert C \vert \geq \sum_{i=1}^{n-1} d(v_i,u_i)$. Hrany mezi $v_i$ a $u_i$ jsou totiž navzájem +různé a každá z~nich je součástí řezu~$C$. Ukážeme, že pravá strana je alespoň $d(v_n)$: $$\eqalign{ -\sum_{i=1}^{n-1} d(v_i,u_i) &= \sum_{i=1}^{n-1} d(\{v_1\ldots 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 \ldots v_{i-1}\},u_{i-1}) = \cr -&= 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 = d(v_n).\cr +\sum_{i=1}^{n-1} d(v_i,u_i) &= \sum_{i=1}^{n-1} d(\{v_1\ldots v_i\},u_i) - d(\{v_1 \ldots v_{i-1}\},u_i) \geq \cr +&\geq \sum_{i=1}^{n-1} d(\{v_1 \ldots v_i\},u_i) - d(\{v_1 \ldots v_{i-1}\},u_{i-1}) = \cr +&= d(\{v_1 \ldots v_{n-1}\},u_{n-1}) - d(\{v_1 \ldots v_0\},u_0) = \cr +&=d(\{v_1 \ldots v_{n-1}\},v_n) - 0 = d(v_n).\cr }$$ \qed -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'$ (sestrojíme nové LU atd.). 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$. - -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$. - -Zde se hodí datová struktura, která doká¾e rychle hledat maxima a zvy¹ovat hodnoty prvkù, -napøíklad Fibonacciho halda. Ta zvládne \ v~èase $\O(\log n)$ a \ v~$\O(1)$ -amortizovanì. Celkem pak ná¹ algoritmus bude mít slo¾itost $\O(n(m+n\log n))$ pro obecné kapacity. - -Pokud jsou kapacity malá celá èísla, mù¾eme vyu¾ít pøíhrádkové struktury. Budeme -si udr¾ovat obousmìrný seznam zatím pou¾itých hodnot $z_v$, ka¾dý prvek takového -seznamu bude obsahovat v¹echny vrcholy se spoleènou hodnotou $z_v$. Kdy¾ budeme -mít seznam seøazený, vybrání minimálního prvku znamená pouze podívat se na -první prvek seznamu a z nìj odebrat jeden vrchol, pøípadnì celý prvek ze seznamu -odstranit. Operace \ poté bude reprezentovat pouze pøesunutí vrcholu o -malý poèet pøihrádek, pøípadnì zalo¾ení nové pøihrádky na správném místì. -\ i \ pak budou mít slo¾itost $\O(1)$ a celý algoritmus $\O(mn)$. +Dokázali jsme, že libovolný řez separující $v_{n-1}$ a $v_n$ je alespoň tak velký jako 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$ kontrahujeme. Rekurzivně najdeme minimální +řez v $G'$ (sestrojíme nové LU atd.). 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$. + +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$, prohlásíme ho za $v_i$ a přepočítáme~$z_v$. + +Zde se hodí datová struktura, která dokáže rychle hledat maxima a zvyÅ¡ovat hodnoty prvků, +například Fibonacciho halda. Ta zvládne \ v~čase $\O(\log n)$ a \ v~$\O(1)$ +amortizovaně. Celkem pak náš algoritmus bude mít složitost $\O(n(m+n\log n))$ pro obecné kapacity. + +Pokud jsou kapacity malá celá čísla, můžeme využít přihrádkové struktury. Budeme +si udržovat obousměrný seznam zatím použitých hodnot $z_v$, každý prvek takového +seznamu bude obsahovat vÅ¡echny vrcholy se společnou hodnotou $z_v$. Když budeme +mít seznam seřazený, vybrání minimálního prvku bude znamenat pouze podívat se na +první prvek seznamu a z něj odebrat jeden vrchol, případně celý prvek ze seznamu +odstranit. Operace \ poté bude reprezentovat pouze přesunutí vrcholu o +malý počet přihrádek, případně založení nové přihrádky na správném místě. +\ proto bude mít složitost $\O(1)$, vÅ¡echny \ dohromady $\O(m)$, +jelikož za~každou hranu přeskakujeme maximálně jednu přihrádku, a celý algoritmus $\O(mn)$. \references \bye