OFFSET
0,4
MAPLE
CL := (f, x) -> PolynomialTools:-CoefficientList(f, x):
A327002 := proc(n) local P, Q;
P := proc(n) option remember; if n = 0 then return 1 fi;
add(binomial(3*n, 3*k+3)*P(n-k-1)*x, k=0..n-1) end:
Q := proc(n) option remember; if n = 0 then return 1 fi;
add(binomial(3*n-1, 3*k)*Q(k)*Q(n-1-k), k=0..n-1) end:
Q(n) - add(CL(P(n), x)[k+1]/k!, k=0..n) end:
seq(A327002(n), n=0..15);
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Aug 13 2019
STATUS
approved