A000327 := proc(n) local a, x1, x2 ; a := 0 ; for x1 from 1 to floor(n^(3/2)) do x2 := (n-x1^(2/3))^(3/2) ; if floor(x2) >= x1+1 then a := a+floor(x2-x1) ; fi; od: a ; end: seq(A000327(n), n=3..80) ; # R. J. Mathar, Sep 29 2009
MATHEMATICA
A000327[n_] := Module[{a, x1, x2 }, a = 0; For[x1 = 1, x1 <= Floor[ n^(3/2)], x1++, x2 = (n - x1^(2/3))^(3/2); If[Floor[x2] >= x1+1, a = a + Floor[x2 - x1]]]; a ]; Table[A000327[n], {n, 3, 80}] (* Jean-François Alcover, Feb 07 2016, after R. J. Mathar *)