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 12 results. Next

A000537 Sum of first n cubes; or n-th triangular number squared.

Original entry on oeis.org

0, 1, 9, 36, 100, 225, 441, 784, 1296, 2025, 3025, 4356, 6084, 8281, 11025, 14400, 18496, 23409, 29241, 36100, 44100, 53361, 64009, 76176, 90000, 105625, 123201, 142884, 164836, 189225, 216225, 246016, 278784, 314721, 354025, 396900, 443556, 494209, 549081
Offset: 0

Views

Author

Keywords

Comments

Number of parallelograms in an n X n rhombus. - Matti De Craene (Matti.DeCraene(AT)rug.ac.be), May 14 2000
Or, number of orthogonal rectangles in an n X n checkerboard, or rectangles in an n X n array of squares. - Jud McCranie, Feb 28 2003. Compare A085582.
Also number of 2-dimensional cage assemblies (cf. A059827, A059860).
The n-th triangular number T(n) = Sum_{r=1..n} r = n(n+1)/2 satisfies the relations: (i) T(n) + T(n-1) = n^2 and (ii) T(n) - T(n-1) = n by definition, so that n^2*n = n^3 = {T(n)}^2 - {T(n-1)}^2 and by summing on n we have Sum_{ r = 1..n } r^3 = {T(n)}^2 = (1+2+3+...+n)^2 = (n*(n+1)/2)^2. - Lekraj Beedassy, May 14 2004
Number of 4-tuples of integers from {0,1,...,n}, without repetition, whose last component is strictly bigger than the others. Number of 4-tuples of integers from {1,...,n}, with repetition, whose last component is greater than or equal to the others.
Number of ordered pairs of two-element subsets of {0,1,...,n} without repetition.
Number of ordered pairs of 2-element multisubsets of {1,...,n} with repetition.
1^3 + 2^3 + 3^3 + ... + n^3 = (1 + 2 + 3 + ... + n)^2.
a(n) is the number of parameters needed in general to know the Riemannian metric g of an n-dimensional Riemannian manifold (M,g), by knowing all its second derivatives; even though to know the curvature tensor R requires (due to symmetries) (n^2)*(n^2-1)/12 parameters, a smaller number (and a 4-dimensional pyramidal number). - Jonathan Vos Post, May 05 2006
Also number of hexagons with vertices in an hexagonal grid with n points in each side. - Ignacio Larrosa Cañestro, Oct 15 2006
Number of permutations of n distinct letters (ABCD...) each of which appears twice with 4 and n-4 fixed points. - Zerinvary Lajos, Nov 09 2006
With offset 1 = binomial transform of [1, 8, 19, 18, 6, ...]. - Gary W. Adamson, Dec 03 2008
The sequence is related to A000330 by a(n) = n*A000330(n) - Sum_{i=0..n-1} A000330(i): this is the case d=1 in the identity n*(n*(d*n-d+2)/2) - Sum_{i=0..n-1} i*(d*i-d+2)/2 = n*(n+1)*(2*d*n-2*d+3)/6. - Bruno Berselli, Apr 26 2010, Mar 01 2012
From Wolfdieter Lang, Jan 11 2013: (Start)
For sums of powers of positive integers S(k,n) := Sum_{j=1..n}j^k one has the recurrence S(k,n) = (n+1)*S(k-1,n) - Sum_{l=1..n} S(k-1,l), n >= 1, k >= 1.
This was used for k=4 by Ibn al-Haytham in an attempt to compute the volume of the interior of a paraboloid. See the Strick reference where the trick he used is shown, and the W. Lang link.
This trick generalizes immediately to arbitrary powers k. For k=3: a(n) = (n+1)*A000330(n) - Sum_{l=1..n} A000330(l), which coincides with the formula given in the previous comment by Berselli. (End)
Regarding to the previous contribution, see also Matem@ticamente in Links field and comments on this recurrences in similar sequences (partial sums of n-th powers). - Bruno Berselli, Jun 24 2013
A rectangular prism with sides A000217(n), A000217(n+1), and A000217(n+2) has surface area 6*a(n+1). - J. M. Bergot, Aug 07 2013, edited with corrected indices by Antti Karttunen, Aug 09 2013
A formula for the r-th successive summation of k^3, for k = 1 to n, is (6*n^2+r*(6*n+r-1)*(n+r)!)/((r+3)!*(n-1)!), (H. W. Gould). - Gary Detlefs, Jan 02 2014
Note that this sequence and its formula were known to (and possibly discovered by) Nicomachus, predating Ibn al-Haytham by 800 years. - Charles R Greathouse IV, Apr 23 2014
a(n) is the number of ways to paint the sides of a nonsquare rectangle using at most n colors. Cf. A039623. - Geoffrey Critzer, Jun 18 2014
For n > 0: A256188(a(n)) = A000217(n) and A256188(m) != A000217(n) for m < a(n), i.e., positions of first occurrences of triangular numbers in A256188. - Reinhard Zumkeller, Mar 26 2015
There is no cube in this sequence except 0 and 1. - Altug Alkan, Jul 02 2016
Also the number of chordless cycles in the complete bipartite graph K_{n+1,n+1}. - Eric W. Weisstein, Jan 02 2018
a(n) is the sum of the elements in the multiplication table [0..n] X [0..n]. - Michel Marcus, May 06 2021

Examples

			G.f. = x + 9*x^2 + 36*x^3 + 100*x^4 + 225*x^5 + 441*x^6 + ... - _Michael Somos_, Aug 29 2022
		

References

  • 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.
  • Avner Ash and Robert Gross, Summing it up, Princeton University Press, 2016, p. 62, eq. (6.3) for k=3.
  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, p. 110ff.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 155.
  • John H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, pp. 36, 58.
  • Clifford Pickover, "Wonders of Numbers, Adventures in Mathematics, Mind and Meaning," Oxford University Press, 2001, p. 325.
  • 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).
  • H. K. Strick, Geschichten aus der Mathematik II, Spektrum Spezial 3/11, p. 13.
  • D. Wells, You Are A Mathematician, "Counting rectangles in a rectangle", Problem 8H, pp. 240; 254, Penguin Books 1995.

Crossrefs

Convolution of A000217 and A008458.
Row sums of triangles A094414 and A094415.
Second column of triangle A008459.
Row 3 of array A103438.
Cf. A236770 (see crossrefs).

Programs

  • GAP
    List([0..40],n->(n*(n+1)/2)^2); # Muniru A Asiru, Dec 05 2018
    
  • Haskell
    a000537 = a000290 . a000217  -- Reinhard Zumkeller, Mar 26 2015
    
  • Magma
    [(n*(n+1)/2)^2: n in [0..50]]; // Wesley Ivan Hurt, Jun 06 2014
    
  • Maple
    a:= n-> (n*(n+1)/2)^2:
    seq(a(n), n=0..40);
  • Mathematica
    Accumulate[Range[0, 50]^3] (* Harvey P. Dale, Mar 01 2011 *)
    f[n_] := n^2 (n + 1)^2/4; Array[f, 39, 0] (* Robert G. Wilson v, Nov 16 2012 *)
    Table[CycleIndex[{{1, 2, 3, 4}, {3, 2, 1, 4}, {1, 4, 3, 2}, {3, 4, 1, 2}}, s] /. Table[s[i] -> n, {i, 1, 2}], {n, 0, 30}] (* Geoffrey Critzer, Jun 18 2014 *)
    Accumulate @ Range[0, 50]^2 (* Waldemar Puszkarz, Jan 24 2015 *)
    Binomial[Range[20], 2]^2 (* Eric W. Weisstein, Jan 02 2018 *)
    LinearRecurrence[{5, -10, 10, -5, 1}, {0, 1, 9, 36, 100}, 20] (* Eric W. Weisstein, Jan 02 2018 *)
    CoefficientList[Series[-((x (1 + 4 x + x^2))/(-1 + x)^5), {x, 0, 20}], x] (* Eric W. Weisstein, Jan 02 2018 *)
  • PARI
    a(n)=(n*(n+1)/2)^2
    
  • Python
    def A000537(n): return (n*(n+1)>>1)**2 # Chai Wah Wu, Oct 20 2023

Formula

a(n) = (n*(n+1)/2)^2 = A000217(n)^2 = Sum_{k=1..n} A000578(k), that is, 1^3 + 2^3 + 3^3 + ... + n^3 = (1 + 2 + 3 + ... + n)^2.
G.f.: (x+4*x^2+x^3)/(1-x)^5. - Simon Plouffe in his 1992 dissertation
a(n) = Sum ( Sum ( 1 + Sum (6*n) ) ), rephrasing the formula in A000578. - Xavier Acloque, Jan 21 2003
a(n) = Sum_{i=1..n} Sum_{j=1..n} i*j, row sums of A127777. - Alexander Adamchuk, Oct 24 2004
a(n) = A035287(n)/4. - Zerinvary Lajos, May 09 2007
This sequence could be obtained from the general formula n*(n+1)*(n+2)*(n+3)*...*(n+k)*(n*(n+k) + (k-1)*k/6)/((k+3)!/6) at k=1. - Alexander R. Povolotsky, May 17 2008
G.f.: x*F(3,3;1;x). - Paul Barry, Sep 18 2008
Sum_{k > 0} 1/a(k) = (4/3)*(Pi^2-9). - Jaume Oliver Lafont, Sep 20 2009
a(n) = Sum_{1 <= k <= m <= n} A176271(m,k). - Reinhard Zumkeller, Apr 13 2010
a(n) = Sum_{i=1..n} J_3(i)*floor(n/i), where J_ 3 is A059376. - Enrique Pérez Herrero, Feb 26 2012
a(n) = Sum_{i=1..n} Sum_{j=1..n} Sum_{k=1..n} min(i,j,k). - Enrique Pérez Herrero, Feb 26 2013 [corrected by Ridouane Oudra, Mar 05 2025]
a(n) = 6*C(n+2,4) + C(n+1,2) = 6*A000332(n+2) + A000217(n), (Knuth). - Gary Detlefs, Jan 02 2014
a(n) = -Sum_{j=1..3} j*Stirling1(n+1,n+1-j)*Stirling2(n+3-j,n). - Mircea Merca, Jan 25 2014
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*(3-4*log(2)). - Vaclav Kotesovec, Feb 13 2015
a(n)*((s-2)*(s-3)/2) = P(3, P(s, n+1)) - P(s, P(3, n+1)), where P(s, m) = ((s-2)*m^2-(s-4)*m)/2 is the m-th s-gonal number. For s=7, 10*a(n) = A000217(A000566(n+1)) - A000566(A000217(n+1)). - Bruno Berselli, Aug 04 2015
From Ilya Gutkovskiy, Jul 03 2016: (Start)
E.g.f.: x*(4 + 14*x + 8*x^2 + x^3)*exp(x)/4.
Dirichlet g.f.: (zeta(s-4) + 2*zeta(s-3) + zeta(s-2))/4. (End)
a(n) = (Bernoulli(4, n+1) - Bernoulli(4, 1))/4, n >= 0, with the Bernoulli polynomial B(4, x) from row n=4 of A053382/A053383. See, e.g., the Ash-Gross reference, p. 62, eq. (6.3) for k=3. - Wolfdieter Lang, Mar 12 2017
a(n) = A000217((n+1)^2) - A000217(n+1)^2. - Bruno Berselli, Aug 31 2017
a(n) = n*binomial(n+2, 3) + binomial(n+2, 4) + binomial(n+1, 4). - Tony Foster III, Nov 14 2017
Another identity: ..., a(3) = (1/2)*(1*(2+4+6)+3*(4+6)+5*6) = 36, a(4) = (1/2)*(1*(2+4+6+8)+3*(4+6+8)+5*(6+8)+7*(8)) = 100, a(5) = (1/2)*(1*(2+4+6+8+10)+3*(4+6+8+10)+5*(6+8+10)+7*(8+10)+9*(10)) = 225, ... - J. M. Bergot, Aug 27 2022
Comment from Michael Somos, Aug 28 2022: (Start)
The previous comment expresses a(n) as the sum of all of the n X n multiplication table array entries.
For example, for n = 4:
1 2 3 4
2 4 6 8
3 6 9 12
4 8 12 16
This array sum can be split up as follows:
+---+---------------+
| 0 | 1 2 3 4 | (0+1)*(1+2+3+4)
| +---+-----------+
| 0 | 2 | 4 6 8 | (1+2)*(2+3+4)
| | +---+-------+
| 0 | 3 | 6 | 9 12 | (2+3)*(3+4)
| | | +---+---+
| 0 | 4 | 8 |12 |16 | (3+4)*(4)
+---+---+---+---+---+
This kind of row+column sums was used by Ramanujan and others for summing Lambert series. (End)
a(n) = 6*A000332(n+4) - 12*A000292(n+1) + 7*A000217(n+1) - n - 1. - Adam Mohamed, Sep 05 2024

Extensions

Edited by M. F. Hasler, May 02 2015

A059860 a(n) = binomial(n+1, 2)^5.

Original entry on oeis.org

1, 243, 7776, 100000, 759375, 4084101, 17210368, 60466176, 184528125, 503284375, 1252332576, 2887174368, 6240321451, 12762815625, 24883200000, 46525874176, 83841135993, 146211169851, 247609900000, 408410100000, 657748550151, 1036579476493, 1601568101376
Offset: 1

Views

Author

Jason Earls, Feb 28 2001

Keywords

Comments

Number of 5-dimensional cage assemblies.
See Chap. 61, "Hyperspace Prisons", of C. Pickover's book "Wonders of Numbers" for full explanation of "cage numbers."

References

  • Clifford A. Pickover, Wonders of Numbers, Oxford University Press, 2001, p. 325.

Crossrefs

Cf. A059827.

Programs

  • Maple
    for n from 1 to 100 do printf(`%d,`,((n^5)*(n + 1)^5)/(2^5) ) od:
    with (combinat):seq(mul(stirling2(n+1,n),k=1..5),n=1..21); # Zerinvary Lajos, Dec 14 2007
  • Mathematica
    m = 5; Table[ ( (n^m)(n + 1)^m )/(2^m), {n, 1, 26} ]
    Table[Binomial[n+1,2]^5,{n,20}] (* Harvey P. Dale, May 04 2018 *)
  • PARI
    a(n) = { (n*(n + 1)/2)^5 } \\ Harry J. Smith, Jun 29 2009

Formula

L(n) = ((n^m)(n + 1)^m)/(2^m) where m is the dimension.
G.f.: x * (x^8 +232*x^7 +5158*x^6+ 27664*x^5 +47290*x^4 +27664*x^3 +5158*x^2 +232*x +1) / (1-x)^11. - Colin Barker, Jun 28 2012
From Amiram Eldar, May 15 2022: (Start)
Sum_{n>=1} 1/a(n) = 4032 - 1120*Pi^2/3 - 32*Pi^4/9.
Sum_{n>=1} (-1)^(n+1)/a(n) = 4480*log(2) + 720*zeta(3) + 60*zeta(5) - 4032. (End)

Extensions

More terms from James Sellers, Feb 28 2001
Better definition from Zerinvary Lajos, May 23 2006
Corrected the definition from binomial(n+2,2)^5 to binomial(n+1,2)^5. - Harry J. Smith, Jun 29 2009

A059977 a(n) = binomial(n+2, 2)^4.

Original entry on oeis.org

1, 81, 1296, 10000, 50625, 194481, 614656, 1679616, 4100625, 9150625, 18974736, 37015056, 68574961, 121550625, 207360000, 342102016, 547981281, 855036081, 1303210000, 1944810000, 2847396321, 4097152081, 5802782976, 8100000000, 11156640625, 15178486401
Offset: 0

Views

Author

Robert G. Wilson v, Mar 06 2001

Keywords

Comments

Number of 4-dimensional cage assemblies.
See Chap. 61, "Hyperspace Prisons", of C. Pickover's book "Wonders of Numbers" for full explanation of "cage numbers."

Examples

			1 = (1 + 1)/2, 81 = (33 + 129)/2, 1296 = (276 + 2316)/2, 10000 = (1300 + 18700)/2, ... - _Philippe Deléham_, May 25 2015
		

References

  • Clifford A. Pickover, "Wonders of Numbers, Adventures in Mathematics, Mind and Meaning," Oxford University Press, 2001, p. 325.

Crossrefs

Cf. A168364 (first differences).

Programs

  • Maple
    with (combinat):seq(mul(stirling2(n+1,n),k=1..4),n=1..24); # Zerinvary Lajos, Dec 16 2007
  • Mathematica
    m = 4; Table[ ( (n^m)(n + 1)^m )/(2^m), {n, 1, 30} ]
  • PARI
    a(n) = { ((n + 1)*(n + 2)/2)^4 } \\ Harry J. Smith, Jun 30 2009
  • Sage
    [stirling_number2(n+1,n)^4for n in range(1,25)] # Zerinvary Lajos, Mar 14 2009
    

Formula

L(n) = ((n^m)(n + 1)^m)/(2^m) where m is the dimension, which in this case is 4.
O.g.f.: -(1+72*x+603*x^2+1168*x^3+603*x^4+72*x^5+x^6)/(-1+x)^9. - R. J. Mathar, Mar 31 2008
a(n) = A000217(n+1)^4. - R. J. Mathar, Dec 13 2011
a(n) = (A000539(n+1) + A000541(n+1))/2. - Philippe Deléham, May 25 2015
From Amiram Eldar, May 15 2022: (Start)
Sum_{n>=0} 1/a(n) = 160*Pi^2/3 + 16*Pi^4/45 - 560.
Sum_{n>=0} (-1)^n/a(n) = 560 - 640*log(2) - 96*zeta(3). (End)

Extensions

Better definition from Zerinvary Lajos, May 23 2006

A059978 a(n) = binomial(n+2,n)^6.

Original entry on oeis.org

1, 729, 46656, 1000000, 11390625, 85766121, 481890304, 2176782336, 8303765625, 27680640625, 82653950016, 225199600704, 567869252041, 1340095640625, 2985984000000, 6327518887936, 12827693806929, 25002110044521, 47045881000000, 85766121000000, 151939915084881
Offset: 0

Views

Author

Robert G. Wilson v, Mar 06 2001

Keywords

Comments

Number of 6-dimensional cage assemblies.

References

  • Clifford A. Pickover, "Wonders of Numbers, Adventures in Mathematics, Mind and Meaning," Oxford University Press, 2001, p. 325.

Crossrefs

Programs

  • Maple
    with (combinat):seq(mul(stirling2(n+1,n),k=1..6),n=1..18); # Zerinvary Lajos, Dec 14 2007
  • Mathematica
    m = 6; Table[n^m (n + 1)^m/2^m, {n, 1, 24}]

Formula

G.f.: (x^10 + 716*x^9 + 37257*x^8 + 450048*x^7 + 1822014*x^6 +2864328*x^5 + 1822014*x^4 + 450048*x^3 + 37257*x^2 + 716*x + 1)/(1-x)^13. - Colin Barker, Jul 09 2012
G.f.: 6F5([3,3,3,3,3,3], [1,1,1,1,1], z). - Benedict W. J. Irwin, Mar 14 2016
a(n) = (1/16)*( 3*S(7,n+1) + 10*S(9,n+1) + 3*S(11,n+1) ), where S(r,n) = Sum_{k = 1..n} k^r. Cf. A059977 and A059980. - Peter Bala, Jul 02 2019
From Amiram Eldar, May 15 2022: (Start)
Sum_{n>=0} 1/a(n) = 2688*Pi^2 + 448*Pi^4/15 + 128*Pi^6/945 - 29568.
Sum_{n>=0} (-1)^n/a(n) = 29568 - 32256*log(2) - 5376*zeta(3) - 720*zeta(5). (End)

Extensions

Better definition from Zerinvary Lajos, May 23 2006

A128629 A triangular array generated by moving Pascal sequences to prime positions and embedding new sequences at the nonprime locations. (cf. A007318 and A000040).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 6, 4, 1, 1, 4, 9, 10, 5, 1, 1, 6, 10, 16, 15, 6, 1, 1, 5, 18, 20, 25, 21, 7, 1, 1, 8, 15, 40, 35, 36, 28, 8, 1, 1, 9, 27, 35, 75, 56, 49, 36, 9, 1
Offset: 1

Views

Author

Alford Arnold, Mar 29 2007

Keywords

Comments

The array can be constructed by beginning with A007318 (Pascal's triangle) placing each diagonal on a prime row. The other rows are filled in by mapping the prime factorization of the row number to the known sequences on the prime rows and multiplying term by term.

Examples

			Row six begins 1 6 18 40 75 126 ... because rows two and three are
1 2 3 4 5 6 ...
1 3 6 10 15 21 ...
The array begins
1 1 1 1 1 1 1 1 1 A000012
1 2 3 4 5 6 7 8 9 A000027
1 3 6 10 15 21 28 36 45 A000217
1 4 9 16 25 36 49 64 81 A000290
1 4 10 20 35 56 84 120 165 A000292
1 6 18 40 75 126 196 288 405 A002411
1 5 15 35 70 126 210 330 495 A000332
1 8 27 64 125 216 343 512 729 A000578
1 9 36 100 225 441 784 1296 2025 A000537
1 8 30 80 175 336 588 960 1485 A002417
1 6 21 56 126 252 462 792 1287 A000389
1 12 54 160 375 756 1372 2304 3645 A019582
1 7 28 84 210 462 924 1716 3003 A000579
1 10 45 140 350 756 1470 2640 4455 A027800
1 12 60 200 525 1176 2352 4320 7425 A004302
1 16 81 256 625 1296 2401 4096 6561 A000583
1 8 36 120 330 792 1716 3432 6435 A000580
1 18 108 400 1125 2646 5488 10368 18225 A019584
1 9 45 165 495 1287 3003 6435 12870 A000581
1 16 90 320 875 2016 4116 7680 13365 A119771
1 15 90 350 1050 2646 5880 11880 22275 A001297
1 12 63 224 630 1512 3234 6336 11583 A027810
1 10 55 220 715 2002 5005 11440 24310 A000582
1 24 162 640 1875 4536 9604 18432 32805 A019583
1 16 100 400 1225 3136 7056 14400 27225 A001249
1 14 84 336 1050 2772 6468 13728 27027 A027818
1 27 216 1000 3375 9261 21952 46656 91125 A059827
1 20 135 560 1750 4536 10290 21120 40095 A085284
		

Crossrefs

Cf. A064553 (second diagonal), A080688 (second diagonal resorted).

Programs

  • Maple
    A128629 := proc(n,m) if n = 1 then 1; elif isprime(n) then p := numtheory[pi](n) ; binomial(p+m-1,p) ; else a := 1 ; for p in ifactors(n)[2] do a := a* procname(op(1,p),m)^ op(2,p) ; od: fi; end: # R. J. Mathar, Sep 09 2009

Extensions

A-number added to each row of the examples by R. J. Mathar, Sep 09 2009

A070333 Expansion of (1 + x)*(1 - x + x^2)/((1 - x)^4*(1 + x + x^2)).

Original entry on oeis.org

1, 3, 6, 12, 21, 33, 50, 72, 99, 133, 174, 222, 279, 345, 420, 506, 603, 711, 832, 966, 1113, 1275, 1452, 1644, 1853, 2079, 2322, 2584, 2865, 3165, 3486, 3828, 4191, 4577, 4986, 5418, 5875, 6357, 6864, 7398, 7959, 8547, 9164, 9810, 10485
Offset: 0

Views

Author

N. J. A. Sloane, May 11 2002

Keywords

Comments

a(n) is the number of 3 X 3 matrices with nonnegative integer entries such that every row sum, column sum and the trace of the matrix is n. - Sharon Sela (sharonsela(AT)hotmail.com), May 20 2002

Crossrefs

Programs

  • Magma
    [Round((2*n+3)*(n^2+3*n+6)/18): n in [0..50]]; // Vincenzo Librandi, Jun 25 2011
    
  • Maple
    A049347 := proc(n) op(1+(n mod 3),[1,-1,0]) ; end proc:
    A070333 := proc(n) 1+7*n/6+n^2/2+n^3/9+2*A049347(n-1)/9 ; end proc: # R. J. Mathar, Dec 03 2010
  • Mathematica
    CoefficientList[ Series[(1 + x^3)/(1 - 3*x + 3*x^2 - 2*x^3 + 3*x^4 - 3*x^5 + x^6), {x, 0, 45}], x]
  • PARI
    Vec((1+x)*(1-x+x^2)/((1-x)^4*(1+x+x^2))+O(x^99)) \\ Charles R Greathouse IV, Sep 27 2012

Formula

a(n) = 1 + 7*n/6 + n^2/2 + n^3/9 + 2*A049347(n-1)/9. - R. J. Mathar, Dec 03 2010
From Mircea Merca, Dec 03 2010: (Start)
a(n) = round((2*n + 3)*(n^2 + 3*n + 6)/18).
a(n) = floor((n + 2)*(2*n^2 + 5*n + 11)/18).
a(n) = ceiling((n + 1)*(2*n^2 + 7*n + 14)/18).
a(n) = round((n + 1)*(2*n^2 + 7*n + 14)/18).
a(n) = a(n-3) + n^2 + 2 for n > 2. (End)
E.g.f.: exp(x)*(1 + x*(32 + x*(15 + 2*x))/18) + 4*exp(-x/2)*sin(sqrt(3)*x/2)/(9*sqrt(3)). - Stefano Spezia, Oct 28 2022

A357178 First differences of cubes of triangular numbers.

Original entry on oeis.org

0, 1, 26, 189, 784, 2375, 5886, 12691, 24704, 44469, 75250, 121121, 187056, 279019, 404054, 570375, 787456, 1066121, 1418634, 1858789, 2402000, 3065391, 3867886, 4830299, 5975424, 7328125, 8915426, 10766601, 12913264, 15389459, 18231750, 21479311, 25174016, 29360529
Offset: 0

Views

Author

Kelvin Voskuijl, Sep 16 2022

Keywords

Comments

Row sums of centered hexagonal numbers A003215 treated as a regular triangle.

Crossrefs

Cf. A059827 (cubes of triangular numbers).
Cf. A000578 (for squares) and A168364 (for fourth powers) of triangular numbers.
Cf. A000217 (triangular numbers), A003215.

Programs

Formula

a(n) = (n^3 + 3*n^5)/4.
G.f.: x*(1 + 20*x + 48*x^2 + 20*x^3 + x^4)/(1 - x)^6. - Stefano Spezia, Sep 19 2022

A059979 Number of 7-dimensional cage assemblies.

Original entry on oeis.org

1, 2187, 279936, 10000000, 170859375, 1801088541, 13492928512, 78364164096, 373669453125, 1522435234375, 5455160701056, 17565568854912, 51676101935731, 140710042265625, 358318080000000, 860542568759296, 1962637152460137, 4275360817613091, 8938717390000000
Offset: 1

Views

Author

Robert G. Wilson v, Mar 06 2001

Keywords

References

  • Clifford A. Pickover, "Wonders of Numbers, Adventures in Mathematics, Mind and Meaning," Oxford University Press, 2001, p. 325.

Crossrefs

Programs

  • Mathematica
    m = 7; Table[ ( (n^m)(n + 1)^m )/(2^m), {n, 1, 20} ]
    (Times@@@Partition[Range[20]^7,2,1])/2^7 (* Harvey P. Dale, Aug 20 2017 *)

Formula

G.f.: -x*(x^12 +2172*x^11 +247236*x^10+ 6030140*x^9 +49258935*x^8 +163809288*x^7 +242384856*x^6 +163809288*x^5 +49258935*x^4 +6030140*x^3 +247236*x^2 +2172*x +1)/(x-1)^15. - Colin Barker, Jul 09 2012
From Benedict W. J. Irwin, Mar 14 2016: (Start)
G.f.: z*7F6([3,3,3,3,3,3,3], [1,1,1,1,1,1], z).
a(n) = n^7*(1+n)^7/128.
(End)
a(n) = binomial(n+1, 2)^7. - Alejandro Rodriguez, Oct 20 2020
From Amiram Eldar, May 15 2022: (Start)
Sum_{n>=1} 1/a(n) = 219648 - 19712*Pi^2 - 3584*Pi^4/15 - 256*Pi^6/135.
Sum_{n>=1} (-1)^(n+1)/a(n) = 236544*log(2) + 40320*zeta(3) + 6720*zeta(5) + 252*zeta(7) - 219648. (End)

A059980 Number of 8-dimensional cage assemblies.

Original entry on oeis.org

1, 6561, 1679616, 100000000, 2562890625, 37822859361, 377801998336, 2821109907456, 16815125390625, 83733937890625, 360040606269696, 1370114370683136, 4702525276151521, 14774554437890625, 42998169600000000, 117033789351264256, 300283484326400961
Offset: 1

Views

Author

Robert G. Wilson v, Mar 06 2001

Keywords

References

  • Clifford A. Pickover, "Wonders of Numbers, Adventures in Mathematics, Mind and Meaning," Oxford University Press, 2001, p. 325.

Crossrefs

Programs

  • Mathematica
    m = 8; Table[n^m (n + 1)^m/2^m, {n, 1, 18}]

Formula

G.f.: -x*(x^14 +6544*x^13 +1568215*x^12 +72338144*x^11 +1086859301*x^10 +6727188848*x^9 +19323413187*x^8 +27306899520*x^7 +19323413187*x^6 +6727188848*x^5 +1086859301*x^4 +72338144*x^3 +1568215*x^2 +6544*x +1)/(x-1)^17. - Colin Barker, Jul 09 2012
From Peter Bala, Jul 02 2019 (Start)
a(n) = (n*(n + 1)/2)^8.
a(n) = (1/16)*( S(9,n) + 7*S(11,n) + 7*S(13,n) + S(15,n) ), where S(r,n) = Sum_{k = 1..n} k^r. Cf. A059977 and A059978. (End)
From Amiram Eldar, May 15 2022: (Start)
Sum_{n>=1} 1/a(n) = 146432*Pi^2 + 5632*Pi^4/3 + 2048*Pi^6/105 + 256*Pi^8/4725 - 1647360.
Sum_{n>=1} (-1)^(n+1)/a(n) = 1647360 - 1757184*log(2) - 304128*zeta(3) - 57600*zeta(5) - 4032*zeta(7). (End)

A071910 a(n) = t(n)*t(n+1)*t(n+2), where t() are the triangular numbers.

Original entry on oeis.org

0, 18, 180, 900, 3150, 8820, 21168, 45360, 89100, 163350, 283140, 468468, 745290, 1146600, 1713600, 2496960, 3558168, 4970970, 6822900, 9216900, 12273030, 16130268, 20948400, 26910000, 34222500, 43120350, 53867268, 66758580, 82123650, 100328400, 121777920
Offset: 0

Views

Author

N. J. A. Sloane, Jun 13 2002

Keywords

Comments

a(n) is also the number of three-dimensional cage assemblies such that the assembly is not a cube. See also A052149 for the two-dimensional version and to A059827 for the non-exclusive version. - Alejandro Rodriguez, Oct 20 2020

Crossrefs

Cf. A006542, (first differences of a(n) /18) A006414, (second differences of a(n) /18) A006322, (third differences of a(n) /18) A004068, (fourth differences of a(n) /18) A005891, (fifth differences of a(n) /18) A008706.

Programs

  • Mathematica
    Join[{0},Times@@@Partition[Accumulate[Range[40]],3,1]] (* or *) LinearRecurrence[{7,-21,35,-35,21,-7,1},{0,18,180,900,3150,8820,21168},40] (* Harvey P. Dale, Aug 08 2025 *)
  • PARI
    t(n) = n*(n+1)/2;
    a(n) = t(n)*t(n+1)*t(n+2); \\ Michel Marcus, Oct 21 2015

Formula

a(n) = 18*A006542(n+3). - Vladeta Jovovic, Jun 14 2002
G.f.: 18*x*(1+3*x+x^2)/(1-x)^7. - Vladeta Jovovic, Jun 14 2002
a(n) = ((n+1)*(n+2))^3/8 - Sum_{i=1..n+1} i^3. - Jon Perry, Feb 13 2004
a(n) = C(2+n, n)*C(3+n, 1+n)*C(4+n, 2+n). - Zerinvary Lajos, Jul 29 2005
a(n) = A059827(n+1) - A000537(n+1). - Michel Marcus, Oct 21 2015
Showing 1-10 of 12 results. Next