Taddeus Kroes 13 ani în urmă
părinte
comite
acdb419a1e
1 a modificat fișierele cu 23 adăugiri și 0 ștergeri
  1. 23 0
      problem50.py

+ 23 - 0
problem50.py

@@ -0,0 +1,23 @@
+#!/usr/bin/env python
+from math import sqrt
+from utils import primes_until
+
+def maxprime(n):
+    primes = list(primes_until(n))
+    maxp = maxlen = 2
+
+    for start in xrange(len(primes)):
+        for end in xrange(start + maxlen + 1, start + 800):
+            if end - start >= maxlen:
+                conseq = primes[start:end]
+                s = sum(conseq)
+
+                if len(conseq) > 1 and s in primes[start:]:
+                    maxlen = end - start
+                    maxp = s
+                    print maxp
+                    #print '%s = %d' % (' + '.join(map(str, primes[start:end])), s)
+
+    return maxp
+
+print maxprime(1000000)