|
%I #61 Jan 15 2024 00:53:04
%S 1,3,21,219,3045,52923,1103781,26857659,746870565,23365498683,
%T 812198635941,31055758599099,1295419975298085,58538439796931643,
%U 2848763394161128101,148537065755389540539,8261178848690959117605,488177936257344615487803,30544839926043868901604261
%N Stirling transform of A032031.
%C Also "AIJ" (ordered, indistinct, labeled) transform of 3,3,3,3...
%C Third row of array A094416 (generalized ordered Bell numbers).
%H Vincenzo Librandi, <a href="/A032033/b032033.txt">Table of n, a(n) for n = 0..200</a>
%H Paul Barry, <a href="https://arxiv.org/abs/1803.06408">Three Études on a sequence transformation pipeline</a>, arXiv:1803.06408 [math.CO], 2018.
%H P. Blasiak, K. A. Penson and A. I. Solomon, <a href="http://arXiv.org/abs/quant-ph/0303030">Dobinski-type relations and the log-normal distribution</a>, arXiv:quant-ph/0303030, J. Phys. A.: Math. Gen 36 (2003) L273.
%H C. G. Bower, <a href="/transforms2.html">Transforms (2)</a>
%H Jacob Sprittulla, <a href="https://arxiv.org/abs/2008.09984">On Colored Factorizations</a>, arXiv:2008.09984 [math.CO], 2020.
%F E.g.f.: 1/(4-3*exp(x)).
%F a(n) = 3 * A050352(n), n > 0.
%F a(n) = Sum_{k=0..n} Stirling2(n, k)*(3^k)*k!.
%F a(n) = (1/4)*Sum_{k>=0} k^n*(3/4)^k. - _Karol A. Penson_, Jan 25 2002
%F a(n) = Sum_{k=0..n} A131689(n,k)*3^k. - _Philippe Deléham_, Nov 03 2008
%F G.f. A(x)=B(x)/x, where B(x)=x+3*x^2+21*x^3+... = Sum_{n>=1} b(n)*x^n satisfies 4*B(x)-x = 3*B(x/(1-x)), and b(n)=3*Sum_{k=1..n-1} binomial(n-1,k-1)*b(k), b(1)=1. - _Vladimir Kruchinin_, Jan 27 2011
%F a(n) = log(4/3)*Integral_{x = 0..inf} (floor(x))^n * (4/3)^(-x) dx. - _Peter Bala_, Feb 14 2015
%F a(0) = 1; a(n) = 3 * Sum_{k=1..n} binomial(n,k) * a(n-k). - _Ilya Gutkovskiy_, Jan 17 2020
%F a(0) = 1; a(n) = 3*a(n-1) - 4*Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - _Seiichi Manyama_, Nov 16 2023
%p b:= proc(n, m) option remember;
%p `if`(n=0, 3^m*m!, m*b(n-1, m)+b(n-1, m+1))
%p end:
%p a:= n-> b(n, 0):
%p seq(a(n), n=0..20); # _Alois P. Heinz_, Aug 04 2021
%t a[n_] := PolyLog[-n, 3/4]/4; a[0] = 1; Table[a[n], {n, 0, 16}] (* _Jean-François Alcover_, Nov 14 2011 *)
%t t = 30; Range[0, t]! CoefficientList[Series[1/(4 - 3 Exp[x]), {x, 0, t}], x] (* _Vincenzo Librandi_, Mar 16 2014 *)
%o (PARI) a(n)=ceil(polylog(-n,3/4)/4) \\ _Charles R Greathouse IV_, Jul 14 2014
%o (PARI) my(N=25,x='x+O('x^N)); Vec(serlaplace(1/(4 - 3*exp(x)))) \\ _Joerg Arndt_, Jan 15 2024
%Y Cf. A032031, A094416, A094418, A094419.
%K nonn,easy
%O 0,2
%A _Christian G. Bower_
|
