def MI():
    return map(int, input().split())

n,d=MI()
upper=2*(n**2)
xydic=[0]*(upper+1)
wzdic=[0]*(upper+1)
for i in range(1,n+1):
    for j in range(1,n+1):
        r=i**2-j**2+d
        l=i**2+j**2
        if 0<r<=upper:
            wzdic[r]+=1
        xydic[l]+=1

ans=0
for i in range(1,upper+1):
    ans+=wzdic[i]*xydic[i]
print(ans)