Goldbergův algoritmus

Je dána síť ${\displaystyle G(V,E)}$ s kapacitou z vrcholu u do vrcholu v danou jako ${\displaystyle c(u,v)}$, zdroj z a spotřebič s. Chceme najít maximální tok, který můžeme poslat přes síť ze z do s. Na vrcholech se provádějí dva typy operací: přemístění přebytku a zvedání vrcholu. Během algoritmu dále udržujeme:

• pro každou hranu hodnotu ${\displaystyle f(u,v)}$, reprezentující současný tok hranou z u do v.
• pro každou hranu volnou kapacitu, danou výrazem ${\displaystyle c(u,v)-f(u,v)}$.
• pro každý vrchol hodnotu ${\displaystyle vyska(u)}$. Operaci přemístění přebytku z vrcholu u do vrcholu v provádíme, pouze když ${\displaystyle vyska(u)>vyska(v)}$.
• pro každý vrchol hodnotu ${\displaystyle prebytek(u)}$, dánu jakou rozdíl součtu toků z u minus součet toků do u (případně pouze jako součet toků do u při zavedení pomocného značení, viz níže).

Z důvodů zjednodušení zápisu popisu algoritmu a implementace přidáme do grafu dočasně pro každou hranu opačně orientovanou hranu s nulovou kapacitou. Během výpočtu pak budeme vyžadovat, aby byl pro každou hranu vždy splněn vztah ${\displaystyle f(u,v)=-f(v,u)}$. Přidání pomocných hran lze samozřejmě udělat pouze v případě, že graf neobsahuje orientované smyčky. V případě existence smyček bychom museli postupovat o trochu opatrněji, rozmyšlení tohoto zřídkakdy se vyskytujícího případu ponecháme čtenáři za domácí cvičení :-). Na hodnoty f(u,v) dále klademe během výpočtu následující podmínky:

• ${\displaystyle \ f(u,v)\leq c(u,v)}$. Tok mezi ${\displaystyle u}$ a ${\displaystyle v}$ nepřevyšuje kapacitu.
• ${\displaystyle \ f(u,v)=-f(v,u)}$. Udržujeme korektní tok (vlnu) sítě.
• ${\displaystyle \ \sum _{v}f(u,v)=prebytek(u)\geq 0}$ pro všechny vrcholy ${\displaystyle u\neq z}$. Pouze zdroj může generovat tok.

Zobrazení ${\displaystyle f}$ splňující tyto podmínky nenazýváme přípustným tokem, tak jak je používán například ve Ford-Fulkensonově metodě, ale pouze vlnou. Rozdíl je v tom, že vlna povoluje kladný přebytek. Vlnu používáme pouze během výpočtu, výstupem algoritmu je přípustný tok (speciální případ vlny). Všimněte si, že poslední podmínka pro přebytek je odvozená z korespondující podmínky pro přípustný tok v běžné síti.

Pozorujeme, že nejdelší možná cesta ze z do s je ${\displaystyle |V|}$ vrcholů dlouhá. Proto pro každý přípustný tok existuje přiřazení výšky vrcholům splňující, že ${\displaystyle vyska(z)=|V|}$ a ${\displaystyle vyska(t)=0}$, a je-li tok hranou z u do v kladný, pak ${\displaystyle vyska(u)>vyska(v)}$. Během přemisťování přebytku vrcholů se vlna chová stejně jako vodní vlna v přírodě—přebytek se přelévá pouze z vrcholů s větší výškou do níže položených. Na rozdíl od Ford-Fulkersonova algoritmu, tok (vlna) v síti nemusí nutně být přípustným tokem během výkonu algoritmu.

Zjednodušeně, nejprve naplníme všechny hrany vedoucí ze zdroje. Dále během výpočtu zvyšujeme výšky vrcholů (kromě z a s) a přebytky se jako vlna šíří mezi vrcholy, dokud veškeré přebytky nedosáhnou spotřebiče. Poté pokračujeme ve zvyšování výšky vnitřních vrcholů sítě, až dokud se přebytky, které nedosáhly spotřebiče, nevrátí zpět do zdroje. Vrchol může dosáhnout nejvýše výšky ${\displaystyle 2|V|-1}$, protože nejdelší možná cesta z libovolného vrcholu kromě s zpátky do z je ${\displaystyle |V|-1}$ dlouhá a ${\displaystyle vyska(z)=|V|}$. Výška s se udržuje na 0.

Implementace v programovacím jazyku Python verze 2:

def goldberg(C, zdroj, spotrebic):
    n = len(C) # C je matice kapacit
    F = [[0] * n for _ in xrange(n)]
    # zbytková kapacita z u do v je C[u][v] - F[u][v]

    vyska = [0] * n # výška vrcholu
    prebytek = [0] * n # rozdíl toku do a z vrcholu
    vyzkousen = [0] * n # sousedi vyzkoušení od posledního zvedání vrcholu
    # "seznam" vrcholů
    seznam   = [i for i in xrange(n) if i != zdroj and i != spotrebic]

    def zvedni(u, v):
        posli = min(prebytek[u], C[u][v] - F[u][v])
        F[u][v] += posli
        F[v][u] -= posli
        prebytek[u] -= posli
        prebytek[v] += posli

    def zvedni(u):
        # najdi nejmenší novou výšku a udělej přemístění přebytku možným,
        # jestli je vůbec možné něco zvednout
        min_vyska = vyska[u]
        for v in xrange(n):
            if C[u][v] - F[u][v] > 0:
                min_vyska = min(min_vyska, vyska[v])
                vyska[u] = min_vyska + 1

    def vyprazdni(u):
        while prebytek[u] > 0:
            if vyzkousen[u] < n: # zkontroluj dalšího souseda
                v = vyzkousen[u]
                if C[u][v] - F[u][v] > 0 and vyska[u] > vyska[v]:
                    zvedni(u, v)
                else:
                    vyzkousen[u] += 1
            else: # zkontrolovali jsme všechny sousedy. musíme ho zvednout
                zvedni(u)
                vyzkousen[u] = 0

    vyska[zdroj] = n   # nejdelší cesta ze zdroje do spotřebiče je kratší než n
    prebytek[zdroj] = Inf # pošli tolik toku, kolik se jen dá sousedovi zdroje
    for v in xrange(n):
        zvednout(zdroj, v)

    p = 0
    while p < len(seznam):
        u = seznam[p]
        old_vyska = vyska[u]
        vyprazdni(u)
        if vyska[u] > old_vyska:
            seznam.insert(0, seznam.pop(p)) # přesuň na začátek seznamu
            p = 0 # začni ze začátku seznamu
        p += 1

    return sum([F[zdroj][i] for i in xrange(n)])

Všimněte si, že tato implementace není moc efektivní. Je pomalejší než Edmonds-Karpův algoritmus dokonce pro hodně husté grafy. Pro zrychlení můžete udělat nejméně dvě věci:

1. Udělejte seznam sousedů pro každý vrchol a nechť je index vyzkousen[u] iterátor, místo intervalu ${\displaystyle 0..n-1}$.
2. Použijte heuristiku mezer. Jestli existuje ${\displaystyle k}$, že pro žádný vrchol ${\displaystyle u}$ neplatí ${\displaystyle vyska(u)=k}$, můžete přiřadit ${\displaystyle vyska(u)=\max(vyska(u),vyska(z)+1)}$ pro všechny vrcholy kromě ${\displaystyle z}$, pro které ${\displaystyle vyska(u)>k}$, protože takové ${\displaystyle k}$ reprezentuje minimální řez grafu a mezi vrcholy ${\displaystyle Z=\{u|vyska(u)>k\}}$ a ${\displaystyle T=\{v|vyska(v)<k\}}$ už nepřeteče víc. Kdyby ${\displaystyle (Z,S)}$ nebyl minimální řez, existovala by hrana ${\displaystyle (u,v)}$, pro niž by ${\displaystyle u\in Z}$, ${\displaystyle v\in S}$ a ${\displaystyle c(u,v)-f(u,v)>0}$. Ale pak by ${\displaystyle vyska(u)}$ nebyla zvětšena víc než ${\displaystyle vyska(v)+1}$, což je v rozporu s předpoklady ${\displaystyle vyska(u)>k}$ a ${\displaystyle vyska(v)<k}$.

• Tok v síti
• Fordův–Fulkersonův algoritmus
• Dinicův algoritmus

• Thomas H. Cormen, Charles E. Leiserson, Ronald L. Rivest, and Clifford Stein. Introduction to Algorithms, Second Edition. MIT Press and McGraw-Hill, 2001. ISBN 0-262-03293-7. Section 26.4: Push-relabel algorithms, and section 26.5: The relabel-to-front-algorithm.
• Jakub Černý: Základní grafové algoritmy
• Martin Mareš, Tomáš Valla: Průvodce labyrintem algoritmů