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

su = [0] * (8*10**6 + 1)
for a in range(1, N+1):
    for b in range(1, N+1):
        su[a*a + b*b] += 1


di = [0] * (8*10**6 + 1)
for c in range(1, N+1):
    for d in range(1, N+1):
        k = c*c - d*d + D
        if k > 0:
            di[c*c - d*d + D] += 1

print(sum([i * j for i, j in zip(su, di)]))