aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorWalt Mankowski <waltman@pobox.com>2022-07-04 18:50:27 -0400
committerWalt Mankowski <waltman@pobox.com>2022-07-04 18:50:27 -0400
commitc5a37622868919d3af2d4d3011df7bfe6eebc689 (patch)
treeaafcb3b4b89a0ad464696b30c95e4e1fbd1d7a0f
parentd75934a5e11be1b9a56286bd169e9a1970d26b42 (diff)
downloadperlweeklychallenge-club-c5a37622868919d3af2d4d3011df7bfe6eebc689.tar.gz
perlweeklychallenge-club-c5a37622868919d3af2d4d3011df7bfe6eebc689.tar.bz2
perlweeklychallenge-club-c5a37622868919d3af2d4d3011df7bfe6eebc689.zip
python code to find primes
-rw-r--r--challenge-172/walt-mankowski/python/primes.py15
1 files changed, 15 insertions, 0 deletions
diff --git a/challenge-172/walt-mankowski/python/primes.py b/challenge-172/walt-mankowski/python/primes.py
new file mode 100644
index 0000000000..686143de0f
--- /dev/null
+++ b/challenge-172/walt-mankowski/python/primes.py
@@ -0,0 +1,15 @@
+from math import sqrt
+
+# find the primes up to n using the sieve of Eratosthenes and return
+# them as list
+def primes_to(n):
+ is_prime = [True] * (n+1)
+ is_prime[0] = is_prime[1] = False
+ for i in range(2, int(sqrt(n))+1):
+ if is_prime[i]:
+ for j in range(i+i, n+1, i):
+ is_prime[j] = False
+
+ return [i for i,val in enumerate(is_prime) if val]
+
+