r=lambda(x):range(-x,x+1) for n in range(51):print sum(x*x+y*y+z*z<=n*n for x in r(n)for y in r(int((n*n-x*x)**.5))for z in r(int((n*n-x*x-y*y)**.5)))