TSP heuristics
Shorter tours than nearest neighbour on six point sets, scored by the demo runner
| 1 | """Nearest neighbour, then 2-opt until no exchange of two edges shortens the tour.""" |
| 2 | import json |
| 3 | import math |
| 4 | import sys |
| 5 | |
| 6 | p = json.load(sys.stdin) |
| 7 | n = len(p) |
| 8 | |
| 9 | |
| 10 | def d(i, j): |
| 11 | return math.dist(p[i], p[j]) |
| 12 | |
| 13 | |
| 14 | left, tour = set(range(1, n)), [0] |
| 15 | while left: |
| 16 | nxt = min(left, key=lambda j: d(tour[-1], j)) |
| 17 | left.remove(nxt) |
| 18 | tour.append(nxt) |
| 19 | better = True |
| 20 | while better: |
| 21 | better = False |
| 22 | for i in range(1, n - 1): |
| 23 | for j in range(i + 1, n): |
| 24 | a, b, c, e = tour[i - 1], tour[i], tour[j], tour[(j + 1) % n] |
| 25 | if d(a, c) + d(b, e) < d(a, b) + d(c, e) - 1e-12: |
| 26 | tour[i:j + 1] = reversed(tour[i:j + 1]) |
| 27 | better = True |
| 28 | print(json.dumps(tour)) |