n, d = map(int, input().split())
c1, c2 = [0] * (2 * n ** 2 + 10 ** 6), [0] * (2 * n ** 2 + 10 ** 6)
for i in range(1, n + 1):
    for j in range(1, n + 1):
        c1[i ** 2 + j ** 2 - 1] += 1
for i in range(1, n + 1):
    for j in range(1, n + 1):
        if i ** 2 - j ** 2 + d - 1 >= 0:
            c2[i ** 2 - j ** 2 + d - 1] += 1
ans = 0
for i in range(2 * n ** 2 + 1):
    ans += c1[i] * c2[i]
print(ans)