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.

A351133 a(n) = Sum_{k=0..n} k! * k^(2*n) * Stirling1(n,k).

This page as a plain text file.
%I A351133 #27 Feb 06 2022 02:15:25
%S A351133 1,1,31,3992,1342294,932514674,1161340476698,2356863300156504,
%T A351133 7278091701243797640,32477694155566998880608,
%U A351133 201155980661221409458717152,1674230688936725338278370413264,18235249164492209082483584810706528
%N A351133 a(n) = Sum_{k=0..n} k! * k^(2*n) * Stirling1(n,k).
%H A351133 Seiichi Manyama, <a href="/A351133/b351133.txt">Table of n, a(n) for n = 0..162</a>
%F A351133 E.g.f.: Sum_{k>=0} log(1 + k^2*x)^k.
%F A351133 a(n) ~ c * d^n * n^(3*n + 1/2), where d = 0.3417329834649268103028466896966197580428514873775849996969994420891... and c = 2.92355271092039591960355156784704285135358... - _Vaclav Kotesovec_, Feb 03 2022
%t A351133 a[0] = 1; a[n_] := Sum[k! * k^(2*n) * StirlingS1[n, k], {k, 1, n}]; Array[a, 13, 0] (* _Amiram Eldar_, Feb 02 2022 *)
%o A351133 (PARI) a(n) = sum(k=0, n, k!*k^(2*n)*stirling(n, k, 1));
%o A351133 (PARI) first(n)=my(x='x+O('x^(n+1))); Vec(serlaplace(sum(k=0, n, log(1+k^2*x)^k)))
%Y A351133 Cf. A006252, A320083, A351134.
%Y A351133 Cf. A229260, A351135, A351136.
%K A351133 nonn
%O A351133 0,3
%A A351133 _Seiichi Manyama_, Feb 02 2022