OFFSET
0,3
LINKS
Muniru A Asiru, Table of n, a(n) for n = 0..550
FORMULA
E.g.f.: (x+x^2+x^2 exp(x))exp(exp(x)+x-1) which is x*A'(x) where A(x) is the e.g.f. for A070071.
MATHEMATICA
nn=20; a=Exp[Exp[x]-1]; Range[0, nn]!CoefficientList[Series[x D[x D[a, x], x], {x, 0, nn}], x]
Table[n^2 BellB[n], {n, 0, 30}] (* Harvey P. Dale, Jul 01 2022 *)
PROG
(GAP) List([0..22], n->n^2*Bell(n)); # Muniru A Asiru, Apr 20 2019
(Python)
from itertools import count, accumulate, islice
def A204618_gen(): # generator of terms
yield 0
blist, b = (1, ), 1
for n in count(1):
blist = list(accumulate(blist, initial=(b:=blist[-1])))
yield b*n**2
CROSSREFS
KEYWORD
nonn
AUTHOR
Geoffrey Critzer, Jan 17 2012
STATUS
approved