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

X = [0] * (8 * 10 ** 6 + 1)
Y = [0] * (8 * 10 ** 6 + 1)
ans = 0

for i in range(1, N + 1):
    
    for j in range(1, N + 1):
        
        X[i ** 2 + j ** 2] += 1
        if i ** 2 - j ** 2 + D > 0:
            Y[i ** 2 - j ** 2 + D] += 1
        
for i in range(8 * 10 ** 6 + 1):
    ans += X[i] * Y[i]
print(ans)