OFFSET
5,1
COMMENTS
a(n) counts the solutions to the inequality x_1^(1/2)+x_2^(1/2)+x_3^(1/2)<=n for any three distinct integers 1<=x_1<x_2<x_3. - R. J. Mathar, Jul 03 2009
REFERENCES
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
B. K. Agarwala, F. C. Auluck, Statistical mechanics and partitions into non-integral powers of integers, Proc. Camb. Phil. Soc., 47 (1951), 207-216.
B. K. Agarwala and F. C. Auluck, Statistical mechanics and partitions into non-integral powers of integers, Proc. Camb. Phil. Soc., 47 (1951), 207-216. [Annotated scanned copy]
MAPLE
A000397 := proc(n) local a, x1, x2, x3 ; a := 0 ; for x1 from 1 to n^2 do for x2 from x1+1 to floor( (n-x1^(1/2))^2 ) do x3 := (n-x1^(1/2)-x2^(1/2))^2 ; if floor(x3) >= x2+1 then a := a+floor(x3-x2) ; fi; od: od: a ; end: for n from 5 do printf("%d, \n", A000397(n)) ; od: # R. J. Mathar, Sep 29 2009
MATHEMATICA
A000397[n_] := Module[{a, x1, x2, x3}, a = 0; For[x1 = 1, x1 <= n^2, x1++, For[x2 = x1+1, x2 <= Floor[(n-x1^(1/2))^2], x2++, x3 = (n-x1^(1/2) - x2^(1/2))^2 ; If[Floor[x3] >= x2+1, a = a + Floor[x3-x2]]]]; a]; Reap[ For[n = 5, n <= 40, n++, Print[an = A000397[n]; Sow[an]]]][[2, 1]] (* Jean-François Alcover, Feb 08 2016, after R. J. Mathar *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from R. J. Mathar, Sep 29 2009
More terms from Sean A. Irvine, Nov 14 2010
STATUS
approved