powered by
Uses an obvious sieve method and some care, working with logical and integers to be quite fast.
Primes(n)
max(n)
n
Factorize
(p1 <- Primes(100)) system.time(p1k <- Primes(1000)) # still lightning .. stopifnot(length(p1k) == 168)
Run the code above in your browser using DataLab