cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A360727 Expansion of Sum_{k>=0} (k * x * (1 + x^2))^k.

This page as a plain text file.
%I A360727 #17 Feb 18 2023 19:00:40
%S A360727 1,1,4,28,264,3206,47684,839249,17058688,393216567,10134918592,
%T A360727 288815780665,9016571143680,306027510946208,11219450971161024,
%U A360727 441846991480590475,18602901833071633792,833832341625621777368,39642569136740054367808
%N A360727 Expansion of Sum_{k>=0} (k * x * (1 + x^2))^k.
%H A360727 Winston de Greef, <a href="/A360727/b360727.txt">Table of n, a(n) for n = 0..385</a>
%F A360727 a(n) = Sum_{k=0..floor(n/3)} (n-2*k)^(n-2*k) * binomial(n-2*k,k).
%F A360727 a(n) ~ n^n * (1 + exp(-2)/n + exp(-4)/(2*n^2)). - _Vaclav Kotesovec_, Feb 18 2023
%t A360727 nmax = 20; CoefficientList[1 + Series[Sum[(k*x*(1 + x^2))^k, {k, 1, nmax}], {x, 0, nmax}], x] (* _Vaclav Kotesovec_, Feb 18 2023 *)
%o A360727 (PARI) my(N=20, x='x+O('x^N)); Vec(sum(k=0, N, (k*x*(1+x^2))^k))
%o A360727 (PARI) a(n) = sum(k=0, n\3, (n-2*k)^(n-2*k)*binomial(n-2*k, k));
%Y A360727 Cf. A360611, A360728.
%Y A360727 Cf. A360730.
%K A360727 nonn
%O A360727 0,3
%A A360727 _Seiichi Manyama_, Feb 18 2023