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

0, 1, 640, 32805, 557056, 5078125, 31072896, 144120025, 545259520, 1764915561, 5050000000, 13075891741, 31173672960, 69337111285, 145365222016, 289606640625, 551903297536, 1011484828945, 1790743593600, 3074024910421, 5132800000000, 8358851918781

Offset

0,3

Links

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

Index entries for linear recurrences with constant coefficients, signature (11,-55,165,-330,462,-462,330,-165,55,-11,1).

Formula

From Colin Barker, Dec 05 2017: (Start)

G.f.: x*(1 + x)*(1 + 628*x + 25192*x^2 + 206044*x^3 + 443470*x^4 + 206044*x^5 + 25192*x^6 + 628*x^7 + x^8) / (1 - x)^11.

a(n) = 11*a(n-1) - 55*a(n-2) + 165*a(n-3) - 330*a(n-4) + 462*a(n-5) - 462*a(n-6) + 330*a(n-7) - 165*a(n-8) + 55*a(n-9) - 11*a(n-10) + a(n-11) for n>10.

(End)

Mathematica

Table[n^8*(n^2+1)/2, {n, 0, 50}] (* G. C. Greubel, Dec 05 2017 *)

Prog

(Magma) [n^8*(n^2+1)/2: n in [0..30]]; // Vincenzo Librandi, Aug 26 2011
(PARI) for(n=0, 30, print1(n^8*(n^2+1)/2, ", ")) \\ G. C. Greubel, Dec 05 2017
(PARI) concat(0, Vec(x*(1 + x)*(1 + 628*x + 25192*x^2 + 206044*x^3 + 443470*x^4 + 206044*x^5 + 25192*x^6 + 628*x^7 + x^8) / (1 - x)^11 + O(x^40))) \\ Colin Barker, Dec 05 2017

Crossrefs

Sequence in context: A097105 A233911 A178975 * A290029 A268875 A252426

Adjacent sequences: A170771 A170772 A170773 * A170775 A170776 A170777

Keyword

nonn,easy

Author

N. J. A. Sloane, Dec 11 2009

Status

approved