# Algoritmo de Floyd-Warshall

En informática, el algoritmo de Floyd-Warshall, descrito en 1959 por Bernard Roy, es un algoritmo de análisis sobre grafos para encontrar el camino mínimo en grafos dirigidos ponderados. El algoritmo encuentra el camino entre todos los pares de vértices en una única ejecución. El algoritmo de Floyd-Warshall es un ejemplo de programación dinámica.

## El algoritmo de Warshall

El algoritmo de Warshall es un ejemplo de algoritmo booleano. A partir de una tabla inicial compuesta de 0`s (no hay correspondencia inicial en el grafo) y 1`s (hay una correspondencia, llamase “flecha”, entre nodos), obtiene una nueva matriz denominada “Matriz de Clausura Transitiva” en la que se muestran todas las posibles uniones entre nodos, directa o indirectamente. Es decir, si de “A” a “B” no hay una “flecha”, es posible que si haya de “A” a “C” y luego de “C” a “B”. Luego, este resultado se verá volcado en la matriz final.

## El algoritmo de Floyd

El algoritmo de Floyd es muy similar, pero trabaja con grafos ponderados. Es decir, el valor de la “flecha” que representamos en la matriz puede ser cualquier número real o infinito. Infinito marca que no existe unión entre los nodos. Esta vez, el resultado será una matriz donde estarán representadas las distancias mínimas entre nodos, seleccionando los caminos más convenientes según su ponderación (“peso”). Por ejemplo, si de “A” a “B” hay 36 (km), pero de “A” a “C” hay 2(km) y de “C” a “B” hay 10 (km), el algoritmo nos devolverá finalmente que de “A” a “B” hay 12 (km).
Los pasos a dar en la aplicación del algoritmo de Floyd son los siguientes:
* Formar las matrices iniciales C y D, donde C es la matriz de adyacencia, y D es una matriz del mismo tamaño cargada con valores iniciales Dij = i.
* Se toma k=1.
* Se selecciona la fila y la columna k de la matriz C y entonces, para i y j, con i≠k, j≠k e i≠j, hacemos:
Si (Cik + Ckj) < Cij → Dij = Dkj y Cij = Cik + Ckj
En caso contrario, dejamos las matrices como están.
* Si k ≤ n, aumentamos k en una unidad y repetimos el paso anterior, en caso contrario paramos las iteraciones.
* La matriz final C contiene los costes óptimos para ir de un vértice a otro, mientras que la matriz D contiene los penúltimos vértices de los caminos óptimos que unen dos vértices, lo cual permite reconstruir cualquier camino óptimo para ir de un vértice a otro.

## Algoritmo

El algoritmo de Floyd-Warshall compara todos los posibles caminos a través del grafo entre cada par de vértices. El algoritmo es capaz de hacer esto con sólo $V^{3}$ comparaciones (esto es notable considerando que puede haber hasta $V^{2}$ aristas en el grafo, y que cada combinación de aristas se prueba). Lo hace mejorando paulatinamente una estimación del camino más corto entre dos vértices, hasta que se sabe que la estimación es óptima.
Sea un grafo $G$ con conjunto de vértices $V$, numerados de 1 a N. Sea además una función ${\textrm {caminoMinimo}}(i,j,k)$ que devuelve el camino mínimo de $i$ a $j$ usando únicamente los vértices de 1 a $k$ como puntos intermedios en el camino. Ahora, dada esta función, nuestro objetivo es encontrar el camino mínimo desde cada $i$ a cada $j$ usando únicamente los vértices de 1 hasta $k+1$.
Hay dos candidatos para este camino: un camino mínimo, que utiliza únicamente los vértices del conjunto $(1...k)$; o bien existe un camino que va desde $i$ hasta $k+1$, y de $k+1$ hasta $j$, que es mejor. Sabemos que el camino óptimo de $i$ a $j$ que únicamente utiliza los vértices de 1 hasta $k$ está definido por ${\textrm {caminoMinimo}}(i,j,k)$, y está claro que si hubiera un camino mejor de $i$ a $k+1$ a $j$, la longitud de este camino sería la concatenación del camino mínimo de $i$ a $k+1$ (utilizando vértices de $(1...k)$) y el camino mínimo de $k+1$ a $j$ (que también utiliza los vértices en $(1...k)$).
Por lo tanto, podemos definir $caminoMinimo(i,j,k)$ de forma recursiva:
${\textrm {caminoMinimo}}(i,j,k)=\min({\textrm {caminoMinimo}}(i,j,k-1),{\textrm {caminoMinimo}}(i,k,k-1)\,+\,{\textrm {caminoMinimo}}(k,j,k-1));\,\!$
${\textrm {caminoMinimo}}(i,j,0)={\textrm {pesoArista}}(i,j);\,\!$
Esta fórmula es la base del algoritmo Floyd-Warshall. Funciona ejecutando primero $caminoMinimo(i,j,1)$ para todos los pares $(i,j)$, usándolos para después hallar $caminoMinimo(i,j,2)$ para todos los pares $(i,j)$... Este proceso continúa hasta que $k=n$, y habremos encontrado el camino más corto para todos los pares de vértices $(i,j)$ usando algún vértice intermedio.

## Pseudocodigo

Convenientemente, cuando calculamos el k-esimo caso, se puede sobreescribir la información salvada en la computación k -1. Esto significa que el algoritmo usa memoria cuadrática. Hay que cuidar la inicialización de las condiciones:

## Comportamiento con ciclos negativos

Para que haya coherencia numérica, Floyd-Warshall supone que no hay ciclos negativos (de hecho, entre cualquier pareja de vértices que forme parte de un ciclo negativo, el camino mínimo no está bien definido porque el camino puede ser infinitamente pequeño). No obstante, si hay ciclos negativos, Floyd-Warshall puede ser usado para detectarlos. Si ejecutamos el algoritmo una vez más, algunos caminos pueden decrementarse pero no garantiza que, entre todos los vértices, caminos entre los cuales puedan ser infinitamente pequeños, el camino se reduzca. Si los números de la diagonal de la matriz de caminos son negativos, es condición necesaria y suficiente para que este vértice pertenezca a un ciclo negativo.

## Ejemplo

Hallar el camino mínimo desde el vértice 3 hasta 4 en el grafo con la siguiente matriz de distancias:
    $D={\begin{bmatrix}0&3&5&1&\infty &\infty \\3&0&\infty &\infty &9&\infty \\5&\infty &0&7&7&1\\1&\infty &7&0&\infty &4\\\infty &9&7&\infty &0&\infty \\\infty &\infty &1&4&\infty &0\end{bmatrix}}$
Aplicamos el algoritmo de Floyd-Warshall, y para ello en cada iteración fijamos un vértice intermedio.
1ª Iteración: nodo intermedio = 1
La matriz es simétrica, por lo que solamente hará falta calcular el triángulo superior de las distancias.
d$_{23}$ = min(d$_{23}$, d$_{21}$ + d$_{13}$) = 8
d$_{24}$ = min(d$_{24}$, d$_{21}$ + d$_{14}$) = 4
d$_{25}$ = min(d$_{25}$, d$_{21}$ + d$_{15}$) = 9
d$_{26}$ = min(d$_{26}$, d$_{21}$ + d$_{16}$) = $\infty$
d$_{32}$ = min(d$_{32}$, d$_{31}$ + d$_{12}$) = 8
d$_{34}$ = min(d$_{34}$, d$_{31}$ + d$_{14}$) = 6
d$_{35}$ = min(d$_{35}$, d$_{31}$ + d$_{15}$) = 7
d$_{36}$ = min(d$_{36}$, d$_{31}$ + d$_{16}$) = 1
d$_{45}$ = min(d$_{45}$, d$_{41}$ + d$_{15}$) = $\infty$
d$_{46}$ = min(d$_{46}$, d$_{41}$ + d$_{16}$) = 4
d$_{56}$ = min(d$_{56}$, d$_{51}$ + d$_{16}$) = $\infty$
La matriz de distancia después de esta iteración es:
    ${\mathcal {W}}_{1}={\begin{bmatrix}0&3&5&1&\infty &\infty \\3&0&8&4&9&\infty \\5&8&0&6&7&1\\1&4&6&0&\infty &4\\\infty &9&7&\infty &0&\infty \\\infty &\infty &1&4&\infty &0\end{bmatrix}}$
2ª Iteración: nodo intermedio = 2
d$_{13}$ = min(d$_{13}$, d$_{12}$ + d$_{23}$) = 5
d$_{14}$ = min(d$_{14}$, d$_{12}$ + d$_{24}$) = 1
d$_{15}$ = min(d$_{15}$, d$_{12}$ + d$_{25}$) = 12
d$_{16}$ = min(d$_{16}$, d$_{12}$ + d$_{26}$) = $\infty$
d$_{34}$ = min(d$_{34}$, d$_{32}$ + d$_{24}$) = 6
d$_{35}$ = min(d$_{35}$, d$_{32}$ + d$_{25}$) = 7
d$_{36}$ = min(d$_{36}$, d$_{32}$ + d$_{26}$) = 1
d$_{45}$ = min(d$_{45}$, d$_{42}$ + d$_{25}$) = 13
d$_{46}$ = min(d$_{46}$, d$_{42}$ + d$_{26}$) = 4
d$_{56}$ = min(d$_{56}$, d$_{52}$ + d$_{26}$) = $\infty$
La matriz de distancia después de esta iteración es:
    ${\mathcal {W}}_{2}={\begin{bmatrix}0&3&5&1&12&\infty \\3&0&8&4&9&\infty \\5&8&0&6&7&1\\1&4&6&0&13&4\\12&9&7&13&0&\infty \\\infty &\infty &1&4&\infty &0\end{bmatrix}}$
3ª Iteración: nodo intermedio = 3
d$_{12}$ = min(d$_{12}$, d$_{13}$ + d$_{32}$) = 3
d$_{14}$ = min(d$_{14}$, d$_{13}$ + d$_{34}$) = 1
d$_{15}$ = min(d$_{15}$, d$_{13}$ + d$_{35}$) = 12
d$_{16}$ = min(d$_{16}$, d$_{13}$ + d$_{36}$) = 6
d$_{24}$ = min(d$_{24}$, d$_{23}$ + d$_{34}$) = 4
d$_{25}$ = min(d$_{25}$, d$_{23}$ + d$_{35}$) = 9
d$_{26}$ = min(d$_{26}$, d$_{23}$ + d$_{36}$) = 9
d$_{45}$ = min(d$_{45}$, d$_{43}$ + d$_{35}$) = 13
d$_{46}$ = min(d$_{46}$, d$_{43}$ + d$_{36}$) = 4
d$_{56}$ = min(d$_{56}$, d$_{53}$ + d$_{36}$) = 8
La matriz de distancia después de esta iteración es:
    ${\mathcal {W}}_{3}={\begin{bmatrix}0&3&5&1&12&6\\3&0&8&4&9&9\\5&8&0&6&7&1\\1&4&6&0&13&4\\12&9&7&13&0&8\\6&9&1&4&8&0\end{bmatrix}}$
4ª Iteración: nodo intermedio = 4
d$_{12}$ = min(d$_{12}$, d$_{14}$ + d$_{42}$) = 3
d$_{13}$ = min(d$_{13}$, d$_{14}$ + d$_{43}$) = 5
d$_{15}$ = min(d$_{15}$, d$_{14}$ + d$_{45}$) = 12
d$_{16}$ = min(d$_{16}$, d$_{14}$ + d$_{46}$) = 5
d$_{23}$ = min(d$_{23}$, d$_{24}$ + d$_{43}$) = 8
d$_{25}$ = min(d$_{25}$, d$_{24}$ + d$_{45}$) = 9
d$_{26}$ = min(d$_{26}$, d$_{24}$ + d$_{46}$) = 8
d$_{35}$ = min(d$_{35}$, d$_{34}$ + d$_{45}$) = 7
d$_{36}$ = min(d$_{36}$, d$_{34}$ + d$_{46}$) = 1
d$_{56}$ = min(d$_{56}$, d$_{54}$ + d$_{46}$) = 8
La matriz de distancia después de esta iteración es:
    ${\mathcal {W}}_{4}={\begin{bmatrix}0&3&5&1&12&5\\3&0&8&4&9&8\\5&8&0&6&7&1\\1&4&6&0&13&4\\12&9&7&13&0&8\\5&8&1&4&8&0\end{bmatrix}}$
5ª Iteración: nodo intermedio = 5
d$_{12}$ = min(d$_{12}$, d$_{15}$ + d$_{52}$) = 3
d$_{13}$ = min(d$_{13}$, d$_{15}$ + d$_{53}$) = 5
d$_{14}$ = min(d$_{14}$, d$_{15}$ + d$_{54}$) = 1
d$_{16}$ = min(d$_{16}$, d$_{15}$ + d$_{56}$) = 5
d$_{23}$ = min(d$_{23}$, d$_{25}$ + d$_{53}$) = 8
d$_{24}$ = min(d$_{24}$, d$_{25}$ + d$_{54}$) = 4
d$_{26}$ = min(d$_{26}$, d$_{25}$ + d$_{56}$) = 8
d$_{34}$ = min(d$_{34}$, d$_{35}$ + d$_{54}$) = 6
d$_{36}$ = min(d$_{36}$, d$_{35}$ + d$_{56}$) = 1
d$_{46}$ = min(d$_{46}$, d$_{45}$ + d$_{56}$) = 4
La matriz de distancia después de esta iteración es:
    ${\mathcal {W}}_{5}={\mathcal {W}}_{4}={\begin{bmatrix}0&3&5&1&12&5\\3&0&8&4&9&8\\5&8&0&6&7&1\\1&4&6&0&13&4\\12&9&7&13&0&8\\5&8&1&4&8&0\end{bmatrix}}$
6ª Iteración: nodo intermedio = 6
d$_{12}$ = min(d$_{12}$, d$_{16}$ + d$_{62}$) = 3
d$_{13}$ = min(d$_{13}$, d$_{16}$ + d$_{63}$) = 5
d$_{14}$ = min(d$_{14}$, d$_{16}$ + d$_{64}$) = 1
d$_{15}$ = min(d$_{15}$, d$_{16}$ + d$_{65}$) = 12
d$_{23}$ = min(d$_{23}$, d$_{26}$ + d$_{63}$) = 8
d$_{24}$ = min(d$_{24}$, d$_{26}$ + d$_{64}$) = 4
d$_{25}$ = min(d$_{25}$, d$_{26}$ + d$_{65}$) = 9
d$_{34}$ = min(d$_{34}$, d$_{36}$ + d$_{64}$) = 5
d$_{35}$ = min(d$_{35}$, d$_{36}$ + d$_{65}$) = 7
d$_{45}$ = min(d$_{45}$, d$_{46}$ + d$_{65}$) = 12
La matriz de distancia después de esta iteración es:
    ${\mathcal {W}}_{6}={\begin{bmatrix}0&3&5&1&12&5\\3&0&8&4&9&8\\5&8&0&5&7&1\\1&4&5&0&12&4\\12&9&7&12&0&8\\5&8&1&4&8&0\end{bmatrix}}$
Ya se han hecho todas las iteraciones posibles. Por tanto, el camino mínimo entre 2 vértices cualesquiera del grafo será el obtenido en la matriz final. En este caso, el camino mínimo entre 3 y 4 vale 5.

## Análisis

Si utilizamos matrices booleanas, para encontrar todos los ${\mathit {n}}^{2}$ de ${\mathcal {W}}_{k}$ desde ${\mathcal {W}}_{{\mathit {k}}-1}$ se necesita hacer $2{\mathit {n}}^{2}$ operaciones binarias. Debido a que empezamos con ${\mathcal {W}}_{0}={\mathcal {W}}_{\mathcal {R}}$ y computamos la secuencia de ${\mathit {n}}$ matrices booleanas ${\mathcal {W}}_{1}$, ${\mathcal {W}}_{2}$, $...$, ${\mathcal {W}}_{\mathit {n}}={\mathcal {M}}_{{\mathcal {R}}^{*}}$, el número total de operaciones binarias es de ${\mathit {n}}\times 2{\mathit {n}}^{2}=2{\mathit {n}}^{3}$. Por lo tanto, la complejidad del algoritmo es $\Theta ({n}^{3})$ y puede ser resuelto por una máquina determinista de Turing en tiempo polinómico.

## Aplicaciones y generalizaciones

El algoritmo de Floyd-Warshall puede ser utilizado para resolver los siguientes problemas, entre otros:
* Camino mínimo en grafos dirigidos (algoritmo de Floyd).
* Cierre transitivo en grafos dirigidos (algoritmo de Warshall). Es la formulación original del algoritmo de Warshall. El grafo es un grafo no ponderado y representado por una matriz booleana de adyacencia. Entonces la operación de adición es reemplazada por la conjunción lógica(AND) y la operación menor por la disyunción lógica (OR).
* Encontrar una expresión regular dada por un lenguaje regular aceptado por un autómata finito (algoritmo de Kleene).
* Inversión de matrices de números reales (algoritmo de Gauss-Jordan).
* Ruta óptima. En esta aplicación es interesante encontrar el camino del flujo máximo entre 2 vértices. Esto significa que en lugar de tomar los mínimos con el pseudocodigo anterior, se coge el máximo. Los pesos de las aristas representan las limitaciones del flujo. Los pesos de los caminos representan cuellos de botella; por ello, la operación de adición anterior es reemplazada por la operación mínimo.
* Comprobar si un grafo no dirigido es bipartito.

## Implementación del algoritmo de Floyd-Warshall

* Implementación en C - joshuarobinson.net Archivado el 7 de septiembre de 2011 en Wayback Machine. (en inglés). (Este enlace ya no es válido)
* Implementación en PHP - julmis.julmajanne.com Archivado el 29 de septiembre de 2011 en Wayback Machine. (gracias a Janne Mikkonen).
* Implementación en Java (explicación paso a paso) - explicación y applet disponible en pms.ifi.lmu.de (en inglés)