A351795 a(n) = n! * Sum_{k=0..n} (k * (n-k))^k/k!.
1, 1, 4, 30, 396, 8360, 256470, 10619952, 564959528, 37370475648, 3001942868490, 287388158562560, 32278318416029532, 4197544986996581376, 625014083479647028622, 105554855135062180485120, 20053957030647088382195280, 4255329207190209023134564352
Offset: 0
Keywords
Links
- Vaclav Kotesovec, Table of n, a(n) for n = 0..266
- Vaclav Kotesovec, Graph - the asymptotic ratio (10000 terms)
Programs
-
Mathematica
a[n_] := n!*(1 + Sum[(k*(n - k))^k/k!, {k, 1, n}]); Array[a, 18, 0] (* Amiram Eldar, Feb 19 2022 *) -
PARI
a(n) = n!*sum(k=0, n, (k*(n-k))^k/k!);
Formula
a(n) ~ sqrt(2*Pi) * n^(2*n + 1/2) / (sqrt(LambertW(exp(2)*n)^2 - 1) * exp(n*(1 - 1/LambertW(exp(2)*n))) * LambertW(exp(2)*n)^n). - Vaclav Kotesovec, Feb 20 2022