OFFSET
0,2
COMMENTS
FORMULA
EXAMPLE
G.f.: A(x) = 1 + 2*x + 14*x^2 + 320*x^3 + 21036*x^4 + 4248736*x^5 + 2753284608*x^6 + 5889659239296*x^7 + 42571528094271584*x^8 +...
such that
A(x) = Sum_{n>=0} log(1 + 1^n*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!.
Equivalently,
A(x) = 1 + log(1 + 1*2*x + 2*2^2*x^2 + 3*2^3*x^3 + 4*2^4*x^4 +...) +
log(1 + 1^2*2^2*x + 2^2*2^4*x^2 + 3^2*2^6*x^3 + 4^2*2^8*x^4 +...)^2/2! +
log(1 + 1^3*2^3*x + 2^3*2^6*x^2 + 3^3*2^9*x^3 + 4^3*2^12*x^4 +...)^3/3! +
log(1 + 1^4*2^4*x + 2^4*2^8*x^2 + 3^4*2^12*x^3 + 4^4*2^16*x^4 +...)^4/4! +
...
The g.f. can be written using the Eulerian numbers like so:
A(x) = 1 + log(1 + 2*x/(1-2*x)^2) + log(1 + 2^2*(x + 2^2*x^2)/(1-2^2*x)^3)^2/2! + log(1 + 2^3*(x + 4*2^3*x^2 + 2^6*x^3)/(1-2^3*x)^4)^3/3! + log(1 + 2^4*(x + 11*2^4*x^2 + 11*2^8*x^3 + 2^24*x^4)/(1-2^4*x)^5)^4/4! + log(1 + 2^5*(x + 26*2^5*x^2 + 66*2^10*x^3 + 26*2^15*x^4 + 2^20*x^5)/(1-2^5*x)^6)^5/5! + log(1 + 2^6*(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)/(1-2^6*x)^7)^6/6! +...+ log(1 + Sum_{k=1..n} A008292(n,k) * 2^(n*k) * x^k / (1 - 2^n*x)^(n+1) )^n/n! +...
PROG
(PARI) {a(n) = my(A=1, Oxn=x*O(x^n));
A = sum(m=0, n+1, log(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 = sum(m=0, n+1, log(1 + sum(k=1, m+1, 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
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Oct 08 2016
STATUS
approved