OFFSET
1,2
COMMENTS
The old name was: a(n)+a(n+1) = cube: 0+8=8=2^3; 8+19=27=3^3; 19+45=64=4^3; ...
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
FORMULA
From R. J. Mathar, Dec 19 2008: (Start)
G.f.: x^2*(8-5*x+4*x^2-x^3)/((1+x)*(1-x)^4).
a(n)+a(n+1) = A000578(n+1).
a(n)= -1/8+3*n^2/4+n^3/2+9*(-1)^n/8. (End)
MATHEMATICA
a=1; lst={}; Do[a=n^3-a; AppendTo[lst, a], {n, 1, 5!}]; lst
CoefficientList[Series[x (8 - 5 x + 4 x^2 - x^3)/((1 + x) (1 - x)^4), {x, 0, 50}], x] (* Vincenzo Librandi, May 19 2014 *)
PROG
(Magma) [n le 1 select n-1 else n^3-Self(n-1): n in [1..50]]; // Vincenzo Librandi, May 19 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vladimir Joseph Stephan Orlovsky, Dec 17 2008
EXTENSIONS
Name improved by Alex Ratushnyak, Aug 05 2012
STATUS
approved