OFFSET
0,3
MATHEMATICA
nmax = 27;
A = Exp[x] + O[x]^(nmax - 1);
B = Exp[1 - Integrate[A, x]]/E;
c = Exp[1 - Integrate[B, x]]/E;
CoefficientList[c, x] Range[0, nmax]! (* Jean-François Alcover, Jul 12 2019, from PARI *)
PROG
(Sage) # uses[bell_transform from A264428]
def A265023_list(len):
uno = [1]*len
complementary_bell_numbers = [sum((-1)^n*b for (n, b) in enumerate (bell_transform(n, uno))) for n in range(len)]
complementary_bell_numbers2 = [sum((-1)^n*b for (n, b) in enumerate (bell_transform(n, complementary_bell_numbers))) for n in range(len)]
return complementary_bell_numbers2
print(A265023_list(28))
(PARI)
\\ For n>28 precision has to be adapted as needed!
A = exp('x + O('x^33) );
B = exp(1 - intformal(A) )/exp(1);
C = exp(1 - intformal(B) )/exp(1);
round(Vec(serlaplace(C)))
CROSSREFS
KEYWORD
sign
AUTHOR
Peter Luschny, Dec 03 2015
STATUS
approved