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-4 of 4 results.

A000326 Pentagonal numbers: a(n) = n*(3*n-1)/2.

Original entry on oeis.org

0, 1, 5, 12, 22, 35, 51, 70, 92, 117, 145, 176, 210, 247, 287, 330, 376, 425, 477, 532, 590, 651, 715, 782, 852, 925, 1001, 1080, 1162, 1247, 1335, 1426, 1520, 1617, 1717, 1820, 1926, 2035, 2147, 2262, 2380, 2501, 2625, 2752, 2882, 3015, 3151
Offset: 0

Views

Author

Keywords

Comments

The average of the first n (n > 0) pentagonal numbers is the n-th triangular number. - Mario Catalani (mario.catalani(AT)unito.it), Apr 10 2003
a(n) is the sum of n integers starting from n, i.e., 1, 2 + 3, 3 + 4 + 5, 4 + 5 + 6 + 7, etc. - Jon Perry, Jan 15 2004
Partial sums of 1, 4, 7, 10, 13, 16, ... (1 mod 3), a(2k) = k(6k-1), a(2k-1) = (2k-1)(3k-2). - Jon Perry, Sep 10 2004
Starting with offset 1 = binomial transform of [1, 4, 3, 0, 0, 0, ...]. Also, A004736 * [1, 3, 3, 3, ...]. - Gary W. Adamson, Oct 25 2007
If Y is a 3-subset of an n-set X then, for n >= 4, a(n-3) is the number of 4-subsets of X having at least two elements in common with Y. - Milan Janjic, Nov 23 2007
Solutions to the duplication formula 2*a(n) = a(k) are given by the index pairs (n, k) = (5,7), (5577, 7887), (6435661, 9101399), etc. The indices are integer solutions to the pair of equations 2(6n-1)^2 = 1 + y^2, k = (1+y)/6, so these n can be generated from the subset of numbers [1+A001653(i)]/6, any i, where these are integers, confined to the cases where the associated k=[1+A002315(i)]/6 are also integers. - R. J. Mathar, Feb 01 2008
a(n) is a binomial coefficient C(n,4) (A000332) if and only if n is a generalized pentagonal number (A001318). Also see A145920. - Matthew Vandermast, Oct 28 2008
Even octagonal numbers divided by 8. - Omar E. Pol, Aug 18 2011
Sequence found by reading the line from 0, in the direction 0, 5, ... and the line from 1, in the direction 1, 12, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. - Omar E. Pol, Sep 08 2011
The hyper-Wiener index of the star-tree with n edges (see A196060, example). - Emeric Deutsch, Sep 30 2011
More generally the n-th k-gonal number is equal to n + (k-2)*A000217(n-1), n >= 1, k >= 3. In this case k = 5. - Omar E. Pol, Apr 06 2013
Note that both Euler's pentagonal theorem for the partition numbers and Euler's pentagonal theorem for the sum of divisors refer more exactly to the generalized pentagonal numbers, not this sequence. For more information see A001318, A175003, A238442. - Omar E. Pol, Mar 01 2014
The Fuss-Catalan numbers are Cat(d,k)= [1/(k*(d-1)+1)]*binomial(k*d,k) and enumerate the number of (d+1)-gon partitions of a (k*(d-1)+2)-gon (cf. Schuetz and Whieldon link). a(n)= Cat(n,3), so enumerates the number of (n+1)-gon partitions of a (3*(n-1)+2)-gon. Analogous sequences are A100157 (k=4) and A234043 (k=5). - Tom Copeland, Oct 05 2014
Binomial transform of (0, 1, 3, 0, 0, 0, ...) (A169585 with offset 1) and second partial sum of (0, 1, 3, 3, 3, ...). - Gary W. Adamson, Oct 05 2015
For n > 0, a(n) is the number of compositions of n+8 into n parts avoiding parts 2 and 3. - Milan Janjic, Jan 07 2016
a(n) is also the number of edges in the Mycielskian of the complete graph K[n]. Indeed, K[n] has n vertices and n(n-1)/2 edges. Then its Mycielskian has n + 3n(n-1)/2 = n(3n-1)/2. See p. 205 of the West reference. - Emeric Deutsch, Nov 04 2016
Sum of the numbers from n to 2n-1. - Wesley Ivan Hurt, Dec 03 2016
Also the number of maximal cliques in the n-Andrásfai graph. - Eric W. Weisstein, Dec 01 2017
Coefficients in the hypergeometric series identity 1 - 5*(x - 1)/(2*x + 1) + 12*(x - 1)*(x - 2)/((2*x + 1)*(2*x + 2)) - 22*(x - 1)*(x - 2)*(x - 3)/((2*x + 1)*(2*x + 2)*(2*x + 3)) + ... = 0, valid for Re(x) > 1. Cf. A002412 and A002418. Column 2 of A103450. - Peter Bala, Mar 14 2019
A generalization of the Comment dated Apr 10 2003 follows. (k-3)*A000292(n-2) plus the average of the first n (2k-1)-gonal numbers is the n-th k-gonal number. - Charlie Marion, Nov 01 2020
a(n+1) is the number of Dyck paths of size (3,3n+1); i.e., the number of NE lattice paths from (0,0) to (3,3n+1) which stay above the line connecting these points. - Harry Richman, Jul 13 2021
a(n) is the largest sum of n positive integers x_1, ..., x_n such that x_i | x_(i+1)+1 for each 1 <= i <= n, where x_(n+1) = x_1. - Yifan Xie, Feb 21 2025

Examples

			Illustration of initial terms:
.
.                                       o
.                                     o o
.                          o        o o o
.                        o o      o o o o
.                o     o o o    o o o o o
.              o o   o o o o    o o o o o
.        o   o o o   o o o o    o o o o o
.      o o   o o o   o o o o    o o o o o
.  o   o o   o o o   o o o o    o o o o o
.
.  1    5     12       22           35
- _Philippe Deléham_, Mar 30 2013
		

References

  • Tom M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, pages 2 and 311.
  • Raymond Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; p. 129.
  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 38, 40.
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6.
  • 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. 1.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.6 Figurate Numbers, p. 291.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 284.
  • Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 64.
  • Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 52-53, 129-130, 132.
  • 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).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 7-10.
  • André Weil, Number theory: an approach through history; from Hammurapi to Legendre, Birkhäuser, Boston, 1984; see p. 186.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, 1987, pp. 98-100.
  • Douglas B. West, Introduction to Graph Theory, 2nd ed., Prentice-Hall, NJ, 2001.

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.
Cf. A001318 (generalized pentagonal numbers), A049452, A033570, A010815, A034856, A051340, A004736, A033568, A049453, A002411 (partial sums), A033579.
See A220083 for a list of numbers of the form n*P(s,n)-(n-1)*P(s,n-1), where P(s,n) is the n-th polygonal number with s sides.
Cf. A240137: sum of n consecutive cubes starting from n^3.
Cf. similar sequences listed in A022288.
Partial sums of A016777.

Programs

  • GAP
    List([0..50],n->n*(3*n-1)/2); # Muniru A Asiru, Mar 18 2019
    
  • Haskell
    a000326 n = n * (3 * n - 1) `div` 2  -- Reinhard Zumkeller, Jul 07 2012
    
  • Magma
    [n*(3*n-1)/2 : n in [0..100]]; // Wesley Ivan Hurt, Oct 15 2015
    
  • Maple
    A000326 := n->n*(3*n-1)/2: seq(A000326(n), n=0..100);
    A000326:=-(1+2*z)/(z-1)**3; # Simon Plouffe in his 1992 dissertation
    a[0]:=0:a[1]:=1:for n from 2 to 50 do a[n]:=2*a[n-1]-a[n-2]+3 od: seq(a[n], n=0..50); # Miklos Kristof, Zerinvary Lajos, Feb 18 2008
  • Mathematica
    Table[n (3 n - 1)/2, {n, 0, 60}] (* Stefan Steinerberger, Apr 01 2006 *)
    Array[# (3 # - 1)/2 &, 47, 0] (* Zerinvary Lajos, Jul 10 2009 *)
    LinearRecurrence[{3, -3, 1}, {0, 1, 5}, 61] (* Harvey P. Dale, Dec 27 2011 *)
    pentQ[n_] := IntegerQ[(1 + Sqrt[24 n + 1])/6]; pentQ[0] = True; Select[Range[0, 3200], pentQ@# &] (* Robert G. Wilson v, Mar 31 2014 *)
    Join[{0}, Accumulate[Range[1, 312, 3]]] (* Harvey P. Dale, Mar 26 2016 *)
    (* For Mathematica 10.4+ *) Table[PolygonalNumber[RegularPolygon[5], n], {n, 0, 46}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
    CoefficientList[Series[x (-1 - 2 x)/(-1 + x)^3, {x, 0, 20}], x] (* Eric W. Weisstein, Dec 01 2017 *)
    PolygonalNumber[5, Range[0, 20]] (* Eric W. Weisstein, Dec 01 2017 *)
  • PARI
    a(n)=n*(3*n-1)/2
    
  • PARI
    vector(100, n, n--; binomial(3*n, 2)/3) \\ Altug Alkan, Oct 06 2015
    
  • PARI
    is_a000326(n) = my(s); n==0 || (issquare (24*n+1, &s) && s%6==5); \\ Hugo Pfoertner, Aug 03 2023
    
  • Python
    # Intended to compute the initial segment of the sequence, not isolated terms.
    def aList():
         x, y = 1, 1
         yield 0
         while True:
             yield x
             x, y = x + y + 3, y + 3
    A000326 = aList()
    print([next(A000326) for i in range(47)]) # Peter Luschny, Aug 04 2019

Formula

Product_{m > 0} (1 - q^m) = Sum_{k} (-1)^k*x^a(k). - Paul Barry, Jul 20 2003
G.f.: x*(1+2*x)/(1-x)^3.
E.g.f.: exp(x)*(x+3*x^2/2).
a(n) = n*(3*n-1)/2.
a(-n) = A005449(n).
a(n) = binomial(3*n, 2)/3. - Paul Barry, Jul 20 2003
a(n) = A000290(n) + A000217(n-1). - Lekraj Beedassy, Jun 07 2004
a(0) = 0, a(1) = 1; for n >= 2, a(n) = 2*a(n-1) - a(n-2) + 3. - Miklos Kristof, Mar 09 2005
a(n) = Sum_{k=1..n} (2*n - k). - Paul Barry, Aug 19 2005
a(n) = 3*A000217(n) - 2*n. - Lekraj Beedassy, Sep 26 2006
a(n) = A126890(n, n-1) for n > 0. - Reinhard Zumkeller, Dec 30 2006
a(n) = A049452(n) - A022266(n) = A033991(n) - A005476(n). - Zerinvary Lajos, Jun 12 2007
Equals A034856(n) + (n - 1)^2. Also equals A051340 * [1,2,3,...]. - Gary W. Adamson, Jul 27 2007
a(n) = binomial(n+1, 2) + 2*binomial(n, 2).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0) = 0, a(1) = 1, a(2) = 5. - Jaume Oliver Lafont, Dec 02 2008
a(n) = a(n-1) + 3*n-2 with n > 0, a(0)=0. - Vincenzo Librandi, Nov 20 2010
a(n) = A000217(n) + 2*A000217(n-1). - Vincenzo Librandi, Nov 20 2010
a(n) = A014642(n)/8. - Omar E. Pol, Aug 18 2011
a(n) = A142150(n) + A191967(n). - Reinhard Zumkeller, Jul 07 2012
a(n) = (A000290(n) + A000384(n))/2 = (A000217(n) + A000566(n))/2 = A049450(n)/2. - Omar E. Pol, Jan 11 2013
a(n) = n*A000217(n) - (n-1)*A000217(n-1). - Bruno Berselli, Jan 18 2013
a(n) = A005449(n) - n. - Philippe Deléham, Mar 30 2013
From Oskar Wieland, Apr 10 2013: (Start)
a(n) = a(n+1) - A016777(n),
a(n) = a(n+2) - A016969(n),
a(n) = a(n+3) - A016777(n)*3 = a(n+3) - A017197(n),
a(n) = a(n+4) - A016969(n)*2 = a(n+4) - A017641(n),
a(n) = a(n+5) - A016777(n)*5,
a(n) = a(n+6) - A016969(n)*3,
a(n) = a(n+7) - A016777(n)*7,
a(n) = a(n+8) - A016969(n)*4,
a(n) = a(n+9) - A016777(n)*9. (End)
a(n) = A000217(2n-1) - A000217(n-1), for n > 0. - Ivan N. Ianakiev, Apr 17 2013
a(n) = A002411(n) - A002411(n-1). - J. M. Bergot, Jun 12 2013
Sum_{n>=1} a(n)/n! = 2.5*exp(1). - Richard R. Forberg, Jul 15 2013
a(n) = floor(n/(exp(2/(3*n)) - 1)), for n > 0. - Richard R. Forberg, Jul 27 2013
From Vladimir Shevelev, Jan 24 2014: (Start)
a(3*a(n) + 4*n + 1) = a(3*a(n) + 4*n) + a(3*n+1).
A generalization. Let {G_k(n)}_(n >= 0) be sequence of k-gonal numbers (k >= 3). Then the following identity holds: G_k((k-2)*G_k(n) + c(k-3)*n + 1) = G_k((k-2)*G_k(n) + c(k-3)*n) + G_k((k-2)*n + 1), where c = A000124. (End)
A242357(a(n)) = 1 for n > 0. - Reinhard Zumkeller, May 11 2014
Sum_{n>=1} 1/a(n)= (1/3)*(9*log(3) - sqrt(3)*Pi). - Enrique Pérez Herrero, Dec 02 2014. See the decimal expansion A244641.
a(n) = (A000292(6*n+k-1)-A000292(k))/(6*n-1)-A000217(3*n+k), for any k >= 0. - Manfred Arens, Apr 26 2015 [minor edits from Wolfdieter Lang, May 10 2015]
a(n) = A258708(3*n-1,1) for n > 0. - Reinhard Zumkeller, Jun 23 2015
a(n) = A007584(n) - A245301(n-1), for n > 0. - Manfred Arens, Jan 31 2016
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*(sqrt(3)*Pi - 6*log(2))/3 = 0.85501000622865446... - Ilya Gutkovskiy, Jul 28 2016
a(m+n) = a(m) + a(n) + 3*m*n. - Etienne Dupuis, Feb 16 2017
In general, let P(k,n) be the n-th k-gonal number. Then P(k,m+n) = P(k,m) + (k-2)mn + P(k,n). - Charlie Marion, Apr 16 2017
a(n) = A023855(2*n-1) - A023855(2*n-2). - Luc Rousseau, Feb 24 2018
a(n) = binomial(n,2) + n^2. - Pedro Caceres, Jul 28 2019
Product_{n>=2} (1 - 1/a(n)) = 3/5. - Amiram Eldar, Jan 21 2021
(n+1)*(a(n^2) + a(n^2+1) + ... + a(n^2+n)) = n*(a(n^2+n+1) + ... + a(n^2+2n)). - Charlie Marion, Apr 28 2024
a(n) = Sum_{k = 0..3*n} (-1)^(n+k+1) * binomial(k, 2)*binomial(3*n+k-1, 2*k). - Peter Bala, Nov 04 2024

Extensions

Incorrect example removed by Joerg Arndt, Mar 11 2010

A000332 Binomial coefficient binomial(n,4) = n*(n-1)*(n-2)*(n-3)/24.

Original entry on oeis.org

0, 0, 0, 0, 1, 5, 15, 35, 70, 126, 210, 330, 495, 715, 1001, 1365, 1820, 2380, 3060, 3876, 4845, 5985, 7315, 8855, 10626, 12650, 14950, 17550, 20475, 23751, 27405, 31465, 35960, 40920, 46376, 52360, 58905, 66045, 73815, 82251, 91390, 101270, 111930, 123410
Offset: 0

Views

Author

Keywords

Comments

Number of intersection points of diagonals of convex n-gon where no more than two diagonals intersect at any point in the interior.
Also the number of equilateral triangles with vertices in an equilateral triangular array of points with n rows (offset 1), with any orientation. - Ignacio Larrosa Cañestro, Apr 09 2002. [See Les Reid link for proof. - N. J. A. Sloane, Apr 02 2016] [See Peter Kagey link for alternate proof. - Sameer Gauria, Jul 29 2025]
Start from cubane and attach amino acids according to the reaction scheme that describes the reaction between the active sites. See the hyperlink on chemistry. - Robert G. Wilson v, Aug 02 2002
For n>0, a(n) = (-1/8)*(coefficient of x in Zagier's polynomial P_(2n,n)). (Zagier's polynomials are used by PARI/GP for acceleration of alternating or positive series.)
Figurate numbers based on the 4-dimensional regular convex polytope called the regular 4-simplex, pentachoron, 5-cell, pentatope or 4-hypertetrahedron with Schlaefli symbol {3,3,3}. a(n)=((n*(n-1)*(n-2)*(n-3))/4!). - Michael J. Welch (mjw1(AT)ntlworld.com), Apr 01 2004, R. J. Mathar, Jul 07 2009
Maximal number of crossings that can be created by connecting n vertices with straight lines. - Cameron Redsell-Montgomerie (credsell(AT)uoguelph.ca), Jan 30 2007
If X is an n-set and Y a fixed (n-1)-subset of X then a(n) is equal to the number of 4-subsets of X intersecting Y. - Milan Janjic, Aug 15 2007
Product of four consecutive numbers divided by 24. - Artur Jasinski, Dec 02 2007
The only prime in this sequence is 5. - Artur Jasinski, Dec 02 2007
For strings consisting entirely of 0's and 1's, the number of distinct arrangements of four 1's such that 1's are not adjacent. The shortest possible string is 7 characters, of which there is only one solution: 1010101, corresponding to a(5). An eight-character string has 5 solutions, nine has 15, ten has 35 and so on, congruent to A000332. - Gil Broussard, Mar 19 2008
For a(n)>0, a(n) is pentagonal if and only if 3 does not divide n. All terms belong to the generalized pentagonal sequence (A001318). Cf. A000326, A145919, A145920. - Matthew Vandermast, Oct 28 2008
Nonzero terms = row sums of triangle A158824. - Gary W. Adamson, Mar 28 2009
Except for the 4 initial 0's, is equivalent to the partial sums of the tetrahedral numbers A000292. - Jeremy Cahill (jcahill(AT)inbox.com), Apr 15 2009
If the first 3 zeros are disregarded, that is, if one looks at binomial(n+3, 4) with n>=0, then it becomes a 'Matryoshka doll' sequence with alpha=0: seq(add(add(add(i,i=alpha..k),k=alpha..n),n=alpha..m),m=alpha..50). - Peter Luschny, Jul 14 2009
For n>=1, a(n) is the number of n-digit numbers the binary expansion of which contains two runs of 0's. - Vladimir Shevelev, Jul 30 2010
For n>0, a(n) is the number of crossing set partitions of {1,2,..,n} into n-2 blocks. - Peter Luschny, Apr 29 2011
The Kn3, Ca3 and Gi3 triangle sums of A139600 are related to the sequence given above, e.g., Gi3(n) = 2*A000332(n+3) - A000332(n+2) + 7*A000332(n+1). For the definitions of these triangle sums, see A180662. - Johannes W. Meijer, Apr 29 2011
For n > 3, a(n) is the hyper-Wiener index of the path graph on n-2 vertices. - Emeric Deutsch, Feb 15 2012
Except for the four initial zeros, number of all possible tetrahedra of any size, having the same orientation as the original regular tetrahedron, formed when intersecting the latter by planes parallel to its sides and dividing its edges into n equal parts. - V.J. Pohjola, Aug 31 2012
a(n+3) is the number of different ways to color the faces (or the vertices) of a regular tetrahedron with n colors if we count mirror images as the same.
a(n) = fallfac(n,4)/4! is also the number of independent components of an antisymmetric tensor of rank 4 and dimension n >= 1. Here fallfac is the falling factorial. - Wolfdieter Lang, Dec 10 2015
Does not satisfy Benford's law [Ross, 2012] - N. J. A. Sloane, Feb 12 2017
Number of chiral pairs of colorings of the vertices (or faces) of a regular tetrahedron with n available colors. Chiral colorings come in pairs, each the reflection of the other. - Robert A. Russell, Jan 22 2020
From Mircea Dan Rus, Aug 26 2020: (Start)
a(n+3) is the number of lattice rectangles (squares included) in a staircase of order n; this is obtained by stacking n rows of consecutive unit lattice squares, aligned either to the left or to the right, which consist of 1, 2, 3, ..., n squares and which are stacked either in the increasing or in the decreasing order of their lengths. Below, there is a staircase or order 4 which contains a(7) = 35 rectangles. [See the Teofil Bogdan and Mircea Dan Rus link, problem 3, under A004320]
_
||
|||_
|||_|_
|||_|_|
(End)
a(n+4) is the number of strings of length n on an ordered alphabet of 5 letters where the characters in the word are in nondecreasing order. E.g., number of length-2 words is 15: aa,ab,ac,ad,ae,bb,bc,bd,be,cc,cd,ce,dd,de,ee. - Jim Nastos, Jan 18 2021
From Tom Copeland, Jun 07 2021: (Start)
Aside from the zeros, this is the fifth diagonal of the Pascal matrix A007318, the only nonvanishing diagonal (fifth) of the matrix representation IM = (A132440)^4/4! of the differential operator D^4/4!, when acting on the row vector of coefficients of an o.g.f., or power series.
M = e^{IM} is the matrix of coefficients of the Appell sequence p_n(x) = e^{D^4/4!} x^n = e^{b. D} x^n = (b. + x)^n = Sum_{k=0..n} binomial(n,k) b_n x^{n-k}, where the (b.)^n = b_n have the e.g.f. e^{b.t} = e^{t^4/4!}, which is that for A025036 aerated with triple zeros, the first column of M.
See A099174 and A000292 for analogous relationships for the third and fourth diagonals of the Pascal matrix. (End)
For integer m and positive integer r >= 3, 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) = (3 - r*m)/2 in the complex plane. - Peter Bala, Jun 02 2024

Examples

			a(5) = 5 from the five independent components of an antisymmetric tensor A of rank 4 and dimension 5, namely A(1,2,3,4), A(1,2,3,5), A(1,2,4,5), A(1,3,4,5) and A(2,3,4,5). 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. Comtet, Advanced Combinatorics, Reidel, 1974, p. 74, Problem 8.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 70.
  • 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.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.6 Figurate Numbers, p. 294.
  • 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. 53, #191.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 127.

Crossrefs

binomial(n, k): A161680 (k = 2), A000389 (k = 5), A000579 (k = 6), A000580 (k = 7), A000581 (k = 8), A000582 (k = 9).
Cf. A000217, A000292, A007318 (column k = 4).
Cf. A158824.
Cf. A006008 (Number of ways to color the faces (or vertices) of a regular tetrahedron with n colors when mirror images are counted as two).
Cf. A104712 (third column, k=4).
See A269747 for a 3-D analog.
Cf. A006008 (oriented), A006003 (achiral) tetrahedron colorings.
Row 3 of A325000, col. 4 of A007318.

Programs

  • GAP
    A000332 := List([1..10^2], n -> Binomial(n, 4)); # Muniru A Asiru, Oct 16 2017
    
  • Magma
    [Binomial(n,4): n in [0..50]]; // Vincenzo Librandi, Nov 23 2014
    
  • Maple
    A000332 := n->binomial(n,4); [seq(binomial(n,4), n=0..100)];
  • Mathematica
    Table[ Binomial[n, 4], {n, 0, 45} ] (* corrected by Harvey P. Dale, Aug 22 2011 *)
    Table[(n-4)(n-3)(n-2)(n-1)/24, {n, 100}] (* Artur Jasinski, Dec 02 2007 *)
    LinearRecurrence[{5,-10,10,-5,1}, {0,0,0,0,1}, 45] (* Harvey P. Dale, Aug 22 2011 *)
    CoefficientList[Series[x^4 / (1 - x)^5, {x, 0, 40}], x] (* Vincenzo Librandi, Nov 23 2014 *)
  • PARI
    a(n)=binomial(n,4);
    
  • Python
    # Starts at a(3), i.e. computes n*(n+1)*(n+2)*(n+3)/24
    # which is more in line with A000217 and A000292.
    def A000332():
        x, y, z, u = 1, 1, 1, 1
        yield 0
        while True:
            yield x
            x, y, z, u = x + y + z + u + 1, y + z + u + 1, z + u + 1, u + 1
    a = A000332(); print([next(a) for i in range(41)]) # Peter Luschny, Aug 03 2019
    
  • Python
    print([n*(n-1)*(n-2)*(n-3)//24 for n in range(50)])
    # Gennady Eremin, Feb 06 2022

Formula

a(n) = n*(n-1)*(n-2)*(n-3)/24.
G.f.: x^4/(1-x)^5. - Simon Plouffe in his 1992 dissertation
a(n) = n*a(n-1)/(n-4). - Benoit Cloitre, Apr 26 2003, R. J. Mathar, Jul 07 2009
a(n) = Sum_{k=1..n-3} Sum_{i=1..k} i*(i+1)/2. - Benoit Cloitre, Jun 15 2003
Convolution of natural numbers {1, 2, 3, 4, ...} and A000217, the triangular numbers {1, 3, 6, 10, ...}. - Jon Perry, Jun 25 2003
a(n) = A110555(n+1,4). - Reinhard Zumkeller, Jul 27 2005
a(n+1) = ((n^5-(n-1)^5) - (n^3-(n-1)^3))/24 - (n^5-(n-1)^5-1)/30; a(n) = A006322(n-2)-A006325(n-1). - Xavier Acloque, Oct 20 2003; R. J. Mathar, Jul 07 2009
a(4*n+2) = Pyr(n+4, 4*n+2) where the polygonal pyramidal numbers are defined for integers A>2 and B>=0 by Pyr(A, B) = B-th A-gonal pyramid number = ((A-2)*B^3 + 3*B^2 - (A-5)*B)/6; For all positive integers i and the pentagonal number function P(x) = x*(3*x-1)/2: a(3*i-2) = P(P(i)) and a(3*i-1) = P(P(i) + i); 1 + 24*a(n) = (n^2 + 3*n + 1)^2. - Jonathan Vos Post, Nov 15 2004
First differences of A000389(n). - Alexander Adamchuk, Dec 19 2004
For n > 3, the sum of the first n-2 tetrahedral numbers (A000292). - Martin Steven McCormick (mathseq(AT)wazer.net), Apr 06 2005 [Corrected by Doug Bell, Jun 25 2017]
Starting (1, 5, 15, 35, ...), = binomial transform of [1, 4, 6, 4, 1, 0, 0, 0, ...]. - Gary W. Adamson, Dec 28 2007
Sum_{n>=4} 1/a(n) = 4/3, from the Taylor expansion of (1-x)^3*log(1-x) in the limit x->1. - R. J. Mathar, Jan 27 2009
A034263(n) = (n+1)*a(n+4) - Sum_{i=0..n+3} a(i). Also A132458(n) = a(n)^2 - a(n-1)^2 for n>0. - Bruno Berselli, Dec 29 2010
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5); a(0)=0, a(1)=0, a(2)=0, a(3)=0, a(4)=1. - Harvey P. Dale, Aug 22 2011
a(n) = (binomial(n-1,2)^2 - binomial(n-1,2))/6. - Gary Detlefs, Nov 20 2011
a(n) = Sum_{k=1..n-2} Sum_{i=1..k} i*(n-k-2). - Wesley Ivan Hurt, Sep 25 2013
a(n) = (A000217(A000217(n-2) - 1))/3 = ((((n-2)^2 + (n-2))/2)^2 - (((n-2)^2 + (n-2))/2))/(2*3). - Raphie Frank, Jan 16 2014
Sum_{n>=0} a(n)/n! = e/24. Sum_{n>=3} a(n)/(n-3)! = 73*e/24. See A067764 regarding the second ratio. - Richard R. Forberg, Dec 26 2013
Sum_{n>=4} (-1)^(n+1)/a(n) = 32*log(2) - 64/3 = A242023 = 0.847376444589... . - Richard R. Forberg, Aug 11 2014
4/(Sum_{n>=m} 1/a(n)) = A027480(m-3), for m>=4. - Richard R. Forberg, Aug 12 2014
E.g.f.: x^4*exp(x)/24. - Robert Israel, Nov 23 2014
a(n+3) = C(n,1) + 3*C(n,2) + 3*C(n,3) + C(n,4). Each term indicates the number of ways to use n colors to color a tetrahedron with exactly 1, 2, 3, or 4 colors.
a(n) = A080852(1,n-4). - R. J. Mathar, Jul 28 2016
From Gary W. Adamson, Feb 06 2017: (Start)
G.f.: Starting (1, 5, 14, ...), x/(1-x)^5 can be written
as (x * r(x) * r(x^2) * r(x^4) * r(x^8) * ...) where r(x) = (1+x)^5;
as (x * r(x) * r(x^3) * r(x^9) * r(x^27) * ...) where r(x) = (1+x+x^2)^5;
as (x * r(x) * r(x^4) * r(x^16) * r(x^64) * ...) where r(x) = (1+x+x^2+x^3)^5;
... (as a conjectured infinite set). (End)
From Robert A. Russell, Jan 22 2020: (Start)
a(n) = A006008(n) - a(n+3) = (A006008(n) - A006003(n)) / 2 = a(n+3) - A006003(n).
a(n+3) = A006008(n) - a(n) = (A006008(n) + A006003(n)) / 2 = a(n) + A006003(n).
a(n) = A007318(n,4).
a(n+3) = A325000(3,n). (End)
Product_{n>=5} (1 - 1/a(n)) = cosh(sqrt(15)*Pi/2)/(100*Pi). - Amiram Eldar, Jan 21 2021

Extensions

Some formulas that referred to another offset corrected by R. J. Mathar, Jul 07 2009

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

A145919 A000332(n) = a(n)*(3*a(n) - 1)/2.

Original entry on oeis.org

0, 0, 0, 0, 1, 2, -3, 5, 7, -9, 12, 15, -18, 22, 26, -30, 35, 40, -45, 51, 57, -63, 70, 77, -84, 92, 100, -108, 117, 126, -135, 145, 155, -165, 176, 187, -198, 210, 222, -234, 247, 260, -273, 287, 301, -315, 330, 345, -360, 376, 392, -408, 425, 442, -459, 477
Offset: 0

Views

Author

Matthew Vandermast, Oct 28 2008

Keywords

Comments

As the formula in the description shows, all members of A000332 belong to the generalized pentagonal sequence (A001318). A001318 also lists all nonnegative numbers that belong to A145919.

Examples

			a(6) = -3 and A000332(6) = (-3)(-10)/2 = 15.
		

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[(-x^4*(x^4 + 2*x^3 - 3*x^2 + 2*x + 1))/((x - 1)^3*(1 + x^2 + x)^3), {x,0,50}], x] (* G. C. Greubel, Jun 13 2017 *)
    LinearRecurrence[{0,0,3,0,0,-3,0,0,1},{0,0,0,0,1,2,-3,5,7},60] (* Harvey P. Dale, Feb 13 2023 *)
  • PARI
    my(x='x+O('x^50)); concat([0, 0, 0, 0], Vec((-x^4*(x^4 +2*x^3 -3*x^2 +2*x +1))/((x-1)^3*(1+x^2+x)^3))) \\ G. C. Greubel, Jun 13 2017

Formula

a(n+3) = A001840(n) when 3 does not divide n, A001840(n)*-1 otherwise.
After first two zeros, this sequence consists of all values of A001318(n) and A045943(n)*(-1), n>=0, sorted in order of increasing absolute value.
G.f.: (-x^4*(x^4+2*x^3-3*x^2+2*x+1))/((x-1)^3*(1+x^2+x)^3). - Maksym Voznyy (voznyy(AT)mail.ru), Jul 27 2009
a(n) = ((n^2-3*n+1)*(1-4*cos(2*Pi*n/3))+3)/18. - Natalia L. Skirrow, Apr 14 2025
Showing 1-4 of 4 results.