N, D = map(int, input().split())

A = [0 for _ in range(5000000)]
B = [0 for _ in range(5000000)]

for x in range(1, N + 1):
    for y in range(1, N + 1):
        s = x**2 + y**2
        if 0 <= s < 5000000:
            A[s] += 1

for z in range(1, N + 1):
    for w in range(1, N + 1):
        t = D + w**2 - z**2
        if 0 <= t < 5000000:
            B[t] += 1

res = 0

for i in range(5000000):
    res += A[i] * B[i]

print(res)