N,D = map(int,input().split())
#LEFTはx**2 + y**2に相当
s = min(N,int(D**(0.5)))
t = (D**(0.5))//3
left = {}
right = {}
for i in range(t,s+1):
    for j in range(t,s+1):
         f = i**2 + j**2
         if(f not in left):
            left[f] = 1
         else:
            left[f] += 1

for i in range(t,s+1):
    for j in range(t,s+1):
        f = i**2 - j**2
        if(f not in right):
            right[f] = 1
        else:
            right[f] += 1

ans = 0
for x in left:
    if(D - x in right):
       ans += left[x]*right[D - x]
print(ans)