main=mapM putStrLn[init$tail$show(n-1,n+1)|n<-[4..2000],all((1/=).mod(n^2))[2..n-2]]