cp's OEIS Frontend

A094417 Generalized ordered Bell numbers Bo(4,n).

%I A094417 #63 Jun 02 2025 05:29:11
%S A094417 1,4,36,484,8676,194404,5227236,163978084,5878837476,237109864804,
%T A094417 10625889182436,523809809059684,28168941794178276,1641079211868751204,
%U A094417 102961115527874385636,6921180217049667005284,496267460209336700111076,37807710659221213027893604
%N A094417 Generalized ordered Bell numbers Bo(4,n).
%C A094417 Fourth row of array A094416, which has more information.
%H A094417 Vincenzo Librandi, <a href="/A094417/b094417.txt">Table of n, a(n) for n = 0..200</a>
%H A094417 Paul Barry, <a href="https://arxiv.org/abs/1803.06408">Three Études on a sequence transformation pipeline</a>, arXiv:1803.06408 [math.CO], 2018.
%F A094417 E.g.f.: 1/(5 - 4*exp(x)).
%F A094417 a(n) = 4 * A050353(n) for n>0.
%F A094417 a(n) = Sum_{k=0..n} A131689(n,k) * 4^k. - _Philippe Deléham_, Nov 03 2008
%F A094417 E.g.f.: A(x) with A_n = 4 * Sum_{k=0..n-1} C(n,k) * A_k; A_0 = 1. - _Vladimir Kruchinin_, Jan 27 2011
%F A094417 G.f.: 2/G(0), where G(k)= 1 + 1/(1 - 8*x*(k+1)/(8*x*(k+1) - 1 + 10*x*(k+1)/G(k+1))); (continued fraction). - _Sergei N. Gladkovskii_, May 30 2013
%F A094417 a(n) = log(5/4)*int {x = 0..inf} (floor(x))^n * (5/4)^(-x) dx. - _Peter Bala_, Feb 14 2015
%F A094417 a(0) = 1; a(n) = 4 * a(n-1) - 5 * Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - _Seiichi Manyama_, Nov 16 2023
%F A094417 From _Seiichi Manyama_, Jun 01 2025: (Start)
%F A094417 a(n) = (-1)^(n+1)/5 * Li_{-n}(5/4), where Li_{n}(x) is the polylogarithm function.
%F A094417 a(n) = (1/5) * Sum_{k>=0} k^n * (4/5)^k.
%F A094417 a(n) = (4/5) * Sum_{k=0..n} 5^k * (-1)^(n-k) * A131689(n,k) for n > 0. (End)
%p A094417 a:= proc(n) option remember;
%p A094417       `if`(n=0, 1, 4* add(binomial(n, k) *a(k), k=0..n-1))
%p A094417     end:
%p A094417 seq(a(n), n=0..20);
%t A094417 max = 16; f[x_] := 1/(5-4*E^x); CoefficientList[Series[f[x], {x, 0, max}], x]*Range[0, max]! (* _Jean-François Alcover_, Nov 14 2011, after g.f. *)
%o A094417 (Magma) m:=20; R<x>:=LaurentSeriesRing(RationalField(), m); b:=Coefficients(R!(1/(5 - 4*Exp(x)))); [Factorial(n-1)*b[n]: n in [1..m]]; // _Bruno Berselli_, Mar 17 2014
%o A094417 (SageMath)
%o A094417 def A094416(n,k): return sum(factorial(j)*n^j*stirling_number2(k,j) for j in range(k+1)) # array
%o A094417 def A094417(k): return A094416(4,k)
%o A094417 [A094417(n) for n in range(31)] # _G. C. Greubel_, Jan 12 2024
%o A094417 (PARI) my(N=25,x='x+O('x^N)); Vec(serlaplace(1/(5 - 4*exp(x)))) \\ _Joerg Arndt_, Jan 15 2024
%Y A094417 Cf. A000670, A004123, A032033, A050353, A094418, A094419, A131689.
%Y A094417 Cf. A346983, A354242, A365567.
%K A094417 nonn
%O A094417 0,2
%A A094417 _Ralf Stephan_, May 02 2004