n,d = map(int,input().split())

yz = [0]*(2*(n**2)+1)

for i in range(1,n+1):
    for j in range(1,n+1):
        yz[i**2+j**2] += 1

ans = 0

for w in range(1,n+1):
    for x in range(1,n+1):
        tmp = w**2+d-x**2
        if tmp > 0 and tmp <= (2*(n**2)):
            ans += yz[tmp]

print(ans)