xyrolle

import numpy as np
import networkx as nx

def warshall(M):
    n = M.shape[0]
    W = M
    for k in range(n):
        for i in range(n):
            for j in range(n):
                W[i, j] = W[i, j] or (W[i, k] and W[k, j])