Solution to https://twitter.com/anniek_p/status/1245351315213074432

permutahedron.py

from itertools import permutations
import networkx as nx
import matplotlib.pyplot as plt

def flips(x):
    for idx, (a, b) in enumerate(zip(x[:-1], x[1:])):
        yield x[:idx] + (b,) + (a,) + x[idx+2:]

G = nx.Graph()
for perm in permutations(range(4)):
    for flipped in flips(perm):
        G.add_edge(perm, flipped)

nx.draw(G, with_labels=True)
plt.show()