N, D = map(int, input().split())
r1 = [0] * 8000001
r2 = [0] * 8000001
ans = 0
for x in range(1, N + 1):
    for y in range(1, N + 1):
        r1[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 > 1:
            r2[w**2 - z**2 + D] += 1
for i in range(8000001):
    ans += r1[i] * r2[i]
print(ans)