{-# OPTIONS_GHC -O2 -fno-cse #-}
primesTME = 2 : gaps 3 (joinT [[p*p, p*p+2*p..] | p <- primes']) 
  where
    primes' = 3 : gaps 5 (joinT -- [[p*p, p*p+2*p..] | p <- primes'])
                  [map (p*) $ scanl (+) p (cycle[2]) | p <- primes'])  
 
joinT ((x:xs):t) = x : union xs (joinT (pairs t))  
pairs ((x:xs):ys:t) = (x : union xs ys) : pairs t 
gaps k s@(x:xs) | k<x  = k:gaps (k+2) s    -- equivalent to
                | True =   gaps (k+2) xs   --  [k,k+2..]`minus`s, k<=head s
 
main = mapM_ print 
         . take 10 . map fst 
         . tail . iterate (splitAt 10.snd)
         $ ([],take 100 primesTME)
 
union (x:xs) (y:ys) = case (compare x y) of 
           LT -> x : union  xs  (y:ys)
           EQ -> x : union  xs     ys 
           GT -> y : union (x:xs)  ys

ey0jIE9QVElPTlNfR0hDIC1PMiAtZm5vLWNzZSAjLX0KcHJpbWVzVE1FID0gMiA6IGdhcHMgMyAoam9pblQgW1twKnAsIHAqcCsyKnAuLl0gfCBwIDwtIHByaW1lcyddKSAKICB3aGVyZQogICAgcHJpbWVzJyA9IDMgOiBnYXBzIDUgKGpvaW5UIC0tIFtbcCpwLCBwKnArMipwLi5dIHwgcCA8LSBwcmltZXMnXSkKICAgICAgICAgICAgICAgICAgW21hcCAocCopICQgc2NhbmwgKCspIHAgKGN5Y2xlWzJdKSB8IHAgPC0gcHJpbWVzJ10pICAKIApqb2luVCAoKHg6eHMpOnQpID0geCA6IHVuaW9uIHhzIChqb2luVCAocGFpcnMgdCkpICAKcGFpcnMgKCh4OnhzKTp5czp0KSA9ICh4IDogdW5pb24geHMgeXMpIDogcGFpcnMgdCAKZ2FwcyBrIHNAKHg6eHMpIHwgazx4ICA9IGs6Z2FwcyAoaysyKSBzICAgIC0tIGVxdWl2YWxlbnQgdG8KICAgICAgICAgICAgICAgIHwgVHJ1ZSA9ICAgZ2FwcyAoaysyKSB4cyAgIC0tICBbayxrKzIuLl1gbWludXNgcywgazw9aGVhZCBzCgptYWluID0gbWFwTV8gcHJpbnQgCiAgICAgICAgIC4gdGFrZSAxMCAuIG1hcCBmc3QgCiAgICAgICAgIC4gdGFpbCAuIGl0ZXJhdGUgKHNwbGl0QXQgMTAuc25kKQogICAgICAgICAkIChbXSx0YWtlIDEwMCBwcmltZXNUTUUpCgp1bmlvbiAoeDp4cykgKHk6eXMpID0gY2FzZSAoY29tcGFyZSB4IHkpIG9mIAogICAgICAgICAgIExUIC0+IHggOiB1bmlvbiAgeHMgICh5OnlzKQogICAgICAgICAgIEVRIC0+IHggOiB1bmlvbiAgeHMgICAgIHlzIAogICAgICAgICAgIEdUIC0+IHkgOiB1bmlvbiAoeDp4cykgIHlz