User:WillNess: Difference between revisions

From HaskellWiki
No edit summary
(var rename for uniformity)
Line 22: Line 22:
The math formula for Sieve of Eratosthenes,
The math formula for Sieve of Eratosthenes,


::::<math>\textstyle\mathbb{S} = \mathbb{N}_{2} \setminus \bigcup_{p\in \mathbb{S}} \{n p:n \in \mathbb{N}_{p}\}</math>  
::::<math>\textstyle\mathbb{S} = \mathbb{N}_{2} \setminus \bigcup_{p\in \mathbb{S}} \{p\,q:q \in \mathbb{N}_{p}\}</math>  
where  
where  
::::<math>\textstyle\mathbb{N}_{k} = \{ n \in \mathbb{N} : n \geq k \}</math> &emsp; . . . or, &ensp;<math>\textstyle\mathbb{N}_{k} = \{k\} \bigcup \mathbb{N}_{k+1}</math> &emsp; :)&emsp;:) .
::::<math>\textstyle\mathbb{N}_{k} = \{ n \in \mathbb{N} : n \geq k \}</math> &emsp; . . . or, &ensp;<math>\textstyle\mathbb{N}_{k} = \{k\} \bigcup \mathbb{N}_{k+1}</math> &emsp; :)&emsp;:) .

Revision as of 08:44, 12 September 2011

I'm interested in Haskell.

I like this:

--   inifinte folding idea due to Richard Bird
--   double staged production idea due to Melissa O'Neill
--   tree folding idea Dave Bayer / simplified formulation Will Ness
primes = 2 : g (fix g) 
  where                
    g xs = 3 : gaps 5 (foldi (\(c:cs) -> (c:) . union cs) 
                             [[x*x, x*x+2*x..] | x <- xs])
    gaps k s@(c:t)                                        
       | k < c = k : gaps (k+2) s     -- minus [k,k+2..] (c:t), k<=c
       | True  =     gaps (k+2) t     --   fused to avoid a space leak

fix g = xs where xs = g xs            -- global defn to avoid space leak

foldi is on Tree-like folds page. union and more at Prime numbers.

The math formula for Sieve of Eratosthenes,

S=N2pS{pq:qNp}

where

Nk={nN:nk}   . . . or,  Nk={k}Nk+1   :) :) .

Trial division sieve:

T={nN2:(pT)(2pn¬(pn))}

If you're put off by self-referentiality, just replace S or T on the right-hand side of equations with N2.