OFFSET
0,3
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Marcos MariƱo and Claudia Rella, On the structure of wave functions in complex Chern-Simons theory, arXiv:2312.00624 [hep-th], 2023. See p. 23.
FORMULA
a(n) = Sum_{k=0..floor(n/2)} ( C(n-k, k) + C(n-k-1, k-1) )^2.
Ignoring initial term, equals the logarithmic derivative of A093128, which gives the number of dissections of a polygon using strictly disjoint diagonals. - Paul D. Hanna, Nov 09 2013
From Vaclav Kotesovec, Feb 28 2014: (Start)
Recurrence (for n>=5): (n-3)*n*a(n) = (2*n^2 - 7*n + 4)*a(n-1) + (n-4)*n*a(n-2) + (2*n^2 - 9*n + 8)*a(n-3) - (n-4)*(n-1)*a(n-4).
G.f.: (2-x+2*x^2)/sqrt((x^2+x+1)*(x^2-3*x+1))-1.
a(n) ~ 5^(3/4) * ((3+sqrt(5))/2)^n / (2*sqrt(Pi*n)).
(End)
MATHEMATICA
Flatten[{1, Table[Sum[(Binomial[n-k, k] + Binomial[n-k-1, k-1])^2, {k, 0, Floor[n/2]}], {n, 1, 20}]}] (* Vaclav Kotesovec, Feb 28 2014 *)
PROG
(PARI) {a(n)=sum(k=0, n\2, (binomial(n-k, k)+binomial(n-k-1, k-1))^2)}
for(n=0, 30, print1(a(n), ", "))
(PARI) /* squared sums of negative powers of Catalan series: */
{a(n)=local(Catalan=2/(1+sqrt(1-4*x +x*O(x^n)))); sum(k=0, n\2, polcoeff(Catalan^-n, k)^2)}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Aug 21 2007
STATUS
approved