OFFSET
0,3
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..760
FORMULA
a(n) = Sum_{k=0..floor(n/2)} binomial(n-k,k) * (n-k)^k.
a(n) ~ exp(exp(1/2)*sqrt(n/2) - 3*exp(1)/8) * n^(n/2) / 2^(n/2 + 1) * (1 + ((exp(1/2) + exp(-1/2))/2^(5/2) + 11*exp(3/2)/2^(9/2))/sqrt(n)).
MAPLE
N:= 40:
S:= series(add((1+k*x)^k*x^k, k=0..N), x, N+1):
seq(coeff(S, x, k), k=0..N); # Robert Israel, Feb 13 2023
MATHEMATICA
nmax = 30; CoefficientList[Series[Sum[(1 + k*x)^k * x^k, {k, 0, nmax}], {x, 0, nmax}], x]
Flatten[{1, Table[Sum[Binomial[n-k, k] * (n-k)^k, {k, 0, n/2}], {n, 1, 30}]}]
PROG
(PARI) {a(n) = polcoeff(sum(m=0, n, (1 + m*x)^m * x^m + x*O(x^n)), n)};
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Feb 13 2023
STATUS
approved