OFFSET
0,2
COMMENTS
Binomial transform is A014915. In general, ((1+x)/(1-r*x))^2 expands to a(n) = ((r+1)*r^n*((r+1)*n + r - 1) + 0^n)/r^2, which is also a(n) = Sum_{k=0..n} C(n,k)*Sum_{j=0..k} (j+1)*(r+1)^j. This is the self-convolution of the coordination sequence for the infinite tree with valency r.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (4,-4).
FORMULA
G.f.: (1+x)^2/(1-2x)^2;
a(n) = 3*2^n(3n+1)/4 + 0^n/4;
a(n) = Sum_{k=0..n} C(n, k)*Sum_{j=0..k} (j+1)*3^j.
a(n) = 4*a(n-1) - 4*a(n-2); a(0)=1, a(1)=6, a(2)=21. - Harvey P. Dale, May 20 2011
Sum_{n>=0} (-1)^n/a(n) = (2^(4/3)*(log(1+2^(1/3)+2^(2/3)) + 2*sqrt(3)*arccot((2^(4/3)-1)/sqrt(3))) - 3)/9. - Amiram Eldar, Feb 04 2026
MATHEMATICA
Join[{1}, LinearRecurrence[{4, -4}, {6, 21}, 30]] (* or *) CoefficientList[ Series[((1+x)/(1-2x))^2, {x, 0, 30}], x] (* Harvey P. Dale, May 20 2011 *)
PROG
(Magma) [3*2^n*(3*n+1)/4+0^n/4: n in [0..30]]; // Vincenzo Librandi, May 21 2011
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Paul Barry, Oct 14 2005
STATUS
approved
