N, D = map(int, input().split())
cnt1 = [0] * (2*N*N+1)
cnt2 = [0] * (2*N*N+1)

for x in range(1, N+1):
    for y in range(1, N+1):
        cnt1[x**2 + y**2] += 1

for w in range(1, N+1):
    for z in range(1, N+1):
        if 0 <= w**2 - z**2 + D <= 2*N*N:
            cnt2[w**2 - z**2 + D] += 1

ans = sum(cnt1_i * cnt2_i for cnt1_i, cnt2_i in zip(cnt1, cnt2))
print(ans)