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

X = []
for i in range(N + 1):
    X.append(i * i)

L = 2 * X[N] + 1
Y = [0] * (L)


for i in range(1, N + 1):
    for j in range(1, N + 1):
        Y[X[i] + X[j]] += 1
ans = 0
for i in range(1, N + 1):
    for j in range(1, N + 1):
        t = D + X[i] - X[j]
        if 0 <= t < L:
            ans += Y[t]

print(ans)