Fordův–Fulkersonův algoritmus

boolean NajdiCestu(Node s) {
 for ( Node u in V(G) )
  stav[u] = FRESH;
 p[s] = +s;
 d[s] = ∞;
 stav[s] = OPEN;
 do  { 
  u = libovolný otevřený uzel;
  stav[u] = CLOSED;
  for (Node v in Nasl(u)) {
   if ((stav[v] == FRESH) && (f(u,v) < q(u,v))) { 
    stav[v] = OPEN;
    p[v] = +u;
    d[v] = min(d[u], q(u,v) – f(u,v));    
   } 
  }
  for (Node v in předch(u)) {
   if (stav[v] == FRESH) && (f(v,u) > 0)) { 
    stav[v] = OPEN;
    p[v] = -u;
    d[v] = min(d[u], f(v,u));
   }
  } 
 } while ((existuje otevřený uzel)  && (u != t));
 return (u == t);
}

void ZvyšTok(Node s) {
 x = t;
 delta = d[t];
 do {
  v = x;
  sgn = p[v];
  u = abs(sgn);
  if (sgn > 0)
   f(u,v) += delta;
  else
   f(v,u) -= delta;
  x = u;
 } while ( v != s );
}

Tento algoritmus začíná tak, že všem hranám přiřadí tok hodnoty nula, neboli nic sítí na počátku neprotéká. Dále pak začne testovat nalezení zlepšující cesty, jestliže tato cesta bude nalezena tak se tok zvýší. Naopak jestliže cesta nalezená už nebude, to znamená že zlepšující cesta už neexistuje (viz obecná teorie k algoritmu) je nalezen maximální tok.

Funkce NajdiCestu(S) na počátku své implementace přiřadí všem uzlům grafu jako stav uzlu hodnotu FRESH. To znamená, že tyto uzly jsou čisté neboli ještě nepoužité. Uzlu S se přiřadí p[s] směr kterým jde hrana a poté ještě delta s jako nekonečno jelikož je to uzel předaný parametrem a tento uzel se otevře. Poté startuje cyklus který na svém začátku vybere uzel který je otevřený (stav uzlu == otevřeno ) a uzavře ho a poté pro všechny jeho následníky zjistí že jestliže je následník, uzel fresh a jeho tok je nižší než kapacita hrany jež spojuje uzel s tímto následníkem, tak ho otevře, přiřadí mu směr kterým je orientovaná hrana a vypočte pro něj delta. Delta se vypočítává právě pro to aby bylo možné zjistit o kolik se dá navýšit tok v hranách orientovaných směrem od zdroje k spotřebiči a snížit tok na hraně směřující od spotřebiče ke zdroji (bilance zůstává stejná). Když se takto projedou všechny následovníci tak se začnou procházet všichni předchůdci tohoto uzlu. U nich se zjišťuje zda jsou FRESH a zda hrana spojující tyto dva uzly má vyšší tok než nula. Když je tato podmínka splněna nastaví se předchůdce na stav OPEN nastaví se mu orientace hrany a určí delta. A toto celé se opakuje dokud existuje alespoň nějaký otevřený uzel či dokud testovaný uzel u se nebude rovnat spotřebiči. A také na základě této naposledy zmiňované rovnosti bude funkce vracet návratovou hodnotu true nebo false.

Odpověď na otázku jak je možné, že tok může "téci" i proti směru orientace hrany v síti?

Vysvětleme si to například na vodovodním potrubí. Existuje jedna trubka, která bude mít orientaci proti směru vodovodu to znamená od domácnosti do vodárny. Jak je možné, ze by touto trubkou tekla voda proti její orientaci? Fyzicky to možné není, ale hodnota jejího toku tomu může ukazovat. Tuto hodnotu trubka má jen proto, že jestliže my potřebujeme touto opačně orientovanou trubkou poslat například 3 litry směrem vodárna→ domácnost, tak jestliže například 10 litru vody jí proteče ve směru domácnost→vodárna, tak ona svůj průtok zmenši třeba na 7 litrů a ty tři litry nechá protéct někde jinde..jinou trubkou..tím pádem navenek vypadá, že průtok směrem vodárna→ domácnost je 3.

V této variantě je pro každou hranu ${\displaystyle e\in E(G)}$ zadáno nejen horní omezení ${\displaystyle q(e)}$, ale i minimální požadavek ${\displaystyle r(e)\in R_{0}^{+}}$ na tok. Tok ${\displaystyle f\colon E(G)\to R_{0}^{+}}$ je přípustný, pokud každá hrana splňuje ${\displaystyle r(e)\leq f(e)\leq q(e)}$. Úlohou pak může být nalezení maximálního, minimálního, nebo libovolného přípustného toku sítí.

Nejprve se budeme zabývat hledáním libovolného přípustného toku. Úlohu převedeme na hledání maximálního toku v síti ${\displaystyle G'}$ bez minimálních požadavků, která ze zadané sítě ${\displaystyle G}$ vznikne následujícími úpravami:

• přidáme hranu ze stoku do zdroje s nekonečnou kapacitou
• přidáme dva nové uzly, fiktivní zdroj ${\displaystyle s'}$ a fiktivní stok ${\displaystyle t'}$, které budou v síti ${\displaystyle G'}$ zdrojem a stokem
• za každou hranu ${\displaystyle e=(u,v)}$ sítě ${\displaystyle G}$ s kladným minimálním požadavkem ${\displaystyle r(e)}$ přidáme do ${\displaystyle G'}$ hrany ${\displaystyle (u,t')}$ a ${\displaystyle (s',v)}$, kterým nastavíme kapacity na ${\displaystyle r(e)}$; minimální požadavek na hraně ${\displaystyle e}$ zrušíme a nastavíme jí kapacitu ${\displaystyle q'(e)=q(e)-r(e)}$
• pokud nyní ze zdroje ${\displaystyle s'}$ vychází více hran do jednoho vrcholu, sloučíme tyto hrany do jedné o kapacitě rovné součtu kapacit původních hran; to samé uděláme, pokud do ${\displaystyle t'}$ vchází více hran ze stejného vrcholu

V síti ${\displaystyle G'}$ nalezneme maximální tok ${\displaystyle f'}$. Nejprve uvažujme případ, kdy každá hrana vycházející z ${\displaystyle s'}$ je nasycena, a tedy i každá hrana vcházející do ${\displaystyle t'}$ je nasycena. To znamená, že za každou hranu ${\displaystyle e=(u,v)}$ původní sítě ${\displaystyle G}$ s kladným ${\displaystyle r(e)}$ teče tok velikosti ${\displaystyle r(e)}$ z ${\displaystyle u}$ do ${\displaystyle t'}$ a z ${\displaystyle s'}$ do ${\displaystyle v}$. Tento tok přesměrujeme, aby tekl hranou ${\displaystyle e}$. Hranou ${\displaystyle e}$ tak teče tok velikosti ${\displaystyle f'(e)+r(e)}$, což nepřekračuje její kapacitu, protože ${\displaystyle f'(e)\leq q'(e)=q(e)-r(e)}$. Po následném odstranění uzlů a hran nově přidaných do ${\displaystyle G'}$ dostaneme přípustný tok sítí ${\displaystyle G}$. V případě, kdy některá hrana vycházející z ${\displaystyle s'}$ nasycena není, znamená to, že v ${\displaystyle G}$ neexistuje žádný přípustný tok. Pokud by totiž existoval, mohli bychom ho přeměnit na tok sítí ${\displaystyle G'}$, kde by každá hrana vycházející z ${\displaystyle s'}$ nasycena byla a tento tok by byl větší než ${\displaystyle f'}$.

Tento přípustný tok můžeme následně zlepšovat přidáváním, nebo ubíráním toku po zlepšujících cestách. V případě hledání maximálního toku je cesta ${\displaystyle P=\langle s=u_{0},u_{1},\ldots ,u_{n}=t\rangle }$, zlepšující, pokud pro každou dvojici ${\displaystyle (u_{i},u_{i+1})}$ platí, že ${\displaystyle (u_{i},u_{i+1})}$ je hranou ${\displaystyle G}$ a ${\displaystyle f((u_{i},u_{i+1}))<q((u_{i},u_{i+1}))}$, nebo ${\displaystyle (u_{i+1},u_{i})}$ je hranou ${\displaystyle G}$ a ${\displaystyle f((u_{i+1},u_{i}))>r((u_{i+1},u_{i}))}$. Po této cestě můžeme tok zvětšit. V případě hledání minimálního toku je ${\displaystyle P=\langle s=u_{0},u_{1},\ldots ,u_{n}=t\rangle }$, zlepšující, pokud pro každou dvojici ${\displaystyle (u_{i},u_{i+1})}$ platí, že ${\displaystyle (u_{i},u_{i+1})}$ je hranou ${\displaystyle G}$ a ${\displaystyle f((u_{i},u_{i+1}))>r((u_{i},u_{i+1}))}$, nebo ${\displaystyle (u_{i+1},u_{i})}$ je hranou ${\displaystyle G}$ a ${\displaystyle f((u_{i+1},u_{i}))<q((u_{i+1},u_{i}))}$. Po této cestě můžeme tok zmenšit.