login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A326271 E.g.f.: Sum_{n>=0} 3^n * (exp(n*x) - 1)^n / n!. 4
1, 3, 39, 948, 34869, 1757163, 114320118, 9226773993, 897658726215, 103005144933870, 13705015429716807, 2085418048857405375, 358813807291982519184, 69146346672687725451039, 14803634157756603606592167, 3496440993213535696041009924, 905508769623362527769907535857, 255762146004426658313683324505247, 78413243604482526944814375388910526, 25983968388767783226046397856822603645 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

More generally, the following sums are equal:

(1) Sum_{n>=0} (p + q^n)^n * r^n/n!,

(2) Sum_{n>=0} q^(n^2) * exp(p*q^n*r) * r^n/n!;

here, q = exp(x) with p = -1, r = 3.

In general, let F(x) be a formal power series in x such that F(0)=1, then

Sum_{n>=0} m^n * F(q^n*r)^b * log( F(q^n*r) )^n / n! =

Sum_{n>=0} r^n * [y^n] F(y)^(m*q^n + p);

here, F(x) = exp(x), q = exp(x), p = -1, r = 3, m = 1.

LINKS

Paul D. Hanna, Table of n, a(n) for n = 0..300

FORMULA

E.g.f.: Sum_{n>=0} 3^n * (exp(n*x) - 1)^n / n!.

E.g.f.: Sum_{n>=0} 3^n * exp(n^2*x) * exp( -3*exp(n*x) ) / n!.

O.g.f.: Sum_{n>=0} 3^n * n^n * x^n / Product_{k=1..n} (1 - n*k*x).

a(n) = Sum_{k=0..n} 3^k * k^n * Stirling2(n,k).

EXAMPLE

E.g.f.: A(x) = 1 + 3*x + 39*x^2/2! + 948*x^3/3! + 34869*x^4/4! + 1757163*x^5/5! + 114320118*x^6/6! + 9226773993*x^7/7! + 897658726215*x^8/8! + 103005144933870*x^9/9! + ...

such that

A(x) = 1 + 3*(exp(x) - 1) + 3^2*(exp(2*x) - 1)^2/2! + 3^3*(exp(3*x) - 1)^3/3! + 3^4*(exp(4*x) - 1)^4/4! + 3^5*(exp(5*x) - 1)^5/5! + 3^6*(exp(6*x) - 1)^6/6! + ...

also

A(x) = exp(-3) + 3*exp(x)*exp(-3*exp(x)) + 3^2*exp(4*x)*exp(-3*exp(2*x))/2! + 3^3*exp(9*x)*exp(-3*exp(3*x))/3! + 3^4*exp(16*x)*exp(-3*exp(4*x))/4! + 3^5*exp(25*x)*exp(-3*exp(5*x))/5! + 3^6*exp(36*x)*exp(-3*exp(6*x))/6! + ...

ORDINARY GENERATING FUNCTION.

O.g.f.: B(x) = 1 + 3*x + 39*x^2 + 948*x^3 + 34869*x^4 + 1757163*x^5 + 114320118*x^6 + 9226773993*x^7 + 897658726215*x^8 + 103005144933870*x^9 + ...

such that

B(x) = 1 + 3*x/(1-x) + 3^2*2^2*x^2/((1-2*x)*(1-4*x)) + 3^3*3^3*x^3/((1-3*x)*(1-6*x)*(1-9*x)) + 3^4*4^4*x^4/((1-4*x)*(1-8*x)*(1-12*x)*(1-16*x)) + 3^5*5^5*x^5/((1-5*x)*(1-10*x)*(1-15*x)*(1-20*x)*(1-25*x)) + ...

PROG

(PARI) {a(n) = sum(k=0, n, 3^k * k^n * stirling(n, k, 2) )}

for(n=0, 30, print1(a(n), ", "))

(PARI) /* E.g.f.: Sum_{n>=0} 3^n * (exp(n*x) - 1)^n / n! */

{a(n) = n! * polcoeff(sum(m=0, n, 3^m * (exp(m*x +x*O(x^n)) - 1)^m / m!), n)}

for(n=0, 30, print1(a(n), ", "))

(PARI) /* O.g.f.: Sum_{n>=0} 3^n * n^n * x^n / Product_{k=1..n} (1 - n*k*x) */

{a(n) = polcoeff(sum(m=0, n, 3^m * m^m * x^m / prod(k=1, m, 1-m*k*x +x*O(x^n))), n)}

for(n=0, 30, print1(a(n), ", "))

CROSSREFS

Cf. A108459, A326270, A326288.

Sequence in context: A092610 A301585 A258923 * A276964 A274573 A278750

Adjacent sequences:  A326268 A326269 A326270 * A326272 A326273 A326274

KEYWORD

nonn

AUTHOR

Paul D. Hanna, Jun 28 2019

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 19 03:54 EST 2020. Contains 331031 sequences. (Running on oeis4.)