A213788 a(n) = Sum_{1<=iA000129(m).

0, 0, 0, 10, 214, 3491, 52001, 748788, 10636260, 150248190, 2117562834, 29816257390, 419662506490, 5905775317025, 83104503504515, 1169392060102440, 16454728773220584, 231536384221100316, 3257968708458764196, 45843125116860034258, 645061876629223784830, 9076710308820189950975

Offset: 0

Author

N. J. A. Sloane, Jun 20 2012

Keywords

• nonn
• easy

Links

• Alois P. Heinz, Table of n, a(n) for n = 0..275
• Index entries for linear recurrences with constant coefficients, signature (21,-98,-34,616,-532,-62,98,-7,-1).

Crossrefs

Cf. A000129, A213785.

Programs

Maple

a:= n-> (Matrix(9, (i, j)-> `if`(i=j-1, 1, `if`(i=9,
         [-1, -7, 98, -62, -532, 616, -34, -98, 21][j], 0)))^(n+3).
         <<5, -1, 0, 0, 0, 0, 10, 214, 3491>>)[1, 1]:
seq (a(n), n=0..30);  # Alois P. Heinz, Jun 20 2012

Mathematica

LinearRecurrence[{21, -98, -34, 616, -532, -62, 98, -7, -1}, {0, 0, 0, 10, 214, 3491, 52001, 748788, 10636260}, 30] (* Jean-François Alcover, Feb 17 2016 *)

PARI

concat(vector(3), Vec((x^4+2*x^3-23*x^2+4*x+10)*x^3 / ((x-1) * (x^2+14*x-1) * (x^2-2*x-1) * (x^2+2*x-1) * (x^2-6*x+1)) + O(x^30))) \\ Colin Barker, Feb 17 2016

Formula

G.f.: (x^4+2*x^3-23*x^2+4*x+10)*x^3 / ((x-1) * (x^2+14*x-1) * (x^2-2*x-1) * (x^2+2*x-1) * (x^2-6*x+1)). - Alois P. Heinz, Jun 20 2012