X-Git-Url: http://mj.ucw.cz/gitweb/?a=blobdiff_plain;f=3-bipcon%2F3-bipcon.tex;h=362bd7b75a79cb904c7c8009f4040da281cb03ca;hb=b2484e9e3f714201f2c3543a2a265d015e6e5915;hp=2b41a5416bef5b66de1f906e0aa5c17374632149;hpb=9caf01a8298a65d9ff4e7e40c9dd61cb6fb60c1d;p=ga.git diff --git a/3-bipcon/3-bipcon.tex b/3-bipcon/3-bipcon.tex index 2b41a54..362bd7b 100644 --- a/3-bipcon/3-bipcon.tex +++ b/3-bipcon/3-bipcon.tex @@ -1,118 +1,120 @@ \input ../sgr.tex -\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ì souvilosti -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 +\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é. -[Tedy $v_j$ je nejbli¾¹í vrchol na~druhé stranì øezu.] +\:$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ï 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=u_{i-1}$, a pak je nerovnost splnìna jako rovnost, nebo je $u_i=v_j$, $j>i$ a -nerovnost plyne z~legálnosti uspoøádá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)$: +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 \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 @@ -121,28 +123,28 @@ $$\eqalign{ }$$ \qed -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)$. +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