Add solution to 2025-12-09

Author: fuglede

2025/day09/solutions.py

@@ -0,0 +1,26 @@
+from itertools import combinations
+from shapely import Polygon, box
+
+with open("input") as f:
+    ls = f.read().strip().split()
+
+ns = [tuple(map(int, line.split(","))) for line in ls]
+poly = Polygon(ns)
+
+answer1 = answer2 = 0
+
+for n1, n2 in combinations(ns, 2):
+    x1, y1 = n1
+    x2, y2 = n2
+    rect_size = (abs(x1 - x2) + 1) * (abs(y1 - y2) + 1)
+    answer1 = max(answer1, rect_size)
+
+    rect = box(min(x1, x2), min(y1, y2), max(x1, x2), max(y1, y2))
+    if poly.contains(rect):
+        answer2 = max(answer2, rect_size)
+
+# Part 1
+print(answer1)
+
+# Part 2
+print(answer2)