OFFSET
0,1
LINKS
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).
FORMULA
G.f.: (3 + 14*x - x^2)/(1 - x)^4.
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4).
a(n) = Sum_{k=0..n} (8*k^2 + 12*k + 3).
Sum_{n>=0} 1/a(n) = 3*(80*log(2) + 5*Pi - 48)/175 = 0.397024075075621559...
EXAMPLE
a(0) = 0*2 + 1*3 = 3;
a(1) = 0*2 + 1*3 + 2*4 + 3*5 = 26;
a(2) = 0*2 + 1*3 + 2*4 + 3*5 + 4*6 + 5*7 = 85;
a(3) = 0*2 + 1*3 + 2*4 + 3*5 + 4*6 + 5*7 + 6*8 + 7*9 = 196;
a(4) = 0*2 + 1*3 + 2*4 + 3*5 + 4*6 + 5*7 + 6*8 + 7*9 + 8*10 + 9*11 = 375, etc.
MATHEMATICA
Table[(n + 1) (2 n + 1) (4 n + 9)/3, {n, 0, 38}]
LinearRecurrence[{4, -6, 4, -1}, {3, 26, 85, 196}, 39]
Table[Sum[8 k^2 + 12 k + 3, {k, 0, n}], {n, 0, 38}]
PROG
(PARI) Vec((3 + 14*x - x^2)/(1 - x)^4 + O(x^50)) \\ Michel Marcus, Feb 25 2016
(Magma) [(n+1)*(2*n+1)*(4*n+9)/3: n in [0..50]]; // Vincenzo Librandi, Feb 25 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Ilya Gutkovskiy, Feb 24 2016
STATUS
approved