f=lambda r,a=4,m=1:r%4>2or r>=0and m*f(r+1)+f(r-a,a+8,2) for n in range(51):print(f(4*n*n))
Note that non-ascii characters in the above source code will be escaped (such as \x9f).
download
return to the top page