OFFSET
1,2
LINKS
Luciano Ancora, Table of n, a(n) for n = 1..1000
Luciano Ancora, Partial sums of m-th powers with Faulhaber polynomials
Index entries for linear recurrences with constant coefficients, signature (13,-78,286,-715,1287,-1716,1716,-1287,715,-286,78,-13,1).
FORMULA
G.f.: (- x - 120*x^2 - 1191*x^3 - 2416*x^4 - 1191*x^5 - 120*x^6 - x^7)/(- 1 + x)^13.
a(n) = n*(1 + n)*(2 + n)*(3 + n)*(4 + n)*(5 + n)*(-3 + 5*n + n^2)*(-2 + 5*n + n^2)*(5 + 5*n + n^2)/95040.
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5) + n^7.
EXAMPLE
First differences: 1, 127, 2060, 14324, 63801, ... (A152726)
----------------------------------------------------------------------
The seventh powers: 1, 128, 2187, 16384, 78125, ... (A001015)
----------------------------------------------------------------------
First partial sums: 1, 129, 2316, 18700, 96825, ... (A000541)
Second partial sums: 1, 130, 2446, 21146, 117971, ... (A250212)
Third partial sums: 1, 131, 2577, 23723, 141694, ... (A254641)
Fourth partial sums: 1, 132, 2709, 26432, 168126, ... (A254646)
Fifth partial sums: 1, 133, 2842, 29274, 197400, ... (this sequence)
MATHEMATICA
Table[n (1 + n) (2 + n) (3 + n) (4 + n) (5 + n) (- 3 + 5 n + n^2) (- 2 + 5 n + n^2) (5 + 5 n + n^2)/95040, {n, 21}] (* or *)
CoefficientList[Series[(- 1 - 120 x - 1191 x^2 - 2416 x^3 - 1191 x^4 - 120 x^5 - x^6)/(-1 + x)^13, {x, 0, 20}], x]
PROG
(PARI) a(n)=n*(1+n)*(2+n)*(3+n)*(4+n)*(5+n)*(-3+5*n+n^2)*(-2+5*n+n^2)*(5+5*n+n^2)/95040 \\ Charles R Greathouse IV, Oct 07 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Luciano Ancora, Feb 12 2015
STATUS
approved