n = int(eval(input()))

a = [0] * (n + 1)

for x in range(1, 101):

    for y in range(1, 101):

        for z in range(1, 101):

            c = x * x + y * y + z * z + x * y + y * z + z * x

            if c <= n:

                a[c] +=1

            else:

                break

for i in a[1:]:

    print(i)

from sys import stdin

inp = lambda : stdin.readline().strip()



n = int(inp())

a = [0]*10**5

for i in range(1,101):

    for j in range(1,101):

        for k in range(1,101):

            a[i**2 + j**2 + k**2 + i*j + i*k + k*j] += 1

for i in a[1:n+1]:

    print(i) 
n = int(eval(input()))
a = [0] * (n + 1)
for x in range(1, 101):
    for y in range(1, 101):
        for z in range(1, 101):
            c = x * x + y * y + z * z + x * y + y * z + z * x
            if c <= n:
                a[c] += 1
            else:
                break
for i in a[1:]:
    print(i)

from sys import stdin

inp = lambda: stdin.readline().strip()
n = int(inp())
a = [0] * 10**5
for i in range(1, 101):
    for j in range(1, 101):
        for k in range(1, 101):
            a[i**2 + j**2 + k**2 + i * j + i * k + k * j] += 1
for i in a[1 : n + 1]:
    print(i)

-n = int(eval(input())) -a = [0] * (n + 1) -for x in range(1, 101): -    for y in range(1, 101): -        for z in range(1, 101): -            c = x * x + y * y + z * z + x * y + y * z + z * x -            if c <= n: -                a[c] += 1 -            else: -                break -for i in a[1:]: +from sys import stdin + +inp = lambda: stdin.readline().strip() +n = int(inp()) +a = [0] * 10**5 +for i in range(1, 101): +    for j in range(1, 101): +        for k in range(1, 101): +            a[i**2 + j**2 + k**2 + i * j + i * k + k * j] += 1 +for i in a[1 : n + 1]: