cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 13 results. Next

A101089 Second partial sums of fourth powers (A000583).

Original entry on oeis.org

1, 18, 116, 470, 1449, 3724, 8400, 17172, 32505, 57838, 97812, 158522, 247793, 375480, 553792, 797640, 1125009, 1557354, 2120020, 2842686, 3759833, 4911236, 6342480, 8105500, 10259145, 12869766, 16011828, 19768546, 24232545, 29506544
Offset: 1

Views

Author

Cecilia Rossiter, Dec 14 2004

Keywords

Comments

a(n) is the n-th antidiagonal sum of the convolution array A213553. - Clark Kimberling, Jun 17 2012
a(n-1)/n^5 is the "retention" of water on a 3 X 3 random surface of n levels - see Knecht et al., 2012, Schrenk et al., 2014. - Robert M. Ziff, Mar 08 2014
The general formula for the second partial sums of m-th powers is: b(n,m) = (n+1)*F(m) - F(m+1), where F(m) is the m-th Faulhaber’s polynomial. - Luciano Ancora, Jan 26 2015

Examples

			a(7) = 8400 = 1*(8-1)^4 + 2*(8-2)^4 + 3*(8-3)^4 + 4*(8-4)^4 + 5*(8-5)^4 + 6*(8-6)^4 + 7*(8-7)^4. - _Bruno Berselli_, Jan 31 2014
		

Crossrefs

Partial sums of A000538.

Programs

  • GAP
    List([1..40], n-> (n+1)^2*(2*(n+1)^4-5*(n+1)^2+3)/60); # G. C. Greubel, Jul 31 2019
  • Magma
    [(1/60)*n*(n+1)^2*(n+2)*(2*n*(n+2)-1): n in [1..40]]; // Vincenzo Librandi, Mar 24 2014
    
  • Maple
    f:=n->(2*n^6-5*n^4+3*n^2)/60;
    [seq(f(n),n=0..50)]; # N. J. A. Sloane, Mar 23 2014
  • Mathematica
    a[n_] := n(n+1)^2(n+2)(2n(n+2) -1)/60; Table[a[n], {n, 40}]
    CoefficientList[Series[(1+x)*(1+10*x+x^2)/(1-x)^7, {x,0,40}], x] (* Vincenzo Librandi, Mar 24 2014 *)
    Nest[Accumulate[#]&,Range[30]^4,2] (* Harvey P. Dale, Aug 13 2024 *)
  • PARI
    a(n)=n*(n+1)^2*(n+2)*(2*n*(n+2)-1)/60 \\ Charles R Greathouse IV, Mar 18 2014
    
  • Sage
    [n*(n+1)^2*(n+2)*(2*n*(n+2)-1)/60 for n in range(1,40)] # Danny Rorabaugh, Apr 20 2015
    

Formula

a(n) = (1/60)*n*(n+1)^2*(n+2)*(2*n*(n+2)-1).
G.f.: x*(1+x)*(1+10*x+x^2)/(1-x)^7. - Colin Barker, Apr 04 2012
a(n) = Sum_{i=1..n} i*(n+1-i)^4, by the definition. - Bruno Berselli, Jan 31 2014
a(n) = 2*a(n-1) - a(n-2) + n^4. - Luciano Ancora, Jan 08 2015
Sum_{n>=1} 1/a(n) = 85/3 + 10*Pi^2/3 - 20*sqrt(2/3)*Pi*cot(sqrt(3/2)*Pi). - Amiram Eldar, Jan 26 2022
a(n) = (1/2)*Sum_{1 <= i, j <= n+1} (i - j)^4 - Peter Bala, Jun 11 2024

Extensions

Edited by Ralf Stephan, Dec 16 2004

A254647 Fourth partial sums of eighth powers (A001016).

Original entry on oeis.org

1, 260, 7595, 94360, 723534, 4037712, 17944290, 67127880, 219319815, 642251428, 1718012933, 4258676240, 9892043980, 21721707840, 45414150132, 90930820464, 175208925885, 326205634020, 588861675535, 1033717781096, 1769137540730, 2958360418000, 4842936861750, 7774492635000
Offset: 1

Views

Author

Luciano Ancora, Feb 05 2015

Keywords

Examples

			The eighth powers:   1, 256, 6561, 65536, 390625, ... (A001016)
First partial sums:  1, 257, 6818, 72354, 462979, ... (A000542)
Second partial sums: 1, 258, 7076, 79430, 542409, ... (A253636)
Third partial sums:  1, 259, 7335, 86765, 629174, ... (A254642)
Fourth partial sums: 1, 260, 7595, 94360, 723534, ... (this sequence)
		

Crossrefs

Programs

  • GAP
    List([1..30], n-> Binomial(n+4,5)*(n+2)*((n+2)^2-3)*(2*(n+2)^4 -28*(n+2)^2 +101)/198); # G. C. Greubel, Aug 28 2019
  • Magma
    [Binomial(n+4,5)*(n+2)*((n+2)^2-3)*(2*(n+2)^4 -28*(n+2)^2 +101)/198: n in [1..30]]; // G. C. Greubel, Aug 28 2019
    
  • Maple
    seq(binomial(n+4,5)*(n+2)*((n+2)^2-3)*(2*(n+2)^4 -28*(n+2)^2 +101)/198, n=1..30); # G. C. Greubel, Aug 28 2019
  • Mathematica
    Table[n(1+n)(2+n)^2(3+n)(4+n)(1+4n+n^2)(21 -48n +20n^2 +16n^3 +2n^4 )/23760, {n,20}] (* or *)
    Accumulate[Accumulate[Accumulate[Accumulate[Range[20]^8]]]] (* or *)
    CoefficientList[Series[(1 +247x +4293x^2 +15619x^3 +15619x^4 +4293x^5 + 247x^6 +x^7)/(1-x)^13, {x,0,19}], x]
  • PARI
    a(n)=n*(1+n)*(2+n)^2*(3+n)*(4+n)*(1+4*n+n^2)*(21-48*n+20*n^2 +16*n^3+2*n^4)/23760 \\ Charles R Greathouse IV, Sep 08 2015
    
  • PARI
    vector(30, n, m=n+2; binomial(m+2,5)*m*(m^2-3)*(2*m^4-28*m^2 +101)/198)
    
  • Sage
    [binomial(n+4,5)*(n+2)*((n+2)^2-3)*(2*(n+2)^4 -28*(n+2)^2 +101)/198 for n in (1..30)] # G. C. Greubel, Aug 28 2019
    

Formula

G.f.: x*(1 +247*x +4293*x^2 +15619*x^3 +15619*x^4 +4293*x^5 +247*x^6 +x^7)/(1-x)^13.
a(n) = n*(1+n)*(2+n)^2*(3+n)*(4+n)*(1 +4*n +n^2)*(21 -48*n +20*n^2 + 16*n^3 +2*n^4)/23760.
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) + n^8.

A101092 Second partial sums of fifth powers (A000584).

Original entry on oeis.org

1, 34, 310, 1610, 6035, 18236, 47244, 109020, 229845, 450670, 832546, 1463254, 2465255, 4005080, 6304280, 9652056, 14419689, 21076890, 30210190, 42543490, 58960891, 80531924, 108539300, 144509300, 190244925, 247861926, 319827834, 409004110, 518691535
Offset: 1

Views

Author

Cecilia Rossiter, Dec 14 2004

Keywords

Crossrefs

Programs

  • Magma
    [(n*(1+n)*(2+n)*(-1+n*(2+n))*(1+2*n*(2+n)))/84: n in [1..40]]; // Vincenzo Librandi, Mar 24 2014
    
  • Maple
    f:=n->(2*n^7-7*n^5+7*n^3-2*n)/84;
    [seq(f(n),n=0..50)];  # N. J. A. Sloane, Mar 23 2014
  • Mathematica
    CoefficientList[Series[(1 + 26 x + 66 x^2 + 26 x^3 + x^4)/, {x, 0, 40}], x] (* Vincenzo Librandi, Mar 24 2014 *)
    Nest[Accumulate,Range[30]^5,2] (* or *) LinearRecurrence[{8,-28,56,-70,56,-28,8,-1},{1,34,310,1610,6035,18236,47244,109020},30] (* Harvey P. Dale, Aug 22 2022 *)
  • PARI
    a(n)=n*(n+1)*(n+2)*(n*(n+2)-1)*(2*n*(2 + n)+1)/84 \\ Charles R Greathouse IV, Apr 21 2015
    
  • Python
    def A101092(n): return n*(n*(n*(n*(n*(n*(n+7<<1)+35)+35)+7)-7)-2)//84 # Chai Wah Wu, Oct 03 2024
  • Sage
    [n*(1+n)*(2+n)*(n*(2+n)-1)*(1+2*n*(2+n))/84 for n in range(1,30)] # Danny Rorabaugh, Apr 21 2015
    

Formula

a(n) = (n*(1 + n)*(2 + n)*(-1 + n*(2 + n))*(1 + 2*n*(2 + n)))/84.
G.f.: x*(1 + 26*x + 66*x^2 + 26*x^3 + x^4)/(1-x)^8. - Colin Barker, Apr 16 2012
a(n) = Sum_{i=1..n} i*(n+1-i)^5, by the definition. - Bruno Berselli, Jan 31 2014
a(n) = 2*a(n-1) - a(n-2) + n^5. - Luciano Ancora, Jan 08 2015
E.g.f.: exp(x)*x*(84 + 1344*x + 2954*x^2 + 1995*x^3 + 525*x^4 + 56*x^5 + 2*x^6)/84. - Stefano Spezia, May 04 2024

Extensions

Edited by Ralf Stephan, Dec 16 2004

A101093 Second partial sums of sixth powers (A001014).

Original entry on oeis.org

1, 66, 860, 5750, 26265, 93436, 278256, 725220, 1703625, 3682030, 7431996, 14167946, 25730705, 44823000, 75305920, 122566056, 193963761, 299373690, 451829500, 668285310, 970507241, 1386109076, 1949746800, 2704487500
Offset: 1

Views

Author

Cecilia Rossiter (cecilia(AT)noticingnumbers.net), Dec 15 2004

Keywords

Crossrefs

Cf. A000540.

Programs

  • GAP
    List([2..30], n-> n^2*(3*n^6 -14*n^4 +21*n^2 -10)/168); # G. C. Greubel, Aug 28 2019
  • Magma
    [n*(1+n)^2*(2+n)*(-1+n*(2+n))*(-2+3*n*(2+n))/168: n in [1..40]]; // Vincenzo Librandi, Mar 24 2014
    
  • Maple
    f:=n->(3*n^8-14*n^6+21*n^4-10*n^2)/168;
    [seq(f(n),n=0..50)];  # N. J. A. Sloane, Mar 23 2014
  • Mathematica
    CoefficientList[Series[(x+1)(x^4+56x^3+246x^2+56x+1)/(1-x)^9, {x,0,40}], x] (* Vincenzo Librandi, Mar 24 2014 *)
    Nest[Accumulate,Range[30]^6,2] (* or *) LinearRecurrence[{9,-36,84,-126, 126,-84,36,-9,1},{1,66,860,5750,26265,93436,278256,725220,1703625},30] (* Harvey P. Dale, Jun 05 2019 *)
  • PARI
    vector(30, n, m=n+1; m^2*(3*m^6 -14*m^4 +21*m^2 -10)/168) \\ G. C. Greubel, Aug 28 2019
    
  • Sage
    [n^2*(3*n^6 -14*n^4 +21*n^2 -10)/168 for n in (2..30)] # G. C. Greubel, Aug 28 2019
    

Formula

a(n) = n*(1 + n)^2*(2 + n)*(-1 + n*(2 + n))*(-2 + 3*n*(2 + n))/168.
G.f.: x*(1+x)*(1 + 56*x + 246*x^2 + 56*x^3 + x^4)/(1-x)^9. - Colin Barker, Dec 18 2012
a(n) = Sum_{i=1..n} i*(n+1-i)^6, by the definition. - Bruno Berselli, Jan 31 2014
a(n) = 2*a(n-1) - a(n-2) + n^6. - Luciano Ancora, Jan 08 2015

Extensions

Edited by Ralf Stephan, Dec 16 2004

A101095 Fourth difference of fifth powers (A000584).

Original entry on oeis.org

1, 28, 121, 240, 360, 480, 600, 720, 840, 960, 1080, 1200, 1320, 1440, 1560, 1680, 1800, 1920, 2040, 2160, 2280, 2400, 2520, 2640, 2760, 2880, 3000, 3120, 3240, 3360, 3480, 3600, 3720, 3840, 3960, 4080, 4200, 4320, 4440, 4560, 4680, 4800, 4920, 5040, 5160, 5280
Offset: 1

Views

Author

Cecilia Rossiter, Dec 15 2004

Keywords

Comments

Original Name: Shells (nexus numbers) of shells of shells of shells of the power of 5.
The (Worpitzky/Euler/Pascal Cube) "MagicNKZ" algorithm is: MagicNKZ(n,k,z) = Sum_{j=0..k+1} (-1)^j*binomial(n + 1 - z, j)*(k - j + 1)^n, with k>=0, n>=1, z>=0. MagicNKZ is used to generate the n-th accumulation sequence of the z-th row of the Euler Triangle (A008292). For example, MagicNKZ(3,k,0) is the 3rd row of the Euler Triangle (followed by zeros) and MagicNKZ(10,k,1) is the partial sums of the 10th row of the Euler Triangle. This sequence is MagicNKZ(5,k-1,2).

Crossrefs

Fourth differences of A000584, third differences of A022521, second differences of A101098, and first differences of A101096.
For other sequences based upon MagicNKZ(n,k,z):
...... | n = 1 | n = 2 | n = 3 | n = 4 | n = 5 | n = 6 | n = 7 | n = 8
--------------------------------------------------------------------------------------
z = 0 | A000007 | A019590 | ....... MagicNKZ(n,k,0) = T(n,k+1) from A008292 .......
z = 1 | A000012 | A040000 | A101101 | A101104 | A101100 | ....... | ....... | .......
z = 2 | A000027 | A005408 | A008458 | A101103 | thisSeq | ....... | ....... | .......
z = 3 | A000217 | A000290 | A003215 | A005914 | A101096 | ....... | ....... | .......
z = 4 | A000292 | A000330 | A000578 | A005917 | A101098 | ....... | ....... | .......
z = 5 | A000332 | A002415 | A000537 | A000583 | A022521 | ....... | A255181 | .......
z = 12 | A001288 | A057788 | ....... | A254870 | A254471 | A254683 | A254646 | A254642
z = 13 | A010965 | ....... | ....... | ....... | A254871 | A254472 | A254684 | A254647
z = 14 | A010966 | ....... | ....... | ....... | ....... | A254872 | ....... | .......
--------------------------------------------------------------------------------------
Cf. A047969.

Programs

  • Magma
    I:=[1,28,121,240,360]; [n le 5 select I[n] else 2*Self(n-1)-Self(n-2): n in [1..50]]; // Vincenzo Librandi, May 07 2015
    
  • Mathematica
    MagicNKZ=Sum[(-1)^j*Binomial[n+1-z, j]*(k-j+1)^n, {j, 0, k+1}];Table[MagicNKZ, {n, 5, 5}, {z, 2, 2}, {k, 0, 34}]
    CoefficientList[Series[(1 + 26 x + 66 x^2 + 26 x^3 + x^4)/(1 - x)^2, {x, 0, 50}], x] (* Vincenzo Librandi, May 07 2015 *)
    Join[{1,28,121,240},Differences[Range[50]^5,4]] (* or *) LinearRecurrence[{2,-1},{1,28,121,240,360},50] (* Harvey P. Dale, Jun 11 2016 *)
  • PARI
    a(n)=if(n>3, 120*n-240, 33*n^2-72*n+40) \\ Charles R Greathouse IV, Oct 11 2015
  • Sage
    [1,28,121]+[120*(k-2) for k in range(4,36)] # Danny Rorabaugh, Apr 23 2015
    

Formula

a(k+1) = Sum_{j=0..k+1} (-1)^j*binomial(n + 1 - z, j)*(k - j + 1)^n; n = 5, z = 2.
For k>3, a(k) = Sum_{j=0..4} (-1)^j*binomial(4, j)*(k - j)^5 = 120*(k - 2).
a(n) = 2*a(n-1) - a(n-2), n>5. G.f.: x*(1+26*x+66*x^2+26*x^3+x^4) / (1-x)^2. - Colin Barker, Mar 01 2012

Extensions

MagicNKZ material edited, Crossrefs table added, SeriesAtLevelR material removed by Danny Rorabaugh, Apr 23 2015
Name changed and keyword 'uned' removed by Danny Rorabaugh, May 06 2015

A254642 Third partial sums of eighth powers (A001016).

Original entry on oeis.org

1, 259, 7335, 86765, 629174, 3314178, 13906578, 49183590, 152191935, 422931613, 1075761505, 2540663307, 5633367740, 11829663860, 23692442292, 45516670332, 84278105421, 150996708135, 262656041515, 444856105561, 735419759634, 1189222877270
Offset: 1

Views

Author

Luciano Ancora, Feb 05 2015

Keywords

Examples

			First differences:   1, 255, 6305, 58975, 325089, ...(A022524)
--------------------------------------------------------------------
The eighth powers:   1, 256, 6561, 65536, 390625, ...(A001016)
--------------------------------------------------------------------
First partial sums:  1, 257, 6818, 72354, 462979, ...(A000542)
Second partial sums: 1, 258, 7076, 79430, 542409, ...(A253636)
Third partial sums:  1, 259, 7335, 86765, 629174, ...(this sequence)
		

Crossrefs

Programs

  • Mathematica
    Table[n (1 + n) (2 + n) (3 + n) (3 + 2 n) (1 + 36 n - 69 n^2 + 45 n^4 + 18 n^5 + 2 n^6)/3960, {n, 22}]
    Accumulate[Accumulate[Accumulate[Range[22]^8]]]
    CoefficientList[Series[(1 + 247 x + 4293 x^2 + 15619 x^3 + 15619 x^4 + 4293 x^5 + 247 x^6 + x^7)/(- 1 + x)^12, {x, 0, 22}], x]
  • PARI
    a(n)=n*(1+n)*(2+n)*(3+n)*(3+2*n)*(1+36*n-69*n^2+45*n^4+18*n^5+2*n^6)/3960 \\ Charles R Greathouse IV, Oct 07 2015

Formula

G.f.: (x + 247*x^2 + 4293*x^3 + 15619*x^4 + 15619*x^5 + 4293*x^6 + 247*x^7 + x^8)/(- 1 + x)^12.
a(n) = n*(1 + n)*(2 + n)*(3 + n)*(3 + 2*n)*(1 + 36*n - 69*n^2 + 45*n^4 + 18*n^5 + 2*n^6)/3960.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + n^8.

A250212 Second partial sums of seventh powers (A001015).

Original entry on oeis.org

1, 130, 2446, 21146, 117971, 494732, 1695036, 4992492, 13072917, 31153342, 68720938, 142120342, 278268263, 519829688, 932250488, 1613106744, 2704301673, 4407716634, 7005003334, 10882290034, 16560665275, 24733398404, 36310956980, 52474986980, 74742532605, 105041888406
Offset: 1

Views

Author

Luciano Ancora, Jan 18 2015

Keywords

Comments

The general formula for the second partial sums of m-th powers is: b(n,m) = (n+1)*F(m) - F(m+1), where F(m) is the m-th Faulhaber’s polynomial.

Crossrefs

Cf. A239094 (same sequence, shifted by 1).

Programs

  • GAP
    List([1..30], n-> Binomial(n+2, 3)*(5*(n+1)^6 -25*(n+1)^4 +38*(n+ 1)^2 -12)/60); # G. C. Greubel, Aug 28 2019
  • Magma
    [(n*(n + 1)*(n + 2)*(5*n^6 + 30*n^5 + 50*n^4 -37*n^2 + 6*n + 6) / 360): n in [1..30]]; // Vincenzo Librandi, Jan 22 2015
    
  • Maple
    seq(binomial(n+2, 3)*(5*(n+1)^6 -25*(n+1)^4 +38*(n+1)^2 -12)/60, n=1..30); # G. C. Greubel, Aug 28 2019
  • Mathematica
    Accumulate[Accumulate[Range[25]^7]] (* Robert G. Wilson v, Jan 21 2015 *)
    Table[(n(n+1)(n+2)(5n^6+30n^5+50n^4-37n^2+6n+6)/360), {n, 30}] (* Vincenzo Librandi, Jan 22 2015 *)
    RecurrenceTable[{a[n]==2a[n-1]-a[n-2]+n^7, a[1]==1,a[2]==130}, a, {n,30}] (* Bruno Berselli, Jan 22 2015 *)
    LinearRecurrence[{10,-45,120,-210,252,-210,120,-45,10,-1},{1,130,2446,21146,117971,494732,1695036,4992492,13072917,31153342},30] (* Harvey P. Dale, Jan 19 2020 *)
  • PARI
    vector(50, n, n*(n+1)*(n+2)*(5*n^6 + 30*n^5 + 50*n^4 - 37*n^2 + 6*n + 6)/360) \\ Michel Marcus, Jan 21 2015
    
  • Sage
    [binomial(n+2, 3)*(5*(n+1)^6 -25*(n+1)^4 +38*(n+1)^2 -12)/60 for n in (1..30)] # G. C. Greubel, Aug 28 2019
    

Formula

a(n) = n*(n+1)*(n+2)*(5*n^6 + 30*n^5 + 50*n^4 - 37*n^2 + 6*n + 6)/360.
a(n) = 2*a(n-1) - a(n-2) + n^7.
G.f.: x*(1 +120*x +1191*x^2 +2416*x^3 +1191*x^4 +120*x^5 +x^6)/(1-x)^10. - Georg Fischer, May 24 2019
a(n) = A239094(n+1). - Danny Rorabaugh, Apr 22 2015

A255178 Second differences of eighth powers (A001016).

Original entry on oeis.org

1, 254, 6050, 52670, 266114, 963902, 2796194, 6927230, 15257090, 30683774, 57405602, 101263934, 170126210, 274309310, 427043234, 644975102, 948713474, 1363412990, 1919399330, 2652834494, 3606422402, 4830154814, 6382097570, 8329217150, 10748247554
Offset: 0

Views

Author

Luciano Ancora, Feb 21 2015

Keywords

Examples

			Second differences:  1, 254, 6050, 52670, 266114, ... (this sequence)
First differences:   1, 255, 6305, 58975, 325089, ... (A022524)
----------------------------------------------------------------------
The eighth powers:   1, 256, 6561, 65536, 390625, ... (A001016)
----------------------------------------------------------------------
First partial sums:  1, 257, 6818, 72354, 462979, ... (A000542)
Second partial sums: 1, 258, 7076, 79430, 542409, ... (A253636)
Third partial sums:  1, 259, 7335, 86765, 629174, ... (A254642)
Fourth partial sums: 1, 260, 7595, 94360, 723534, ... (A254647)
		

Crossrefs

Programs

  • Magma
    [n eq 0 select 1 else 2*(28*n^6+70*n^4+28*n^2+1): n in [0..30]]; // Vincenzo Librandi, Mar 12 2015
  • Mathematica
    Join[{1}, Table[2 (28 n^6 + 70 n^4 + 28 n^2 + 1), {n, 1, 30}]]
    Join[{1},Differences[Range[0,30]^8,2]] (* Harvey P. Dale, Aug 26 2024 *)

Formula

G.f.: (1 + x)*(1 + 246*x + 4047*x^2 + 11572*x^3 + 4047*x^4 + 246*x^5 + x^6)/(1 - x)^7.
a(n) = 2*(28*n^6 + 70*n^4 + 28*n^2 + 1) for n>0, a(0)=1.

Extensions

Edited by Bruno Berselli, Mar 19 2015

A253637 Second partial sums of ninth powers (A001017).

Original entry on oeis.org

1, 514, 20710, 303050, 2538515, 14851676, 67518444, 254402940, 828707925, 2403012910, 6335265586, 15427298614, 35123831015, 75481410200, 154282348760, 301802764056, 567911055849, 1032378638010, 1819533917950, 3118689197890, 5212124524411, 8512829068724, 13614686274500, 21358351020500, 32916713032125, 49904578722726
Offset: 1

Views

Author

Luciano Ancora, Jan 07 2015

Keywords

Comments

The formula for the second partial sums of m-th powers is: b(n,m) = (n+1)*F(m) - F(m+1), where F(m) are the m-th Faulhaber's formulas.

Crossrefs

Cf. A001017.

Programs

  • GAP
    List([1..30], n-> n*(n+1)*(n+2)*(n^2+n-1)*(n^2+3*n+1)*(6*n^4+24*n^3 +5*n^2-38*n+ 25)/660 ); # G. C. Greubel, Aug 28 2019
  • Magma
    [n*(n+1)*(n+2)*(n^2+n-1)*(n^2+3*n+1)*(6*n^4+24*n^3+5*n^2-38*n+ 25)/660: n in [1..30]]; // G. C. Greubel, Aug 28 2019
    
  • Maple
    seq(n*(n+1)*(n+2)*(n^2+n-1)*(n^2+3*n+1)*(6*n^4+24*n^3+5*n^2-38*n+ 25)/660, n=1..30); # G. C. Greubel, Aug 28 2019
  • Mathematica
    CoefficientList[Series[(1 +502x +14608x^2 +88234x^3 +156190x^4 +88234x^5 +14608x^6 +502x^7 +x^8)/(1-x)^12, {x, 0, 30}], x] (* Vincenzo Librandi, Jan 19 2015 *)
    Nest[Accumulate,Range[30]^9,2] (* Harvey P. Dale, Apr 18 2021 *)
  • PARI
    a(n) = (6*n^11 + 66*n^10 + 275*n^9 + 495*n^8 + 198*n^7 - 462*n^6 - 330*n^5 + 330*n^4 + 231*n^3 - 99*n^2 - 50*n)/660; \\ Michel Marcus, Jan 08 2015
    
  • Sage
    [n*(n+1)*(n+2)*(n^2+n-1)*(n^2+3*n+1)*(6*n^4+24*n^3+5*n^2-38*n+ 25)/660 for n in (1..30)] # G. C. Greubel, Aug 28 2019
    

Formula

a(n) = n*(n+1)*(n+2)*(n^2+n-1)*(n^2+3*n+1)*(6*n^4 + 24*n^3 + 5*n^2 - 38*n + 25)/660.
a(n) = 2*a(n-1) - a(n-2) + n^9.
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)^12. - Vincenzo Librandi, Jan 19 2015

A253710 Second partial sums of tenth powers (A008454).

Original entry on oeis.org

1, 1026, 61100, 1169750, 12044025, 83384476, 437200176, 1864757700, 6779099625, 21693441550, 62545208076, 165314338826, 405941961425, 935824239000, 2042356907200, 4248401203176, 8470439399601, 16262944822650, 30186516503500, 54350088184350, 95193540843401, 162596916293876, 271426802958000, 443660070587500
Offset: 1

Views

Author

Luciano Ancora, Jan 10 2015

Keywords

Comments

The formula for the second partial sums of m-th powers is: b(n,m) = (n+1)*F(m) - F(m+1), where F(m) are the m-th Faulhaber's formulas.

Programs

  • Mathematica
    a253710[n_] := Block[{f}, f[1] = 1; f[2] = 1026; f[x_] := 2*f[x - 1] - f[x - 2] + x^10; Array[f, n]]; a253710[21] (* Michael De Vlieger, Jan 11 2015 *)
    CoefficientList[Series[(1 + 1013 x + 47840 x^2 + 455192 x^3 + 1310354 x^4 + 1310354 x^5 + 455192 x^6 + 47840 x^7 + 1013 x^8 + x^9) / (1 - x)^13, {x, 0, 40}], x] (* Vincenzo Librandi, Jan 19 2015 *)
    Nest[Accumulate,Range[30]^10,2] (* Harvey P. Dale, May 10 2019 *)

Formula

a(n) = n*(n+1)^2*(n+2)*(n^2 + 2*n - 2)*(2*n^6 + 12*n^5 + 16*n^4 - 16*n^3 - 17*n^2 + 30*n - 5)/264.
a(n) = 2*a(n-1)-a(n-2)+n^10.
G.f.: x*(1 + 1013*x + 47840*x^2 + 455192*x^3 + 1310354*x^4 + 1310354*x^5 + 455192*x^6 + 47840*x^7 + 1013*x^8 + x^9)/(1-x)^13. - Vincenzo Librandi, Jan 19 2015
Showing 1-10 of 13 results. Next