cp's OEIS Frontend

A225116 a(n) = 3^n*A_{n, 1/3}(-1) where A_{n, k}(x) are the generalized Eulerian polynomials.

%I A225116 #19 Mar 19 2020 16:38:13
%S A225116 1,5,24,110,480,2000,8064,32240,130560,531200,2095104,8030720,
%T A225116 33546240,156569600,536838144,243660800,8589803520,244224819200,
%U A225116 137438429184,-28539130347520,2199021158400,4960294141952000,35184363700224,-1015283149035274240,562949919866880
%N A225116 a(n) = 3^n*A_{n, 1/3}(-1) where A_{n, k}(x) are the generalized Eulerian polynomials.
%H A225116 Peter Luschny, <a href="http://oeis.org/wiki/User:Peter_Luschny/EulerianPolynomialsGeneralized">Generalized Eulerian polynomials</a>.
%F A225116 a(n) = 2^(1+n)*(3^n+sum_{j=0..n}(binomial(n,j)*Li_{-j}(-1)*3^(n-j))).
%F A225116 a(n) = 2^(t+1)*(zeta(-t)*(1-2^(t+1))+(2^t-1)). - _Peter Luschny_, Jul 20 2013
%p A225116 EulerianPolynomial := proc(n,k,x) local j; if x = 1 then k^n*n! else (1-x)^(1+n)*(1+add(binomial(n,j)* polylog(-j,x)*k^j, j = 0..n)) fi end:
%p A225116 A225116 := n -> 3^n*EulerianPolynomial(n, 1/3, -1);
%p A225116 seq(round(evalf(A225116(i), 24)), i = 0..24);
%t A225116 Table[2^(t+1)*(Zeta[-t]*(1-2^(t+1))+(2^t-1)), {t,0,24}] (* _Peter Luschny_, Jul 20 2013 *)
%t A225116 Table[EulerE[n, 3] 2^n , {n, 0, 20}] (* _Vladimir Reshetnikov_, Oct 21 2015 *)
%o A225116 (Sage)
%o A225116 from mpmath import mp, polylog
%o A225116 mp.dps = 32; mp.pretty = True
%o A225116 def A225116(n): return 2^(1+n)*(3^n+add(binomial(n,j)*polylog(-j,-1) *3^(n-j) for j in (0..n)))
%o A225116 [int(A225116(n)) for n in (0..24)]
%Y A225116 Cf. A155585(n) = 1^n*A_{n, 1/1}(-1), A119881(n) = 2^n*A_{n, 1/2}(-1).
%K A225116 sign
%O A225116 0,2
%A A225116 _Peter Luschny_, Apr 29 2013