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