Taddeus Kroes 13 жил өмнө
parent
commit
df4f504153
1 өөрчлөгдсөн 20 нэмэгдсэн , 0 устгасан
  1. 20 0
      problem44.py

+ 20 - 0
problem44.py

@@ -0,0 +1,20 @@
+#!/usr/bin/env python
+from math import sqrt
+from sys import exit
+
+def P(n):
+    return n * (3 * n - 1) / 2
+
+def is_pent(n):
+    return ((1 + sqrt(24 * n + 1)) / 6.).is_integer()
+
+for k in xrange(2, 10000):
+    Pk = P(k)
+
+    for j in range(k - 1, 0, -1):
+        Pj = P(j)
+        D = Pk - Pj
+
+        if is_pent(D) and is_pent(Pk + Pj):
+            print D
+            exit()