lazy erathostenes algorithm

daedalus · Oct 16, 2023 · 02087e3 · 02087e3
1 parent f7a424a
commit 02087e3
Show file tree

Hide file tree

Showing 2 changed files with 14 additions and 1 deletion.
diff --git a/FormalDerivative.py b/FormalDerivative.py
@@ -6,7 +6,7 @@ def polyrep(Fx):
   s = ["%dx ^ %d" % x for x in Fx]
   return " + ".join(s)
 
-fx = [(1,6),(1,0)]
+fx = [(-1,6),(1,0)]
 print(polyrep(fx))
 fx = FormalDerivative(fx)
 print(polyrep(fx))
diff --git a/lazy_erathostenes.py b/lazy_erathostenes.py
@@ -0,0 +1,13 @@
+from itertools import islice
+
+def nats(n):
+  yield n
+  yield from nats(n+1)
+
+def sieve(s):
+  #print("sieve",s)
+  n = next(s)
+  yield n
+  yield from sieve(i for i in s if i % n != 0)
+
+print(list(islice(sieve(nats(2)), 30)))