OFFSET
0,2
LINKS
Alois P. Heinz and Vaclav Kotesovec, Table of n, a(n) for n = 0..422 (terms n=51..98 from Vaclav Kotesovec)
FORMULA
From Vaclav Kotesovec, Nov 18 2016: (Start)
Recurrence: 2*(n+2)^3*(n+3)*(3*n + 4)*(3*n + 5)*(21546*n^6 - 45513*n^5 - 3699*n^4 + 13101*n^3 + 6745*n^2 + 2032*n + 228)*a(n) = (n+2)*(100641366*n^11 + 49968333*n^10 - 359008281*n^9 - 154345878*n^8 + 328506480*n^7 - 50516019*n^6 - 255412581*n^5 + 41482508*n^4 + 116397064*n^3 + 50875696*n^2 + 13946672*n + 1423680)*a(n-1) - 2*(173359116*n^12 + 510099768*n^11 - 127771911*n^10 - 1779706188*n^9 - 1320942528*n^8 + 1504634418*n^7 + 2094287811*n^6 + 74489810*n^5 - 783633768*n^4 - 274711936*n^3 + 13818000*n^2 + 12588128*n + 862080)*a(n-2) - 16*(2*n - 3)*(2*n + 1)*(3*n - 2)*(3*n - 1)*(4*n + 1)*(4*n + 3)*(21546*n^6 + 83763*n^5 + 91926*n^4 - 25905*n^3 - 108086*n^2 - 58260*n - 5560)*a(n-3).
a(n) ~ 2^(8*n+51/2) / (3^11 * Pi^(3/2) * n^(15/2)).
(End)
EXAMPLE
a(1) = 8: [(4),(0)], [(4),(1),(0)], [(4),(2),(0)], [(4),(2),(1),(0)], [(4),(3),(0)], [(4),(3),(1),(0)], [(4),(3),(2),(0)], [(4),(3),(2),(1),(0)].
MAPLE
b:= proc(l) option remember; `if`(l[-1]=0, 1, add(add(b(subsop(
i=j, l)), j=`if`(i=1, 0, l[i-1])..l[i]-1), i=1..nops(l)))
end:
a:= n-> `if`(n=0, 1, b([4$n])):
seq(a(n), n=0..16);
MATHEMATICA
b[l_] := b[l] = If[l[[-1]] == 0, 1, Sum[Sum[b[ReplacePart[l, i -> j]], {j, If[i == 1, 0, l[[i - 1]]], l[[i]] - 1}], {i, 1, Length[l]}]];
a[n_] := If[n == 0, 1, b[Array[4&, n]]];
a /@ Range[0, 16] (* Jean-François Alcover, Dec 20 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Jul 16 2013
STATUS
approved