turingprimes.hs

import System.Environment
import Control.Parallel.Strategies

import Math.NumberTheory.Primes.Sieve (primes)
import Math.NumberTheory.Primes.Testing (isPrime)

-- Converts 42d -> 101010d
binread :: Integer -> Integer
binread 0 = 0
binread n = let (q, r) = n `divMod` 2
            in  r + 10 * binread q

turingPrimes :: [Integer]
turingPrimes = (filter (isPrime.binread) primes)
                   `using` (parListChunk 4096 rdeepseq)

main = do
    args <- getArgs
    let bitCount = read $ head args
    let ps       = takeWhile (< (2^bitCount)) turingPrimes
    mapM_ print $ map (\ n -> (n, binread n)) ps
    print $ length ps