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 A016036 #47 Jan 20 2025 09:02:54
%S A016036 1,4,31,361,5626,109951,2585269,71066626,2236441141,79289379361,
%T A016036 3127129674736,135802922499949,6439320471558781,331026965612789356,
%U A016036 18338413238239145731,1089132347371148170381,69033182553940825258594,4651256393180943757676371
%N A016036 Row sums of triangle A000369.
%H A016036 Vincenzo Librandi, <a href="/A016036/b016036.txt">Table of n, a(n) for n = 1..200</a>
%H A016036 Wolfdieter Lang, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL3/LANG/lang.html">On generalizations of Stirling number triangles</a>, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
%F A016036 E.g.f.: exp(1 - (1-4*x)^(1/4)) - 1.
%F A016036 a(n) = 6*(2*n-5)*a(n-1) - 3*(16*n^2-96*n+145)*a(n-2) + 2*(4*n-15)*(2*n-7)*(4*n-13)*a(n-3) + a(n-4), n >= 4; a(0) = 1, a(1) = 1, a(2) = 4, a(3) = 31.
%F A016036 a(n) = 1 + (n-1)!*Sum_{m=1..n-1} ( Sum_{k=1..n-m} binomial(n+k-1,n-1) * ( Sum_{j=0..k} binomial(j,n-m-3*k+2*j)*binomial(k,j)*3^(-n+m+3*k-j)*2^(n-m-5*k+3*j)*(-1)^(n-m-k) ) )/(m-1)!. - _Vladimir Kruchinin_, Oct 18 2011
%F A016036 a(n) = D^n(exp(x)) evaluated at x = 0, where D is the operator 1/(1-x)^3*d/dx. Cf. A001515, A015735 and A028575. - _Peter Bala_, Nov 25 2011
%F A016036 a(n) ~ 2^(2*n-3/2)*n^(n-3/4)*exp(1-n)*sqrt(Pi)/Gamma(3/4) * (1 - Gamma(3/4)/(n^(1/4)*sqrt(Pi)) + Gamma(3/4)^2/(4*sqrt(n/2)*Pi)). - _Vaclav Kotesovec_, Aug 10 2013
%F A016036 From _Seiichi Manyama_, Jan 20 2025: (Start)
%F A016036 a(n) = Sum_{k=0..n} (-1)^k * 4^(n-k) * |Stirling1(n,k)| * A000587(k).
%F A016036 a(n) = e * (-4)^n * n! * Sum_{k>=0} (-1)^k * binomial(k/4,n)/k!. (End)
%t A016036 a[n_, m_] /; (n>= m>= 1):= a[n, m]= (4*(n-1)-m)*a[n-1,m] + a[n-1,m-1]; a[n_, m_] /; n<m=0; a[_,0]= 0; a[1,1] = 1; a[n_]:= Sum[a[n,m], {m, n}]; Table[a[n], {n,20}] (* _Jean-François Alcover_, Feb 28 2013 *)
%t A016036 With[{nn=20},CoefficientList[Series[Exp[1-Surd[1-4x,4]]-1,{x,0,nn}],x] Range[0,nn]!]//Rest (* _Harvey P. Dale_, Apr 20 2016 *)
%o A016036 (Maxima)
%o A016036 a(n):=((n-1)!*sum((sum(binomial(n+k-1,n-1)*sum(binomial(j,n-m-3*k+2*j)*binomial(k,j)*3^(-n+m+3*k-j)*2^(n-m-5*k+3*j)*(-1)^(n-m-k),j,0,k),k,1,n-m))/(m-1)!,m,1,n-1))+1; /* _Vladimir Kruchinin_, Oct 18 2011 */
%o A016036 (Magma) R<x>:=PowerSeriesRing(Rationals(), 30); Coefficients(R!(Laplace( Exp(1-(1-4*x)^(1/4)) -1 ))); // _G. C. Greubel_, Oct 02 2023
%o A016036 (SageMath)
%o A016036 def A016036_list(prec):
%o A016036 P.<x> = PowerSeriesRing(QQ, prec)
%o A016036 return P( exp(1-(1-4*x)^(1/4)) -1 ).egf_to_ogf().list()
%o A016036 a=A016036_list(40); a[1:] # _G. C. Greubel_, Oct 02 2023
%Y A016036 Sequences with e.g.f. exp(1-(1-m*x)^(1/m)) - 1: A000012 (m=1), A001515 (m=2), A015735 (m=3), this sequence (m=4), A028575 (m=5), A028844 (m=6).
%Y A016036 Cf. A000369.
%K A016036 nonn
%O A016036 1,2
%A A016036 _Wolfdieter Lang_