N,D=map(int,input().split())
v=[n*n for n in range(1,N+1)]
p=[0]*8000001
for x in v:
    for y in v:
        p[x+y]+=1
sm=0
for z in v:
    for w in v:
        n=w-z+D
        if n>0: sm+=p[n]
print(sm)