N,D=map(int,input().split())
c=[0 for i in range(2*(4000000)+1)]
for i in range(1,1+N):
    for j in range(1,1+N):
        a=i**2+j**2
        c[a]+=1
m=dict()
for i in range(1,1+N):
    for j in range(1,1+N):
        a=i**2-j**2
        if a in m:
            m[a]+=1
        else:
            m[a]=1
ans=0
for i in m:
    k=D-i
    if k>=0:
      ans+=m[i]*c[k]
print(ans)