main=mapM print$filter((==1).f.f.f.f.f.f)[1..2000] f 0=0 f x=f(x`div`10)+(x`mod`10)^2