f r=print$sum[1+4*sum[div(r^2-z^2)i*(2-mod i 4)|i<-[1,3..6^5]]|z<-[-r..r]] main=mapM f[0..50]

Note that non-ascii characters in the above source code will be escaped (such as \x9f).

download

return to the top page