%I #9 Jan 04 2024 14:04:14
%S 0,0,0,1,4,10,35,140,476,1624,6070,22495,81455,301301,1131494,4230681,
%T 15852396,59881956,226877648,860447129,3273728234,12493453344,
%U 47760610689,182905145214,701883651799,2697952583635,10385325566785,40033903418860,154534663044346
%N a(n) = (1/6) * 2^(-n) * n! * [x^n] Exp(2*x, 1)*(Exp(2*x, 3) - 1), where Exp(x, m) = Sum_{k>=0} (x^k / k!)^m.
%F a(n) = A344854(n) / 2^n.
%F a(n) = (1/6)*(hypergeom([-n/3, (1 - n)/3, (2 - n)/3], [1, 1], -27) - 1).
%p Exp := (x, m) -> sum((x^k / k!)^m, k=0..infinity):
%p gf := Exp(2*x, 1)*(Exp(2*x, 3) - 1): ser := series(gf, x, 34):
%p seq((1/6)*2^(-n)*n!*simplify(coeff(ser, x, n)), n = 0..28);
%t a[n_] := (1/6) (HypergeometricPFQ[{-n/3, (1 - n)/3, (2 - n)/3}, {1, 1}, -27] - 1);
%t Table[a[n], {n, 0, 28}]
%o (Python)
%o from sympy import hyperexpand, Rational
%o from sympy.functions import hyper
%o def A344559(n): return (hyperexpand(hyper((Rational(-n,3),Rational(1-n,3),Rational(2-n,3)),(1,1),-27))-1)//6 # _Chai Wah Wu_, Jan 04 2024
%Y Cf. A344854.
%K nonn
%O 0,5
%A _Peter Luschny_, Jun 01 2021
|