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 10 results.

A101090 Third partial sums of fourth powers (A000583).

Original entry on oeis.org

1, 19, 135, 605, 2054, 5778, 14178, 31350, 63855, 121693, 219505, 378027, 625820, 1001300, 1555092, 2352732, 3477741, 5035095, 7155115, 9997801, 13757634, 18668870, 25011350, 33116850, 43375995, 56245761, 72257589, 92026135, 116258680, 145765224, 181469288
Offset: 1

Views

Author

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

Keywords

Comments

In general, the r-th successive summation of the fourth powers from 1 to n = (2*n+r)*(12*n^2+12*n*r+r^2-5*r)*(r+n)!/((r+4)!*(n-1)!). Here r = 3. - Gary Detlefs, Mar 01 2013

Crossrefs

Programs

  • Mathematica
    Nest[Accumulate, Range[50]^4, 3] (* Paolo Xausa, Jun 17 2024 *)

Formula

a(n) = (n*(1+n)*(2+n)*(3+n)*(3+2*n)*(-1+2*n*(3+n)))/840.
G.f.: x*(1+x)*(1+10*x+x^2)/(1-x)^8. - Colin Barker, Apr 04 2012
a(n) = (2*n+3)*(12*n^2+36*n-6)*(n+3)!/(5040*(n-1)!). - Gary Detlefs, Mar 01 2013

Extensions

Edited by Ralf Stephan, Dec 16 2004

A254681 Fifth partial sums of fourth powers (A000583).

Original entry on oeis.org

1, 21, 176, 936, 3750, 12342, 35112, 89232, 207207, 446875, 906048, 1743248, 3206268, 5670588, 9690000, 16062144, 25912029, 40797009, 62837104, 94875000, 140670530, 205134930, 294610680, 417203280, 583171875, 805386231
Offset: 1

Views

Author

Luciano Ancora, Feb 12 2015

Keywords

Examples

			Fourth differences:  1, 12,  23,  24, (repeat 24)  ...   (A101104)
Third differences:   1, 13,  36,  60,   84,   108, ...   (A101103)
Second differences:  1, 14,  50, 110,  194,   302, ...   (A005914)
First differences:   1, 15,  65, 175,  369,   671, ...   (A005917)
-------------------------------------------------------------------------
The fourth powers:   1, 16,  81, 256,  625,  1296, ...   (A000583)
-------------------------------------------------------------------------
First partial sums:  1, 17,  98, 354,  979,  2275, ...   (A000538)
Second partial sums: 1, 18, 116, 470, 1449,  3724, ...   (A101089)
Third partial sums:  1, 19, 135, 605, 2054,  5778, ...   (A101090)
Fourth partial sums: 1, 20, 155, 760, 2814,  8592, ...   (A101091)
Fifth partial sums:  1, 21, 176, 936, 3750, 12342, ...   (this sequence)
		

Crossrefs

Programs

  • Magma
    [Binomial(n+5,6)*n*(n+5)*(2*n+5)/42: n in [1..30]]; // G. C. Greubel, Dec 01 2018
    
  • Maple
    seq(coeff(series((x+11*x^2+11*x^3+x^4)/(1-x)^10,x,n+1), x, n), n = 1 .. 30); # Muniru A Asiru, Dec 02 2018
  • Mathematica
    Table[n^2(1+n)(2+n)(3+n)(4+n)(5+n)^2(5+2n)/30240, {n,26}] (* or *)
    CoefficientList[Series[(1 + 11 x + 11 x^2 + x^3)/(1-x)^10, {x,0,25}], x]
    CoefficientList[Series[(1/30240)E^x (30240 + 604800 x + 2041200 x^2 + 2368800 x^3 + 1233540 x^4 + 326592 x^5 + 46410 x^6 + 3540 x^7 + 135 x^8 + 2 x^9), {x, 0, 50}], x]*Table[n!, {n, 0, 50}] (* Stefano Spezia, Dec 02 2018 *)
    Nest[Accumulate[#]&,Range[30]^4,5] (* Harvey P. Dale, Jan 03 2022 *)
  • PARI
    my(x='x+O('x^30)); Vec((x+11*x^2+11*x^3+x^4)/(1-x)^10) \\ G. C. Greubel, Dec 01 2018
    
  • Sage
    [binomial(n+5,6)*n*(n+5)*(2*n+5)/42 for n in (1..30)] # G. C. Greubel, Dec 01 2018

Formula

G.f.:(x + 11*x^2 + 11*x^3 + x^4)/(1 - x)^10.
a(n) = n^2*(1 + n)*(2 + n)*(3 + n)*(4 + n)*(5 + n)^2*(5 + 2*n)/30240.
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5) + n^4.
E.g.f.: (1/30240)*exp(x)*(30240 + 604800*x + 2041200*x^2 + 2368800*x^3 + 1233540*x^4 + 326592*x^5 + 46410*x^6 + 3540*x^7 + 135*x^8 + 2*x^9). - Stefano Spezia, Dec 02 2018
From Amiram Eldar, Jan 26 2022: (Start)
Sum_{n>=1} 1/a(n) = 172032*log(2)/125 - 2382233/2500.
Sum_{n>=1} (-1)^(n+1)/a(n) = 42*Pi^2/25 - 43008*Pi/125 + 2663213/2500. (End)

A101104 a(1)=1, a(2)=12, a(3)=23, and a(n)=24 for n>=4.

Original entry on oeis.org

1, 12, 23, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24
Offset: 1

Views

Author

Cecilia Rossiter, Dec 15 2004

Keywords

Comments

Original name: The first summation of row 4 of Euler's triangle - a row that will recursively accumulate to the power of 4.

Crossrefs

For other sequences based upon MagicNKZ(n,k,z):
..... | n = 1 | n = 2 | n = 3 | n = 4 | n = 5 | n = 6 | n = 7
---------------------------------------------------------------------------
z = 0 | A000007 | A019590 | .......MagicNKZ(n,k,0) = A008292(n,k+1) .......
z = 1 | A000012 | A040000 | A101101 | thisSeq | A101100 | ....... | .......
z = 2 | A000027 | A005408 | A008458 | A101103 | A101095 | ....... | .......
z = 3 | A000217 | A000290 | A003215 | A005914 | A101096 | ....... | .......
z = 4 | A000292 | A000330 | A000578 | A005917 | A101098 | ....... | .......
z = 5 | A000332 | A002415 | A000537 | A000583 | A022521 | ....... | A255181
Cf. A101095 for an expanded table and more about MagicNKZ.

Programs

  • Mathematica
    MagicNKZ = Sum[(-1)^j*Binomial[n+1-z, j]*(k-j+1)^n, {j, 0, k+1}];Table[MagicNKZ, {n, 4, 4}, {z, 1, 1}, {k, 0, 34}]
    Join[{1, 12, 23},LinearRecurrence[{1},{24},56]] (* Ray Chandler, Sep 23 2015 *)

Formula

a(k) = MagicNKZ(4,k,1) where MagicNKZ(n,k,z) = Sum_{j=0..k+1} (-1)^j*binomial(n+1-z,j)*(k-j+1)^n (cf. A101095). That is, a(k) = Sum_{j=0..k+1} (-1)^j*binomial(4, j)*(k-j+1)^4.
a(1)=1, a(2)=12, a(3)=23, and a(n)=24 for n>=4. - Joerg Arndt, Nov 30 2014
G.f.: x*(1+11*x+11*x^2+x^3)/(1-x). - Colin Barker, Apr 16 2012

Extensions

New name from Joerg Arndt, Nov 30 2014
Original Formula edited and Crossrefs table added by Danny Rorabaugh, Apr 22 2015

A254644 Fourth partial sums of fifth powers (A000584).

Original entry on oeis.org

1, 36, 381, 2336, 10326, 36552, 110022, 292512, 704847, 1567852, 3263403, 6422208, 12046268, 21675408, 37608828, 63194304, 103199469, 164281524, 255573769, 389409504, 582206130, 855534680, 1237402530, 1763779680, 2480401755, 3444885756, 4729197591, 6422513536, 8634521016, 11499207456
Offset: 1

Views

Author

Luciano Ancora, Feb 05 2015

Keywords

Examples

			Fifth differences:   1, 27,  93,  119,   120, (repeat 120) (A101100)
Fourth differences:  1, 28, 121,  240,   360,   480, ...   (A101095)
Third differences:   1, 29, 150,  390,   750,  1230, ...   (A101096)
Second differences:  1, 30, 180,  570,  1320,  2550, ...   (A101098)
First differences:   1, 31, 211,  781,  2101,  4651, ...   (A022521)
-------------------------------------------------------------------------
The fifth powers:    1, 32, 243, 1024,  3125,  7776, ...   (A000584)
-------------------------------------------------------------------------
First partial sums:  1, 33, 276, 1300,  4425, 12201, ...   (A000539)
Second partial sums: 1, 34, 310, 1610,  6035, 18236, ...   (A101092)
Third partial sums:  1, 35, 345, 1955,  7990, 26226, ...   (A101099)
Fourth partial sums: 1, 36, 381, 2336, 10326, 36552, ...   (this sequence)
		

Crossrefs

Cf. A101091 (fourth partial sums of fourth powers).

Programs

  • GAP
    List([1..30], n-> Binomial(n+4,5)*(5*(n+2)^4 -35*(n+2)^2 +36)/126); # G. C. Greubel, Aug 28 2019
  • Magma
    [Binomial(n+4,5)*(5*(n+2)^4 -35*(n+2)^2 +36)/126: n in [1..30]]; // G. C. Greubel, Aug 28 2019
    
  • Maple
    seq(binomial(n+4,5)*(5*(n+2)^4 -35*(n+2)^2 +36)/126, n=1..30); # G. C. Greubel, Aug 28 2019
  • Mathematica
    Table[n(1+n)(2+n)(3+n)(4+n)(-24 +20n +85n^2 +40n^3 +5n^4)/15120, {n, 30}] (* or *) Accumulate[Accumulate[Accumulate[Accumulate[Range[24]^5]]]] (* or *) CoefficientList[Series[(1 +26x +66x^2 +26x^3 +x^4)/(1-x)^10, {x, 0, 30}], x]
    Nest[Accumulate,Range[30]^5,4] (* or *) LinearRecurrence[{10,-45,120, -210,252,-210,120,-45,10,-1}, {1,36,381,2336,10326,36552,110022,292512, 704847,1567852},30] (* Harvey P. Dale, May 08 2016 *)
  • PARI
    vector(30, n, m=n+2; binomial(m+2,5)*(5*m^4 -35*m^2 +36)/126) \\ G. C. Greubel, Aug 28 2019
    
  • Sage
    [binomial(n+4,5)*(5*(n+2)^4 -35*(n+2)^2 +36)/126 for n in (1..30)] # G. C. Greubel, Aug 28 2019
    

Formula

G.f.: x*(1 + 26*x + 66*x^2 + 26*x^3 + x^4)/(1 - x)^10.
a(n) = n*(1 + n)*(2 + n)*(3 + n)*(4 + n)*(-24 + 20*n + 85*n^2 + 40*n^3 + 5*n^4)/15120.
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) + n^5.

Extensions

Edited by Bruno Berselli, Feb 10 2015

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

A254470 Sixth partial sums of fourth powers (A000583).

Original entry on oeis.org

1, 22, 198, 1134, 4884, 17226, 52338, 141570, 348777, 795652, 1701700, 3444948, 6651216, 12321804, 22011804, 38073948, 63985977, 104782986, 167620090, 262495090, 403165620, 608300550, 902911230, 1320114510, 1903286385, 2708672616, 3808530792, 5294887048
Offset: 1

Views

Author

Luciano Ancora, Feb 15 2015

Keywords

Examples

			First differences:   1, 15,  65, 175,  369,   671, ... (A005917)
-------------------------------------------------------------------------
The fourth powers:   1, 16,  81, 256,  625,  1296, ... (A000583)
-------------------------------------------------------------------------
First partial sums:  1, 17,  98, 354,  979,  2275, ... (A000538)
Second partial sums: 1, 18, 116, 470, 1449,  3724, ... (A101089)
Third partial sums:  1, 19, 135, 605, 2054,  5778, ... (A101090)
Fourth partial sums: 1, 20, 155, 760, 2814,  8592, ... (A101091)
Fifth partial sums:  1, 21, 176, 936, 3750, 12342, ... (A254681)
Sixth partial sums:  1, 22, 198,1134, 4884, 17226, ... (this sequence)
		

Crossrefs

Programs

  • Magma
    [n*(1+n)*(2+n)*(3+n)^2*(4+n)*(5+n)*(6+n)*(1+12*n+ 2*n^2)/302400: n in [1..30]]; // Vincenzo Librandi, Feb 15 2015
    
  • Mathematica
    Table[n (1 + n) (2 + n) (3 + n)^2 (4 + n) (5 + n) (6 + n) (1 + 12 n + 2 n^2)/302400,{n, 25}] (* or *) CoefficientList[Series[(- 1 - 11 x - 11 x^2 - x^3)/(- 1 + x)^11, {x, 0, 24}], x]
    Nest[Accumulate,Range[30]^4,6] (* or *) LinearRecurrence[{11,-55,165,-330,462,-462,330,-165,55,-11,1},{1,22,198,1134,4884,17226,52338,141570,348777,795652,1701700},30] (* Harvey P. Dale, Apr 23 2016 *)
  • PARI
    vector(50,n,n*(1 + n)*(2 + n)*(3 + n)^2*(4 + n)*(5 + n)*(6 + n)*(1 + 12*n + 2*n^2)/302400) \\ Derek Orr, Feb 19 2015

Formula

G.f.: (-x - 11*x^2 - 11*x^3 - x^4)/(- 1 + x)^11.
a(n) = n*(1 + n)*(2 + n)*(3 + n)^2*(4 + n)*(5 + n)*(6 + n)*(1 + 12*n + 2*n^2)/302400.
a(n) = 6*a(n-1) - 15*a(n-2) + 20*a(n-3) - 15*a(n-4) + 6*a(n-5) - a(n-6) + n^4.
Sum_{n>=1} 1/a(n) = 3320303/2601 + 1400*Pi^2/17 + (8960/17)*sqrt(2/17)*Pi*cot(sqrt(17/2)*Pi). - Amiram Eldar, Jan 26 2022

A254870 Seventh partial sums of fourth powers (A000583).

Original entry on oeis.org

1, 23, 221, 1355, 6239, 23465, 75803, 217373, 566150, 1361802, 3063502, 6508450, 13159666, 25481470, 47493274, 85567222, 149553199, 254336185, 421956275, 684451365, 1087616985, 1695917535, 2598828765, 3918943275, 5822229660, 8530902276, 12339433068
Offset: 1

Views

Author

Luciano Ancora, Feb 17 2015

Keywords

Examples

			Second differences:   2, 14,  50,  110,  194,   302, ...   A120328(2k+1)
First differences:    1, 15,  65,  175,  369,   671, ...   A005917
--------------------------------------------------------------------------
The fourth powers:    1, 16,  81,  256,  625,  1296, ...   A000583
--------------------------------------------------------------------------
First partial sums:   1, 17,  98,  354,  979,  2275, ...   A000538
Second partial sums:  1, 18, 116,  470, 1449,  3724, ...   A101089
Third partial sums:   1, 19, 135,  605, 2054,  5778, ...   A101090
Fourth partial sums:  1, 20, 155,  760, 2814,  8592, ...   A101091
Fifth partial sums:   1, 21, 176,  936, 3750, 12342, ...   A254681
Sixth partial sums:   1, 22, 198, 1134, 4884, 17226, ...   A254470
Seventh partial sums: 1, 23, 221, 1355, 6239, 23465, ...   (this sequence)
		

Crossrefs

Programs

  • Magma
    [n*(1+n)*(2+n)*(3+n)*(4+n)*(5+n)*(6+n)*(7+n)*(7+2*n)*(7 +42*n+6*n^2)/19958400: n in [1..30]]; // Vincenzo Librandi, Feb 19 2015
  • Mathematica
    Table[n (1 + n) (2 + n) (3 + n) (4 + n) (5 + n) (6 + n) (7 + n) (7 + 2 n)((7 + 42 n + 6 n^2)/19958400), {n, 24}] (* or *)
    CoefficientList[Series[(1 + 11 x + 11 x^2 + x^3)/(- 1 + x)^12, {x, 0, 23}], x]
  • PARI
    vector(50,n,n*(1 + n)*(2 + n)*(3 + n)*(4 + n)*(5 + n)*(6 + n)*(7 + n)*(7 + 2*n)*(7 + 42*n + 6*n^2)/19958400) \\ Derek Orr, Feb 19 2015
    

Formula

G.f.: (x + 11*x^2 + 11*x^3 + x^4)/(- 1 + x)^12.
a(n) = n*(1 + n)*(2 + n)*(3 + n)*(4 + n)*(5 + n)*(6 + n)*(7 + n)*(7 + 2*n)*(7 + 42*n + 6*n^2)/19958400.
a(n) = 7*a(n-1) - 21*a(n-2) + 35*a(n-3) - 35*a(n-4) + 21*a(n-5) - 7*a(n-6) + a(n-7) + n^4.

A302352 a(n) = Sum_{k=0..n} k^4*binomial(2*n-k,n).

Original entry on oeis.org

0, 1, 19, 155, 936, 4884, 23465, 107107, 472600, 2036838, 8631206, 36119798, 149724940, 616104450, 2520629685, 10265200035, 41650094640, 168481778790, 679847488650, 2737640364810, 11005139655744, 44176226269728, 177114113623194, 709364594864910, 2838599638596176, 11350436081373340
Offset: 0

Views

Author

Ilya Gutkovskiy, Apr 06 2018

Keywords

Comments

Main diagonal of iterated partial sums array of fourth powers (starting with the first partial sums). For nonnegative integers see A002054, for squares see A265612, for cubes see A293550.

Crossrefs

Programs

  • Mathematica
    Table[Sum[k^4 Binomial[2 n - k, n], {k, 0, n}], {n, 0, 25}]
    Table[SeriesCoefficient[x (1 + 11 x + 11 x^2 + x^3)/(1 - x)^(n + 6), {x, 0, n}], {n, 0, 25}]
    Table[2^(2 n + 1) n (75 n^3 + 52 n^2 - 3 n - 4) Gamma[n + 3/2]/(Sqrt[Pi] Gamma[n + 6]), {n, 0, 25}]
    CoefficientList[Series[(24 - 180 x + 410 x^2 - 285 x^3 + 31 x^4 + Sqrt[1 - 4 x] (-24 + 132 x - 194 x^2 + 65 x^3 - x^4))/(2 Sqrt[1 - 4 x] x^5), {x, 0, 25}], x]
    CoefficientList[Series[E^(2 x) (-576 + 360 x - 244 x^2 + 75 x^3) BesselI[0, 2 x]/x^3 + E^(2 x) (576 - 360 x + 532 x^2 - 255 x^3 + 75 x^4) BesselI[1, 2 x]/x^4, {x, 0, 25}], x]* Range[0, 25]!
  • PARI
    a(n) = sum(k=0, n, k^4*binomial(2*n-k,n)); \\ Michel Marcus, Apr 07 2018

Formula

a(n) = [x^n] x*(1 + 11*x + 11*x^2 + x^3)/(1 - x)^(n+6).
a(n) = 2^(2*n+1)*n*(75*n^3 + 52*n^2 - 3*n - 4)*Gamma(n+3/2)/(sqrt(Pi)*Gamma(n+6)).
a(n) ~ 75*2^(2*n+1)/sqrt(Pi*n).

A257449 a(n) = 75*(2^n - 1) - 4*n^3 - 18*n^2 - 52*n.

Original entry on oeis.org

1, 17, 99, 373, 1115, 2901, 6907, 15509, 33483, 70405, 145451, 296997, 601819, 1213493, 2439195, 4893301, 9804587, 19630629, 39286603, 78602885, 157240251, 314520277, 629086139, 1258224213, 2516507275, 5033080901, 10066236267, 20132555749, 40265204123
Offset: 1

Views

Author

Luciano Ancora, Apr 23 2015

Keywords

Comments

See the first comment of A257448.

Examples

			This sequence provides the antidiagonal sums of the array:
1, 16,  81, 256,  625,  1296, ...   A000583
1, 17,  98, 354,  979,  2275, ...   A000538
1, 18, 116, 470, 1449,  3724, ...   A101089
1, 19, 135, 605, 2054,  5778, ...   A101090
1, 20, 155, 760, 2814,  8592, ...   A101091
1, 21, 176, 936, 3750, 12342, ...   A254681
...
See also A254681 (Example field).
		

Crossrefs

Programs

  • Magma
    [75*(2^n-1)-4*n^3-18*n^2-52*n: n in [1..30]]; // Vincenzo Librandi, Apr 24 2015
  • Mathematica
    Table[75 (2^n - 1) - 4 n^3 - 18 n^2 - 52 n, {n, 30}]

Formula

G.f.: -x*(1 + x)*(1 + 10*x + x^2)/((-1 + x)^4*(-1 + 2*x)).
a(n) = 6*a(n-1) -14*a(n-2) +16*a(n-3) -9*a(n-4) +2*a(n-5) for n>5.

A302353 a(n) = Sum_{k=0..n} k^n*binomial(2*n-k,n).

Original entry on oeis.org

1, 1, 7, 69, 936, 16290, 345857, 8666413, 250355800, 8191830942, 299452606190, 12095028921250, 534924268768540, 25710497506696860, 1334410348734174285, 74379234152676275325, 4431350132232658244400, 281020603194039519937590, 18900157831016574533520330, 1343698678390575915132318870
Offset: 0

Views

Author

Ilya Gutkovskiy, Apr 06 2018

Keywords

Comments

a(n) is the n-th term of the main diagonal of iterated partial sums array of n-th powers (starting with the first partial sums).

Examples

			For n = 4 we have:
------------------------
0   1    2    3    [4]
------------------------
0,  1,  17,   98,  354,  ... A000538 (partial sums of fourth powers)
0,  1,  18,  116,  470,  ... A101089 (partial sums of A000538)
0,  1,  19,  135,  605,  ... A101090 (partial sums of A101089)
0,  1,  20,  155,  760,  ... A101091 (partial sums of A101090)
0,  1,  21,  176, [936], ... A254681 (partial sums of A101091)
------------------------
therefore a(4) = 936.
		

Crossrefs

Programs

  • Mathematica
    Join[{1}, Table[Sum[k^n Binomial[2 n - k, n], {k, 0, n}], {n, 19}]]
    Table[SeriesCoefficient[HurwitzLerchPhi[x, -n, 0]/(1 - x)^(n + 1), {x, 0, n}], {n, 0, 19}]

Formula

a(n) ~ c * (r * (2-r)^(2-r) / (1-r)^(1-r))^n * n^n, where r = 0.69176629470097668698335106516328398961170464277337300459988208658267146... is the root of the equation (2-r) = (1-r) * exp(1/r) and c = 0.96374921279011282619632879505754646526289414675402231447188230355850496... - Vaclav Kotesovec, Apr 08 2018
Showing 1-10 of 10 results.