|
%I #16 Mar 28 2020 05:28:16
%S 1,-1,2,-4,9,-22,54,-139,372,-948,2607,-7388,16058,-58957,174854,
%T 210448,4345025,-2008714,-165872030,-1756557123,-6144936528,
%U 60244093040,1164910003567,8228177887688,-10562519450714,-967088274083133,-11322641425582454,-37483806372774364
%N Second order complementary Bell numbers.
%t nmax = 27;
%t A = Exp[x] + O[x]^(nmax - 1);
%t B = Exp[1 - Integrate[A, x]]/E;
%t c = Exp[1 - Integrate[B, x]]/E;
%t CoefficientList[c, x] Range[0, nmax]! (* _Jean-François Alcover_, Jul 12 2019, from PARI *)
%o (Sage) # uses[bell_transform from A264428]
%o def A265023_list(len):
%o uno = [1]*len
%o complementary_bell_numbers = [sum((-1)^n*b for (n, b) in enumerate (bell_transform(n, uno))) for n in range(len)]
%o complementary_bell_numbers2 = [sum((-1)^n*b for (n, b) in enumerate (bell_transform(n, complementary_bell_numbers))) for n in range(len)]
%o return complementary_bell_numbers2
%o print(A265023_list(28))
%o (PARI)
%o \\ For n>28 precision has to be adapted as needed!
%o A = exp('x + O('x^33) );
%o B = exp(1 - intformal(A) )/exp(1);
%o C = exp(1 - intformal(B) )/exp(1);
%o round(Vec(serlaplace(C)))
%Y Cf. A000587 (complementary Bell numbers), A264428.
%K sign
%O 0,3
%A _Peter Luschny_, Dec 03 2015
| 