a=[0]*1001
p=1
while True:
    mi = 1001
    for q in xrange(0, p):
        for r in xrange(0, p):
            if q + r > p: break
            for s in xrange(0, p):
                if r + s>p: break
                if s + q>p: break
                v = p*p - q*q - r*r - s*s
                if v>0 and v<=1000: a[v]+=1
                if mi>v: mi = v
    if mi>1000:break
    p+=1
for n in xrange(1, 1001): print [n, a[n]]