OFFSET
0,3
COMMENTS
If n is divisible by 10, then a(n) is divisible by 1000.
LINKS
Index entries for linear recurrences with constant coefficients, signature (6,-15,20,-15,6,-1).
FORMULA
From Alois P. Heinz, Feb 10 2023: (Start)
a(n) = Sum_{k=0..n} k^2 * A061579(n,k).
G.f.: x*(x+1)*(9*x+1)/(x-1)^6. (End)
From Amiram Eldar, Feb 13 2023: (Start)
Sum_{n>=1} 1/a(n) = 39/8 - 3*Pi^2/4 + 3*zeta(3).
Sum_{n>=1} (-1)^(n+1)/a(n) = 12*log(2) - 51/8 - 3*Pi^2/8 + 9*zeta(3)/4. (End)
EXAMPLE
a(25) = n^3*(n+1)*(n+2)/6 = 25^3*(25+1)*(25+2)/6 = 15625*26*27/6 = 15625*13*9 = 1828125.
MAPLE
with(combinat):a:=n->sum(sum(sum(binomial(n+2, 2)/3, j=1..n), k=1..n), m=1..n): seq(a(n), n=0..31); # Zerinvary Lajos, May 30 2007
MATHEMATICA
a[n_] := n^3*(n+1)*(n+2)/6; Array[a, 35, 0] (* Amiram Eldar, Feb 13 2023 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Brandon Ang (xyz1236(AT)verizon.net), Jun 28 2006
STATUS
approved