from bisect import bisect_left, bisect_right

def main():
    n, d = map(int, input().split())

    L = []
    for w in range(1, n + 1):
        for z in range(1, n + 1):
            L.append(w**2 - z**2)
    L.sort()
    ans = 0
    for x in range(1, n + 1):
        for y in range(1, n + 1):
            res = x**2 + y**2 - d
            ans += bisect_right(L, res) - bisect_left(L, res)
    print(ans)

if __name__ == '__main__':
    main()