import sys

input = sys.stdin.readline

import math

def main():

    N = int(eval(input()))

    n = int(math.sqrt(N)) + 1

    ans = [0]*10001

    for x in range(1,n):

        for y in range(1,n):

            for z in range(1,n):

                v = x**2 + y**2 + z**2 + x*y + y*z + z*x

                if v <= 10000:

                    ans[v] += 1 

    for i in range(1,N+1):

        print((ans[i]))



if __name__ == '__main__':

    main()
import sys

input = sys.stdin.readline

from math import sqrt

def main():

    N = int(eval(input()))

    A = [0]*(N+1)

    n = int(sqrt(N)) + 1

    for x in range(1, n):

        for y in range(1, n):

            for z in range(1, n):

                tmp = x*x + y*y + z*z + x*y + y*z + z*x

                if tmp <= N:

                    A[tmp] += 1

    for i in range(1,N+1):

        print((A[i]))



if __name__ == '__main__':

    main()
import sys

input = sys.stdin.readline
import math


def main():
    N = int(eval(input()))
    n = int(math.sqrt(N)) + 1
    ans = [0] * 10001
    for x in range(1, n):
        for y in range(1, n):
            for z in range(1, n):
                v = x**2 + y**2 + z**2 + x * y + y * z + z * x
                if v <= 10000:
                    ans[v] += 1
    for i in range(1, N + 1):
        print((ans[i]))


if __name__ == "__main__":
    main()

import sys

input = sys.stdin.readline
from math import sqrt


def main():
    N = int(eval(input()))
    A = [0] * (N + 1)
    n = int(sqrt(N)) + 1
    for x in range(1, n):
        for y in range(1, n):
            for z in range(1, n):
                tmp = x * x + y * y + z * z + x * y + y * z + z * x
                if tmp <= N:
                    A[tmp] += 1
    for i in range(1, N + 1):
        print((A[i]))


if __name__ == "__main__":
    main()

-import math +from math import sqrt -    n = int(math.sqrt(N)) + 1 -    ans = [0] * 10001 +    A = [0] * (N + 1) +    n = int(sqrt(N)) + 1 -                v = x**2 + y**2 + z**2 + x * y + y * z + z * x -                if v <= 10000: -                    ans[v] += 1 +                tmp = x * x + y * y + z * z + x * y + y * z + z * x +                if tmp <= N: +                    A[tmp] += 1 -        print((ans[i])) +        print((A[i]))