A007487 -id:A007487 - cp's OEIS frontend

A000538 Sum of fourth powers: 0^4 + 1^4 + ... + n^4.

0, 1, 17, 98, 354, 979, 2275, 4676, 8772, 15333, 25333, 39974, 60710, 89271, 127687, 178312, 243848, 327369, 432345, 562666, 722666, 917147, 1151403, 1431244, 1763020, 2153645, 2610621, 3142062, 3756718, 4463999, 5273999, 6197520, 7246096, 8432017, 9768353

Offset: 0

N. J. A. Sloane

This sequence is related to A000537 by the transform a(n) = n*A000537(n) - Sum_{i=0..n-1} A000537(i). - Bruno Berselli, Apr 26 2010
A formula for the r-th successive summation of k^4, for k = 1 to n, is ((12*n^2+(12*n-5)*r+r^2)*(2*n+r)*(n+r)!)/((r+4)!*(n-1)!), (H. W. Gould). - Gary Detlefs, Jan 02 2014
The number of four dimensional hypercubes in a 4D grid with side lengths n. This applies in general to k dimensions. That is, the number of k-dimensional hypercubes in a k-dimensional grid with side lengths n is equal to the sum of 1^k + 2^k + ... + n^k. - Alejandro Rodriguez, Oct 20 2020

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 813.
A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 222.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 155.
R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1991, p. 275.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 0..1000
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
J. L. Bailey, Jr., A table to facilitate the fitting of certain logistic curves, Annals Math. Stat., 2 (1931), 355-359.
J. L. Bailey, A table to facilitate the fitting of certain logistic curves, Annals Math. Stat., 2 (1931), 355-359. [Annotated scanned copy]
Bruno Berselli, A description of the transform in Comments lines: website Matem@ticamente (in Italian).
Stefano Capparelli, Notes on Discrete Math, Società Editrice Esculapio SRL (2019) 3-4.
Feihu Liu, Guoce Xin, and Chen Zhang, Ehrhart Polynomials of Order Polytopes: Interpreting Combinatorial Sequences on the OEIS, arXiv:2412.18744 [math.CO], 2024. See p. 13.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Eric Weisstein, MathWorld: Faulhaber's Formula
Wikipedia, Faulhaber's formula
Index entries for linear recurrences with constant coefficients, signature (6,-15,20,-15,6,-1).

Cf. A000217, A000330, A000537, A000539, A000540, A000541, A000542, A007487, A023002, A064538, A101089.
Row 4 of array A103438.
Cf. A000583.
Cf. A254640.

a000538 n = (3 * n * (n + 1) - 1) * (2 * n + 1) * (n + 1) * n `div` 30
-- Reinhard Zumkeller, Nov 11 2012

[n*(1+n)*(1+2*n)*(-1+3*n+3*n^2)/30: n in [0..35]]; // Vincenzo Librandi, Apr 04 2015

A000538 := n-> n*(n+1)*(2*n+1)*(3*n^2+3*n-1)/30;

Accumulate[Range[0,40]^4] (* Harvey P. Dale, Jan 13 2011 *)
CoefficientList[Series[x (1 + 11 x + 11 x^2 + x^3)/(1 - x)^6, {x, 0, 40}], x] (* Vincenzo Librandi, Dec 07 2015 *)
LinearRecurrence[{6, -15, 20, -15, 6, -1}, {0, 1, 17, 98, 354, 979}, 35] (* Jean-François Alcover, Feb 09 2016 *)
Table[x^5/5+x^4/2+x^3/3-x/30,{x,40}] (* Harvey P. Dale, Jun 06 2021 *)

A000538(n):=n*(n+1)*(2*n+1)*(3*n^2+3*n-1)/30$
makelist(A000538(n),n,0,30); /* Martin Ettl, Nov 12 2012 */

a(n) = n*(1+n)*(1+2*n)*(-1+3*n+3*n^2)/30 \\ Charles R Greathouse IV, Nov 20 2012

concat(0, Vec(x*(1+11*x+11*x^2+x^3)/(1-x)^6 + O(x^100))) \\ Altug Alkan, Dec 07 2015

A000538_list, m = [0], [24, -36, 14, -1, 0, 0]
for _ in range(10**2):
    for i in range(5):
        m[i+1] += m[i]
    A000538_list.append(m[-1]) # Chai Wah Wu, Nov 05 2014

def A000538(n): return n*(n**2*(n*(6*n+15)+10)-1)//30 # Chai Wah Wu, Oct 03 2024

[bernoulli_polynomial(n,5)/5 for n in range(1, 35)] # Zerinvary Lajos, May 17 2009

a(n) = n*(1+n)*(1+2*n)*(-1+3*n+3*n^2)/30.
The preceding formula is due to al-Kachi (1394-1437). - Juri-Stepan Gerasimov, Jul 12 2009
G.f.: x*(x+1)*(1+10*x+x^2)/(1-x)^6. Simon Plouffe in his 1992 dissertation. More generally, the o.g.f. for Sum_{k=0..n} k^m is x*E(m, x)/(1-x)^(m+2), where E(m, x) is the Eulerian polynomial of degree m (cf. A008292). The e.g.f. for these o.g.f.s is: x/(1-x)^2*(exp(y/(1-x))-exp(x*y/(1-x)))/(exp(x*y/(1-x))-x*exp(y/(1-x))). - Vladeta Jovovic, May 08 2002
a(n) = Sum_{i = 1..n} J_4(i)*floor(n/i), where J_4 is A059377. - Enrique Pérez Herrero, Feb 26 2012
a(n) = 5*a(n-1) - 10* a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5) + 24. - Ant King, Sep 23 2013
a(n) = -Sum_{j=1..4} j*Stirling1(n+1,n+1-j)*Stirling2(n+4-j,n). - Mircea Merca, Jan 25 2014
Sum_{n>=1} (-1)^(n+1)/a(n) = -30*(4 + 3/cos(sqrt(7/3)*Pi/2))*Pi/7. - Vaclav Kotesovec, Feb 13 2015
a(n) = (n + 1)*(n + 1/2)*n*(n + 1/2 + sqrt(7/12))*(n + 1/2  - sqrt(7/12))/5, see the Graham et al. reference, p. 275. - Wolfdieter Lang, Apr 02 2015

The general V. Jovovic formula has been slightly changed after his approval by Wolfdieter Lang, Nov 03 2011

A103438 Square array T(m,n) read by antidiagonals: Sum_{k=1..n} k^m.

Original entry on oeis.org

0, 0, 1, 0, 1, 2, 0, 1, 3, 3, 0, 1, 5, 6, 4, 0, 1, 9, 14, 10, 5, 0, 1, 17, 36, 30, 15, 6, 0, 1, 33, 98, 100, 55, 21, 7, 0, 1, 65, 276, 354, 225, 91, 28, 8, 0, 1, 129, 794, 1300, 979, 441, 140, 36, 9, 0, 1, 257, 2316, 4890, 4425, 2275, 784, 204, 45, 10

Offset: 0

Ralf Stephan, Feb 11 2005

For the o.g.f.s of the column sequences for this array, see A196837 and the link given there. - Wolfdieter Lang, Oct 15 2011
T(m,n)/n is the m-th moment of the discrete uniform distribution on {1,2,...,n}. - Geoffrey Critzer, Dec 31 2018
T(1,n) divides T(m,n) for odd m. - Franz Vrabec, Dec 23 2020

			Square array begins:
  0, 1,  2,   3,    4,     5,     6,      7,      8,      9, ... A001477;
  0, 1,  3,   6,   10,    15,    21,     28,     36,     45, ... A000217;
  0, 1,  5,  14,   30,    55,    91,    140,    204,    285, ... A000330;
  0, 1,  9,  36,  100,   225,   441,    784,   1296,   2025, ... A000537;
  0, 1, 17,  98,  354,   979,  2275,   4676,   8772,  15333, ... A000538;
  0, 1, 33, 276, 1300,  4425, 12201,  29008,  61776, 120825, ... A000539;
  0, 1, 65, 794, 4890, 20515, 67171, 184820, 446964, 978405, ... A000540;
Antidiagonal triangle begins as:
  0;
  0, 1;
  0, 1,  2;
  0, 1,  3,  3;
  0, 1,  5,  6,  4;
  0, 1,  9, 14, 10,  5;
  0, 1, 17, 36, 30, 15, 6;

J. Faulhaber, Academia Algebrae, Darinnen die miraculosische inventiones zu den höchsten Cossen weiters continuirt und profitirt werden, Augspurg, bey Johann Ulrich Schönigs, 1631.

G. C. Greubel, Antidiagonals n = 0..50, flattened
José L. Cereceda, Sums of powers of integers and hyperharmonic numbers, arXiv:2005.03407 [math.NT], 2020.
T. A. Gulliver, Divisibility of sums of powers of odd integers, Int. Math. For. 5 (2010) 3059-3066.
T. A. Gulliver, Sums of Powers of Integers Divisible by Three, Int. J. Contemp. Math. Sciences, Vol. 7, 2012, no. 38, pp. 1895-1901. - From _N. J. A. Sloane_, Dec 22 2012
V. J. W. Guo and J. Zeng, A q-analogue of Faulhaber's formula for sums of powers, arXiv:math/0501441 [math.CO], 2005.
H. Helfgott and I. M. Gessel, Enumeration of tilings of diamonds and hexagons with defects, arXiv:math/9810143 [math.CO], 1998.
T. Kim, q-analogues of the sums of powers of consecutive integers, arXiv:math/0502113 [math.NT], 2005.
D. E. Knuth, Johann Faulhaber and sums of powers, Math. Comp. 61 (1993), no. 203, 277-294.
Eric Weisstein's World of Mathematics, Discrete Uniform Distribution.
Wikipedia, Faulhaber's formula

Rows include A000027, A000217, A000330, A000537, A000538, A000539, A000540, A000541, A000542, A007487, A023002.
Columns include A000051, A001550, A001551, A001552, A001553, A001554, A001555, A001556, A001557.
Diagonals include A076015 and A031971.
Antidiagonal sums are in A103439.
Antidiagonals are the rows of triangle A192001.

T:= func< n,k | n eq 0 select k else (&+[j^n: j in [0..k]]) >;
[T(n-k,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Dec 22 2021

seq(print(seq(Zeta(0,-k,1)-Zeta(0,-k,n+1),n=0..9)),k=0..6);
# (Produces the square array from the example.) Peter Luschny, Nov 16 2008
# alternative
A103438 := proc(m,n)
    (bernoulli(m+1,n+1)-bernoulli(m+1))/(m+1) ;
    if m = 0 then
        %-1 ;
    else
        % ;
    end if;
end proc: # R. J. Mathar, May 10 2013
# simpler:
A103438 := proc(m,n)
    (bernoulli(m+1,n+1)-bernoulli(m+1,1))/(m+1) ;
end proc: # Peter Luschny, Mar 20 2024

T[m_, n_]:= HarmonicNumber[m, -n]; Flatten[Table[T[m-n, n], {m, 0, 11}, {n, m, 0, -1}]] (* Jean-François Alcover, May 11 2012 *)

PARI
```
T(m,n)=sum(k=0,n,k^m)
```

from itertools import count, islice
from math import comb
from fractions import Fraction
from sympy import bernoulli
def A103438_T(m,n): return sum(k**m for k in range(1,n+1)) if n<=m else int(sum(comb(m+1,i)*(bernoulli(i) if i!=1 else Fraction(1,2))*n**(m-i+1) for i in range(m+1))/(m+1))
def A103438_gen(): # generator of terms
    for m in count(0):
        for n in range(m+1):
            yield A103438_T(m-n,n)
A103438_list = list(islice(A103438_gen(),100)) # Chai Wah Wu, Oct 23 2024

def T(n,k): return (bernoulli_polynomial(k+1, n+1) - bernoulli_polynomial(1, n+1)) /(n+1)
flatten([[T(n-k,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Dec 22 2021

E.g.f.: e^x*(e^(x*y)-1)/(e^x-1).
T(m, n) = Zeta(-n, 1) - Zeta(-n, m + 1), for m >= 0 and n >= 0, where Zeta(z, v) is the Hurwitz zeta function. - Peter Luschny, Nov 16 2008
T(m, n) = HarmonicNumber(m, -n). - Jean-François Alcover, May 11 2012
T(m, n) = (Bernoulli(m + 1, n + 1) - Bernoulli(m + 1, 1)) / (m + 1). - Peter Luschny, Mar 20 2024
T(m, n) = Sum_{k=0...m-n} B(k)*(-1)^k*binomial(m-n,k)*n^(m-n-k+1)/(m-n-k+1), where B(k) = Bernoulli number A027641(k) / A027642(k). - Robert B Fowler, Aug 20 2024
T(m, n) = Sum_{i=1..n} J_m(i)*floor(n/i), where J_m is the m-th Jordan totient function. - Ridouane Oudra, Jul 19 2025

A023002 Sum of 10th powers.

Original entry on oeis.org

0, 1, 1025, 60074, 1108650, 10874275, 71340451, 353815700, 1427557524, 4914341925, 14914341925, 40851766526, 102769130750, 240627622599, 529882277575, 1106532668200, 2206044295976, 4222038196425, 7792505423049, 13923571680850

Offset: 0

David W. Wilson

T. D. Noe, Table of n, a(n) for n = 0..1000
Bruno Berselli, A description of the recursive method in Formula lines (second formula): website Matem@ticamente (in Italian).
Feihu Liu, Guoce Xin, and Chen Zhang, Ehrhart Polynomials of Order Polytopes: Interpreting Combinatorial Sequences on the OEIS, arXiv:2412.18744 [math.CO], 2024. See p. 13.
Eric Weisstein's World of Mathematics, Power Sum.
Index entries for linear recurrences with constant coefficients, signature (12,-66,220,-495,792,-924,792,-495,220,-66,12,-1).

Sequences of the form Sum_{j=0..n} j^m : A000217 (m=1), A000330 (m=2), A000537 (m=3), A000538 (m=4), A000539 (m=5), A000540 (m=6), A000541 (m=7), A000542 (m=8), A007487 (m=9), this sequence (m=10), A123095 (m=11), A123094 (m=12), A181134 (m=13).
Row 10 of array A103438.
Cf. A215083, A069095.

[&+[n^10: n in [0..m]]: m in [0..19]];  // Bruno Berselli, Aug 23 2011

A023002:= n-> bernoulli(11, n+1)/11; seq(A023002(n), n=0..30); # G. C. Greubel, Jul 21 2021

Table[Sum[k^10, {k, n}], {n, 0, 30}] (* Vladimir Joseph Stephan Orlovsky, Aug 14 2008 *)
Accumulate[Range[0,20]^10] (* Harvey P. Dale, Aug 23 2011 *)

a(n)=(6*x^11+33*x^10+55*x^9-66*x^7+66*x^5-33*x^3+5*x)/66 \\ Charles R Greathouse IV, Aug 23 2011

a(n)=sum(i=0,10,binomial(11,i)*bernfrac(i)*n^(11-i))/11+n^10 \\ Charles R Greathouse IV, Aug 23 2011

A023002_list, m = [0], [3628800, -16329600, 30240000, -29635200, 16435440, -5103000, 818520, -55980, 1022, -1, 0 , 0]
for _ in range(20):
    for i in range(11):
        m[i+1]+= m[i]
    A023002_list.append(m[-1])
print(A023002_list) # Chai Wah Wu, Nov 05 2014

[bernoulli_polynomial(n,11)/11 for n in range(2, 21)]# Zerinvary Lajos, May 17 2009

a(n) = n*(n+1)*(2*n+1)*(n^2+n-1)(3*n^6 +9*n^5 +2*n^4 -11*n^3 +3*n^2 +10*n -5)/66 (see MathWorld, Power Sum, formula 40). - Bruno Berselli, Apr 26 2010
a(n) = n*A007487(n) - Sum_{i=0..n-1} A007487(i). - Bruno Berselli, Apr 27 2010
From Bruno Berselli, Aug 23 2011: (Start)
a(n) = -a(-n-1).
G.f.: x*(1+x)*(1 +1012*x +46828*x^2 +408364*x^3 +901990*x^4 +408364*x^5 +46828*x^6 +1012*x^7 +x^8)/(1-x)^12. (End)
a(n) = (-1)*Sum_{j=1..10} j*Stirling1(n+1,n+1-j)*Stirling2(n+10-j,n). - Mircea Merca, Jan 25 2014
a(n) = Sum_{i=1..n} J_10(i)*floor(n/i), where J_10 is A069095. - Ridouane Oudra, Jul 17 2025

A123094 Sum of first n 12th powers.

Original entry on oeis.org

0, 1, 4097, 535538, 17312754, 261453379, 2438235715, 16279522916, 84998999652, 367428536133, 1367428536133, 4505856912854, 13421957361110, 36720042483591, 93413954858887, 223160292749512, 504635269460168, 1087257506689929, 2244088888116105, 4457403807182266

Offset: 0

Zerinvary Lajos, Sep 27 2006

T. D. Noe, Table of n, a(n) for n = 0..1000
Bruno Berselli, A description of the recursive method in Formula lines (first formula): website Matem@ticamente (in Italian).

Sequences of the form Sum_{j=0..n} j^m : A000217 (m=1), A000330 (m=2), A000537 (m=3), A000538 (m=4), A000539 (m=5), A000540 (m=6), A000541 (m=7), A000542 (m=8), A007487 (m=9), A023002 (m=10), A123095 (m=11), this sequence (m=12), A181134 (m=13).
Cf. A008456, A215083.
Cf. A008275, A008277.

[(&+[j^12: j in [0..n]]): j in [0..30]]; // G. C. Greubel, Jul 21 2021

Maple
```
[seq(add(i^12, i=1..n), n=0..18)];
```

Table[Sum[k^12, {k, n}], {n, 0, 30}] (* Vladimir Joseph Stephan Orlovsky, Aug 14 2008 *)
Accumulate[Range[0,30]^12]  (* Harvey P. Dale, Apr 26 2011 *)

A123094_list, m = [0], [479001600, -2634508800, 6187104000, -8083152000, 6411968640, -3162075840, 953029440, -165528000, 14676024, -519156, 4094, -1, 0 , 0]
for _ in range(10**2):
    for i in range(13):
        m[i+1]+= m[i]
    A123094_list.append(m[-1]) # Chai Wah Wu, Nov 05 2014

[bernoulli_polynomial(n,13)/13 for n in range(1, 30)] # Zerinvary Lajos, May 17 2009

a(n) = n*A123095(n) - Sum_{i=0..n-1} A123095(i). - Bruno Berselli, Apr 27 2010
a(n) = n * (n+1) * (2*n+1) * (105*n^10 +525*n^9 +525*n^8 -1050*n^7 -1190*n^6 +2310*n^5 +1420*n^4 -3285*n^3 -287*n^2 +2073*n -691)/2730. - Bruno Berselli, Oct 03 2010
a(n) = (-1)*Sum_{j=1..12} j*Stirling1(n+1,n+1-j)*Stirling2(n+12-j,n). - Mircea Merca, Jan 25 2014

A123095 Sum of first n 11th powers.

Original entry on oeis.org

0, 1, 2049, 179196, 4373500, 53201625, 415998681, 2393325424, 10983260016, 42364319625, 142364319625, 427675990236, 1170684360924, 2962844754961, 7012409924625, 15662165784000, 33254351828416, 67526248136049, 131794658215281, 248284917113500, 453084917113500

Offset: 0

Zerinvary Lajos, Sep 27 2006

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Bruno Berselli, A description of the recursive method in Comments lines: website Matem@ticamente (in Italian).
Index entries for linear recurrences with constant coefficients, signature (13,-78,286,-715,1287,-1716,1716,-1287,715,-286,78,-13,1).

Sequences of the form Sum_{j=0..n} j^m : A000217 (m=1), A000330 (m=2), A000537 (m=3), A000538 (m=4), A000539 (m=5), A000540 (m=6), A000541 (m=7), A000542 (m=8), A007487 (m=9), A023002 (m=10), this sequence (m=11), A123094 (m=12), A181134 (m=13).
Cf. A008455, A215083.
Cf. A008275, A008277.

[(&+[j^11: j in [0..n]]): n in [0..30]]; // G. C. Greubel, Jul 21 2021

[seq(add(i^11, i=1..n), n=0..20)];
a[0]:=0:a[1]:=1:for n from 2 to 50 do a[n]:=a[n-1]+n^11 od: seq(a[n], n=0..13); # Zerinvary Lajos, Feb 22 2008

Table[Sum[k^11, {k, n}], {n, 0, 30}] (* Vladimir Joseph Stephan Orlovsky, Aug 14 2008 *)
Accumulate[Range[0,20]^11] (* Harvey P. Dale, Sep 17 2021 *)

A123095_list, m = [0], [39916800, -199584000, 419126400, -479001600, 322494480, -129230640, 29607600, -3498000, 171006, -2046, 1, 0 , 0]
for _ in range(10**2):
    for i in range(12):
        m[i+1]+= m[i]
    A123095_list.append(m[-1]) # Chai Wah Wu, Nov 05 2014

[(bernoulli_polynomial(n+1, 12) - bernoulli(12))/12  for n in (0..30)] # G. C. Greubel, Jul 21 2021

a(n) = n*A023002(n) - Sum_{i=0..n-1} A023002(i). - Bruno Berselli, Apr 27 2010
a(n) = n^2*(n+1)^2*(2*n^8 +8*n^7 +4*n^6 -16*n^5 -5*n^4 +26*n^3 -3*n^2 -20*n +10)/24. - Bruno Berselli, Oct 03 2010
G.f.: x*(x^10 +2036*x^9 +152637*x^8 +2203488*x^7 +9738114*x^6 +15724248*x^5 +9738114*x^4 +2203488*x^3 +152637*x^2 +2036*x +1)/(1-x)^13. - Colin Barker, May 27 2012
a(n) = (-1)*Sum_{j=1..11} j*Stirling1(n+1,n+1-j)*Stirling2(n+11-j,n). - Mircea Merca, Jan 25 2014
a(n) = 1728*A006542(n+2)^2 + 216*A288876(n-2) + 96*A006542(n+2) + A000537(n). - Yasser Arath Chavez Reyes, May 25 2024

A181134 Sum of 13th powers: a(n) = Sum_{j=0..n} j^13.

Original entry on oeis.org

0, 1, 8193, 1602516, 68711380, 1289414505, 14350108521, 111239118928, 660994932816, 3202860761145, 13202860761145, 47725572905076, 154718778284148, 457593884876401, 1251308658130545, 3197503726489920

Offset: 0

Bruno Berselli, Oct 05 2010 - Oct 18 2010

This form of recurrence is a general property of the array in A103438 (sums of the first n-th powers).

Bruno Berselli, Table of n, a(n) for n = 0..10000.
Bruno Berselli, A description of the recursive method in Comments lines: website Matem@ticamente (in Italian), 2008.
Index entries for linear recurrences with constant coefficients, signature (15,-105,455,-1365,3003,-5005,6435,-6435,5005,-3003,1365,-455,105,-15,1).

Cf. A010801.
Sequences of the form Sum_{j=0..n} j^m : A000217 (m=1), A000330 (m=2), A000537 (m=3), A000538 (m=4), A000539 (m=5), A000540 (m=6), A000541 (m=7), A000542 (m=8), A007487 (m=9), A023002 (m=10), A123095 (m=11), A123094 (m=12), A181134 (m=13).
Cf. A215083, A253712.

[(&+[j^13: j in [0..n]]): n in [0..30]]; // G. C. Greubel, Jul 21 2021

A181134 := proc(n) (bernoulli(14,n+1) - bernoulli(14))/14 ; end proc: seq(A181134(n), n=0..10); # R. J. Mathar, Oct 14 2010

Accumulate[Range[0,20]^13] (* Harvey P. Dale, Oct 30 2017 *)

A181134_list, m = [0], [6227020800, -37362124800, 97037740800, -142702560000, 130456085760, -76592355840, 28805736960, -6711344640, 901020120, -60780720, 1569750, -8190, 1, 0 , 0]
for _ in range(10**2):
    for i in range(14):
        m[i+1]+= m[i]
    A181134_list.append(m[-1]) # Chai Wah Wu, Nov 06 2014

[(bernoulli_polynomial(n+1, 14) - bernoulli(14))/14  for n in (0..30)] # G. C. Greubel, Jul 21 2021

For n>0, a(n) = n*A123094(n) - Sum_{i=0..n-1} A123094(i), where Sum_{i=0..n-1} A123094(i) = A253712(n-1) = (n-1)*n^2*(n+1)*(30*n^10 - 425*n^8 + 2578*n^6 - 8147*n^4 + 12874*n^2 - 7601)/5460.
a(n) = a(-n-1) = (n*(n + 1))^2*(30*n^10 + 150*n^9 + 125*n^8 - 400*n^7 - 326*n^6 + 1052*n^5 + 367*n^4 - 1786*n^3 + 202*n^2 + 1382*n - 691)/420.
G.f.: see comment of Vladeta Jovovic in A000538.
a(n) = -Sum_{j=1..13} j*Stirling1(n+1,n+1-j)*Stirling2(n+13-j,n). - Mircea Merca, Jan 25 2014

A254643 Third partial sums of ninth powers (A001017).

Original entry on oeis.org

1, 515, 21225, 324275, 2862790, 17714466, 85232910, 339635850, 1168343775, 3571356685, 9906622271, 25333920885, 60457751900, 135939162100, 290221510860, 592024274916, 1159935330765, 2192313968775, 4011847886725, 7130537084615

Offset: 1

Luciano Ancora, Feb 05 2015

			First differences:   1, 511, 19171, 242461, 1690981, ... (A022525)
------------------------------------------------------------------------
The ninth powers:    1, 512, 19683, 262144, 1953125, ... (A001017)
------------------------------------------------------------------------
First partial sums:  1, 513, 20196, 282340, 2235465, ... (A007487)
Second partial sums: 1, 514, 20710, 303050, 2538515, ... (A253637)
Third partial sums:  1, 515, 21225, 324275, 2862790, ... (this sequence)

Luciano Ancora, Table of n, a(n) for n = 1..1000
Luciano Ancora, Partial sums of m-th powers with Faulhaber polynomials
Luciano Ancora, Pascal’s triangle and recurrence relations for partial sums of m-th powers
Index entries for linear recurrences with constant coefficients, signature (13,-78,286,-715,1287,-1716,1716,-1287,715,-286,78,-13,1).

Cf. A001017, A007487, A022525, A253637.

List([1..30], n-> Binomial(n+3,4)*(2*n^8 +24*n^7 +98*n^6 +126*n^5 -97*n^4 -203*n^3 +127*n^2 +84*n -50)/110); # G. C. Greubel, Aug 28 2019

[Binomial(n+3,4)*(2*n^8 +24*n^7 +98*n^6 +126*n^5 -97*n^4 -203*n^3 +127*n^2 +84*n -50)/110: n in [1..30]]; // G. C. Greubel, Aug 28 2019

seq(binomial(n+3,4)*(2*n^8 +24*n^7 +98*n^6 +126*n^5 -97*n^4 -203*n^3 +127*n^2 +84*n -50)/110, n=1..30); # G. C. Greubel, Aug 28 2019

Table[n(1+n)(2+n)(3+n)(-50 +84n +127n^2 -204n^3 -97n^4 +126n^5 +98n^6 +24n^7 +2n^8)/2640, {n, 20}] (* or *)
CoefficientList[Series[(1 +502x +14608x^2 +88234x^3 +156190x^4 +88234x^5 +14608x^6 +502x^7 +x^8)/(1-x)^13, {x, 0, 19}], x] (* Ancora *)
Accumulate[Accumulate[Accumulate[Range[10]^9]]] (* Alonso del Arte, Feb 09 2015 *)

vector(30, n, m=n+3; binomial(m,4)*(2*(n*m)^4 -10*(n*m)^3 +11*(n*m)^2 +28*(n*m) -50)/110) \\ G. C. Greubel, Aug 28 2019

[binomial(n+3,4)*(2*n^8 +24*n^7 +98*n^6 +126*n^5 -97*n^4 -203*n^3 +127*n^2 +84*n -50)/110 for n in (1..30)] # G. C. Greubel, Aug 28 2019

G.f.: x*(1 +502*x +14608*x^2 +88234*x^3 +156190*x^4 +88234*x^5 +14608*x^6 +502*x^7 +x^8)/(1-x)^13.
a(n) = n*(1+n)*(2+n)*(3+n)*(-50 +84*n +127*n^2 -204*n^3 -97*n^4 +126*n^5 +98*n^6 +24*n^7 +2*n^8)/2640.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + n^9.

Edited by Alonso del Arte and Bruno Berselli, Feb 10 2015

A255179 Second differences of ninth powers (A001017).

Original entry on oeis.org

1, 510, 18660, 223290, 1448520, 6433590, 22151340, 63588210, 159338640, 359376750, 745368180, 1443884970, 2642886360, 4611828390, 7725765180, 12493804770, 19592282400, 29903014110, 44556993540, 64983894810, 92967744360, 130709124630, 180894272460

Offset: 0

Luciano Ancora, Feb 21 2015

			Second differences:  1, 510, 18660, 223290, 1448520, ... (this sequence)
First differences:   1, 511, 19171, 242461, 1690981, ... (A022525)
------------------------------------------------------------------------
The ninth powers:    1, 512, 19683, 262144, 1953125, ... (A001017)
------------------------------------------------------------------------
First partial sums:  1, 513, 20196, 282340, 2235465, ... (A007487)
Second partial sums: 1, 514, 20710, 303050, 2538515, ... (A253637)
Third partial sums:  1, 515, 21225, 324275, 2862790, ... (A254643)

Luciano Ancora, Table of n, a(n) for n = 0..1000
Luciano Ancora, Sums of powers of positive integers and their recurrence relations, section 0.5.
Index entries for linear recurrences with constant coefficients, signature (8,-28,56,-70,56,-28,8,-1).

Cf. A001017, A007487, A022525, A253637, A254643, A255177, A255178.

[1] cat [6*n*(3+28*n^2+42*n^4+12*n^6): n in [1..30]]; // Vincenzo Librandi, Mar 12 2015

Join[{1}, Table[6 n (3 + 28 n^2 + 42 n^4 + 12 n^6), {n, 1, 30}]]
Join[{1},Differences[Range[0,30]^9,2]] (* or *) LinearRecurrence[{8,-28,56,-70,56,-28,8,-1},{1,510,18660,223290,1448520,6433590,22151340,63588210,159338640},30] (* Harvey P. Dale, Jan 26 2019 *)

G.f.: (1 + 502*x + 14608*x^2 + 88234*x^3 + 156190*x^4 + 88234*x^5 + 14608*x^6 + 502*x^7 + x^8)/(1 - x)^8.

a(n) = 6*n*(3 + 28*n^2 + 42*n^4 + 12*n^6) for n>0, a(0)=1.

Corrected g.f. by Bruno Berselli, Feb 25 2015

Offset changed by Bruno Berselli, Mar 20 2015

A081175 Numbers of the form Sum_{i=1..k} i^j, j >= 1, k >= 1.

Original entry on oeis.org

1, 3, 5, 6, 9, 10, 14, 15, 17, 21, 28, 30, 33, 36, 45, 55, 65, 66, 78, 91, 98, 100, 105, 120, 129, 136, 140, 153, 171, 190, 204, 210, 225, 231, 253, 257, 276, 285, 300, 325, 351, 354, 378, 385, 406, 435, 441, 465, 496, 506, 513, 528, 561, 595, 630, 650, 666, 703

Offset: 1

N. J. A. Sloane, Apr 18 2003

Union of sums of k-th powers, for k >= 1.

			30 is in the set because 30 = 1^2 + 2^2 + 3^2 + 4^2 (j=2, k=4).

Robert Israel, Table of n, a(n) for n = 1..10000
Michael Penn, an excruciatingly deep dive into the power sum., YouTube video, 2022.

Cf. A000217, A000330, A000537, A000538, A000539, A000540, A000541, A000542, A007487, A023002.

For primes in this sequence see A164307.

N:= 1000: # to get all terms <= N
A:=select(`<=`,{1, seq(seq(sum(i^k,i=1..m), m=2..floor((N*(k+1))^(1/(k+1)))),k = 1 ..ilog2(N-1))},N):
sort(convert(A,list)); # Robert Israel, Jan 26 2015

Take[ Union[ Flatten[ Table[ Sum[ i^j, {i, 1, n}], {j, 1, 9}, {n, 1, 40}]]], 60]

Corrected and extended by Robert G. Wilson v, May 08 2003

A255183 Third differences of ninth powers (A001017).

Original entry on oeis.org

1, 509, 18150, 204630, 1225230, 4985070, 15717750, 41436870, 95750430, 200038110, 385991430, 698516790, 1199001390, 1968942030, 3113936790, 4768039590, 7098477630, 10310731710, 14653979430, 20426901270

Offset: 0

Luciano Ancora, Mar 18 2015

			Third differences:   1, 509, 18150, 204630, 1225230, ...  (this sequence)
Second differences:  1, 510, 18660, 223290, 1448520, ...  (A255179)
First differences:   1, 511, 19171, 242461, 1690981, ...  (A022525)
---------------------------------------------------------------------
The ninth powers:    1, 512, 19683, 262144, 1953125, ...  (A001017)
---------------------------------------------------------------------

Luciano Ancora, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (7,-21,35,-35,21,-7,1).

Cf. A001017, A007487, A022525, A253637, A255179, A255182.

[1,509] cat [6*(84*n^6-252*n^5+630*n^4-840*n^3+756*n^2-378*n+85): n in [2..30]]; // Vincenzo Librandi, Mar 18 2015

Join[{1, 509}, Table[6 (84 n^6 - 252 n^5 + 630 n^4 - 840 n^3 + 756 n^2 - 378 n + 85), {n, 2, 30}]]
Join[{1,509},Differences[Range[0,20]^9,3]] (* Harvey P. Dale, Apr 24 2015 *)

G.f.: (1 + 502*x + 14608*x^2 + 88234*x^3 + 156190*x^4 + 88234*x^5 + 14608*x^6 + 502*x^7 + x^8)/(1 - x)^7.

a(n) = 6*(84*n^6 - 252*n^5 + 630*n^4 - 840*n^3 + 756*n^2 - 378*n + 85) for n>1, a(0)=1, a(1)=509.

Edited by Bruno Berselli, Mar 20 2015

cp's OEIS Frontend

A000538 Sum of fourth powers: 0^4 + 1^4 + ... + n^4.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

Maxima

PARI

PARI

Python

Python

Sage

Formula

Extensions

A103438 Square array T(m,n) read by antidiagonals: Sum_{k=1..n} k^m.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Python

SageMath

Formula

A023002 Sum of 10th powers.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Python

Sage

Formula

A123094 Sum of first n 12th powers.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Python

Sage

Formula

A123095 Sum of first n 11th powers.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Python

Sage

Formula