primes = [ x | x <- [1..], length (factor x) == 1]
factor n = [x | x <- [2..n+1], n `mod` x == 0]