-
Notifications
You must be signed in to change notification settings - Fork 6
/
Copy pathsieve.enso
70 lines (56 loc) · 2.02 KB
/
sieve.enso
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
import Standard.Base.Any.Any
import Standard.Base.Meta
import Standard.Base.Data.Time.Date_Time.Date_Time
import Standard.Base.Data.Numbers.Integer
from Standard.Base.Data.Boolean import True, False
from Standard.Base import IO
type Gen
Generator a:Integer tail:Any->Gen
next : Gen
next self = self.tail 0
n self = self.a
natural : Gen
natural =
gen = n -> Gen.Generator n (_ -> gen (n + 1))
gen 2
type Filter
Empty
Item n:Integer next:Filter
Head first:Filter set_last:Filter->Nothing
acceptAndAdd : Integer
acceptAndAdd self n =
upto = n.sqrt
iterate f = case f of
Filter.Item p next -> if n % p == 0 then False else
if p > upto then True else @Tail_Call iterate next
Filter.Empty -> True
appendN head set_last =
pair = Meta.atom_with_hole (h -> Filter.Item n h)
newLast = pair.value
set_last newLast
Filter.Head head.first pair.fill
case self of
Filter.Empty ->
pair = Meta.atom_with_hole (h -> Filter.Item n h)
Filter.Head pair.value pair.fill
Filter.Head _ _ -> if iterate self.first then appendN self self.set_last else self
type Primes
Alg generator filter
Primes.next self = case self of
Primes.Alg g f ->
filter = f.acceptAndAdd g.n
new_primes = Primes.Alg g.next filter
if Meta.is_same_object filter f then @Tail_Call new_primes.next else new_primes
Primes.lastPrime self = case self of
Primes.Alg g _ -> g.n - 1
compute_nth_prime : Primes -> Integer -> Any
compute_nth_prime p n = if n == 0 then p else
pn = p.next
@Tail_Call compute_nth_prime pn n-1
main =
at = Date_Time.now.to_unix_epoch_milliseconds
p = Primes.Alg natural Filter.Empty
r = compute_nth_prime p 100000
now = Date_Time.now.to_unix_epoch_milliseconds
IO.println <| "Hundred thousand prime numbers in " + (now - at).to_text + " ms. Last one is " + r.lastPrime.to_text
@Tail_Call main