OFFSET
0,3
FORMULA
G.f.: Sum_{k>=0} x^k / (1 - 3*k*x)^(k+1).
a(n) = Sum_{k=0..n} (3*k)^(n-k) * binomial(n,k).
MAPLE
A356827 := proc(n)
add((3*k)^(n-k) * binomial(n, k), k=0..n) ;
end proc:
seq(A356827(n), n=0..70) ; # R. J. Mathar, Dec 04 2023
PROG
(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(exp(x*exp(3*x))))
(PARI) my(N=20, x='x+O('x^N)); Vec(sum(k=0, N, x^k/(1-3*k*x)^(k+1)))
(PARI) a(n) = sum(k=0, n, (3*k)^(n-k)*binomial(n, k));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Aug 29 2022
STATUS
approved