login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A277037 G.f.: exp( Sum_{n>=1} [Sum_{k>=1} k^n * 2^(n*k) * x^k]^n / n ), a power series in x with integer coefficients. 2
1, 2, 18, 484, 54630, 26638924, 53843811956, 442942117297000, 14725418961500037126, 1971239927985067569365772, 1060292226589575099894174194524, 2288290973515256950275126683431946552, 19795837218795604674370624304477542380054748, 685985356865646724678258830150265065104998427771576, 95174256167264272421248219248338459257647770713814222870312 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

More generally, for fixed integer q, G(x,q) = exp( Sum_{n>=1} [Sum_{k>=1} k^n * q^(n*k) * x^k]^n / n ) is an integer series such that G(x,q) = exp( Sum_{n>=1} q^(n^2)/(1 - q^n*x)^(n^2+n) * [ Sum_{k=1..n} A008292(n,k) * q^(n*k-n) * x^k ]^n / n ), where A008292 are the Eulerian numbers.

LINKS

Table of n, a(n) for n=0..14.

FORMULA

G.f.: exp( Sum_{n>=1} [ Sum_{k=1..n} A008292(n,k) * 2^(n*k) * x^k ]^n / (1 - 2^n*x)^(n*(n+1)) / n ), where A008292 are the Eulerian numbers.

EXAMPLE

G.f.: A(x) = 1 + 2*x + 18*x^2 + 484*x^3 + 54630*x^4 + 26638924*x^5 + 53843811956*x^6 + 442942117297000*x^7 +...

such that the logarithm of g.f. A(x) equals the series:

log(A(x)) = Sum_{n>=1} (2^n*x + 2^n*2^(2*n)*x^2 + 3^n*2^(3*n)*x^3 +...+ k^n*2^(k*n)*x^k +...)^n/n.

This logarithmic series can be written using the Eulerian numbers like so:

log(A(x)) = 2*x/(1-2*x)^2 + 2^4*(x + 2^2*x^2)^2/(1-2^2*x)^6/2 + 2^9*(x + 4*2^3*x^2 + 2^6*x^3)^3/(1-2^3*x)^12/3 + 2^16*(x + 11*2^4*x^2 + 11*2^8*x^3 + 2^24*x^4)^4/(1-2^4*x)^20/4 + 2^25*(x + 26*2^5*x^2 + 66*2^10*x^3 + 26*2^15*x^4 + 2^20*x^5)^5/(1-2^5*x)^30/5 + 2^36*(x + 57*2^6*x^2 + 302*2^12*x^3 + 302*2^18*x^4 + 57*2^24*x^5 + 2^30*x^6)^6/(1-2^6*x)^42/6 +...+ [ Sum_{k=1..n} A008292(n,k) * 2^(n*k) * x^k ]^n / (1 - 2^n*x)^(n*(n+1))/n +...

Explicitly,

log(A(x)) = 2*x + 32*x^2/2 + 1352*x^3/3 + 214272*x^4/4 + 132616992*x^5/5 + 322738100480*x^6/6 + 3099838240135296*x^7/7 + 117796258487089512448*x^8/8 +...

PROG

(PARI) {a(n) = my(A=1, Oxn=x*O(x^n)); A = exp( sum(m=1, n+1, sum(k=1, n+1, k^m * 2^(m*k) * x^k +x*O(x^n) )^m / m )); polcoeff(A, n)}

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

(PARI) {A008292(n, k) = sum(j=0, k, (-1)^j * (k-j)^n * binomial(n+1, j))}

{a(n) = my(A=1, Oxn=x*O(x^n)); A = exp( sum(m=1, n+1, sum(k=1, m, A008292(m, k) * 2^(m*k) * x^k / (1 - 2^m*x +Oxn)^(m+1) )^m / m ) ); polcoeff(A, n)}

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

CROSSREFS

Cf. A277035, A156170, A277036, A008292.

Sequence in context: A053916 A355131 A355134 * A015203 A208885 A255433

Adjacent sequences:  A277034 A277035 A277036 * A277038 A277039 A277040

KEYWORD

nonn

AUTHOR

Paul D. Hanna, Oct 08 2016

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 26 06:36 EDT 2022. Contains 354877 sequences. (Running on oeis4.)