N,D=map(int,input().split())

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

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