Sexy primes: Difference between revisions
Content deleted Content added
Line 695:
printf "Number of unsexy primes: %d\n Last 10: %s\n\n" (length unsexy) (show $ drop (length unsexy - 10) unsexy)
where d = takeWhile (< 1000035) primes
unsexy = filter (\n -> (not $ isPrime (n+6)) && (not $ isPrime (n-6))) d
pairs = foldr (\a b -> if isPrime (a-6) then (a-6,a) : b else b) [] d
triplets = foldr (\a b -> if isPrime (a-12) && isPrime (a-6) then (a-12,a-6,a) : b else b) [] d
| |||