OFFSET
0,3
LINKS
Stanislav Sykora, Table of n, a(n) for n = 0..9999
Stanislav Sýkora, Magnetic Resonance on OEIS, Stan's NMR Blog (Dec 31, 2014), Retrieved Nov 12, 2019.
Index entries for linear recurrences with constant coefficients, signature (8,-27,48,-42,0,42,-48,27,-8,1).
FORMULA
a(n) = Sum_{k=0..floor(n/2)} (n-2k)^7.
From Colin Barker, Dec 22 2015: (Start)
a(n) = 1/96*(6*n^8+48*n^7+112*n^6-224*n^4+256*n^2+51*((-1)^n-1)).
G.f.: x*(1+120*x+1191*x^2+2416*x^3+1191*x^4+120*x^5+x^6) / ((1-x)^9*(1+x)).
(End)
EXAMPLE
a(5) = 5^7 + 3^7 + 1^7 = 80313.
MATHEMATICA
Table[SeriesCoefficient[x (1 + 120 x + 1191 x^2 + 2416 x^3 + 1191 x^4 + 120 x^5 + x^6)/((1 - x)^9 (1 + x)), {x, 0, n}], {n, 0, 26}] (* Michael De Vlieger, Dec 22 2015 *)
LinearRecurrence[{8, -27, 48, -42, 0, 42, -48, 27, -8, 1}, {0, 1, 128, 2188, 16512, 80313, 296448, 903856, 2393600, 5686825}, 30] (* Vincenzo Librandi, Dec 23 2015 *)
nxt[{n_, a_, b_}]:={n+1, b, a+(n+1)^7}; NestList[nxt, {1, 0, 1}, 30][[All, 2]] (* Harvey P. Dale, Jun 16 2022 *)
PROG
(PARI) nmax=100; a = vector(nmax); a[2]=1; for(i=3, #a, a[i]=a[i-2]+(i-1)^7); print(a);
(PARI) concat(0, Vec(x*(1+120*x+1191*x^2+2416*x^3+1191*x^4+120*x^5+x^6)/((1-x)^9*(1+x)) + O(x^50))) \\ Colin Barker, Dec 22 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Stanislav Sykora, Nov 07 2013
STATUS
approved