n,d = map(int, input().split())

leng=max(1+2*n*n, 1+2*n*n + d)

cnt1=[0]*(leng)
for x in range(1,n+1):
    for y in range(1,n+1):
        cnt1[x**2 + y**2]+=1

cnt2=[0]*(leng)
for w in range(1,n+1):
    for z in range(1,n+1):
        ind = w**2 - z**2 + d
        if ind >= 0:
            cnt2[ind]+=1

ans=0
for i in range(len(cnt1)):
    ans+=(cnt1[i]*cnt2[i])
print(ans)