main=mapM putStrLn[init$tail$show(n-1,n+1)|n<-[4..2000],all((1/=).mod(n^2))[2..n-2]]
Note that non-ascii characters in the above source code will be escaped (such as \x9f).
download
return to the top page