login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A350366 a(n) = [x^n] Product_{k=1..n} (1 + k*x)/(1 - k*x). 4
1, 2, 18, 312, 8000, 271770, 11502162, 583036832, 34437042432, 2322677883330, 176137593178250, 14835018315726312, 1373972097646792800, 138787120025382437882, 15184417945878202716450, 1788809909368939651651200, 225755544056485027686459392 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
LINKS
FORMULA
a(n) = Sum_{k=0..n} (-1)^(n-k) * Stirling1(n+1, k+1) * Stirling2(k+n, n).
a(n) ~ c * d^n * (n-1)!, where d = (1+r) / ((-1 + exp(r + LambertW(-1, -exp(-r)*r))) * LambertW(-exp(-1-r)*(1+r))) = 8.406107401279769476199925123910168..., r = 0.7545302104650497245839827141610818561001159135034... is the root of the equation r*(1 + r + LambertW(-exp(-1 - r)*(1 + r))) = -(1 + r)*(r + LambertW(-1, -exp(-r)*r)) and c = 0.281498742412700978029375818376931142913157133987685... - Vaclav Kotesovec, Dec 29 2021
MATHEMATICA
a[n_] := Coefficient[Series[Product[(1 + k*x)/(1 - k*x), {k, 1, n}], {x, 0, n}], x, n]; Array[a, 17, 0] (* Amiram Eldar, Dec 27 2021 *)
Table[Sum[(-1)^(n - k)*StirlingS1[n + 1, k + 1] * StirlingS2[k + n, n], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Dec 29 2021 *)
PROG
(PARI) a(n) = sum(k=0, n, (-1)^(n-k)*stirling(n+1, k+1, 1)*stirling(k+n, n, 2));
CROSSREFS
Sequence in context: A258922 A192555 A370941 * A370056 A179497 A296837
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Dec 27 2021
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 09:04 EDT 2024. Contains 371240 sequences. (Running on oeis4.)