N,D = map(int,input().split())
C = N * N * 2 + 1
dat = [0] * C
for x in range(1,N + 1):
    u = x * x
    for y in range(1,N + 1):
        dat[u + y * y] += 1
ans = 0
for w in range(1,N + 1):
    u = D + w * w
    for z in range(1, N + 1):
        v = u - z * z
        if 1 <= v < C:
            ans += dat[v]
print(ans)