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]