Add solution to 2025-12-08

Author: fuglede

2025/day08/solutions.py

@@ -0,0 +1,19 @@
+from itertools import combinations
+from math import dist, prod
+
+from scipy.cluster.hierarchy import DisjointSet
+
+with open("input") as f:
+    ls = f.read().strip().split("\n")
+
+ns = [tuple(map(int, l.split(","))) for l in ls]
+edges = sorted(combinations(ns, 2), key=lambda pair: dist(*pair))
+target = edges[1000]
+
+ds = DisjointSet(ns)
+while ds.n_subsets > 1:
+    if (edge := edges.pop(0)) == target:
+        print(prod(sorted(map(len, ds.subsets()))[-3:]))
+    ds.merge(*edge)
+
+print(prod(next(zip(*edge))))