A168665 a(n) = n^7*(n^8 + 1)/2.

0, 1, 16448, 7175547, 536879104, 15258828125, 235092632256, 2373781166743, 17592187092992, 102945568438809, 500000005000000, 2088624094451411, 7703510805209088, 25592946538419637, 77784047831612864

Offset

0,3

Links

Vincenzo Librandi, Table of n, a(n) for n = 0..10000

Formula

From G. C. Greubel, Jul 29 2016: (Start)

G.f.: x*(1 + 16432*x + 6912499*x^2 + 424043552*x^3 + 7520619041*x^4 + 51388499408*x^5 + 155693786091*x^6 + 223769429952*x^7 + 155693786091*x^8 + 51388499408*x^9 + 7520619041*x^10 + 424043552*x^11 + 6912499*x^12 + 16432*x^13 + x^14)/(1 - x)^16.

E.g.f.: (1/2)*x*(2 + 16446*x + 2375402*x^2 + 42356300*x^3 + 210767060*x^4 + 420693294*x^5 + 408741334*x^6 + 216627840*x^7 + 67128490*x^8 + 12662650*x^9 + 1479478*x^10 + 106470*x^11 + 4550*x^12 + 105*x^13 + x^14)*exp(x). (End)

Mathematica

Table[n^7*(n^8 + 1)/2, {n, 0, 50}] (* G. C. Greubel, Jul 29 2016 *)

Prog

(Magma) [n^7*(n^8+1)/2: n in [0..20]]; // Vincenzo Librandi, Aug 28 2011

Crossrefs

Sequence in context: A133533 A250844 A076166 * A283027 A031829 A057680

Adjacent sequences: A168662 A168663 A168664 * A168666 A168667 A168668

Keyword

nonn,easy

Author

N. J. A. Sloane, Dec 11 2009

Status

approved