|
%I #18 Aug 20 2021 06:15:04
%S 1,1,1,3,9,21,161,911,3473,48329,406241,2150171,44216921,491897693,
%T 3327845249,90934644711,1257256962081,10352273016081,353351881109313,
%U 5836715156967219,56621346170765481,2319460179075419941,44545835926727113441,497433851743810193983,23782590451590763744689
%N E.g.f.: exp(x-(1-x^3)^(1/3)+1).
%F a(n) ~ n! * (exp(2) + 2*exp(1/2) * cos((4*Pi*n - 3*sqrt(3))/6)) / (3^(2/3)*Gamma(2/3)*n^(4/3)) * (1 - 3^(5/6)*Gamma(2/3)^2 / (2*Pi*n^(1/3))). - _Vaclav Kotesovec_, Jun 08 2021
%t nmax = 30; CoefficientList[Series[E^(x-(1-x^3)^(1/3)+1), {x, 0, nmax}], x] * Range[0, nmax]! (* _Vaclav Kotesovec_, Jun 08 2021 *)
%Y Cf. A259239,A067622.
%K nonn
%O 0,4
%A _Karol A. Penson_ and Katarzyna Gorska, Jun 25 2015
|
