This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A372163 #13 Feb 16 2025 08:34:06
%S A372163 1,1,5,28,321,3636,65947,1154238,28622001,684987400,21513702771,
%T A372163 656705784714,24936869827465,928288327257084,41315505985090443,
%U A372163 1817727059210127286,92749773791662574433,4712674616532693996432,271157923143678988333027
%N A372163 E.g.f. A(x) satisfies A(x) = exp( x * (1+3*x)^(1/3) * A(x) ).
%H A372163 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/LambertW-Function.html">Lambert W-Function</a>.
%F A372163 E.g.f.: A(x) = exp( -LambertW(-x * (1+3*x)^(1/3)) ).
%F A372163 a(n) = n! * Sum_{k=0..n} (k+1)^(k-1) * 3^(n-k) * binomial(k/3,n-k)/k!.
%F A372163 a(n) ~ sqrt(4*r+1) * n^(n-1) / (exp(n - 5/2) * r^(n - 3/2)), where r = 0.29742497866288781360719311656731644994668261137281157848090655000... is the root of the equation r*(1 + 3*r)^(1/3) = exp(-1). - _Vaclav Kotesovec_, Apr 22 2024
%o A372163 (PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(exp(-lambertw(-x*(1+3*x)^(1/3)))))
%o A372163 (PARI) a(n) = n!*sum(k=0, n, (k+1)^(k-1)*3^(n-k)*binomial(k/3, n-k)/k!);
%Y A372163 Cf. A362771, A372162.
%K A372163 nonn
%O A372163 0,3
%A A372163 _Seiichi Manyama_, Apr 20 2024