OFFSET
0,2
COMMENTS
Compare the g.f. to: Sum_{n>=0} A000108(n)*(1-x)^n*x^n = 1/(1-x).
LINKS
Robert Israel, Table of n, a(n) for n = 0..56
FORMULA
a(n) = [x^n] 2/(1 + sqrt(1 - 4*2^n*(x-x^2))).
a(n) = Sum_{k=0..[n/2]} (-1)^k*2^(n(n-k))*C(n-k,k)*A000108(n-k).
EXAMPLE
G.f.: A(x) = 1 + 2*x + 28*x^2 + 2304*x^3 + 856576*x^4 +...
MAPLE
S:= add(2^(n^2)*binomial(2*n, n)/(n+1)*(1-2^n*x)^n*x^n, n=0..30):
seq(coeff(S, x, n), n=0..30); # Robert Israel, Nov 13 2016
MATHEMATICA
Table[Sum[(-1)^k*2^(n(n-k))*Binomial[n-k, k]*Binomial[2*(n-k), (n-k)]/(n-k+1), {k, 0, Floor[n/2]}], {n, 0, 20}] (* G. C. Greubel, Nov 13 2016 *)
PROG
(PARI) {a(n)=polcoeff(sum(m=0, n, 2^(m^2)*binomial(2*m, m)/(m+1)*x^m*(1-2^m*x)^m)+x*O(x^n), n)}
(PARI) {a(n)=polcoeff(2/(1+sqrt(1-4*2^n*(x-x^2) +x*O(x^n))), n)}
(PARI) {a(n)=sum(k=0, n\2, (-1)^k*2^(n*(n-k))*binomial(n-k, k)*binomial(2*n-2*k, n-k)/(n-k+1))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 29 2009
STATUS
approved