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.
%I A293847 #25 Dec 02 2021 19:48:04 %S A293847 1,1,5,49,793,19361,672061,31721425,1963804529,154746407233, %T A293847 15136503333301,1799712380844401,255578390749947145, %U A293847 42713809784784354529,8296411053128532892013,1852797862395580239567121,471358206112272764630500321,135500644700064476406317390465 %N A293847 E.g.f.: exp(Sum_{n>=1} n!*x^n). %H A293847 Seiichi Manyama, <a href="/A293847/b293847.txt">Table of n, a(n) for n = 0..253</a> %F A293847 a(0) = 1 and a(n) = (n-1)! * Sum_{k=1..n} k*k!*a(n-k)/(n-k)! for n > 0. %F A293847 a(n) ~ n!^2. - _Vaclav Kotesovec_, Oct 18 2017 %p A293847 a:= proc(n) option remember; `if`(n=0, 1, add( %p A293847 a(n-i)*binomial(n-1, i-1)*i!^2, i=1..n)) %p A293847 end: %p A293847 seq(a(n), n=0..20); # _Alois P. Heinz_, Dec 02 2021 %t A293847 nmax = 20; CoefficientList[Series[E^Sum[k!*x^k, {k, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]! (* _Vaclav Kotesovec_, Oct 18 2017 *) %o A293847 (PARI) {a(n) = n!*polcoeff(exp(sum(k=1, n, k!*x^k)+x*O(x^n)), n)} %Y A293847 Cf. A158876, A159476. %K A293847 nonn %O A293847 0,3 %A A293847 _Seiichi Manyama_, Oct 17 2017