Svoboda | Graniru | BBC Russia | Golosameriki | Facebook
login
A352981
a(n) = Sum_{k=0..floor(n/2)} k^n.
6
1, 0, 1, 1, 17, 33, 794, 2316, 72354, 282340, 10874275, 53201625, 2438235715, 14350108521, 762963987380, 5249352196144, 317685943157892, 2502137235710736, 169842891165484965, 1506994510201252425, 113394131858832552133, 1119223325228757961465
OFFSET
0,5
LINKS
FORMULA
G.f.: Sum_{k>=0} (k * x)^(2 * k) / (1 - k * x).
a(n) ~ exp((3 + (-1)^n)/2) * (n/2)^n / (exp(2) - 1). - Vaclav Kotesovec, Apr 14 2022
MATHEMATICA
a[0] = 1; a[n_] := Sum[k^n, {k, 0, Floor[n/2]}]; Array[a, 22, 0] (* Amiram Eldar, Apr 13 2022 *)
PROG
(PARI) a(n) = sum(k=0, n\2, k^n);
(PARI) my(N=40, x='x+O('x^N)); Vec(sum(k=0, N, (k*x)^(2*k)/(1-k*x)))
(Magma) [(&+[k^n: k in [0..Floor(n/2)]]): n in [0..40]]; // G. C. Greubel, Nov 01 2022
(SageMath) [sum( k^n for k in range((n//2)+1)) for n in range(41)] # G. C. Greubel, Nov 01 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Seiichi Manyama, Apr 13 2022
STATUS
approved