OFFSET
0,3
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Index entries for linear recurrences with constant coefficients, signature (17,-136,680,-2380,6188,-12376,19448,-24310,24310,-19448,12376,-6188,2380,-680,136,-17,1).
FORMULA
G.f.: (x + 33007*x^2 + 20971930*x^3 + 1786040870*x^4 + 42691533167*x^5 +
391057936193*x^6 + 1603741633836*x^7 + 3191398997460*x^8 + 3191399928015*x^9 + 1603741544321*x^10 + 391057582106*x^11 + 42691705382*x^12 + 1786044385*x^13 + 20960815*x^14 + 32512*x^15)/(1-x)^17. - G. C. Greubel, Dec 06 2017
MATHEMATICA
Table[n^9*(n^7+1)/2, {n, 0, 30}] (* G. C. Greubel, Dec 06 2017 *)
LinearRecurrence[{17, -136, 680, -2380, 6188, -12376, 19448, -24310, 24310, -19448, 12376, -6188, 2380, -680, 136, -17, 1}, {0, 1, 33024, 21533202, 2147614720, 76294921875, 1410559992576, 16616485461604, 140737555464192, 926510288136165, 5000000500000000, 22974866110759926, 92442132027408384, 332708309893839607, 1088976679235208960, 3284204197078125000, 9223372071214514176}, 20] (* Harvey P. Dale, Nov 21 2020 *)
PROG
(Magma) [n^9*(n^7+1)/2: n in [0..20]]; // Vincenzo Librandi, Aug 26 2011
(PARI) for(n=0, 30, print1(n^9*(n^7+1)/2, ", ")) \\ G. C. Greubel, Dec 06 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Dec 11 2009
STATUS
approved