OFFSET
1,3
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..230
FORMULA
a(n) = [q^n] { ([n+1]_q)! * Sum_{m=0..n+1} (-1)^m * q^(m(m-1)/2) / ([m]_q)! }; here, the q-factorial of n is denoted by ([n]_q)! = Product_{j=1..n} (1-q^j)/(1-q).
a(n) ~ c * 4^n / sqrt(Pi*n), where c = A048651^2 = QPochhammer(1/2)^2 = 0.083398563863748521534541808155312... - Vaclav Kotesovec, Aug 30 2023
PROG
(PARI) {a(n)=polcoeff(prod(j=1, n+1, (1-q^j)/(1-q))* sum(k=0, n+1, (-1)^k*q^(k*(k-1)/2)/if(k==0, 1, prod(j=1, k, (1-q^j)/(1-q)))), n, q)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 05 2008
STATUS
approved