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.

Previous Showing 31-40 of 396 results. Next

A001318 Generalized pentagonal numbers: m*(3*m - 1)/2, m = 0, +-1, +-2, +-3, ....

Original entry on oeis.org

0, 1, 2, 5, 7, 12, 15, 22, 26, 35, 40, 51, 57, 70, 77, 92, 100, 117, 126, 145, 155, 176, 187, 210, 222, 247, 260, 287, 301, 330, 345, 376, 392, 425, 442, 477, 495, 532, 551, 590, 610, 651, 672, 715, 737, 782, 805, 852, 876, 925, 950, 1001, 1027, 1080, 1107, 1162, 1190, 1247, 1276, 1335
Offset: 0

Views

Author

Keywords

Comments

Partial sums of A026741. - Jud McCranie; corrected by Omar E. Pol, Jul 05 2012
From R. K. Guy, Dec 28 2005: (Start)
"Conway's relation twixt the triangular and pentagonal numbers: Divide the triangular numbers by 3 (when you can exactly):
0 1 3 6 10 15 21 28 36 45 55 66 78 91 105 120 136 153 ...
0 - 1 2 .- .5 .7 .- 12 15 .- 22 26 .- .35 .40 .- ..51 ...
.....-.-.....+..+.....-..-.....+..+......-...-.......+....
"and you get the pentagonal numbers in pairs, one of positive rank and the other negative.
"Append signs according as the pair have the same (+) or opposite (-) parity.
"Then Euler's pentagonal number theorem is easy to remember:
"p(n-0) - p(n-1) - p(n-2) + p(n-5) + p(n-7) - p(n-12) - p(n-15) ++-- = 0^n
where p(n) is the partition function, the left side terminates before the argument becomes negative and 0^n = 1 if n = 0 and = 0 if n > 0.
"E.g. p(0) = 1, p(7) = p(7-1) + p(7-2) - p(7-5) - p(7-7) + 0^7 = 11 + 7 - 2 - 1 + 0 = 15."
(End)
The sequence may be used in order to compute sigma(n), as described in Euler's article. - Thomas Baruchel, Nov 19 2003
Number of levels in the partitions of n + 1 with parts in {1,2}.
a(n) is the number of 3 X 3 matrices (symmetrical about each diagonal) M = {{a, b, c}, {b, d, b}, {c, b, a}} such that a + b + c = b + d + b = n + 2, a,b,c,d natural numbers; example: a(3) = 5 because (a,b,c,d) = (2,2,1,1), (1,2,2,1), (1,1,3,3), (3,1,1,3), (2,1,2,3). - Philippe Deléham, Apr 11 2007
Also numbers a(n) such that 24*a(n) + 1 = (6*m - 1)^2 are odd squares: 1, 25, 49, 121, 169, 289, 361, ..., m = 0, +-1, +-2, ... . - Zak Seidov, Mar 08 2008
From Matthew Vandermast, Oct 28 2008: (Start)
Numbers n for which A000326(n) is a member of A000332. Cf. A145920.
This sequence contains all members of A000332 and all nonnegative members of A145919. For values of n such that n*(3*n - 1)/2 belongs to A000332, see A145919. (End)
Starting with offset 1 = row sums of triangle A168258. - Gary W. Adamson, Nov 21 2009
Starting with offset 1 = Triangle A101688 * [1, 2, 3, ...]. - Gary W. Adamson, Nov 27 2009
Starting with offset 1 can be considered the first in an infinite set generated from A026741. Refer to the array in A175005. - Gary W. Adamson, Apr 03 2010
Vertex number of a square spiral whose edges have length A026741. The two axes of the spiral forming an "X" are A000326 and A005449. The four semi-axes forming an "X" are A049452, A049453, A033570 and the numbers >= 2 of A033568. - Omar E. Pol, Sep 08 2011
A general formula for the generalized k-gonal numbers is given by n*((k - 2)*n - k + 4)/2, n=0, +-1, +-2, ..., k >= 5. - Omar E. Pol, Sep 15 2011
a(n) is the number of 3-tuples (w,x,y) having all terms in {0,...,n} and 2*w = 2*x + y. - Clark Kimberling, Jun 04 2012
Generalized k-gonal numbers are second k-gonal numbers and positive terms of k-gonal numbers interleaved, k >= 5. - Omar E. Pol, Aug 04 2012
a(n) is the sum of the largest parts of the partitions of n+1 into exactly 2 parts. - Wesley Ivan Hurt, Jan 26 2013
Conway's relation mentioned by R. K. Guy is a relation between triangular numbers and generalized pentagonal numbers, two sequences from different families, but as triangular numbers are also generalized hexagonal numbers in this case we have a relation between two sequences from the same family. - Omar E. Pol, Feb 01 2013
Start with the sequence of all 0's. Add n to each value of a(n) and the next n - 1 terms. The result is the generalized pentagonal numbers. - Wesley Ivan Hurt, Nov 03 2014
(6k + 1) | a(4k). (3k + 1) | a(4k+1). (3k + 2) | a(4k+2). (6k + 5) | a(4k+3). - Jon Perry, Nov 04 2014
Enge, Hart and Johansson proved: "Every generalised pentagonal number c >= 5 is the sum of a smaller one and twice a smaller one, that is, there are generalised pentagonal numbers a, b < c such that c = 2a + b." (see link theorem 5). - Peter Luschny, Aug 26 2016
The Enge, et al. result for c >= 5 also holds for c >= 2 if 0 is included as a generalized pentagonal number. That is, 2 = 2*1 + 0. - Michael Somos, Jun 02 2018
Suggestion for title, where n actually matches the list and b-file: "Generalized pentagonal numbers: k(n)*(3*k(n) - 1)/2, where k(n) = A001057(n) = [0, 1, -1, 2, -2, 3, -3, ...], n >= 0" - Daniel Forgues, Jun 09 2018 & Jun 12 2018
Generalized k-gonal numbers are the partial sums of the sequence formed by the multiples of (k - 4) and the odd numbers (A005408) interleaved, with k >= 5. - Omar E. Pol, Jul 25 2018
The last digits form a symmetric cycle of length 40 [0, 1, 2, 5, ..., 5, 2, 1, 0], i.e., a(n) == a(n + 40) (mod 10) and a(n) == a(40*k - n - 1) (mod 10), 40*k > n. - Alejandro J. Becerra Jr., Aug 14 2018
Only 2, 5, and 7 are prime. All terms are of the form k*(k+1)/6, where 3 | k or 3 | k+1. For k > 6, the value divisible by 3 must have another factor d > 2, which will remain after the division by 6. - Eric Snyder, Jun 03 2022
8*a(n) is the product of two even numbers one of which is n + n mod 2. - Peter Luschny, Jul 15 2022
a(n) is the dot product of [1, 2, 3, ..., n] and repeat[1, 1/2]. a(5) = 12 = [1, 2, 3, 4, 5] dot [1, 1/2, 1, 1/2, 1] = [1 + 1 + 3 + 2 + 5]. - Gary W. Adamson, Dec 10 2022
Every nonnegative number is the sum of four terms of this sequence [S. Realis]. - N. J. A. Sloane, May 07 2023
From Peter Bala, Jan 06 2025: (Start)
The sequence terms are the exponents in the expansions of the following infinite products:
1) Product_{n >= 1} (1 - s(n)*q^n) = 1 + q + q^2 + q^5 + q^7 + q^12 + q^15 + ..., where s(n) = (-1)^(1 + mod(n+1,3)).
2) Product_{n >= 1} (1 - q^(2*n))*(1 - q^(3*n))^2/((1 - q^n)*(1 - q^(6*n))) = 1 + q + q^2 + q^5 + q^7 + q^12 + q^15 + ....
3) Product_{n >= 1} (1 - q^n)*(1 - q^(4*n))*(1 - q^(6*n))^5/((1 - q^(2*n))*(1 - q^(3*n))*(1 - q^(12*n)))^2 = 1 - q + q^2 - q^5 - q^7 + q^12 - q^15 + q^22 + q^26 - q^35 + ....
4) Product_{n >= 1} (1 - q^(2*n))^13/((1 - (-1)^n*q^n)*(1 - q^(4*n)))^5 = 1 - 5*q + 7*q^2 - 11*q^5 + 13*q^7 - 17*q^12 + 19*q^15 - + .... See Oliver, Theorem 1.1. (End)

Examples

			G.f. = x + 2*x^2 + 5*x^3 + 7*x^4 + 12*x^5 + 15*x^6 + 22*x^7 + 26*x^8 + 35*x^9 + ...
		

References

  • Enoch Haga, A strange sequence and a brilliant discovery, chapter 5 of Exploring prime numbers on your PC and the Internet, first revised ed., 2007 (and earlier ed.), pp. 53-70.
  • Ross Honsberger, Ingenuity in Mathematics, Random House, 1970, p. 117.
  • Donald E. Knuth, The Art of Computer Programming, vol. 4A, Combinatorial Algorithms, (to appear), section 7.2.1.4, equation (18).
  • Ivan Niven and Herbert S. Zuckerman, An Introduction to the Theory of Numbers, 2nd ed., Wiley, NY, 1966, p. 231.
  • 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).

Crossrefs

Cf. A080995 (characteristic function), A026741 (first differences), A034828 (partial sums), A165211 (mod 2).
Cf. A000326 (pentagonal numbers), A005449 (second pentagonal numbers), A000217 (triangular numbers).
Indices of nonzero terms of A010815, i.e., the (zero-based) indices of 1-bits of the infinite binary word to which the terms of A068052 converge.
Union of A036498 and A036499.
Sequences of generalized k-gonal numbers: this sequence (k=5), A000217 (k=6), A085787 (k=7), A001082 (k=8), A118277 (k=9), A074377 (k=10), A195160 (k=11), A195162 (k=12), A195313 (k=13), A195818 (k=14), A277082 (k=15), A274978 (k=16), A303305 (k=17), A274979 (k=18), A303813 (k=19), A218864 (k=20), A303298 (k=21), A303299 (k=22), A303303 (k=23), A303814 (k=24), A303304 (k=25), A316724 (k=26), A316725 (k=27), A303812 (k=28), A303815 (k=29), A316729 (k=30).
Column 1 of A195152.
Squares in APs: A221671, A221672.
Quadrisection: A049453(k), A033570(k), A033568(k+1), A049452(k+1), k >= 0.
Cf. A002620.

Programs

  • GAP
    a:=[0,1,2,5];; for n in [5..60] do a[n]:=2*a[n-2]-a[n-4]+3; od; a; # Muniru A Asiru, Aug 16 2018
    
  • Haskell
    a001318 n = a001318_list !! n
    a001318_list = scanl1 (+) a026741_list -- Reinhard Zumkeller, Nov 15 2015
    
  • Magma
    [(6*n^2 + 6*n + 1 - (2*n + 1)*(-1)^n)/16 : n in [0..50]]; // Wesley Ivan Hurt, Nov 03 2014
    
  • Magma
    [(3*n^2 + 2*n + (n mod 2) * (2*n + 1)) div 8: n in [0..70]]; // Vincenzo Librandi, Nov 04 2014
    
  • Maple
    A001318 := -(1+z+z**2)/(z+1)**2/(z-1)**3; # Simon Plouffe in his 1992 dissertation; gives sequence without initial zero
    A001318 := proc(n) (6*n^2+6*n+1)/16-(2*n+1)*(-1)^n/16 ; end proc: # R. J. Mathar, Mar 27 2011
  • Mathematica
    Table[n*(n+1)/6, {n, Select[Range[0, 100], Mod[#, 3] != 1 &]}]
    Select[Accumulate[Range[0,200]]/3,IntegerQ] (* Harvey P. Dale, Oct 12 2014 *)
    CoefficientList[Series[x (1 + x + x^2) / ((1 + x)^2 (1 - x)^3), {x, 0, 70}], x] (* Vincenzo Librandi, Nov 04 2014 *)
    LinearRecurrence[{1,2,-2,-1,1},{0,1,2,5,7},70] (* Harvey P. Dale, Jun 05 2017 *)
    a[ n_] := With[{m = Quotient[n + 1, 2]}, m (3 m + (-1)^n) / 2]; (* Michael Somos, Jun 02 2018 *)
  • PARI
    {a(n) = (3*n^2 + 2*n + (n%2) * (2*n + 1)) / 8}; /* Michael Somos, Mar 24 2011 */
    
  • PARI
    {a(n) = if( n<0, n = -1-n); polcoeff( x * (1 - x^3) / ((1 - x) * (1-x^2))^2 + x * O(x^n), n)}; /* Michael Somos, Mar 24 2011 */
    
  • PARI
    {a(n) = my(m = (n+1) \ 2); m * (3*m + (-1)^n) / 2}; /* Michael Somos, Jun 02 2018 */
    
  • Python
    def a(n):
        p = n % 2
        return (n + p)*(3*n + 2 - p) >> 3
    print([a(n) for n in range(60)])  # Peter Luschny, Jul 15 2022
    
  • Python
    def A001318(n): return n*(n+1)-(m:=n>>1)*(m+1)>>1 # Chai Wah Wu, Nov 23 2024
  • Sage
    @CachedFunction
    def A001318(n):
        if n == 0 : return 0
        inc = n//2 if is_even(n) else n
        return inc + A001318(n-1)
    [A001318(n) for n in (0..59)] # Peter Luschny, Oct 13 2012
    

Formula

Euler: Product_{n>=1} (1 - x^n) = Sum_{n=-oo..oo} (-1)^n*x^(n*(3*n - 1)/2).
A080995(a(n)) = 1: complement of A090864; A000009(a(n)) = A051044(n). - Reinhard Zumkeller, Apr 22 2006
Euler transform of length-3 sequence [2, 2, -1]. - Michael Somos, Mar 24 2011
a(-1 - n) = a(n) for all n in Z. a(2*n) = A005449(n). a(2*n - 1) = A000326(n). - Michael Somos, Mar 24 2011. [The extension of the recurrence to negative indices satisfies the signature (1,2,-2,-1,1), but not the definition of the sequence m*(3*m -1)/2, because there is no m such that a(-1) = 0. - Klaus Purath, Jul 07 2021]
a(n) = 3 + 2*a(n-2) - a(n-4). - Ant King, Aug 23 2011
Product_{k>0} (1 - x^k) = Sum_{k>=0} (-1)^k * x^a(k). - Michael Somos, Mar 24 2011
G.f.: x*(1 + x + x^2)/((1 + x)^2*(1 - x)^3).
a(n) = n*(n + 1)/6 when n runs through numbers == 0 or 2 mod 3. - Barry E. Williams
a(n) = A008805(n-1) + A008805(n-2) + A008805(n-3), n > 2. - Ralf Stephan, Apr 26 2003
Sequence consists of the pentagonal numbers (A000326), followed by A000326(n) + n and then the next pentagonal number. - Jon Perry, Sep 11 2003
a(n) = (6*n^2 + 6*n + 1)/16 - (2*n + 1)*(-1)^n/16; a(n) = A034828(n+1) - A034828(n). - Paul Barry, May 13 2005
a(n) = Sum_{k=1..floor((n+1)/2)} (n - k + 1). - Paul Barry, Sep 07 2005
a(n) = A000217(n) - A000217(floor(n/2)). - Pierre CAMI, Dec 09 2007
If n even a(n) = a(n-1) + n/2 and if n odd a(n) = a(n-1) + n, n >= 2. - Pierre CAMI, Dec 09 2007
a(n)-a(n-1) = A026741(n) and it follows that the difference between consecutive terms is equal to n if n is odd and to n/2 if n is even. Hence this is a self-generating sequence that can be simply constructed from knowledge of the first term alone. - Ant King, Sep 26 2011
a(n) = (1/2)*ceiling(n/2)*ceiling((3*n + 1)/2). - Mircea Merca, Jul 13 2012
a(n) = (A008794(n+1) + A000217(n))/2 = A002378(n) - A085787(n). - Omar E. Pol, Jan 12 2013
a(n) = floor((n + 1)/2)*((n + 1) - (1/2)*floor((n + 1)/2) - 1/2). - Wesley Ivan Hurt, Jan 26 2013
From Oskar Wieland, Apr 10 2013: (Start)
a(n) = a(n+1) - A026741(n),
a(n) = a(n+2) - A001651(n),
a(n) = a(n+3) - A184418(n),
a(n) = a(n+4) - A007310(n),
a(n) = a(n+6) - A001651(n)*3 = a(n+6) - A016051(n),
a(n) = a(n+8) - A007310(n)*2 = a(n+8) - A091999(n),
a(n) = a(n+10)- A001651(n)*5 = a(n+10)- A072703(n),
a(n) = a(n+12)- A007310(n)*3,
a(n) = a(n+14)- A001651(n)*7. (End)
a(n) = (A007310(n+1)^2 - 1)/24. - Richard R. Forberg, May 27 2013; corrected by Zak Seidov, Mar 14 2015; further corrected by Jianing Song, Oct 24 2018
a(n) = Sum_{i = ceiling((n+1)/2)..n} i. - Wesley Ivan Hurt, Jun 08 2013
G.f.: x*G(0), where G(k) = 1 + x*(3*k + 4)/(3*k + 2 - x*(3*k + 2)*(3*k^2 + 11*k + 10)/(x*(3*k^2 + 11*k + 10) + (k + 1)*(3*k + 4)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 16 2013
Sum_{n>=1} 1/a(n) = 6 - 2*Pi/sqrt(3). - Vaclav Kotesovec, Oct 05 2016
a(n) = Sum_{i=1..n} numerator(i/2) = Sum_{i=1..n} denominator(2/i). - Wesley Ivan Hurt, Feb 26 2017
a(n) = A000292(A001651(n))/A001651(n), for n>0. - Ivan N. Ianakiev, May 08 2018
a(n) = ((-5 + (-1)^n - 6n)*(-1 + (-1)^n - 6n))/96. - José de Jesús Camacho Medina, Jun 12 2018
a(n) = Sum_{k=1..n} k/gcd(k,2). - Pedro Caceres, Apr 23 2019
Quadrisection. For r = 0,1,2,3: a(r + 4*k) = 6*k^2 + sqrt(24*a(r) + 1)*k + a(r), for k >= 1, with inputs (k = 0) {0,1,2,5}. These are the sequences A049453(k), A033570(k), A033568(k+1), A049452(k+1), for k >= 0, respectively. - Wolfdieter Lang, Feb 12 2021
a(n) = a(n-4) + sqrt(24*a(n-2) + 1), n >= 4. - Klaus Purath, Jul 07 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = 6*(log(3)-1). - Amiram Eldar, Feb 28 2022
a(n) = A002620(n) + A008805(n-1). Gary W. Adamson, Dec 10 2022
E.g.f.: (x*(7 + 3*x)*cosh(x) + (1 + 5*x + 3*x^2)*sinh(x))/8. - Stefano Spezia, Aug 01 2024

A211795 Number of (w,x,y,z) with all terms in {1,...,n} and w*x < 2*y*z.

Original entry on oeis.org

0, 1, 11, 58, 177, 437, 894, 1659, 2813, 4502, 6836, 10008, 14121, 19449, 26117, 34372, 44422, 56597, 71044, 88160, 108115, 131328, 158074, 188773, 223604, 263172, 307719, 357715, 413493, 475690, 544480, 620632, 704381, 796413
Offset: 0

Views

Author

Clark Kimberling, Apr 27 2012

Keywords

Comments

Each sequence in the following guide counts 4-tuples
(w,x,y,z) such that the indicated relation holds and the four numbers w,x,y,z are in {1,...,n}. The notation "m div" means that m divides every term of the sequence.
A211058 ... wx <= yz
A211787 ... wx <= 2yz
A211795 ... wx < 2yz
A211797 ... wx > 2yz
A211809 ... wx >= 2yz
A211812 ... wx <= 3yz
A211917 ... wx < 3yz
A211918 ... wx > 3yz
A211919 ... wx >= 3yz
A211920 ... 2wx < 3yz
A211921 ... 2wx <= 3yz
A211922 ... 2wx > 3yz
A211923 ... 2wx >= 3yz
A212019 ... wx = 2yz ..... 2 div
A212020 ... wx = 3yz ..... 2 div
A212021 ... 2wx = 3yz .... 2 div
A212047 ... wx = 4yz
A212048 ... 3wx = 4yz .... 2 div
A212049 ... wx = 5yz ..... 2 div
A212050 ... 2wx = 5yz .... 2 div
A212051 ... 3wx = 5yz .... 2 div
A212052 ... 4wx = 5yz .... 2 div
A209978 ... wx = yz + 1 .. 2 div
A212053 ... wx <= yz + 1
A212054 ... wx > yz + 1
A212055 ... wx <= yz + 2
A212056 ... wx > yz + 2
A197168 ... wx = yz + 2 .. 2 div
A061201 ... w = xyz
A212057 ... w < xyz
A212058 ... w >= xyz
A212059 ... w = xyz - 1
A212060 ... w = xyz - 2
A212061 ... wx = (yz)^2
A212062 ... w^2 = xyz
A212063 ... w^2 < xyz
A212064 ... w^2 >= xyz
A212065 ... w^2 <= xyz
A212066 ... w^2 > xyz
A212067 ... w^3 = xyz
A002623 ... w = 2x + y + z
A006918 ... w = 2x + 2y + z
A000601 ... w = x + 2y + 3z (except for initial 0's)
A212068 ... 2w = x + y + z
A212069 ... 3w = x + y + z (w = average{x,y,z})
A212088 ... 3w < x + y + z
A212089 ... 3w >= x + y + z
A212090 ... w < x + y + z
A000332 ... w >= x + y + z
A212145 ... w < 2x + y + z
A001752 ... w >= 2x + y + z
A001400 ... w = 2x +3y + 4z
A005900 ... w = -x + y + z
A192023 ... w = -x + y + z + 2
A212091 ... w^2 = x^2 + y^2 + z^2 ... 3 div
A212087 ... w^2 + x^2 = y^2 + z^2
A212092 ... w^2 < x^2 + y^2 + z^2
A212093 ... w^2 <= x^2 + y^2 + z^2
A212094 ... w^2 > x^2 + y^2 + z^2
A212095 ... w^2 >= x^2 + y^2 + z^2
A212096 ... w^3 = x^3 + y^3 + z^3 ... 6 div
A212097 ... w^3 < x^3 + y^3 + z^3
A212098 ... w^3 <= x^3 + y^3 + z^3
A212099 ... w^3 > x^3 + y^3 + z^3
A212100 ... w^3 >= x^3 + y^3 + z^3
A212101 ... wx^2 = yz^2
A212102 ... 1/w = 1/x + 1/y + 1/z
A212103 ... 3/w = 1/x + 1/y + 1/z; w = h.m. of {x,y,z}
A212104 ... 3/w >= 1/x + 1/y + 1/z; w >= h.m.
A212105 ... 3/w < 1/x + 1/y + 1/z; w < h.m.
A212106 ... 3/w > 1/x + 1/y + 1/z; w > h.m.
A212107 ... 3/w <= 1/x + 1/y + 1/z; w <= h.m.
A212133 ... median(w,x,y,z) = mean(w,x,y,z)
A212134 ... median(w,x,y,z) <= mean(w,x,y,z)
A212135 ... median(w,x,y,z) > mean(w,x,y,z)
A212241 ... wx + yz > n
A212243 ... 2wx + yz = n
A212244 ... w = xyz - n
A212245 ... w = xyz - 2n
A212246 ... 2w = x + y + z - n
A212247 ... 3w = x + y + z + n
A212249 ... 3w < x + y + z + n
A212250 ... 3w >= x + y + z + n
A212251 ... 3w = x + y + z + n + 1
A212252 ... 3w = x + y + z + n + 2
A212254 ... w = x + 2y + 3z - n
A212255 ... w^2 = mean(x^2, y^2, z^2)
A212256 ... 4/w = 1/x + 1/y +1/z + 1/n
In the list above, if the relation in the second column is of the form "w rel ax + by + cz" then the sequence is linearly recurrent. In the list below, the same is true for expressions involving more than one relation.
A000332 ... w < x <= y < z .... C(n,4)
A000914 ... w < x <= y < z .... Stirling 1st kind
A000914 ... w < x <= y >= z ... Stirling 1st kind
A050534 ... w < x < y >= z .... tritriangular
A001296 ... w <= x <= y >= z .. 4-dim pyramidal
A006322 ... x < x > y >= z
A002418 ... w < x >= y < z
A050534 ... w < x >=y >= z
A212415 ... w < x >= y <= z
A001296 ... w < x >= y <= z
A212246 ... w <= x > y <= z
A006322 ... w <= x >= y <= z
A212501 ... w > x < y >= z
A212503 ... w < 2x and y < 2z ..... A (note below)
A212504 ... w < 2x and y > 2z ..... A
A212505 ... w < 2x and y >= 2z .... A
A212506 ... w <= 2x and y <= 2z ... A
A212507 ... w < 2x and y <= 2z .... B
A212508 ... w < 2x and y < 3z ..... C
A212509 ... w < 2x and y <= 3z .... C
A212510 ... w < 2x and y > 3z ..... C
A212511 ... w < 2x and y >= 3z .... C
A212512 ... w <= 2x and y < 3z .... C
A212513 ... w <= 2x and y <= 3z ... C
A212514 ... w <= 2x and y > 3z .... C
A212515 ... w <= 2x and y >= 3z ... C
A212516 ... w > 2x and y < 3z ..... C
A212517 ... w > 2x and y <= 3z .... C
A212518 ... w > 2x and y > 3z ..... C
A212519 ... w > 2x and y >= 3z .... C
A212520 ... w >= 2x and y < 3z .... C
A212521 ... w >= 2x and y <= 3z ... C
A212522 ... w >= 2x and y > 3z .... C
A212523 ... w + x < y + z
A212560 ... w + x <= y + z
A212561 ... w + x = 2y + 2z
A212562 ... w + x < 2y + 2z ....... B
A212563 ... w + x <= 2y + 2z ...... B
A212564 ... w + x > 2y + 2z ....... B
A212565 ... w + x >= 2y + 2z ...... B
A212566 ... w + x = 3y + 3z
A212567 ... 2w + 2x = 3y + 3z
A212570 ... |w - x| = |x - y| + |y - z|
A212571 ... |w - x| < |x - y| + |y - z| ... B ... 4 div
A212572 ... |w - x| <= |x - y| + |y - z| .. B
A212573 ... |w - x| > |x - y| + |y - z| ... B ... 2 div
A212574 ... |w - x| >= |x - y| + |y - z| .. B
A212575 ... 2|w - x| = |x - y| + |y - z|
A212576 ... |w - x| = 2|x - y| + 2|y - z|
A212577 ... |w - x| = 2|x - y| - |y - z|
A212578 ... 2|w - x| = |x - y| - |y - z|
A212579 ... min{|w-x|,|w-y|} = min{|x-y|,|x-z|}
A212692 ... w = |x - y| + |y - z| ............... 2 div
A212568 ... w < |x - y| + |y - z| ............... 2 div
A212573 ... w <= |x - y| + |y - z| .............. 2 div
A212574 ... w > |x - y| + |y - z|
A212575 ... w >= |x - y| + |y - z|
A212676 ... w + x = |x - y| + |y - z| ......... H
A212677 ... w + y = |x - y| + |y - z|
A212678 ... w + x + y = |x - y| + |y - z|
A006918 ... w + x + y + z = |x - y| + |y - z| . H
A212679 ... |x - y| = |y - z| ................. H
A212680 ... |x - y| = |y - z| + 1 ..............H 2 div
A212681 ... |x - y| < |y - z| ................... 2 div
A212682 ... |x - y| >= |y - z|
A212683 ... |x - y| = w + |y - z| ............... 2 div
A212684 ... |x - y| = n - w + |y - z|
A212685 ... |w - x| = w + |y - z|
A186707 ... |w - x| < w + |y - z| ... (Note D)
A212714 ... |w - x| >= w + |y - z| .......... H . 2 div
A212686 ... 2*|w - x| = n + |y - z| ............. 4 div
A212687 ... 2*|w - x| < n + |y - z| ......... B
A212688 ... 2*|w - x| < n + |y - z| ......... B . 2 div
A212689 ... 2*|w - x| > n + |y - z| ......... B . 2 div
A212690 ... 2*|w - x| <= n + |y - z| ........ B
A212691 ... w + |x - y| = |x - z| + |y - z| . E . 2 div
...
In the above lists, all the terms of (w,x,y,z) are in {1,...,n}, but in the next lists they are all in {0,...,n}, and sequences are all linearly recurrent.
R=range{w,x,y,z}=max{w,x,y,z}-min{w,x,y,z}.
A212740 ... max{w,x,y,z} < 2*min{w,x,y,z} .... A
A212741 ... max{w,x,y,z} >= 2*min{w,x,y,z} ... A
A212742 ... max{w,x,y,z} <= 2*min{w,x,y,z} ... A
A212743 ... max{w,x,y,z} > 2*min{w,x,y,z} .... A . 2 div
A212744 ... w=range (=max-min) ............... E
A212745 ... w=max{w,x,y,z} - 2*min{w,x,y,z}
A212746 ... R is in {w,x,y,z} ................ E
A212569 ... R is not in {w,x,y,z} ............ E
A212749 ... w=R or x
A212750 ... w=R or x=R or y
A212751 ... w=R or x=R or y
A212752 ... wR ......... A
A212753 ... wR or z>R ......... D
A212754 ... wR or y>R or z>R ......... D
A002415 ... w = x + R ........................ D
A212755 ... |w - x| = R ...................... D
A212756 ... 2w = x + R
A212757 ... 2w = R
A212758 ... w = floor(R/2)
A002413 ... w = floor((x+y+z/2))
A212759 ... w, x, y are even
A212760 ... w is even and x = y + z .......... E
A212761 ... w is odd and x and y are even .... F . 2 div
A212762 ... w and x are odd y is even ........ F . 2 div
A212763 ... w, x, y are odd .................. F
A212764 ... w, x, y are even and z is odd .... F
A030179 ... w and x are even and y and z odd
A212765 ... w is even and x,y,z are odd ...... F
A212766 ... w is even and x is odd ........... A . 2 div
A212767 ... w and x are even and w+x=y+z ..... E
A212889 ... R is even ........................ A
A212890 ... R is odd ......................... A . 2 div
A212742 ... w-x, x-y, y-z are all even ....... A
A212892 ... w-x, x-y, y-z are all odd ........ A
A212893 ... w-x, x-y, y-z have same parity ... A
A005915 ... min{|w-x|, |x-y|, |y-z|} = 0
A212894 ... min{|w-x|, |x-y|, |y-z|} = 1
A212895 ... min{|w-x|, |x-y|, |y-z|} = 2
A179824 ... min{|w-x|, |x-y|, |y-z|} > 0
A212896 ... min{|w-x|, |x-y|, |y-z|} <= 1
A212897 ... min{|w-x|, |x-y|, |y-z|} > 1
A212898 ... min{|w-x|, |x-y|, |y-z|} <= 2
A212899 ... min{|w-x|, |x-y|, |y-z|} > 2
A212901 ... |w-x| = |x-y| = |y-z|
A212900 ... |w-x|, |x-y|, |y-z| are distinct . G
A212902 ... |w-x| < |x-y| < |y-z| ............ G
A212903 ... |w-x| <= |x-y| <= |y-z| .......... G
A212904 ... |w-x| + |x-y| + |y-z| = n ........ H
A212905 ... |w-x| + |x-y| + |y-z| = 2n ....... H
...
Note A: A212503-A212506 (and others) have these recurrence coefficients: 2,2,-6,0,6,-2,-2,1.
B: 3,-1,-5,5,1,-3,1
C: 0,2,2,-1,-4,0,2,0,-2,0,4,1,-2,-2,0,1
D: 4,-5,0,5,-4,1
E: 1,3,-3,-3,3,1,-1
F: 1,4,-4,-6,6,4,-4,-1,1
G: 2,1,-3,-1,1,3,-1,-2,1
H: 2,1,-4,1,2,-1

Examples

			a(2)=11 counts these (w,x,y,z): (1,1,1,1), (1,1,1,2), (1,1,2,1), (2,1,2,1), (2,1,1,2), (1,2,2,1), (1,2,1,2), (1,1,2,2), (1,2,2,2), (2,1,2,2), (2,2,2,2).
		

References

  • A. Barvinok, Lattice Points and Lattice Polytopes, Chapter 7 in Handbook of Discrete and Computational Geometry, CRC Press, 1997, 133-152.
  • P. Gritzmann and J. M. Wills, Lattice Points, Chapter 3.2 in Handbook of Convex Geometry, vol. B, North-Holland, 1993, 765-797.

Crossrefs

Programs

  • Mathematica
    t = Compile[{{n, _Integer}}, Module[{s = 0},
        (Do[If[w*x < 2 y*z, s = s + 1], {w, 1, #},
          {x, 1, #}, {y, 1, #}, {z, 1, #}] &[n]; s)]];
    Map[t[#] &, Range[0, 40]] (* A211795 *)
    (* Peter J. C. Moses, Apr 13 2012 *)

Formula

a(n) = n^4 - A211809(n).

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

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

A000389 Binomial coefficients C(n,5).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 6, 21, 56, 126, 252, 462, 792, 1287, 2002, 3003, 4368, 6188, 8568, 11628, 15504, 20349, 26334, 33649, 42504, 53130, 65780, 80730, 98280, 118755, 142506, 169911, 201376, 237336, 278256, 324632, 376992, 435897, 501942, 575757, 658008, 749398
Offset: 0

Keywords

Comments

a(n+4) is the number of inequivalent ways of coloring the vertices of a regular 4-dimensional simplex with n colors, under the full symmetric group S_5 of order 120, with cycle index (x1^5 + 10*x1^3*x2 + 20*x1^2*x3 + 15*x1*x2^2 + 30*x1*x4 + 20*x2*x3 + 24*x5)/120.
Figurate numbers based on 5-dimensional regular simplex. According to Hyun Kwang Kim, it appears that every nonnegative integer can be represented as the sum of g = 10 of these 5-simplex(n) numbers (compared with g=3 for triangular numbers, g=5 for tetrahedral numbers and g=8 for pentatope numbers). - Jonathan Vos Post, Nov 28 2004
The convolution of the nonnegative integers (A001477) with the tetrahedral numbers (A000292), which are the convolution of the nonnegative integers with themselves (making appropriate allowances for offsets of all sequences). - Graeme McRae, Jun 07 2006
a(n) is the number of terms in the expansion of (a_1 + a_2 + a_3 + a_4 + a_5 + a_6)^n. - Sergio Falcon, Feb 12 2007
Product of five consecutive numbers divided by 120. - Artur Jasinski, Dec 02 2007
Equals binomial transform of [1, 5, 10, 10, 5, 1, 0, 0, 0, ...]. - Gary W. Adamson, Feb 02 2009
Equals INVERTi transform of A099242 (1, 7, 34, 153, 686, 3088, ...). - Gary W. Adamson, Feb 02 2009
For a team with n basketball players (n>=5), this sequence is the number of possible starting lineups of 5 players, without regard to the positions (center, forward, guard) of the players. - Mohammad K. Azarian, Sep 10 2009
a(n) is the number of different patterns, regardless of order, when throwing (n-5) 6-sided dice. For example, one die can display the 6 numbers 1, 2, ..., 6; two dice can display the 21 digit-pairs 11, 12, ..., 56, 66. - Ian Duff, Nov 16 2009
Sum of the first n pentatope numbers (1, 5, 15, 35, 70, 126, 210, ...), see A000332. - Paul Muljadi, Dec 16 2009
Sum_{n>=0} a(n)/n! = e/120. Sum_{n>=4} a(n)/(n-4)! = 501*e/120. See A067764 regarding the second ratio. - Richard R. Forberg, Dec 26 2013
For a set of integers {1,2,...,n}, a(n) is the sum of the 2 smallest elements of each subset with 4 elements, which is 3*C(n+1,5) (for n>=4), hence a(n) = 3*C(n+1,5) = 3*A000389(n+1). - Serhat Bulut, Mar 11 2015
a(n) = fallfac(n,5)/5! is also the number of independent components of an antisymmetric tensor of rank 5 and dimension n >= 1. Here fallfac is the falling factorial. - Wolfdieter Lang, Dec 10 2015
Number of compositions (ordered partitions) of n+1 into exactly 6 parts. - Juergen Will, Jan 02 2016
Number of weak compositions (ordered weak partitions) of n-5 into exactly 6 parts. - Juergen Will, Jan 02 2016
a(n+3) could be the general number of all geodetic graphs of diameter n>=2 homeomorphic to the Petersen Graph. - Carlos Enrique Frasser, May 24 2018
From Robert A. Russell, Dec 24 2020: (Start)
a(n) is the number of chiral pairs of colorings of the 5 tetrahedral facets (or vertices) of the regular 4-D simplex (5-cell, pentachoron, Schläfli symbol {3,3,3}) using subsets of a set of n colors. Each member of a chiral pair is a reflection but not a rotation of the other.
a(n+4) is the number of unoriented colorings of the 5 tetrahedral facets of the regular 4-D simplex (5-cell, pentachoron) using subsets of a set of n colors. Each chiral pair is counted as one when enumerating unoriented arrangements. (End)
For integer m and positive integer r >= 4, the polynomial a(n) + a(n + m) + a(n + 2*m) + ... + a(n + r*m) in n has its zeros on the vertical line Re(n) = (4 - r*m)/2 in the complex plane. - Peter Bala, Jun 02 2024

Examples

			G.f. = x^5 + 6*x^6 + 21*x^7 + 56*x^8 + 126*x^9 + 252*x^10 + 462*x^11 + ...
For A={1,2,3,4}, the only subset with 4 elements is {1,2,3,4}; sum of 2 minimum elements of this subset: a(4) = 1+2 = 3 = 3*C(4+1,5).
For A={1,2,3,4,5}, the subsets with 4 elements are {1,2,3,4}, {1,2,3,5}, {1,2,4,5}, {1,3,4,5}, {2,3,4,5}; sum of 2 smallest elements of each subset: a(5) = (1+2)+(1+2)+(1+2)+(1+3)+(2+3) = 18 = 3*C(5+1,5). - _Serhat Bulut_, Mar 11 2015
a(6) = 6 from the six independent components of an antisymmetric tensor A of rank 5 and dimension 6: A(1,2,3,4,5), A(1,2,3,4,6), A(1,2,3,5,6), A(1,2,4,5,6), A(1,3,4,5,6), A(2,3,4,5,6). See the Dec 10 2015 comment. - _Wolfdieter Lang_, Dec 10 2015
		

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. 828.
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 196.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 7.
  • Gupta, Hansraj; Partitions of j-partite numbers into twelve or a smaller number of parts. Collection of articles dedicated to Professor P. L. Bhatnagar on his sixtieth birthday. Math. Student 40 (1972), 401-441 (1974).
  • J. C. P. Miller, editor, Table of Binomial Coefficients. Royal Society Mathematical Tables, Vol. 3, Cambridge Univ. Press, 1954.
  • 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).

Crossrefs

Cf. A099242. - Gary W. Adamson, Feb 02 2009
Cf. A242023. A104712 (fourth column, k=5).
5-cell colorings: A337895 (oriented), A132366(n-1) (achiral).
Unoriented colorings: A063843 (5-cell edges, faces), A128767 (8-cell vertices, 16-cell facets), A337957 (16-cell vertices, 8-cell facets), A338949 (24-cell), A338965 (600-cell vertices, 120-cell facets).
Chiral colorings: A331352 (5-cell edges, faces), A337954 (8-cell vertices, 16-cell facets), A234249 (16-cell vertices, 8-cell facets), A338950 (24-cell), A338966 (600-cell vertices, 120-cell facets).

Programs

  • Haskell
    a000389 n = a000389_list !! n
    a000389_list = 0 : 0 : f [] a000217_list where
       f xs (t:ts) = (sum $ zipWith (*) xs a000217_list) : f (t:xs) ts
    -- Reinhard Zumkeller, Mar 03 2015, Apr 13 2012
    
  • Magma
    [Binomial(n, 5): n in [0..40]]; // Vincenzo Librandi, Mar 12 2015
  • Maple
    f:=n->(1/120)*(n^5-10*n^4+35*n^3-50*n^2+24*n): seq(f(n), n=0..60);
    ZL := [S, {S=Prod(B,B,B,B,B,B), B=Set(Z, 1 <= card)}, unlabeled]: seq(combstruct[count](ZL, size=n+1), n=0..42); # Zerinvary Lajos, Mar 13 2007
    A000389:=1/(z-1)**6; # Simon Plouffe, 1992 dissertation
  • Mathematica
    Table[Binomial[n, 5], {n, 5, 50}] (* Stefan Steinerberger, Apr 02 2006 *)
    CoefficientList[Series[x^5 / (1 - x)^6, {x, 0, 40}], x] (* Vincenzo Librandi, Mar 12 2015 *)
    LinearRecurrence[{6,-15,20,-15,6,-1},{0,0,0,0,0,1},50] (* Harvey P. Dale, Jul 17 2016 *)
  • PARI
    (conv(u,v)=local(w); w=vector(length(u),i,sum(j=1,i,u[j]*v[i+1-j])); w);
    (t(n)=n*(n+1)/2); u=vector(10,i,t(i)); conv(u,u)
    

Formula

G.f.: x^5/(1-x)^6.
a(n) = n*(n-1)*(n-2)*(n-3)*(n-4)/120.
a(n) = (n^5-10*n^4+35*n^3-50*n^2+24*n)/120. (Replace all x_i's in the cycle index with n.)
a(n+2) = Sum_{i+j+k=n} i*j*k. - Benoit Cloitre, Nov 01 2002
Convolution of triangular numbers (A000217) with themselves.
Partial sums of A000332. - Alexander Adamchuk, Dec 19 2004
a(n) = -A110555(n+1,5). - Reinhard Zumkeller, Jul 27 2005
a(n+3) = (1/2!)*(d^2/dx^2)S(n,x)|A049310.%20-%20_Wolfdieter%20Lang">{x=2}, n>=2, one half of second derivative of Chebyshev S-polynomials evaluated at x=2. See A049310. - _Wolfdieter Lang, Apr 04 2007
a(n) = A052787(n+5)/120. - Zerinvary Lajos, Apr 26 2007
Sum_{n>=5} 1/a(n) = 5/4. - R. J. Mathar, Jan 27 2009
For n>4, a(n) = 1/(Integral_{x=0..Pi/2} 10*(sin(x))^(2*n-9)*(cos(x))^9). - Francesco Daddi, Aug 02 2011
Sum_{n>=5} (-1)^(n + 1)/a(n) = 80*log(2) - 655/12 = 0.8684411114... - Richard R. Forberg, Aug 11 2014
a(n) = -a(4-n) for all n in Z. - Michael Somos, Oct 07 2014
0 = a(n)*(+a(n+1) + 4*a(n+2)) + a(n+1)*(-6*a(n+1) + a(n+2)) for all n in Z. - Michael Somos, Oct 07 2014
a(n) = 3*C(n+1, 5) = 3*A000389(n+1). - Serhat Bulut, Mar 11 2015
From Ilya Gutkovskiy, Jul 23 2016: (Start)
E.g.f.: x^5*exp(x)/120.
Inverse binomial transform of A054849. (End)
From Robert A. Russell, Dec 24 2020: (Start)
a(n) = A337895(n) - a(n+4) = (A337895(n) - A132366(n-1)) / 2 = a(n+4) - A132366(n-1).
a(n+4) = A337895(n) - a(n) = (A337895(n) + A132366(n-1)) / 2 = a(n) + A132366(n-1).
a(n+4) = 1*C(n,1) + 4*C(n,2) + 6*C(n,3) + 4*C(n,4) + 1*C(n,5), where the coefficient of C(n,k) is the number of unoriented pentachoron colorings using exactly k colors. (End)

Extensions

Corrected formulas that had been based on other offsets. - R. J. Mathar, Jun 16 2009
I changed the offset to 0. This will require some further adjustments to the formulas. - N. J. A. Sloane, Aug 01 2010

A006003 a(n) = n*(n^2 + 1)/2.

Original entry on oeis.org

0, 1, 5, 15, 34, 65, 111, 175, 260, 369, 505, 671, 870, 1105, 1379, 1695, 2056, 2465, 2925, 3439, 4010, 4641, 5335, 6095, 6924, 7825, 8801, 9855, 10990, 12209, 13515, 14911, 16400, 17985, 19669, 21455, 23346, 25345, 27455, 29679, 32020, 34481, 37065, 39775
Offset: 0

Keywords

Comments

Write the natural numbers in groups: 1; 2,3; 4,5,6; 7,8,9,10; ... and add the groups. In other words, "sum of the next n natural numbers". - Felice Russo
Number of rhombi in an n X n rhombus, if 'crossformed' rhombi are allowed. - Matti De Craene (Matti.DeCraene(AT)rug.ac.be), May 14 2000
Also the sum of the integers between T(n-1)+1 and T(n), the n-th triangular number (A000217). Sum of n-th row of A000027 regarded as a triangular array.
Unlike the cubes which have a similar definition, it is possible for 2 terms of this sequence to sum to a third. E.g., a(36) + a(37) = 23346 + 25345 = 48691 = a(46). Might be called 2nd-order triangular numbers, thus defining 3rd-order triangular numbers (A027441) as n(n^3+1)/2, etc. - Jon Perry, Jan 14 2004
Also as a(n)=(1/6)*(3*n^3+3*n), n > 0: structured trigonal diamond numbers (vertex structure 4) (cf. A000330 = alternate vertex; A000447 = structured diamonds; A100145 for more on structured numbers). - James A. Record (james.record(AT)gmail.com), Nov 07 2004
The sequence M(n) of magic constants for n X n magic squares (numbered 1 through n^2) from n=3 begins M(n) = 15, 34, 65, 111, 175, 260, ... - Lekraj Beedassy, Apr 16 2005 [comment corrected by Colin Hall, Sep 11 2009]
The sequence Q(n) of magic constants for the n-queens problem in chess begins 0, 0, 0, 0, 34, 65, 111, 175, 260, ... - Paul Muljadi, Aug 23 2005
Alternate terms of A057587. - Jeremy Gardiner, Apr 10 2005
Also partial differences of A063488(n) = (2*n-1)*(n^2-n+2)/2. a(n) = A063488(n) - A063488(n-1) for n>1. - Alexander Adamchuk, Jun 03 2006
In an n X n grid of numbers from 1 to n^2, select -- in any manner -- one number from each row and column. Sum the selected numbers. The sum is independent of the choices and is equal to the n-th term of this sequence. - F.-J. Papp (fjpapp(AT)umich.edu), Jun 06 2006
Nonnegative X values of solutions to the equation (X-Y)^3 - (X+Y) = 0. To find Y values: b(n) = (n^3-n)/2. - Mohamed Bouhamida, May 16 2006
For the equation: m*(X-Y)^k - (X+Y) = 0 with X >= Y, k >= 2 and m is an odd number the X values are given by the sequence defined by a(n) = (m*n^k+n)/2. The Y values are given by the sequence defined by b(n) = (m*n^k-n)/2. - Mohamed Bouhamida, May 16 2006
If X is an n-set and Y a fixed 3-subset of X then a(n-3) is equal to the number of 4-subsets of X intersecting Y. - Milan Janjic, Jul 30 2007
(m*(2n)^k+n, m*(2n)^k-n) solves the Diophantine equation: 2m*(X-Y)^k - (X+Y) = 0 with X >= Y, k >= 2 where m is a positive integer. - Mohamed Bouhamida, Oct 02 2007
Also c^(1/2) in a^(1/2) + b^(1/2) = c^(1/2) such that a^2 + b = c. - Cino Hilliard, Feb 09 2008
a(n) = n*A000217(n) - Sum_{i=0..n-1} A001477(i). - Bruno Berselli, Apr 25 2010
a(n) is the number of triples (w,x,y) having all terms in {0,...,n} such that at least one of these inequalities fails: x+y < w, y+w < x, w+x < y. - Clark Kimberling, Jun 14 2012
Sum of n-th row of the triangle in A209297. - Reinhard Zumkeller, Jan 19 2013
The sequence starting with "1" is the third partial sum of (1, 2, 3, 3, 3, ...). - Gary W. Adamson, Sep 11 2015
a(n) is the largest eigenvalue of the matrix returned by the MATLAB command magic(n) for n > 0. - Altug Alkan, Nov 10 2015
a(n) is the number of triples (x,y,z) having all terms in {1,...,n} such that all these triangle inequalities are satisfied: x+y > z, y+z > x, z+x > y. - Heinz Dabrock, Jun 03 2016
Shares its digital root with the stella octangula numbers (A007588). See A267017. - Peter M. Chema, Aug 28 2016
Can be proved to be the number of nonnegative solutions of a system of three linear Diophantine equations for n >= 0 even: 2*a_{11} + a_{12} + a_{13} = n, 2*a_{22} + a_{12} + a_{23} = n and 2*a_{33} + a_{13} + a_{23} = n. The number of solutions is f(n) = (1/16)*(n+2)*(n^2 + 4n + 8) and a(n) = n*(n^2 + 1)/2 is obtained by remapping n -> 2*n-2. - Kamil Bradler, Oct 11 2016
For n > 0, a(n) coincides with the trace of the matrix formed by writing the numbers 1...n^2 back and forth along the antidiagonals (proved, see A078475 for the examples of matrix). - Stefano Spezia, Aug 07 2018
The trace of an n X n square matrix where the elements are entered on the ascending antidiagonals. The determinant is A069480. - Robert G. Wilson v, Aug 07 2018
Bisections are A317297 and A005917. - Omar E. Pol, Sep 01 2018
Number of achiral colorings of the vertices (or faces) of a regular tetrahedron with n available colors. An achiral coloring is identical to its reflection. - Robert A. Russell, Jan 22 2020
a(n) is the n-th centered triangular pyramidal number. - Lechoslaw Ratajczak, Nov 02 2021
a(n) is the number of words of length n defined on 4 letters {b,c,d,e} that contain one or no b's, one c or two d's, and any number of e's. For example, a(3) = 15 since the words are (number of permutations in parentheses): bce (6), bdd (3), cee (3), and dde (3). - Enrique Navarrete, Jun 21 2025

Examples

			G.f. = x + 5*x^2 + 15*x^3 + 34*x^4 + 65*x^5 + 111*x^6 + 175*x^7 + 260*x^8 + ...
For a(2)=5, the five tetrahedra have faces AAAA, AAAB, AABB, ABBB, and BBBB with colors A and B. - _Robert A. Russell_, Jan 31 2020
		

References

  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 15, p. 5, Ellipses, Paris 2008.
  • F.-J. Papp, Colloquium Talk, Department of Mathematics, University of Michigan-Dearborn, March 6, 2005.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000330, A000537, A066886, A057587, A027480, A002817 (partial sums).
Cf. A000578 (cubes).
(1/12)*t*(n^3-n)+n for t = 2, 4, 6, ... gives A004006, A006527, this sequence, A005900, A004068, A000578, A004126, A000447, A004188, A004466, A004467, A007588, A062025, A063521, A063522, A063523.
Antidiagonal sums of array in A000027. Row sums of the triangular view of A000027.
Cf. A063488 (sum of two consecutive terms), A005917 (bisection), A317297 (bisection).
Cf. A105374 / 8.
Tetrahedron colorings: A006008 (oriented), A000332(n+3) (unoriented), A000332 (chiral), A037270 (edges).
Other polyhedron colorings: A337898 (cube faces, octahedron vertices), A337897 (octahedron faces, cube vertices), A337962 (dodecahedron faces, icosahedron vertices), A337960 (icosahedron faces, dodecahedron vertices).
Row 3 of A325001 (simplex vertices and facets) and A337886 (simplex faces and peaks).

Programs

  • GAP
    a_n:=List([0..nmax], n->n*(n^2 + 1)/2); # Stefano Spezia, Aug 12 2018
    
  • Haskell
    a006003 n = n * (n ^ 2 + 1) `div` 2
    a006003_list = scanl (+) 0 a005448_list
    -- Reinhard Zumkeller, Jun 20 2013
    
  • MATLAB
    % Also works with FreeMat.
    for(n=0:nmax); tm=n*(n^2 + 1)/2; fprintf('%d\t%0.f\n', n, tm); end
    % Stefano Spezia, Aug 12 2018
    
  • Magma
    [n*(n^2 + 1)/2 : n in [0..50]]; // Wesley Ivan Hurt, Sep 11 2015
    
  • Magma
    [Binomial(n,3)+Binomial(n-1,3)+Binomial(n-2,3): n in [2..60]]; // Vincenzo Librandi, Sep 12 2015
    
  • Mathematica
    Table[ n(n^2 + 1)/2, {n, 0, 45}]
    LinearRecurrence[{4,-6,4,-1}, {0,1,5,15},50] (* Harvey P. Dale, May 16 2012 *)
    CoefficientList[Series[x (1 + x + x^2)/(x - 1)^4, {x, 0, 45}], x] (* Vincenzo Librandi, Sep 12 2015 *)
    With[{n=50},Total/@TakeList[Range[(n(n^2+1))/2],Range[0,n]]] (* Requires Mathematica version 11 or later *) (* Harvey P. Dale, Nov 28 2017 *)
  • Maxima
    a(n):=n*(n^2 + 1)/2$ makelist(a(n), n, 0, nmax); /* Stefano Spezia, Aug 12 2018 */
    
  • PARI
    {a(n) = n * (n^2 + 1) / 2}; /* Michael Somos, Dec 24 2011 */
    
  • PARI
    concat(0, Vec(x*(1+x+x^2)/(x-1)^4 + O(x^20))) \\ Felix Fröhlich, Oct 11 2016
    
  • Python
    def A006003(n): return n*(n**2+1)>>1 # Chai Wah Wu, Mar 25 2024

Formula

a(n) = binomial(n+2, 3) + binomial(n+1, 3) + binomial(n, 3). [corrected by Michel Marcus, Jan 22 2020]
G.f.: x*(1+x+x^2)/(x-1)^4. - Floor van Lamoen, Feb 11 2002
Partial sums of A005448. - Jonathan Vos Post, Mar 16 2006
Binomial transform of [1, 4, 6, 3, 0, 0, 0, ...] = (1, 5, 15, 34, 65, ...). - Gary W. Adamson, Aug 10 2007
a(n) = -a(-n) for all n in Z. - Michael Somos, Dec 24 2011
a(n) = Sum_{k = 1..n} A(k-1, k-1-n) where A(i, j) = i^2 + i*j + j^2 + i + j + 1. - Michael Somos, Jan 02 2012
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4), with a(0)=0, a(1)=1, a(2)=5, a(3)=15. - Harvey P. Dale, May 16 2012
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 3. - Ant King, Jun 13 2012
a(n) = A000217(n) + n*A000217(n-1). - Bruno Berselli, Jun 07 2013
a(n) = A057145(n+3,n). - Luciano Ancora, Apr 10 2015
E.g.f.: (1/2)*(2*x + 3*x^2 + x^3)*exp(x). - G. C. Greubel, Dec 18 2015; corrected by Ilya Gutkovskiy, Oct 12 2016
a(n) = T(n) + T(n-1) + T(n-2), where T means the tetrahedral numbers, A000292. - Heinz Dabrock, Jun 03 2016
From Ilya Gutkovskiy, Oct 11 2016: (Start)
Convolution of A001477 and A008486.
Convolution of A000217 and A158799.
Sum_{n>=1} 1/a(n) = H(-i) + H(i) = 1.343731971048019675756781..., where H(k) is the harmonic number, i is the imaginary unit. (End)
a(n) = A000578(n) - A135503(n). - Miquel Cerda, Dec 25 2016
Euler transform of length 3 sequence [5, 0, -1]. - Michael Somos, Dec 25 2016
a(n) = A037270(n)/n for n > 0. - Kritsada Moomuang, Dec 15 2018
a(n) = 3*A000292(n-1) + n. - Bruce J. Nicholson, Nov 23 2019
a(n) = A011863(n) - A011863(n-2). - Bruce J. Nicholson, Dec 22 2019
From Robert A. Russell, Jan 22 2020: (Start)
a(n) = C(n,1) + 3*C(n,2) + 3*C(n,3), where the coefficient of C(n,k) is the number of tetrahedron colorings using exactly k colors.
a(n) = C(n+3,4) - C(n,4).
a(n) = 2*A000332(n+3) - A006008(n) = A006008(n) - 2*A000332(n) = A000332(n+3) - A000332(n).
a(n) = A325001(3,n). (End)
From Amiram Eldar, Aug 21 2023: (Start)
Sum_{n>=1} 1/a(n) = 2 * (A248177 + A001620).
Product_{n>=2} (1 - 1/a(n)) = cosh(sqrt(7)*Pi/2)*cosech(Pi)/4.
Product_{n>=1} (1 + 1/a(n)) = cosh(sqrt(7)*Pi/2)*cosech(Pi). (End)

Extensions

Better description from Albert Rich (Albert_Rich(AT)msn.com), Mar 1997

A002294 a(n) = binomial(5*n, n)/(4*n + 1).

Original entry on oeis.org

1, 1, 5, 35, 285, 2530, 23751, 231880, 2330445, 23950355, 250543370, 2658968130, 28558343775, 309831575760, 3390416787880, 37377257159280, 414741863546285, 4628362722856425, 51912988256282175, 584909606696793885, 6617078646960613370
Offset: 0

Keywords

Comments

From Wolfdieter Lang, Sep 14 2007: (Start)
a(n), n >= 1, enumerates quintic trees (rooted, ordered, incomplete) with n vertices (including the root).
This is the Pfaff-Fuss-Catalan sequence C^{m}_n for m = 5. See the Graham et al. reference, p. 347. eq. 7.66. See also the Pólya-Szegő reference.
Also 5-Raney sequence. See the Graham et al. reference, pp. 346-347. (End)
a(n) = A258708(3*n, 2*n) for n > 0. - Reinhard Zumkeller, Jun 23 2015
Conjecturally, a(n) is the number of 4-uniform words on the alphabet [n] that avoid the patterns 231 and 221 (see the Defant and Kravitz link). - Colin Defant, Sep 26 2018
From Stillwell (1995), p. 62: "Eisenstein's Theorem. If y^5 + y = x, then y has a power series expansion y = x - x^5 + 10*x^9/2^1 - 15 * 14 * x^13/3! + 20 * 19 * 18*x^17/4! - ...." - Michael Somos, Sep 19 2019
a(n) is the total number of down steps before the first up step in all 4_1-Dyck paths of length 5*n. A 4_1-Dyck path is a lattice path with steps (1, 4), (1, -1) that starts and ends at y = 0 and stays above the line y = -1. - Sarah Selkirk, May 10 2020
Dropping the first 1 (starting from 1, 5, 35, ... with offset 1), the series reversion gives 1, -5, 15, -35, 70, ... (again offset 1), essentially A000332 and row 5 of A027555. - R. J. Mathar, Aug 17 2023
Number of rooted polyominoes composed of n hexagonal cells of the hyperbolic regular tiling with Schläfli symbol {6,oo}. A rooted polyomino has one external edge identified, and chiral pairs are counted as two. A stereographic projection of the {6,oo} tiling on the Poincaré disk can be obtained via the Christensson link. - Robert A. Russell, Jan 27 2024
This is instance k = 5 of the generalized Catalan family {C(k, n)}_{n>=0} given in a comment of A130564. - Wolfdieter Lang, Feb 05 2024

Examples

			There are a(2) = 5 quintic trees (vertex degree <= 5 and 5 possible branchings) with 2 vertices (one of them the root). Adding one more branch (one more vertex) to these five trees yields 5*5 + binomial(5,2) = 35 = a(3) such trees.
G.f. = 1 + x + 5*x^2 + 35*x^3 + 285*x^4 + 2530*x^5 + 23751*x^6 + 231880*x^7 + ...
G.f. = t + t^5 + 5*t^9 + 35*t^13 + 285*t^17 + 2530*t^21 + 23751*t^25 + 231880*t^29 + ...
		

References

  • Archiv der Mathematik u. Physik, Editor's note: "Über die Bestimmung der Anzahl der verschiedenen Arten, auf welche sich ein n-Eck durch Diagonalen in lauter m-Ecke zerlegen laesst, mit Bezug auf einige Abhandlungen der Herren Lame, Rodrigues, Binet, Catalan und Duhamel in dem Journal de Mathematiques pures et appliquees, publie par Joseph Liouville. T. III. IV.", Archiv der Mathematik u. Physik, 1 (1841), pp. 193ff; see especially p. 198.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 23.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, pp. 200, 347.
  • G. Pólya and G. Szegő, Problems and Theorems in Analysis, Springer-Verlag, Heidelberg, New York, 2 vols., 1972, Vol. 1, problem 211, p. 146 with solution on p. 348.
  • Ulrike Sattler, Decidable classes of formal power series with nice closure properties, Diplomarbeit im Fach Informatik, Univ. Erlangen - Nürnberg, Jul 27 1994.
  • 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).

Crossrefs

Cf. A001764, A002296, A258708, A346647 (binomial transform), A346665 (inverse binomial transform).
Fourth column of triangle A062993.
Polyominoes: A221184{n-1} (oriented), A004127 (unoriented), A369473 (chiral), A143546 (achiral), A002293 {5,oo}, A002295 {7,oo}.
Cf. A130564.

Programs

  • GAP
    List([0..22],n->Binomial(5*n,n)/(4*n+1)); # Muniru A Asiru, Nov 01 2018
  • Haskell
    a002294 n = a002294_list !! n
    a002294_list = [a258708 (3 * n) (2 * n) | n <- [1..]]
    -- Reinhard Zumkeller, Jun 23 2015
    
  • Magma
    [ Binomial(5*n,n)/(4*n+1): n in [0..100]]; // Vincenzo Librandi, Mar 24 2011
    
  • Maple
    seq(binomial(5*k+1,k)/(5*k+1),k=0..30); # Robert FERREOL, Apr 03 2015
    n:=30:G:=series(RootOf(g = 1+x*g^5, g),x=0,n+1):seq(coeff(G,x,k),k=0..n); # Robert FERREOL, Apr 03 2015
  • Mathematica
    CoefficientList[InverseSeries[ Series[ y - y^5, {y, 0, 100}], x], x][[Range[2, 100, 4]]]
    Table[Binomial[5n,n]/(4n+1),{n,0,20}] (* Harvey P. Dale, Dec 30 2011 *)
    a[ n_] := SeriesCoefficient[ HypergeometricPFQ[ {1, 2, 3, 4}/5, {2, 3, 5}/4, x 5^5/4^4], {x, 0, n}]; (* Michael Somos, May 06 2015 *)
    a[ n_] := With[{m = 4 n + 1}, SeriesCoefficient[ InverseSeries @ Series[ x - x^5, {x, 0, m}], {x, 0, m}]]; (* Michael Somos, May 06 2015 *)
  • PARI
    {a(n) = binomial( 5 * n, n) / (4*n + 1)}; /* Michael Somos, Mar 17 2011 */
    
  • PARI
    {a(n) = if( n<0, 0, n = 4*n + 1; polcoeff( serreverse( x - x^5 + x * O(x^n) ), n))}; /* Michael Somos, Mar 17 2011 */
    

Formula

For the connection with the solution of the quintic, hypergeometric series, and Lagrange inversion, see Beukers (2014). - N. J. A. Sloane, Mar 12 2014
G.f.: hypergeometric([1, 2, 3, 4] / 5, [2, 3, 5] / 4, x * 5^5 / 4^4). - Michael Somos, Mar 17 2011
O.g.f. A(x) satisfies A(x) = 1 + x * A(x)^5 = 1 / (1 - x * A(x)^4).
Given g.f. A(x) then z = t * A(t^4) satisfies 0 = z^5 - z + t. - Michael Somos, Mar 17 2011
a(n) = binomial(5*n, n - 1)/n, n >= 1, a(0) = 1. From the Lagrange series of the o.g.f. A(x) with its above given implicit equation.
a(n) = upper left term in M^n, M = the production matrix:
1, 1;
4, 4, 1;
10, 10, 4, 1;
20, 20, 10, 4, 1;
...
where (1, 4, 10, 20, ...) is the tetrahedral sequence, A000292. - Gary W. Adamson, Jul 08 2011
D-finite with recurrence: 8*n*(4*n+1)*(2*n-1)*(4*n-1)*a(n) - 5*(5*n-4)*(5*n-3)*(5*n-2)*(5*n-1)*a(n-1) = 0. - R. J. Mathar, Dec 02 2014
a(n) = binomial(5*n + 1, n)/(5*n + 1) = A062993(n+3,3). - Robert FERREOL, Apr 03 2015
a(0) = 1; a(n) = Sum_{i1 + i2 + ... + i5 = n - 1} a(i1) * a(i2) * ... *a(i5) for n >= 1. - Robert FERREOL, Apr 03 2015
From Ilya Gutkovskiy, Jan 15 2017: (Start)
O.g.f.: 5F4([1/5, 2/5, 3/5, 4/5, 1]; [1/2, 3/4, 1, 5/4]; 3125*x/256).[Cancellation of the 1s, see G.f. the above. - Wolfdieter Lang, Feb 05 2024]
E.g.f.: 4F4([1/5, 2/5, 3/5, 4/5]; [1/2, 3/4, 1, 5/4]; 3125*x/256).
a(n) ~ 5^(5*n + 1/2)/(sqrt(Pi) * 2^(8*n + 7/2) * n^(3/2)). (End)
x*A'(x)/A(x) = (A(x) - 1)/(- 4*A(x) + 5) = x + 9*x^2 + 91*x^3 + 969*x^4 + ... is the o.g.f. of A163456. Cf. A001764 and A002293 - A002296. - Peter Bala, Feb 04 2022
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^9). - Seiichi Manyama, Jun 16 2025

Extensions

More terms from Olivier Gérard, Jul 05 2001

A045943 Triangular matchstick numbers: a(n) = 3*n*(n+1)/2.

Original entry on oeis.org

0, 3, 9, 18, 30, 45, 63, 84, 108, 135, 165, 198, 234, 273, 315, 360, 408, 459, 513, 570, 630, 693, 759, 828, 900, 975, 1053, 1134, 1218, 1305, 1395, 1488, 1584, 1683, 1785, 1890, 1998, 2109, 2223, 2340, 2460, 2583, 2709, 2838, 2970, 3105, 3243, 3384, 3528
Offset: 0

Author

Keywords

Comments

Also, 3 times triangular numbers, a(n) = 3*A000217(n).
In the 24-bit RGB color cube, the number of color-lattice-points in r+g+b = n planes at n < 256 equals the triangular numbers. For n = 256, ..., 765 the number of legitimate color partitions is less than A000217(n) because {r,g,b} components cannot exceed 255. For n = 256, ..., 511, the number of non-color partitions are computable with A045943(n-255), while for n = 512, ..., 765, the number of color points in r+g+b planes equals A000217(765-n). - Labos Elemer, Jun 20 2005
If a 3-set Y and an (n-3)-set Z are disjoint subsets of an n-set X then a(n-3) is the number of 3-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007
a(n) is also the smallest number that may be written both as the sum of n-1 consecutive positive integers and n consecutive positive integers. - Claudio Meller, Oct 08 2010
For n >= 3, a(n) equals 4^(2+n)*Pi^(1 - n) times the coefficient of zeta(3) in the following integral with upper bound Pi/4 and lower bound 0: int x^(n+1) tan x dx. - John M. Campbell, Jul 17 2011
The difference a(n)-a(n-1) = 3*n, for n >= 1. - Stephen Balaban, Jul 25 2011 [Comment clarified by N. J. A. Sloane, Aug 01 2024]
Sequence found by reading the line from 0, in the direction 0, 3, ..., and the same line from 0, in the direction 0, 9, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. This is one of the orthogonal axes of the spiral; the other is A032528. - Omar E. Pol, Sep 08 2011
A005449(a(n)) = A000332(3n + 3) = C(3n + 3, 4), a second pentagonal number of triangular matchstick number index number. Additionally, a(n) - 2n is a pentagonal number (A000326). - Raphie Frank, Dec 31 2012
Sum of the numbers from n to 2n. - Wesley Ivan Hurt, Nov 24 2015
Number of orbits of Aut(Z^7) as function of the infinity norm (n+1) of the representative integer lattice point of the orbit, when the cardinality of the orbit is equal to 5376 or 17920 or 20160. - Philippe A.J.G. Chevalier, Dec 28 2015
Also the number of 4-cycles in the (n+4)-triangular honeycomb acute knight graph. - Eric W. Weisstein, Jul 27 2017
Number of terms less than 10^k, k=0,1,2,3,...: 1, 3, 8, 26, 82, 258, 816, 2582, 8165, 25820, 81650, 258199, 816497, 2581989, 8164966, ... - Muniru A Asiru, Jan 24 2018
Numbers of the form 3*m*(2*m + 1) for m = 0, -1, 1, -2, 2, -3, 3, ... - Bruno Berselli, Feb 26 2018
Partial sums of A008585. - Omar E. Pol, Jun 20 2018
Column 1 of A273464. (Number of ways to select a unit lozenge inside an isosceles triangle of side length n; all vertices on a hexagonal lattice.) - R. J. Mathar, Jul 10 2019
Total number of pips in the n-th suit of a double-n domino set. - Ivan N. Ianakiev, Aug 23 2020

Examples

			From _Stephen Balaban_, Jul 25 2011: (Start)
T(n), the triangular numbers = number of nodes,
a(n-1) = number of edges in the T(n) graph:
       o    (T(1) = 1, a(0) = 0)
       o
      / \   (T(2) = 3, a(1) = 3)
     o - o
       o
      / \
     o - o  (T(3) = 6, a(2) = 9)
    / \ / \
   o - o - o
... [Corrected by _N. J. A. Sloane_, Aug 01 2024] (End)
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 543.

Crossrefs

The generalized pentagonal numbers b*n+3*n*(n-1)/2, for b = 1 through 12, form sequences A000326, A005449, A045943, A115067, A140090, A140091, A059845, A140672, A140673, A140674, A140675, A151542.
A diagonal of A010027.
Orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A115067, A008585, A005843, A001477, A000217.
Cf. A027480 (partial sums).
Cf. A002378 (3-cycles in triangular honeycomb acute knight graph), A028896 (5-cycles), A152773 (6-cycles).
This sequence: Sum_{k = n..2*n} k.
Cf. A304993: Sum_{k = n..2*n} k*(k+1)/2.
Cf. A050409: Sum_{k = n..2*n} k^2.
Similar sequences are listed in A316466.

Programs

Formula

a(n) is the sum of n+1 integers starting from n, i.e., 1+2, 2+3+4, 3+4+5+6, 4+5+6+7+8, etc. - Jon Perry, Jan 15 2004
a(n) = A126890(n+1,n-1) for n>1. - Reinhard Zumkeller, Dec 30 2006
a(n) + A145919(3*n+3) = 0. - Matthew Vandermast, Oct 28 2008
a(n) = A000217(2*n) - A000217(n-1); A179213(n) <= a(n). - Reinhard Zumkeller, Jul 05 2010
a(n) = a(n-1)+3*n, n>0. - Vincenzo Librandi, Nov 18 2010
G.f.: 3*x/(1-x)^3. - Bruno Berselli, Jan 21 2011
a(n) = A005448(n+1) - 1. - Omar E. Pol, Oct 03 2011
a(n) = A001477(n)+A000290(n)+A000217(n). - J. M. Bergot, Dec 08 2012
a(n) = 3*a(n-1)-3*a(n-2)+a(n-3) for n>2. - Wesley Ivan Hurt, Nov 24 2015
a(n) = A027480(n)-A027480(n-1). - Peter M. Chema, Jan 18 2017.
2*a(n)+1 = A003215(n). - Miquel Cerda, Jan 22 2018
a(n) = T(2*n) - T(n-1), where T(n) = A000217(n). In general, T(k)*T(n) = Sum_{i=0..k-1} (-1)^i*T((k-i)*(n-i)). - Charlie Marion, Dec 06 2020
E.g.f.: 3*exp(x)*x*(2 + x)/2. - Stefano Spezia, May 19 2021
From Amiram Eldar, Jan 10 2022: (Start)
Sum_{n>=1} 1/a(n) = 2/3.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*(2*log(2)-1)/3. (End)
Product_{n>=1} (1 - 1/a(n)) = -(3/(2*Pi))*cos(sqrt(11/3)*Pi/2). - Amiram Eldar, Feb 21 2023

A000579 Figurate numbers or binomial coefficients C(n,6).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 7, 28, 84, 210, 462, 924, 1716, 3003, 5005, 8008, 12376, 18564, 27132, 38760, 54264, 74613, 100947, 134596, 177100, 230230, 296010, 376740, 475020, 593775, 736281, 906192, 1107568, 1344904, 1623160, 1947792, 2324784, 2760681, 3262623
Offset: 0

Keywords

Comments

Number of triangles (all of whose vertices lie inside the circle) formed when n points in general position on a circle are joined by straight lines - Antreas P. Hatzipolakis (xpolakis(AT)otenet.gr), May 25 2000
Figurate numbers based on 6-dimensional regular simplex. According to Hyun Kwang Kim, it appears that every nonnegative integer can be represented as the sum of g = 13 of these numbers. - Jonathan Vos Post, Nov 28 2004
a(n) = A110555(n+1,6). - Reinhard Zumkeller, Jul 27 2005
a(n) is the number of terms in the expansion of (a_1 + a_2 + a_3 + a_4 + a_5 + a_6 + a_7)^n. - Sergio Falcon, Feb 12 2007
Only prime in this sequence is 7. - Artur Jasinski, Dec 02 2007
6-dimensional triangular numbers, sixth partial sums of binomial transform of [1, 0, 0, 0, ...]. - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009, R. J. Mathar, Jul 07 2009
The number of n-digit numbers the binary expansion of which contains 3 runs of 0's. Generally, the number of n-digit numbers with k runs of 0's is Sum_{i = k..n-k} binomial(i-1, k-1)*binomial(n-i, k) = C(n,2*k) = A034839(n,k) - Vladimir Shevelev, Jul 30 2010
The dimension of the space spanned by a 6-form that couples to M5-brane worldsheets wrapping 6-cycles inside tori (ref. Green,Miller,Vanhove eq. 3.10). - Stephen Crowley, Jan 09 2012
For a set of integers {1,2,...,n}, A253943(n) is the sum of the 2 smallest elements of each subset with 5 elements, which is 3*C(n+1,6) (for n>=5), hence A253943(n) = 3*a(n+1). - Serhat Bulut, Oktay Erkan Temizkan, Mar 13 2015
a(n) = fallfac(n, 6)/6! is also the number of independent components of an antisymmetric tensor of rank 6 and dimension n >= 1. Here fallfac is the falling factorial. - Wolfdieter Lang, Dec 10 2015
Number of orbits of Aut(Z^7) as function of the infinity norm n of the representative integer lattice point of the orbit, when the cardinality of the orbit is equal to 645120. - Philippe A.J.G. Chevalier, Dec 28 2015
Coordination sequence for 6-dimensional cyclotomic lattice Z[zeta_7].

Examples

			a(9) = 84 = (1, 3, 3, 1) dot (1, 6, 15, 20) = (1 + 18 + 45 + 20). - _Gary W. Adamson_, Aug 02 2008
G.f. = x^6 + 7*x^7 + 28*x^8 + 84*x^9 + 210*x^10 + 462*x^11 + 924*x^12 + ...
For A = {1,2,3,4,5,6} subsets with 5 elements are {1,2,3,4,5}, {1,2,3,4,6}, {1,2,3,5,6}, {1,2,4,5,6}, {1,3,4,5,6}, {2,3,4,5,6}. Sum of 2 smallest elements of each subset: a(6) = (1+2) + (1+2) + (1+2) + (1+2) + (1+3) + (2+3) = 21 = 3*C(6+1,6) = 3*A000579(6+1). - _Serhat Bulut_, Oktay Erkan Temizkan, Mar 13 2015
a(7) = 7 from the seven independent components of an antisymmetric tensor A of rank 6 and dimension 7: A(1,2,3,4,5,6), A(1,2,3,4,5,7), A(1,2,3,4,6,7), A(1,2,3,5,6,7) A(1,2,4,5,6,7), A(1,2,3,5,6,7) and A(2,3,4,5,6,7). See a Dec 10 2015 comment. - _Wolfdieter Lang_, Dec 10 2015
		

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. 828.
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 196.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 7.
  • J. C. P. Miller, editor, Table of Binomial Coefficients. Royal Society Mathematical Tables, Vol. 3, Cambridge Univ. Press, 1954.
  • 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).
  • Charles W. Trigg: Mathematical Quickies. New York: Dover Publications, Inc., 1985, p. 11, #32

Crossrefs

Cf. A053135, A053128, A000580 (partial sums), A000581, A000582, A000217, A000292, A000332, A000389 (first differences), A104712 (fifth column, k=6).

Programs

  • Magma
    [Binomial(n,6) : n in [0..50]]; // Wesley Ivan Hurt, Jul 13 2014
    
  • Maple
    A000579 := n->binomial(n,6);
    ZL := [S, {S=Prod(B,B,B,B,B,B,B), B=Set(Z, 1 <= card)}, unlabeled]: seq(combstruct[count](ZL, size=n), n=7..40); # Zerinvary Lajos, Mar 13 2007
    A000579:=-1/(z-1)**7; # Simon Plouffe in his 1992 dissertation, referring to offset 0.
    seq(binomial(n,6),n=0..33); # Zerinvary Lajos, Jun 16 2008
    G(x):=x^6*exp(x): f[0]:=G(x): for n from 1 to 39 do f[n]:=diff(f[n-1],x) od: x:=0: seq(f[n]/6!,n=6..39); # Zerinvary Lajos, Apr 05 2009
  • Mathematica
    Table[Binomial[n, 6], {n, 6, 50}] (* Stefan Steinerberger, Apr 02 2006 *)
    Table[n(n - 1)(n - 2)(n - 3)(n - 4)(n - 5)/720, {n, 0, 100}] (* Artur Jasinski, Dec 02 2007 *)
    LinearRecurrence[{7,-21,35,-35,21,-7,1},{0,0,0,0,0,0,1},50] (* Harvey P. Dale, Dec 30 2012 *)
    CoefficientList[ Series[ -7x^6/(x-1)^7,{x, 0, 35}], x]/7 (* Robert G. Wilson v, Jan 29 2015 *)
  • PARI
    a(n)=binomial(n,6) \\ Charles R Greathouse IV, Nov 20 2012
    
  • Python
    A000579_list, m = [], [1, -5, 10, -10, 5, -1, 0]
    for _ in range(10**2):
        A000579_list.append(m[-1])
        for i in range(6):
            m[i+1] += m[i] # Chai Wah Wu, Jan 24 2016

Formula

G.f.: x^6/(1-x)^7.
E.g.f.: exp(x)*x^6/720.
a(n) = (n^6 - 15*n^5 + 85*n^4 - 225*n^3 + 274*n^2 - 120*n)/720.
Conjecture: a(n+3) = Sum_{0 <= k, L, m <= n; k + L + m <= n} k*L*m. - Ralf Stephan, May 06 2005
Convolution of the nonnegative numbers (A001477) with the hexagonal numbers (A000389). Also convolution of the triangular numbers (A000217) with the tetrahedral numbers (A000292). - Sergio Falcon, Feb 12 2007
a(n) = n*(n - 1)*(n - 2)*(n - 3)*(n - 4)*(n - 5)/720. - Artur Jasinski, Dec 02 2007, R. J. Mathar, Jul 07 2009
Equals binomial transform of [1, 6, 15, 20, 15, 6, 1, 0, 0, 0, ...]. - Gary W. Adamson, Aug 02 2008
a(0) = 0, a(1) = 0, a(2) = 0, a(3) = 0, a(4) = 0, a(5) = 0, a(6) = 1, 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). - Harvey P. Dale, Dec 30 2012
Sum_{n >= 0} a(n)/n! = e/720. Sum_{n >= 5} a(n)/(n-5)! = 4051*e/720. See A067653 regarding the second ratio. - Richard R. Forberg, Dec 26 2013
Sum_{n >= 6} 1/a(n) = 6/5. - Hermann Stamm-Wilbrandt, Jul 13 2014
Sum_{n >= 6} (-1)^(n + 1)/a(n) = 192*log(2) - 661/5 = 0.8842586675... Also see A242023. - Richard R. Forberg, Aug 11 2014
a(n) = a(5-n) for all n in Z. - Michael Somos, Oct 07 2014
0 = a(n)*(+a(n+1) +5*a(n+2)) + a(n+1)*(-7*a(n+1) +a(n+2)) for all n in Z. - Michael Somos, Oct 07 2014
a(n) = 3*C(n+1,6) = 3*A000579(n+1). - Serhat Bulut, Oktay Erkan Temizkan, Mar 13 2015
a(n) = A000292(n-5)*A000292(n-2)/20. - R. J. Mathar, Nov 29 2015

Extensions

Some formulas that referred to other offsets corrected by R. J. Mathar, Jul 07 2009
I changed the offset to 0. This will require some further adjustments to the formulas. - N. J. A. Sloane, Aug 01 2010
Shevelev comment inserted and further adaptations to offset by R. J. Mathar, Aug 03 2010

A002417 4-dimensional figurate numbers: a(n) = n*binomial(n+2, 3).

Original entry on oeis.org

1, 8, 30, 80, 175, 336, 588, 960, 1485, 2200, 3146, 4368, 5915, 7840, 10200, 13056, 16473, 20520, 25270, 30800, 37191, 44528, 52900, 62400, 73125, 85176, 98658, 113680, 130355, 148800, 169136, 191488, 215985, 242760, 271950, 303696, 338143, 375440, 415740
Offset: 1

Keywords

Comments

a(n) is 1/6 the number of colorings of a 2 X 2 hexagonal array with n+2 colors. - R. H. Hardin, Feb 23 2002
a(n) is the sum of all numbers that cannot be written as t*(n+1) + u*(n+2) for nonnegative integers t,u (see Schuh). - Floor van Lamoen, Oct 09 2002
a(n) is the total number of rectangles (including squares) contained in a stepped pyramid of n rows (or of base 2n-1) of squares. A stepped pyramid of squares of base 2*6 - 1 = 11, for instance, has the following vertices:
..........X.X
........X.X.X.X
......X.X.X.X.X.X
....X.X.X.X.X.X.X.X
..X.X.X.X.X.X.X.X.X.X
X.X.X.X.X.X.X.X.X.X.X.X
X.X.X.X.X.X.X.X.X.X.X.X - Lekraj Beedassy, Sep 02 2003
Partial sums of A002412. - Jonathan Vos Post, Mar 16 2006
a(n) equals -1 times the coefficient of x^3 of the characteristic polynomial of the (n + 2) X (n + 2) matrix with 2's along the main diagonal and 1's everywhere else (see Mathematica code below). - John M. Campbell, May 28 2011
a(n) is the n-th antidiagonal sum of the convolution array A213750. - Clark Kimberling, Jun 20 2012
Convolution of A000027 with A000384 (excluding 0). - Bruno Berselli, Dec 06 2012
The sequence is the binomial transform of (1, 7, 15, 13, 4, 0, 0, 0, ...). - Gary W. Adamson, Jul 31 2015
Also the number of 3-cycles in the (n+2)-triangular graph. - Eric W. Weisstein, Aug 14 2017

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 195.
  • K. -W. Lau, Solution to Problem 2495, Journal of Recreational Mathematics 2002-3 31(1) 79-80.
  • Fred. Schuh, Vragen betreffende een onbepaalde vergelijking, Nieuw Tijdschrift voor Wiskunde, 52 (1964-1965) 193-198.
  • 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).

Crossrefs

Bisection of A002624.
a(n) = A093561(n+3, 4).
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers with the k-gonal numbers.
Cf. A151974 (number of 4-cycles in the triangular graph), A290939 (5-cycles), A290940 (6-cycles).

Programs

  • GAP
    List([1..40], n-> n^2*(n+1)*(n+2)/6 ); # G. C. Greubel, Jul 03 2019
  • Magma
    /* A000027 convolved with A000384 (excluding 0): */ A000384:=func; [&+[(n-i+1)*A000384(i): i in [1..n]]: n in [1..40]]; // Bruno Berselli, Dec 06 2012
    
  • Magma
    [n*Binomial(n+2,3):n in [1..40]]; // Vincenzo Librandi, Aug 02 2015
    
  • Maple
    seq(n^2*(n+1)*(n+2)/6, n=1..50);
  • Mathematica
    Table[n Binomial[n+2, 3], {n, 40}]
    Table[-Coefficient[CharacteristicPolynomial[Array[KroneckerDelta[#1, #2] + 1 &, {n+2, n+2}], x], x^3], {n, 40}] (* John M. Campbell, May 28 2011 *)
    Nest[Accumulate, Range[1, 170, 4], 3] (* Vladimir Joseph Stephan Orlovsky, Jan 21 2012 *)
    LinearRecurrence[{5, -10, 10, -5, 1}, {1, 8, 30, 80, 175}, 40] (* Harvey P. Dale, Jan 11 2014 *)
    Table[n Pochhammer[n, 3]/6, {n, 40}] (* or *) CoefficientList[Series[ (1+3x)/(1-x)^5, {x,0,40}], x] (* Eric W. Weisstein, Aug 14 2017 *)
  • PARI
    a(n)=n^2*(n+1)*(n+2)/6 \\ Charles R Greathouse IV, Jun 10 2011
    
  • Sage
    [n*binomial(n+2,3) for n in (1..40)] # G. C. Greubel, Jul 03 2019
    

Formula

a(n) = n^2*(n+1)*(n+2)/6.
G.f.: x*(1+3*x)/(1-x)^5. - Simon Plouffe in his 1992 dissertation
a(n) = C(n+2, 2)*n^2/3. - Paul Barry, Jun 26 2003
a(n) = C(n+3, n)*C(n+1, 1). - Zerinvary Lajos, Apr 27 2005
a(n) = (binomial(n+3,n-1) - binomial(n+2,n-2))*(binomial(n+1,n-1) - binomial(n,n-2)). - Zerinvary Lajos, May 12 2006
a(n) = 5*a(n-1) -10*a(n-2) +10*a(n-3) -5*a(n-4) +a(n-5), n>5. - Wesley Ivan Hurt, Aug 01 2015
G.f.: x*hypergeometric2F1(2,4;1;x). - R. J. Mathar, Aug 09 2015
a(n) = A080852(4,n-1). - R. J. Mathar, Jul 28 2016
Sum_{n>=1} 1/a(n) = Pi^2/2 - 15/4. - Jaume Oliver Lafont, Jul 13 2017
E.g.f.: x*(6 + 18*x + 9*x^2 + x^3)*exp(x)/3!. - G. C. Greubel, Jul 03 2019
Sum_{n>=1} (-1)^(n+1)/a(n) = (Pi^2 + 27 - 48*log(2))/4. - Amiram Eldar, Jun 28 2020
a(n) = A000332(n+3) + 3*A000332(n+2). - Mircea Dan Rus, Jul 29 2020

Extensions

Edited and extended by Floor van Lamoen, Oct 09 2002

A007531 a(n) = n*(n-1)*(n-2) (or n!/(n-3)!).

Original entry on oeis.org

0, 0, 0, 6, 24, 60, 120, 210, 336, 504, 720, 990, 1320, 1716, 2184, 2730, 3360, 4080, 4896, 5814, 6840, 7980, 9240, 10626, 12144, 13800, 15600, 17550, 19656, 21924, 24360, 26970, 29760, 32736, 35904, 39270, 42840, 46620, 50616, 54834, 59280, 63960, 68880
Offset: 0

Keywords

Comments

Ed Pegg Jr conjectures that n^3 - n = k! has a solution if and only if n is 2, 3, 5 or 9 (when k is 3, 4, 5 and 6).
Three-dimensional promic (or oblong) numbers, cf. A002378. - Alexandre Wajnberg, Dec 29 2005
Doubled first differences of tritriangular numbers A050534(n) = (1/8)n(n + 1)(n - 1)(n - 2). a(n) = 2*(A050534(n+1) - A050534(n)). - Alexander Adamchuk, Apr 11 2006
If Y is a 4-subset of an n-set X then, for n >= 6, a(n-4) is the number of (n-5)-subsets of X having exactly two elements in common with Y. - Milan Janjic, Dec 28 2007
Convolution of A005843 with A008585. - Reinhard Zumkeller, Mar 07 2009
a(n) = A000578(n) - A000567(n). - Reinhard Zumkeller, Sep 18 2009
For n > 3: a(n) = A173333(n, n-3). - Reinhard Zumkeller, Feb 19 2010
Let H be the n X n Hilbert matrix H(i, j) = 1/(i+j-1) for 1 <= i, j <= n. Let B be the inverse matrix of H. The sum of the elements in row 2 of B equals (-1)^n a(n+1). - T. D. Noe, May 01 2011
a(n) equals 2^(n-1) times the coefficient of log(3) in 2F1(n-2, n-2, n, -2). - John M. Campbell, Jul 16 2011
For n > 2 a(n) = 1/(Integral_{x = 0..Pi/2} (sin(x))^5*(cos(x))^(2*n-5)). - Francesco Daddi, Aug 02 2011
a(n) is the number of functions f:[3] -> [n] that are injective since there are n choices for f(1), (n-1) choices for f(2), and (n-2) choices for f(3). Also, a(n+1) is the number of functions f:[3] -> [n] that are width-2 restricted (that is, the pre-image under f of any element in [n] is of size 2 or less). See "Width-restricted finite functions" link below. - Dennis P. Walsh, Mar 01 2012
This sequence is produced by three consecutive triangular numbers t(n-1), t(n-2) and t(n-3) in the expression 2*t(n-1)*(t(n-2)-t(n-3)) for n = 0, 1, 2, ... - J. M. Bergot, May 14 2012
For n > 2: A020639(a(n)) = 2; A006530(a(n)) = A093074(n-1). - Reinhard Zumkeller, Jul 04 2012
Number of contact points between equal spheres arranged in a tetrahedron with n - 1 spheres in each edge. - Ignacio Larrosa Cañestro, Jan 07 2013
Also for n >= 3, area of Pythagorean triangle in which one side differs from hypotenuse by two units. Consider any Pythagorean triple (2n, n^2-1, n^2+1) where n > 1. The area of such a Pythagorean triangle is n(n^2-1). For n = 2, 3, 4,.. the areas are 6, 24, 60, .... which are the given terms of the series. - Jayanta Basu, Apr 11 2013
Cf. A130534 for relations to colored forests, disposition of flags on flagpoles, and colorings of the vertices (chromatic polynomial) of the complete graph K_3. - Tom Copeland, Apr 05 2014
Starting with 6, 24, 60, 120, ..., a(n) is the number of permutations of length n>=3 avoiding the partially ordered pattern (POP) {1>2} of length 5. That is, the number of length n permutations having no subsequences of length 5 in which the first element is larger than the second element. - Sergey Kitaev, Dec 11 2020
For integer m and positive integer r >= 2, the polynomial a(n) + a(n + m) + a(n + 2*m) + ... + a(n + r*m) in n has its zeros on the vertical line Re(n) = (2 - r*m)/2 in the complex plane. - Peter Bala, Jun 02 2024

References

  • R. K. Guy, Unsolved Problems in Theory of Numbers, Section D25.
  • L. B. W. Jolley, "Summation of Series", Dover Publications, 1961, p. 40.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

binomial(n, k): A161680 (k = 2), A000332 (k = 4), A000389 (k = 5), A000579 (k = 6), A000580 (k = 7), A000581 (k = 8), A000582 (k = 9).
Cf. A028896.

Programs

  • Haskell
    a007531 n = product [n-2..n]  -- Reinhard Zumkeller, Jul 04 2012
    
  • Magma
    [n*(n-1)*(n-2): n in [0..40]]; // Vincenzo Librandi, May 02 2011
    
  • Maple
    [seq(6*binomial(n,3),n=0..41)]; # Zerinvary Lajos, Nov 24 2006
  • Mathematica
    Table[n^3 - 3n^2 + 2n, {n, 0, 42}]
    Table[FactorialPower[n, 3], {n, 0, 42}] (* Arkadiusz Wesolowski, Oct 29 2012 *)
  • PARI
    a(n)=n*(n-1)*(n-2)
    
  • Sage
    [n*(n-1)*(n-2) for n in range(40)] # G. C. Greubel, Feb 11 2019

Formula

a(n) = 6*A000292(n-2).
a(n) = Sum_{i=1..n} polygorial(3,i) where polygorial(3,i) = A028896(i-1). - Daniel Dockery (peritus(AT)gmail.com), Jun 16 2003
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 6, n > 2. - Zak Seidov, Feb 09 2006
G.f.: 6*x^2/(1-x)^4.
a(-n) = -a(n+2).
1/6 + 3/24 + 5/60 + ... = Sum_{k>=1} (2*k-1)/(k*(k+1)*(k+2)) = 3/4. [Jolley Eq. 213]
a(n+1) = n^3 - n. - Mohammad K. Azarian, Jul 26 2007
E.g.f.: x^3*exp(x). - Geoffrey Critzer, Feb 08 2009
If the first 0 is eliminated, a(n) = floor(n^5/(n^2+1)). - Gary Detlefs, Feb 11 2010
1/6 + 1/24 + 1/60 + ... = Sum_{n>=1} 1/(n*(n+1)*(n+2)) = 1/4. - Mohammad K. Azarian, Dec 29 2010
a(0) = 0, a(n) = a(n-1) + 3*(n-1)*(n-2). - Jean-François Alcover, Jan 08 2013
(a(n+1) - a(n))/6 = A000217(n-2) for n > 0. - J. M. Bergot, Jul 30 2013
Partial sums of A028896. - R. J. Mathar, Aug 28 2014
1/6 + 1/24 + 1/60 + ... + 1/(n*(n+1)*(n+2)) = n*(n+3)/(4*(n+1)*(n+2)). - Christina Steffan, Jul 20 2015
a(n+2)^2 = A005563(n)^3 + A005563(n)^2. - Bruno Berselli, May 03 2018
a(n)*a(n+1) + A000096(n-3)^2 = m^2 (a perfect square), m = ((a(n)+a(n+1))/2)-n. - Ezhilarasu Velayutham, May 21 2019
Sum_{n>=3} (-1)^(n+1)/a(n) = 2*log(2) - 5/4. - Amiram Eldar, Jul 02 2020
For n >= 3, (a(n) + (a(n) + (a(n) + ...)^(1/3))^(1/3))^(1/3) = n - 1. - Paolo Xausa, Apr 09 2022
Previous Showing 31-40 of 396 results. Next