A032033 - OEIS

%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_