OFFSET
0,2
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..50
EXAMPLE
G.f.: A(x) = 1 + 2*x + 16*x^2 + 380*x^3 + 15280*x^4 + 842672*x^5 + 57985144*x^6 + 4735508672*x^7 + 445364211760*x^8 + 47281191656960*x^9 + 5586025249211056*x^10 + ...
MATHEMATICA
a[n_] := Module[{A}, A = {1}; Do[AppendTo[A, 0]; A[[-1]] = Round[ Coefficient[ Sum[(1+x + x*O[x]^Length[A])^(m*(m-1)/2)/(A.x^Range[0, Length[A] - 1])^m*2^m/3^(m + 1), {m, 0, 30 Length[A] + 200}]/2, x, Length[A] - 1]], {i, 1, n}]; Print[A[[n + 1]]]; A[[n + 1]]];
a /@ Range[0, 25] (* Jean-François Alcover, Oct 02 2019 *)
PROG
(PARI) /* Requires adequate precision */
\p500
{a(n) = my(A=[1]); for(i=1, n, A = concat(A, 0); A[#A] = round( polcoeff( sum(m=0, 30*#A+200, (1+x+x*O(x^#A))^(m*(m-1)/2)/Ser(A)^m*2^m/3^(m+1)*1.)/2, #A-1))); A[n+1]}
for(n=0, 25, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Apr 23 2019
EXTENSIONS
Added missing parentheses to definition. - N. J. A. Sloane, Aug 01 2019
STATUS
approved