[Haskell] これで素数が出るなんて！

primes = map head (iterate sieve [2..])
sieve (p:xs) = [x| x<-xs, x `mod` p /=0]

・・・たったこれだけで素数がずらずら出てくる。
すごい。ミラクル。かっこいい。
クイックソートの例は見たことあったけど、これは初めて。
まだ理解できていないけど、凄いパワーだHaskell