%I #26 Jan 12 2023 06:44:54
%S 0,1,5,26,162,1267,12343,145652,2036148,33192789,622384729,
%T 13263528350,318121600694,8517247764135,252725694989611,
%U 8258153081400856,295515712276222952,11523986940937975401,487562536078882116717,22291094729329088403298
%N Expansion of Sum_{k>=0} k^2 * x^k/(1 - k^2 * x).
%H Seiichi Manyama, <a href="/A349882/b349882.txt">Table of n, a(n) for n = 0..339</a>
%F a(n) = Sum_{k=0..n} k^(2*(n-k+1)).
%F a(n) = A234568(n+1) - 1. - _Hugo Pfoertner_, Dec 04 2021
%F a(n) ~ sqrt(Pi) * ((n+1)/LambertW(exp(1)*(n+1)))^(5/2 + 2*n - 2*(n+1)/LambertW(exp(1)*(n+1))) / sqrt(1 + LambertW(exp(1)*(n+1))). - _Vaclav Kotesovec_, Dec 04 2021
%F G.f.: Sum_{k>=1} x^k/(1 - (k+1)^2 * x). - _Seiichi Manyama_, Jan 12 2023
%t a[n_] := Sum[If[k == n - k + 1 == 0, 1, k^(2*(n - k + 1))], {k, 0, n}]; Array[a, 20, 0] (* _Amiram Eldar_, Dec 04 2021 *)
%o (PARI) a(n, s=2, t=2) = sum(k=0, n, k^(t*(n-k)+s));
%o (PARI) my(N=20, x='x+O('x^N)); concat(0, Vec(sum(k=0, N, k^2*x^k/(1-k^2*x))))
%o (PARI) my(N=20, x='x+O('x^N)); concat(0, Vec(sum(k=1, N, x^k/(1-(k+1)^2*x)))) \\ _Seiichi Manyama_, Jan 12 2023
%Y Cf. A003101, A234568, A249459, A349863, A359659.
%K nonn
%O 0,3
%A _Seiichi Manyama_, Dec 03 2021
|