TSP: pyconcorde: https://github.com/jvkersch/pyconcorde

pyconcorde.py

# pyconcorde: https://github.com/jvkersch/pyconcorde
from concorde.tsp import TSPSolver

num=30
X = np.random.random(num)*10**6
Y = np.random.random(num)*10**6

solver = TSPSolver.from_data(X, Y, norm='EUC_2D')
solution = solver.solve()

print(solution.found_tour)
print(solution.optimal_value/(10**6))

path = np.concatenate([solution.tour, [0]])
print(path)