Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
E.g.f.: Sum_{n>=0} (exp(2^n*x) - 1)^n/n!, an analog of the Bell numbers (A000110).
5

%I #6 Oct 08 2016 08:03:57

%S 1,2,20,712,91920,44874784,85939843136,660213878210688,

%T 20540390859740217600,2592165941692975372042752,

%U 1324271564605167892188248409088,2730585827960928853182474922961668096

%N E.g.f.: Sum_{n>=0} (exp(2^n*x) - 1)^n/n!, an analog of the Bell numbers (A000110).

%F a(n) = exp(-2^n) * Sum_{k>=0} (2^k*k)^n/k!.

%F a(n) = [x^n/n! ] B(x)^(2^n), where B(x) = exp(exp(x) - 1) is the e.g.f. of the Bell numbers.

%F a(n) = Sum_{k=0..n} S2(n,k) * 2^(n*k), where S2(n,k) = A008277(n,k) are the Stirling numbers of the second kind.

%F G.f.: A(x) = Sum_{n>=0} 2^(n^2) * x^n / [Product_{k=1..n} (1 - k*2^n*x)].

%e E.g.f.: A(x) = 1 + 2*x + 20*x^2/2! + 712*x^3/3! + 91920*x^4/4! +...

%e A(x) = 1 + (exp(2*x) - 1) + (exp(4*x) - 1)^2/2! + (exp(8*x) - 1)^3/3! +...+ (exp(2^n*x) - 1)^n/n! +...

%e a(n) = coefficient of x^n/n! in Bell(x)^(2^n) where Bell(x) = exp(exp(x)-1):

%e Bell(x) = 1 + x + 2*x^2/2! + 5*x^3/3! + 15*x^4/4! + 52*x^5/5! + 203*x^6/6! +...+ A000110(n)*x^n/n! +...

%o (PARI) {a(n)=local(infnty=n^4+10);round(exp(-2^n)*sum(k=0,infnty,(2^k*k)^n/k!))}

%o (PARI) {a(n)=n!*polcoeff(sum(k=0,n,(exp(2^k*x +x*O(x^n))-1)^k/k!),n)}

%o (PARI) {a(n)=n!*polcoeff(exp(2^n*(exp(x +x*O(x^n))-1)),n)}

%o (PARI) {S2(n,k)=(1/k!)*sum(i=0,k,(-1)^(k-i)*binomial(k,i)*i^n)}

%o {a(n)=sum(k=0,n,S2(n,k)*2^(n*k))}

%o (PARI) {a(n)=polcoeff(sum(k=0,n,(2^k*x)^k/prod(j=1,k,1-j*2^k*x+x*O(x^n))),n)}

%Y Cf. A000110, A168408, A277036, A008277.

%K nonn

%O 0,2

%A _Paul D. Hanna_, Nov 25 2009, Nov 25 2009, Feb 16 2010