INF = 8000001
left = [0]*INF
right = [0]*INF
N, D = map(int, input().split())

for x in range(1, N+1):
    for y in range(1, N+1):
        left[x**2 + y**2] += 1
for z in range(1, N+1):
    for w in range(1, N+1):
        if w**2 - z**2 + D > 0:
            right[w**2 - z**2 + D] += 1
ans = 0;
for i in range(INF):
    ans += left[i] * right[i]
print(ans)