primes :: Integral a => [a]
primes = sieve [2..]
  where sieve (x:xs) = x : sieve [ y | y <- xs, rem y x /= 0 ]