n, d = map(int, input().split())
ans = 0
sahen = [0] * 8000001
uhen = [0] * 5000001
for x in range(1, n+1):
	for y in range(1, n+1):
		sahen[x**2 + y**2] += 1
for z in range(1, n+1):
	for w in range(1, n+1):
		if w**2 + d - z**2 >= 0:
			uhen[w**2 + d - z**2] += 1
ans = 0
for i in range(5000001):
	ans += sahen[i] * uhen[i]
print(ans)