N, D = map(int, input().split())
# x >= y >= z とする
cnt = 0
ans = 0
for w in range(1, N + 1):
    wd = w * w + D
    for x in range(max(1, int((wd / 3) ** 0.5)), N + 1):
        x2 = x * x
        if x2 > wd:
            break
        for y in range(max(1, int(((wd - x2) / 2) ** 0.5)), x + 1):
            y2 = y * y
            if x2 + y2 > wd:
                break
            wdxy = wd - x2 - y2
            tz = int(wdxy ** 0.5)
            for i in range(max(1, tz - 1), tz + 2):
                if i * i == wdxy:
                    if i > N or i > y:
                        continue
                    z = i
                    # print(x, y, z, w, x*x+y*y+z*z,w*w +D)
                    if x == y and y == z:
                        ans += 1
                    elif x == y or y == z:
                        ans += 3
                    else:
                        ans += 6
                    break

print(ans)