|
%I #27 Mar 08 2021 04:51:57
%S 1,0,0,0,0,0,1,1,1,1,1,1,463,1717,4720,11441,25884,56135,2977313,
%T 23524737,125212889,552517341,2183244857,8025931950,124257251233,
%U 1468856487536,12433365625566,85767520652726,518324768774506,2858925345803536,26181976719735061
%N E.g.f.: exp(exp(x) - Sum_{i=0..5} x^i/i!).
%C a(n) is the number of set partitions of [n] into blocks of size > 5.
%H Seiichi Manyama, <a href="/A293025/b293025.txt">Table of n, a(n) for n = 0..595</a>
%F E.g.f.: Product_{i>5} exp(x^i/i!).
%p a:= proc(n) option remember; `if`(n=0, 1, add(
%p a(n-j)*binomial(n-1, j-1), j=6..n))
%p end:
%p seq(a(n), n=0..35); # _Alois P. Heinz_, Sep 28 2017
%t m = 31;
%t Exp[Exp[x] - Sum[x^i/i!, {i, 0, 5}]] + O[x]^m // CoefficientList[#, x]& // (# Range[0, m-1]!)& (* _Jean-François Alcover_, Mar 08 2021 *)
%o (PARI) my(x='x+O('x^66)); Vec(serlaplace(exp(exp(x)-1-x-x^2/2-x^3/6-x^4/24-x^5/120)))
%Y Column k=5 of A293024.
%K nonn
%O 0,13
%A _Seiichi Manyama_, Sep 28 2017
|
