main=mapM print[n|n<-[1..2000],f n!!9<2] f n=n:f(sum[read[x]^2|x<-show n])
Note that non-ascii characters in the above source code will be escaped (such as \x9f).
download
return to the top page