cp's OEIS Frontend

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

Showing 1-10 of 82 results. Next

A008585 a(n) = 3*n.

Original entry on oeis.org

0, 3, 6, 9, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48, 51, 54, 57, 60, 63, 66, 69, 72, 75, 78, 81, 84, 87, 90, 93, 96, 99, 102, 105, 108, 111, 114, 117, 120, 123, 126, 129, 132, 135, 138, 141, 144, 147, 150, 153, 156, 159, 162, 165, 168, 171, 174, 177
Offset: 0

Views

Author

Keywords

Comments

If n != 1 and n^2+2 is prime then n is a member of this sequence. - Cino Hilliard, Mar 19 2007
Multiples of 3. Positive members of this sequence are the third transversal numbers (or 3-transversal numbers): Numbers of the 3rd column of positive numbers in the square array of nonnegative and polygonal numbers A139600. Also, numbers of the 3rd column in the square array A057145. - Omar E. Pol, May 02 2008
Numbers n for which polynomial 27*x^6-2^n is factorizable. - Artur Jasinski, Nov 01 2008
1/7 in base-2 notation = 0.001001001... = 1/2^3 + 1/2^6 + 1/2^9 + ... - Gary W. Adamson, Jan 24 2009
A165330(a(n)) = 153 for n > 0; subsequence of A031179. - Reinhard Zumkeller, Sep 17 2009
A011655(a(n)) = 0. - Reinhard Zumkeller, Nov 30 2009
A215879(a(n)) = 0. - Reinhard Zumkeller, Dec 28 2012
Moser conjectured, and Newman proved, that the terms of this sequence are more likely to have an even number of 1s in binary than an odd number. The excess is an undulating multiple of n^(log 3/log 4). See also Coquet, who refines this result. - Charles R Greathouse IV, Jul 17 2013
Integer areas of medial triangles of integer-sided triangles.
Also integer subset of A188158(n)/4.
A medial triangle MNO is formed by joining the midpoints of the sides of a triangle ABC. The area of a medial triangle is A/4 where A is the area of the initial triangle ABC. - Michel Lagneau, Oct 28 2013
From Derek Orr, Nov 22 2014: (Start)
Let b(0) = 0, and b(n) = the number of distinct terms in the set of pairwise sums {b(0), ... b(n-1)} + {b(0), ... b(n-1)}. Then b(n+1) = a(n), for n > 0.
Example: b(1) = the number of distinct sums of {0} + {0}. The only possible sum is {0} so b(1) = 1. b(2) = the number of distinct sums of {0,1} + {0,1}. The possible sums are {0,1,2} so b(2) = 3. b(3) = the number of distinct sums of {0,1,3} + {0,1,3}. The possible sums are {0, 1, 2, 3, 4, 6} so b(3) = 6. This continues and one can see that b(n+1) = a(n). (End)
Number of partitions of 6n into exactly 2 parts. - Colin Barker, Mar 23 2015
Partial sums are in A045943. - Guenther Schrack, May 18 2017
Number of edges in a maximal planar graph with n+2 vertices, n > 0 (see A008486 comments). - Jonathan Sondow, Mar 03 2018
Also numbers such that when the leftmost digit is moved to the unit's place the result is divisible by 3. - Stefano Spezia, Jul 08 2025

Examples

			G.f.: 3*x + 6*x^2 + 9*x^3 + 12*x^4 + 15*x^5 + 18*x^6 + 21*x^7 + ...
		

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.

Crossrefs

Row / column 3 of A004247 and of A325820.
Cf. A016957, A057145, A139600, A139606, A001651 (complement), A032031 (partial products), A190944 (binary), A061819 (base 4).

Programs

Formula

G.f.: 3*x/(1-x)^2. - R. J. Mathar, Oct 23 2008
a(n) = A008486(n), n > 0. - R. J. Mathar, Oct 28 2008
G.f.: A(x) - 1, where A(x) is the g.f. of A008486. - Gennady Eremin, Feb 20 2021
a(n) = Sum_{k=0..inf} A030308(n,k)*A007283(k). - Philippe Deléham, Oct 17 2011
E.g.f.: 3*x*exp(x). - Ilya Gutkovskiy, May 18 2016
From Guenther Schrack, May 18 2017: (Start)
a(3*k) = a(a(k)) = A008591(n).
a(3*k+1) = a(a(k) + 1) = a(A016777(n)) = A017197(n).
a(3*k+2) = a(a(k) + 2) = a(A016789(n)) = A017233(n). (End)

Extensions

Partially edited by Joerg Arndt, Mar 11 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

Views

Author

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

A002411 Pentagonal pyramidal numbers: a(n) = n^2*(n+1)/2.

Original entry on oeis.org

0, 1, 6, 18, 40, 75, 126, 196, 288, 405, 550, 726, 936, 1183, 1470, 1800, 2176, 2601, 3078, 3610, 4200, 4851, 5566, 6348, 7200, 8125, 9126, 10206, 11368, 12615, 13950, 15376, 16896, 18513, 20230, 22050, 23976, 26011, 28158, 30420, 32800, 35301, 37926, 40678
Offset: 0

Views

Author

Keywords

Comments

a(n) = n^2(n+1)/2 is half the number of colorings of three points on a line with n+1 colors. - R. H. Hardin, Feb 23 2002
Sum of n smallest multiples of n. - Amarnath Murthy, Sep 20 2002
a(n) = number of (n+6)-bit binary sequences with exactly 7 1's none of which is isolated. A 1 is isolated if its immediate neighbor(s) are 0. - David Callan, Jul 15 2004
Also as a(n) = (1/6)*(3*n^3+3*n^2), n > 0: structured trigonal prism numbers (cf. A100177 - structured prisms; A100145 for more on structured numbers). - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Kekulé numbers for certain benzenoids. - Emeric Deutsch, Nov 18 2005
If Y is a 3-subset of an n-set X then, for n >= 5, a(n-4) is the number of 5-subsets of X having at least two elements in common with Y. - Milan Janjic, Nov 23 2007
a(n-1), n >= 2, is the number of ways to have n identical objects in m=2 of altogether n distinguishable boxes (n-2 boxes stay empty). - Wolfdieter Lang, Nov 13 2007
a(n+1) is the convolution of (n+1) and (3n+1). - Paul Barry, Sep 18 2008
The number of 3-character strings from an alphabet of n symbols, if a string and its reversal are considered to be the same.
Partial sums give A001296. - Jonathan Vos Post, Mar 26 2011
a(n-1):=N_1(n), n >= 1, is the number of edges of n planes in generic position in three-dimensional space. See a comment under A000125 for general arrangement. Comment to Arnold's problem 1990-11, see the Arnold reference, p.506. - Wolfdieter Lang, May 27 2011
Partial sums of pentagonal numbers A000326. - Reinhard Zumkeller, Jul 07 2012
From Ant King, Oct 23 2012: (Start)
For n > 0, the digital roots of this sequence A010888(A002411(n)) form the purely periodic 9-cycle {1,6,9,4,3,9,7,9,9}.
For n > 0, the units' digits of this sequence A010879(A002411(n)) form the purely periodic 20-cycle {1,6,8,0,5,6,6,8,5,0,6,6,3,0,0,6,1,8,0,0}.
(End)
a(n) is the number of inequivalent ways to color a path graph having 3 nodes using at most n colors. Note, here there is no restriction on the color of adjacent nodes as in the above comment by R. H. Hardin (Feb 23 2002). Also, here the structures are counted up to graph isomorphism, where as in the above comment the "three points on a line" are considered to be embedded in the plane. - Geoffrey Critzer, Mar 20 2013
After 0, row sums of the triangle in A101468. - Bruno Berselli, Feb 10 2014
Latin Square Towers: Take a Latin square of order n, with symbols from 1 to n, and replace each symbol x with a tower of height x. Then the total number of unit cubes used is a(n). - Arun Giridhar, Mar 29 2015
This is the case k = n+4 of b(n,k) = n*((k-2)*n-(k-4))/2, which is the n-th k-gonal number. Therefore, this is the 3rd upper diagonal of the array in A139600. - Luciano Ancora, Apr 11 2015
For n > 0, a(n) is the number of compositions of n+7 into n parts avoiding the part 2. - Milan Janjic, Jan 07 2016
Also the Wiener index of the n-antiprism graph. - Eric W. Weisstein, Sep 07 2017
For n > 0, a(2n+1) is the number of non-isomorphic 5C_m-snakes, where m = 2n+1 or m = 2n (for n >= 2). A kC_n-snake is a connected graph in which the k >= 2 blocks are isomorphic to the cycle C_n and the block-cutpoint graph is a path. - Christian Barrientos, May 15 2019
For n >= 1, a(n-1) is the number of 0°- and 45°-tilted squares that can be drawn by joining points in an n X n lattice. - Paolo Xausa, Apr 13 2021
a(n) is the number of all possible products of n rolls of a six-sided die. This can be easily seen by the recursive formula a(n) = a(n - 1) + 2 * binomial(n, 2) + binomial(n + 1, 2). - Rafal Walczak, Jun 15 2024
a(n) is the number of all triples consisting of nonnegative integers smaller than n such that the sum of the first two integers is less than n. - Ruediger Jehn, Aug 17 2025

Examples

			a(3)=18 because 4 identical balls can be put into m=2 of n=4 distinguishable boxes in binomial(4,2)*(2!/(1!*1!) + 2!/2!) = 6*(2+1) = 18 ways. The m=2 part partitions of 4, namely (1,3) and (2,2), specify the filling of each of the 6 possible two-box choices. - _Wolfdieter Lang_, Nov 13 2007
		

References

  • V. I. Arnold (ed.), Arnold's Problems, Springer, 2004, comments on Problem 1990-11 (p. 75), pp. 503-510. Numbers N_1.
  • Christian Barrientos, Graceful labelings of cyclic snakes, Ars Combin., Vol. 60 (2001), pp. 85-96.
  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 194.
  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (see p. 166, Table 10.4/I/5).
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93.
  • 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. 2.
  • 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

A006002(n) = -a(-1-n).
a(n) = A093560(n+2, 3), (3, 1)-Pascal column.
A row or column of A132191.
Second column of triangle A103371.
Cf. similar sequences listed in A237616.

Programs

  • GAP
    List([0..45], n->n^2*(n+1)/2); # Muniru A Asiru, Feb 19 2018
  • Haskell
    a002411 n = n * a000217 n  -- Reinhard Zumkeller, Jul 07 2012
    
  • Magma
    [n^2*(n+1)/2: n in [0..40]]; // Wesley Ivan Hurt, May 25 2014
    
  • Maple
    seq(n^2*(n+1)/2, n=0..40);
  • Mathematica
    Table[n^2 (n + 1)/2, {n, 0, 40}]
    LinearRecurrence[{4, -6, 4, -1}, {0, 1, 6, 18}, 50] (* Harvey P. Dale, Oct 20 2011 *)
    Nest[Accumulate, Range[1, 140, 3], 2] (* Vladimir Joseph Stephan Orlovsky, Jan 21 2012 *)
    CoefficientList[Series[x (1 + 2 x) / (1 - x)^4, {x, 0, 45}], x] (* Vincenzo Librandi, Jan 08 2016 *)
  • PARI
    a(n)=n^2*(n+1)/2
    
  • PARI
    concat(0, Vec(x*(1+2*x)/(1-x)^4 + O(x^100))) \\ Altug Alkan, Jan 07 2016
    

Formula

Average of n^2 and n^3.
G.f.: x*(1+2*x)/(1-x)^4. - Simon Plouffe in his 1992 dissertation
a(n) = n*Sum_{k=0..n} (n-k) = n*Sum_{k=0..n} k. - Paul Barry, Jul 21 2003
a(n) = n*A000217(n). - Xavier Acloque, Oct 27 2003
a(n) = (1/2)*Sum_{j=1..n} Sum_{i=1..n} (i+j) = (1/2)*(n^2+n^3) = (1/2)*A011379(n). - Alexander Adamchuk, Apr 13 2006
Row sums of triangle A127739, triangle A132118; and binomial transform of [1, 5, 7, 3, 0, 0, 0, ...] = (1, 6, 18, 40, 75, ...). - Gary W. Adamson, Aug 10 2007
G.f.: x*F(2,3;1;x). - Paul Barry, Sep 18 2008
Sum_{j>=1} 1/a(j) = hypergeom([1, 1, 1], [2, 3], 1) = -2 + 2*zeta(2) = A195055 - 2. - Stephen Crowley, Jun 28 2009
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4); a(0)=0, a(1)=1, a(2)=6, a(3)=18. - Harvey P. Dale, Oct 20 2011
From Ant King, Oct 23 2012: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 3.
a(n) = (n+1)*(2*A000326(n)+n)/6 = A000292(n) + 2*A000292(n-1).
a(n) = A000330(n)+A000292(n-1) = A000217(n) + 3*A000292(n-1).
a(n) = binomial(n+2,3) + 2*binomial(n+1,3).
(End)
a(n) = (A000330(n) + A002412(n))/2 = (A000292(n) + A002413(n))/2. - Omar E. Pol, Jan 11 2013
a(n) = (24/(n+3)!)*Sum_{j=0..n} (-1)^(n-j)*binomial(n,j)*j^(n+3). - Vladimir Kruchinin, Jun 04 2013
Sum_{n>=1} a(n)/n! = (7/2)*exp(1). - Richard R. Forberg, Jul 15 2013
E.g.f.: x*(2 + 4*x + x^2)*exp(x)/2. - Ilya Gutkovskiy, May 31 2016
From R. J. Mathar, Jul 28 2016: (Start)
a(n) = A057145(n+4,n).
a(n) = A080851(3,n-1). (End)
For n >= 1, a(n) = (Sum_{i=1..n} i^2) + Sum_{i=0..n-1} i^2*((i+n) mod 2). - Paolo Xausa, Apr 13 2021
a(n) = Sum_{k=1..n} GCD(k,n) * LCM(k,n). - Vaclav Kotesovec, May 22 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = 2 + Pi^2/6 - 4*log(2). - Amiram Eldar, Jan 03 2022

A003627 Primes of the form 3n-1.

Original entry on oeis.org

2, 5, 11, 17, 23, 29, 41, 47, 53, 59, 71, 83, 89, 101, 107, 113, 131, 137, 149, 167, 173, 179, 191, 197, 227, 233, 239, 251, 257, 263, 269, 281, 293, 311, 317, 347, 353, 359, 383, 389, 401, 419, 431, 443, 449, 461, 467, 479, 491, 503, 509, 521, 557, 563, 569, 587
Offset: 1

Views

Author

Keywords

Comments

Inert rational primes in the field Q(sqrt(-3)). - N. J. A. Sloane, Dec 25 2017
Primes p such that 1+x+x^2 is irreducible over GF(p). - Joerg Arndt, Aug 10 2011
Primes p dividing sum(k=0,p,C(2k,k)) -1 = A006134(p)-1. - Benoit Cloitre, Feb 08 2003
A039701(A049084(a(n))) = 2; A134323(A049084(a(n))) = -1. - Reinhard Zumkeller, Oct 21 2007
The set of primes of the form 3n - 1 is a superset of the set of lesser of twin primes larger than three (A001359). - Paul Muljadi, Jun 05 2008
Primes of this form do not occur in or as divisors of {n^2+n+1}. See A002383 (n^2+n+1 = prime), A162471 (prime divisors of n^2+n+1 not in A002383), and A002061 (numbers of the form n^2-n+1). - Daniel Tisdale, Jul 04 2009
Or, primes not in A007645. A003627 UNION A007645 = A000040. Also, primes of the form 6*k-5/2-+3/2. - Juri-Stepan Gerasimov, Jan 28 2010
Except for first term "2", all these prime numbers are of the form: 6*n-1. - Vladimir Joseph Stephan Orlovsky, Jul 13 2011
A088534(a(n)) = 0. - Reinhard Zumkeller, Oct 30 2011
For n>1: Numbers k such that (k-4)! mod k =(-1)^(floor(k/3)+1)*floor((k+1)/6), k>4. - Gary Detlefs, Jan 02 2012
Binomial(a(n),3)/a(n)= (3*A024893(n)^2+A024893(n))/2, n>1. - Gary Detlefs, May 06 2012
For every prime p in this sequence, 3 is a 9th power mod p. See Williams link. - Michel Marcus, Nov 12 2017
2 adjoined to A007528. - David A. Corneth, Nov 12 2017
For n >= 2 there exists a polygonal number P_s(3) = 3s - 3 = a(n) + 1. These are the only primes p with P_s(k) = p + 1, s >= 3, k >= 3, since P_s(k) - 1 is composite for k > 3. - Ralf Steiner, May 17 2018

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. 870.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Primes of form 3n+1 give A002476.
These are the primes arising in A024893, A087370, A088879. A091177 gives prime index.
Subsequence of A034020.

Programs

  • Haskell
    a003627 n = a003627_list !! (n-1)
    a003627_list = filter ((== 2) . (`mod` 3)) a000040_list
    -- Reinhard Zumkeller, Oct 30 2011
    
  • Magma
    [n: n in PrimesUpTo(720) | n mod 3 eq 2]; // Bruno Berselli, Apr 05 2011
    
  • Maple
    t1 := {}; for n from 0 to 500 do if isprime(3*n+2) then t1 := {op(t1),3*n+2}; fi; od: A003627 := convert(t1,list);
  • Mathematica
    Select[Range[-1, 600, 3], PrimeQ[#] &] (* Vincenzo Librandi, Jun 17 2015 *)
    Select[Prime[Range[200]],Mod[#,3]==2&] (* Harvey P. Dale, Jan 31 2023 *)
  • PARI
    is(n)=n%3==2 && isprime(n) \\ Charles R Greathouse IV, Mar 20 2013

Formula

From R. J. Mathar, Apr 03 2011: (Start)
Sum_{n>=1} 1/a(n)^2 = 0.30792... = A085548 - 1/9 - A175644.
Sum_{n>=1} 1/a(n)^3 = 0.134125... = A085541 - 1/27 - A175645. (End)

A007528 Primes of the form 6k-1.

Original entry on oeis.org

5, 11, 17, 23, 29, 41, 47, 53, 59, 71, 83, 89, 101, 107, 113, 131, 137, 149, 167, 173, 179, 191, 197, 227, 233, 239, 251, 257, 263, 269, 281, 293, 311, 317, 347, 353, 359, 383, 389, 401, 419, 431, 443, 449, 461, 467, 479, 491, 503, 509, 521, 557, 563, 569, 587
Offset: 1

Views

Author

Keywords

Comments

For values of k see A024898.
Also primes p such that p^q - 2 is not prime where q is an odd prime. These numbers cannot be prime because the binomial p^q = (6k-1)^q expands to 6h-1 some h. Then p^q-2 = 6h-1-2 is divisible by 3 thus not prime. - Cino Hilliard, Nov 12 2008
a(n) = A211890(3,n-1) for n <= 4. - Reinhard Zumkeller, Jul 13 2012
There exists a polygonal number P_s(3) = 3s - 3 = a(n) + 1. These are the only primes p with P_s(k) = p + 1, s >= 3, k >= 3, since P_s(k) - 1 is composite for k > 3. - Ralf Steiner, May 17 2018
From Bernard Schott, Feb 14 2019: (Start)
A theorem due to Andrzej Mąkowski: every integer greater than 161 is the sum of distinct primes of the form 6k-1. Examples: 162 = 5 + 11 + 17 + 23 + 47 + 59; 163 = 17 + 23 + 29 + 41 + 53. (See Sierpiński and David Wells.)
{2,3} Union A002476 Union {this sequence} = A000040.
Except for 2 and 3, all Sophie Germain primes are of the form 6k-1.
Except for 3, all the lesser of twin primes are also of the form 6k-1.
Dirichlet's theorem on arithmetic progressions states that this sequence is infinite. (End)
For all elements of this sequence p=6*k-1, there are no (x,y) positive integers such that k=6*x*y-x+y. - Pedro Caceres, Apr 06 2019

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. 870.
  • A. Mąkowski, Partitions into unequal primes, Bull. Acad. Polon. Sci. Sér. Sci. Math. Astr. Phys. 8 (1960), 125-126.
  • Wacław Sierpiński, Elementary Theory of Numbers, p. 144, Warsaw, 1964.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, Revised edition, 1997, p. 127.

Crossrefs

Intersection of A016969 and A000040.
Prime sequences A# (k,r) of the form k*n+r with 0 <= r <= k-1 (i.e., primes == r (mod k), or primes p with p mod k = r) and gcd(r,k)=1: A000040 (1,0), A065091 (2,1), A002476 (3,1), A003627 (3,2), A002144 (4,1), A002145 (4,3), A030430 (5,1), A045380 (5,2), A030431 (5,3), A030433 (5,4), A002476 (6,1), this sequence (6,5), A140444 (7,1), A045392 (7,2), A045437 (7,3), A045471 (7,4), A045458 (7,5), A045473 (7,6), A007519 (8,1), A007520 (8,3), A007521 (8,5), A007522 (8,7), A061237 (9,1), A061238 (9,2), A061239 (9,4), A061240 (9,5), A061241 (9,7), A061242 (9,8), A030430 (10,1), A030431 (10,3), A030432 (10,7), A030433 (10,9), A141849 (11,1), A090187 (11,2), A141850 (11,3), A141851 (11,4), A141852 (11,5), A141853 (11,6), A141854 (11,7), A141855 (11,8), A141856 (11,9), A141857 (11,10), A068228 (12,1), A040117 (12,5), A068229 (12,7), A068231 (12,11).
Cf. A034694 (smallest prime == 1 (mod n)).
Cf. A038700 (smallest prime == n-1 (mod n)).
Cf. A038026 (largest possible value of smallest prime == r (mod n)).
Cf. A001359 (lesser of twin primes), A005384 (Sophie Germain primes).

Programs

  • GAP
    Filtered(List([1..100],n->6*n-1),IsPrime); # Muniru A Asiru, May 19 2018
  • Haskell
    a007528 n = a007528_list !! (n-1)
    a007528_list = [x | k <- [0..], let x = 6 * k + 5, a010051' x == 1]
    -- Reinhard Zumkeller, Jul 13 2012
    
  • Maple
    select(isprime,[seq(6*n-1,n=1..100)]); # Muniru A Asiru, May 19 2018
  • Mathematica
    Select[6 Range[100]-1,PrimeQ]  (* Harvey P. Dale, Feb 14 2011 *)
  • PARI
    forprime(p=2, 1e3, if(p%6==5, print1(p, ", "))) \\ Charles R Greathouse IV, Jul 15 2011
    
  • PARI
    forprimestep(p=5,1000,6, print1(p", ")) \\ Charles R Greathouse IV, Mar 03 2025
    

Formula

A003627 \ {2}. - R. J. Mathar, Oct 28 2008
Conjecture: Product_{n >= 1} ((a(n) - 1) / (a(n) + 1)) * ((A002476(n) + 1) / (A002476(n) - 1)) = 3/4. - Dimitris Valianatos, Feb 11 2020
From Vaclav Kotesovec, May 02 2020: (Start)
Product_{k>=1} (1 - 1/a(k)^2) = 9*A175646/Pi^2 = 1/1.060548293.... =4/(3*A333240).
Product_{k>=1} (1 + 1/a(k)^2) = A334482.
Product_{k>=1} (1 - 1/a(k)^3) = A334480.
Product_{k>=1} (1 + 1/a(k)^3) = A334479. (End)
Legendre symbol (-3, a(n)) = -1 and (-3, A002476(n)) = +1, for n >= 1. For prime 3 one sets (-3, 3) = 0. - Wolfdieter Lang, Mar 03 2021

A017329 a(n) = 10*n + 5.

Original entry on oeis.org

5, 15, 25, 35, 45, 55, 65, 75, 85, 95, 105, 115, 125, 135, 145, 155, 165, 175, 185, 195, 205, 215, 225, 235, 245, 255, 265, 275, 285, 295, 305, 315, 325, 335, 345, 355, 365, 375, 385, 395, 405, 415, 425, 435, 445, 455, 465, 475, 485, 495, 505, 515, 525, 535
Offset: 0

Views

Author

Keywords

Comments

Continued fraction expansion of tanh(1/5). - Benoit Cloitre, Dec 17 2002
n such that 5 divides the numerator of B(2n) where B(2n) = the 2n-th Bernoulli number. - Benoit Cloitre, Jan 01 2004
5 times odd numbers. - Omar E. Pol, May 02 2008
5th transversal numbers (or 5-transversal numbers): Numbers of the 5th column of positive numbers in the square array of nonnegative and polygonal numbers A139600. Also, numbers of the 5th column in the square array A057145. - Omar E. Pol, May 02 2008
Successive sums: 5, 20, 45, 80, 125, ... (see A033429). - Philippe Deléham, Dec 08 2011
3^a(n) + 1 is divisible by 61. - Vincenzo Librandi, Feb 05 2013
If the initial 5 is changed to 1, giving 1,15,25,35,45,..., these are values of m such that A323288(m)/m reaches a new record high value. - N. J. A. Sloane, Jan 23 2019

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189. - From N. J. A. Sloane, Dec 01 2012

Crossrefs

Programs

Formula

a(n) = 5*A005408(n). - Omar E. Pol, Oct 19 2008
a(n) = 20*n - a(n-1) (with a(0)=5). - Vincenzo Librandi, Nov 19 2010
G.f.: 5*(x+1)/(x-1)^2. - Colin Barker, Nov 14 2012
a(n) = A057145(n+2,5). - R. J. Mathar, Jul 28 2016
E.g.f.: 5*exp(x)*(1 + 2*x). - Stefano Spezia, Feb 14 2020
Sum_{n>=0} (-1)^n/a(n) = Pi/20. - Amiram Eldar, Dec 12 2021
From Amiram Eldar, Nov 23 2024: (Start)
Product_{n>=0} (1 - (-1)^n/a(n)) = sqrt(5-sqrt(5))/2 = sqrt(2)*sin(Pi/5) = A182007/A002193.
Product_{n>=0} (1 + (-1)^n/a(n)) = phi/sqrt(2) (A094884). (End)
a(n) = (n+3)^2 - (n-2)^2. - Alexander Yutkin, Mar 16 2025
From Elmo R. Oliveira, Apr 12 2025: (Start)
a(n) = 2*a(n-1) - a(n-2).
a(n) = A008587(2*n+1). (End)

A016957 a(n) = 6*n + 4.

Original entry on oeis.org

4, 10, 16, 22, 28, 34, 40, 46, 52, 58, 64, 70, 76, 82, 88, 94, 100, 106, 112, 118, 124, 130, 136, 142, 148, 154, 160, 166, 172, 178, 184, 190, 196, 202, 208, 214, 220, 226, 232, 238, 244, 250, 256, 262, 268, 274, 280, 286, 292, 298, 304, 310, 316, 322, 328
Offset: 0

Views

Author

Keywords

Comments

Number of 2 X n binary matrices avoiding simultaneously the right-angled numbered polyomino patterns (ranpp) (00;1), (01,1) and (11;0). An occurrence of a ranpp (xy;z) in a matrix A=(a(i,j)) is a triple (a(i1,j1), a(i1,j2), a(i2,j1)) where i1 < i2, j1 < j2 and these elements are in the same relative order as those in the triple (x,y,z). In general, the number of m X n 0-1 matrices in question is given by (n+2)*2^(m-1) + 2*m*(n-1) - 2 for m > 1 and n > 1. - Sergey Kitaev, Nov 12 2004
If Y is a 4-subset of an n-set X then, for n >= 4, a(n-4) is the number of 3-subsets of X having at least two elements in common with Y. - Milan Janjic, Dec 08 2007
4th transversal numbers (or 4-transversal numbers): Numbers of the 4th column of positive numbers in the square array of nonnegative and polygonal numbers A139600. Also, numbers of the 4th column in the square array A057145. - Omar E. Pol, May 02 2008
a(n) is the maximum number such that there exists an edge coloring of the complete graph with a(n) vertices using n colors and every subgraph whose edges are of the same color (subgraph induced by edge color) is planar. - Srikanth K S, Dec 18 2010
Also numbers having two antecedents in the Collatz problem: 12*n+8 and 2*n+1 (respectively A017617(n) and A005408(n)). - Michel Lagneau, Dec 28 2012
a(n) = 6n+4 has three undirected edges e1 = (3n+2, 6n+4), e2 = (6n+4, 12n+8) and e3 = (2n+1, 6n+4) in the Collatz graph of A006370. - Heinz Ebert, Mar 16 2021
Conjecture: this sequence contains some but not all, even numbers with odd abundance A088827. They appear in this sequence at indices A186424(n) - 1. - John Tyler Rascoe, Jul 09 2022

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189. - From N. J. A. Sloane, Dec 01 2012

Crossrefs

Programs

Formula

A008615(a(n)) = n+1. - Reinhard Zumkeller, Feb 27 2008
a(n) = A016789(n)*2. - Omar E. Pol, May 02 2008
A157176(a(n)) = A067412(n+1). - Reinhard Zumkeller, Feb 24 2009
a(n) = sqrt(A016958(n)). - Zerinvary Lajos, Jun 30 2009
a(n) = 2*(6*n+1) - a(n-1) (with a(0)=4). - Vincenzo Librandi, Nov 20 2010
a(n) = floor((sqrt(36*n^2 - 36*n + 1) + 6*n + 1)/2). - Srikanth K S, Dec 18 2010
From Colin Barker, Jan 30 2012: (Start)
G.f.: 2*(2+x)/(1-2*x+x^2).
a(n) = 2*a(n-1) - a(n-2). (End)
A089911(2*a(n)) = 9. - Reinhard Zumkeller, Jul 05 2013
a(n) = 3 * A005408(n) + 1. - Fred Daniel Kline, Oct 24 2015
a(n) = A057145(n+2,4). - R. J. Mathar, Jul 28 2016
a(4*n+2) = 4 * a(n). - Zhandos Mambetaliyev, Sep 22 2018
Sum_{n>=0} (-1)^n/a(n) = sqrt(3)*Pi/18 - log(2)/6. - Amiram Eldar, Dec 10 2021
E.g.f.: 2*exp(x)*(2 + 3*x). - Stefano Spezia, May 29 2024

A139600 Square array T(n,k) = n*(k-1)*k/2+k, of nonnegative numbers together with polygonal numbers, read by antidiagonals upwards.

Original entry on oeis.org

0, 0, 1, 0, 1, 2, 0, 1, 3, 3, 0, 1, 4, 6, 4, 0, 1, 5, 9, 10, 5, 0, 1, 6, 12, 16, 15, 6, 0, 1, 7, 15, 22, 25, 21, 7, 0, 1, 8, 18, 28, 35, 36, 28, 8, 0, 1, 9, 21, 34, 45, 51, 49, 36, 9, 0, 1, 10, 24, 40, 55, 66, 70, 64, 45, 10, 0, 1, 11, 27, 46, 65, 81, 91, 92, 81, 55, 11
Offset: 0

Views

Author

Omar E. Pol, Apr 27 2008

Keywords

Comments

A general formula for polygonal numbers is P(n,k) = (n-2)*(k-1)*k/2 + k, where P(n,k) is the k-th n-gonal number.
The triangle sums, see A180662 for their definitions, link this square array read by antidiagonals with twelve different sequences, see the crossrefs. Most triangle sums are linear sums of shifted combinations of a sequence, see e.g. A189374. - Johannes W. Meijer, Apr 29 2011

Examples

			The square array of nonnegatives together with polygonal numbers begins:
=========================================================
....................... A   A   .   .   A    A    A    A
....................... 0   0   .   .   0    0    1    1
....................... 0   0   .   .   1    1    3    3
....................... 0   0   .   .   6    7    9    9
....................... 0   0   .   .   9    3    6    6
....................... 0   1   .   .   5    2    0    0
....................... 4   2   .   .   7    9    6    7
=========================================================
Nonnegatives . A001477: 0,  1,  2,  3,  4,   5,   6,   7, ...
Triangulars .. A000217: 0,  1,  3,  6, 10,  15,  21,  28, ...
Squares ...... A000290: 0,  1,  4,  9, 16,  25,  36,  49, ...
Pentagonals .. A000326: 0,  1,  5, 12, 22,  35,  51,  70, ...
Hexagonals ... A000384: 0,  1,  6, 15, 28,  45,  66,  91, ...
Heptagonals .. A000566: 0,  1,  7, 18, 34,  55,  81, 112, ...
Octagonals ... A000567: 0,  1,  8, 21, 40,  65,  96, 133, ...
9-gonals ..... A001106: 0,  1,  9, 24, 46,  75, 111, 154, ...
10-gonals .... A001107: 0,  1, 10, 27, 52,  85, 126, 175, ...
11-gonals .... A051682: 0,  1, 11, 30, 58,  95, 141, 196, ...
12-gonals .... A051624: 0,  1, 12, 33, 64, 105, 156, 217, ...
...
=========================================================
The column with the numbers 2, 3, 4, 5, 6, ... is formed by the numbers > 1 of A000027. The column with the numbers 3, 6, 9, 12, 15, ... is formed by the positive members of A008585.
		

Crossrefs

A formal extension negative n is in A326728.
Triangle sums (see the comments): A055795 (Row1), A080956 (Row2; terms doubled), A096338 (Kn11, Kn12, Kn13, Fi1, Ze1), A002624 (Kn21, Kn22, Kn23, Fi2, Ze2), A000332 (Kn3, Ca3, Gi3), A134393 (Kn4), A189374 (Ca1, Ze3), A011779 (Ca2, Ze4), A101357 (Ca4), A189375 (Gi1), A189376 (Gi2), A006484 (Gi4). - Johannes W. Meijer, Apr 29 2011
Sequences of m-gonal numbers: A000217 (m=3), A000290 (m=4), A000326 (m=5), A000384 (m=6), A000566 (m=7), A000567 (m=8), A001106 (m=9), A001107 (m=10), A051682 (m=11), A051624 (m=12), A051865 (m=13), A051866 (m=14), A051867 (m=15), A051868 (m=16), A051869 (m=17), A051870 (m=18), A051871 (m=19), A051872 (m=20), A051873 (m=21), A051874 (m=22), A051875 (m=23), A051876 (m=24), A255184 (m=25), A255185 (m=26), A255186 (m=27), A161935 (m=28), A255187 (m=29), A254474 (m=30).

Programs

  • Magma
    T:= func< n,k | k*(n*(k-1)+2)/2 >;
    A139600:= func< n,k | T(n-k, k) >;
    [A139600(n,k): k in  [0..n], n in [0..12]]; // G. C. Greubel, Jul 12 2024
    
  • Maple
    T:= (n, k)-> n*(k-1)*k/2+k:
    seq(seq(T(d-k, k), k=0..d), d=0..14);  # Alois P. Heinz, Oct 14 2018
  • Mathematica
    T[n_, k_] := (n + 1)*(k - 1)*k/2 + k; Table[T[n - k - 1, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Robert G. Wilson v, Jul 12 2009 *)
  • Python
    def A139600Row(n):
        x, y = 1, 1
        yield 0
        while True:
            yield x
            x, y = x + y + n, y + n
    for n in range(8):
        R = A139600Row(n)
        print([next(R) for  in range(11)]) # _Peter Luschny, Aug 04 2019
    
  • SageMath
    def T(n,k): return k*(n*(k-1)+2)/2
    def A139600(n,k): return T(n-k, k)
    flatten([[A139600(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Jul 12 2024

Formula

T(n,k) = n*(k-1)*k/2+k.
T(n,k) = A057145(n+2,k). - R. J. Mathar, Jul 28 2016
From Stefano Spezia, Apr 12 2024: (Start)
G.f.: y*(1 - x - y + 2*x*y)/((1 - x)^2*(1 - y)^3).
E.g.f.: exp(x+y)*y*(2 + x*y)/2. (End)

Extensions

Edited by Omar E. Pol, Jan 05 2009

A077028 The rascal triangle, read by rows: T(n,k) (n >= 0, 0 <= k <= n) = k(n-k) + 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 5, 4, 1, 1, 5, 7, 7, 5, 1, 1, 6, 9, 10, 9, 6, 1, 1, 7, 11, 13, 13, 11, 7, 1, 1, 8, 13, 16, 17, 16, 13, 8, 1, 1, 9, 15, 19, 21, 21, 19, 15, 9, 1, 1, 10, 17, 22, 25, 26, 25, 22, 17, 10, 1, 1, 11, 19, 25, 29, 31, 31, 29, 25, 19, 11, 1, 1, 12, 21, 28, 33, 36, 37, 36, 33, 28, 21, 12, 1
Offset: 0

Views

Author

Clark Kimberling, Oct 19 2002

Keywords

Comments

Pascal's triangle is formed using the rule South = West + East, whereas the rascal triangle uses the rule South = (West*East+1)/North. [Anggoro et al.]
The n-th diagonal is congruent to 1 mod n-1.
Row sums are the cake numbers, A000125. Alternating sum of row n is 0 if n even and (3-n)/2 if n odd. Rows are symmetric, beginning and ending with 1. The number of occurrences of k in this triangle is the number of divisors of k-1, given by A000005.
The triangle can be generated by numbers of the form k*(n-k) + 1 for k = 0 to n. Conjecture: except for n = 0,1 and 6 every row contains a prime. - Amarnath Murthy, Jul 15 2005
Above conjecture needs more exceptions, rows 30 and 54 do not contain primes. - Alois P. Heinz, Aug 31 2017
From Moshe Shmuel Newman, Apr 06 2008: (Start)
Consider the semigroup of words in x,y,q subject to the relationships: yx = xyq, qx = xq, qy = yq.
Now take words of length n in x and y, with exactly k y's. If there had been no relationships, the number of different words of this type would be n choose k, sequence A007318. Thanks to the relationships, the number of words of this type is the k-th entry in the n-th row of this sequence (read as a triangle, with the first row indexed by zero and likewise the first entry in each row.)
For example: with three letters and one y, we have three possibilities: xxy, xyx = xxyq, yxx = xxyqq. No two of them are equal, so this entry is still 3, as in Pascal's triangle.
With four letters, two y's, we have the first reduction: xyyx = yxxy = xxyyqq and this is the only reduction for 4 letters. So the middle entry of the fourth row is 5 instead of 6, as in the Pascal triangle. (End)
Main diagonals of this triangle sum to polygonal numbers. See A057145. - Raphie Frank, Oct 30 2012
T(n,k) gives the number of distinct sums of k elements in {1,2,...,n}, e.g., T(5,4) = the number of distinct sums of 4 elements in {1,2,3,4,5}, which is (5+4+3+2) - (4+3+2+1) + 1 = 5. - Derek Orr, Nov 26 2014
Conjecture: excluding the starting and ending 1's in each row, those that contain only prime numbers are n = 2, 3, 5, 7, 13, and 17. Tested up to row 10^9. - Rogério Serôdio, Sep 20 2017
The rascal triangle also uses the rule South = (West+East+1)-North. [Abstracts of AMS, Winter 2019, p. 526, 1145-VS-280, refers to Julian Fleron] - Michael Somos, Jan 12 2019
As a square array read by antidiagonals, selecting terms that give a remainder of 1 when divided by a prime gives evenly sized squares. For example, when each term is divided by 2, showing the remainder looks like:
1 1 1 1 1
1 0 1 0 1
1 1 1 1 1
1 0 1 0 1
1 1 1 1 1 - Nathaniel J. Strout, Jan 01 2020
T(n,k) is the number of binary words of length n which contain exactly k 1s and have at most 1 ascent. T(n,k) is also the number of ascent sequences avoiding 001 and 210 with length n+1 and exactly k ascents. - Amelia Gibbs, May 21 2024
T(n,k) represents the first and foundational instance R1 of a new family of Pascal-like triangles called Iterated Rascal triangles; A374378 is triangle R2; A374452 is triangle R3. - Kolosov Petro, Sep 28 2024

Examples

			Third diagonal (1,3,5,7,...) consists of the positive integers congruent to 1 mod 2.
The triangle T(n, k) begins:
  n\k  0  1  2  3  4  5  6  7  8  9 10 ...
  0:   1
  1:   1  1
  2:   1  2  1
  3:   1  3  3  1
  4:   1  4  5  4  1
  5:   1  5  7  7  5  1
  6:   1  6  9 10  9  6  1
  7:   1  7 11 13 13 11  7  1
  8:   1  8 13 16 17 16 13  8  1
  9:   1  9 15 19 21 21 19 15  9  1
 10:   1 10 17 22 25 26 25 22 17 10  1
 ... reformatted. - _Wolfdieter Lang_, Dec 19 2017
As a square array read by antidiagonals, the first rows are:
  1,  1,  1,  1,  1,  1, ...
  1,  2,  3,  4,  5,  6, ...
  1,  3,  5,  7,  9, 11, ...
  1,  4,  7, 10, 13, 16, ...
  1,  5,  9, 13, 17, 21, ...
		

Crossrefs

The maximum value for each antidiagonal is given by sequence A033638.
Equals A004247(n) + 1.

Programs

  • Maple
    A077028 := proc(n,k)
       if n <0 or k<0 or k > n then
           0;
       else
           k*(n-k)+1 ;
       end if;
    end proc: # R. J. Mathar, Jul 28 2016
  • Mathematica
    t[n_, k_] := k (n - k) + 1; t[0, 0] = 1; Table[ t[n, k], {n, 0, 12}, {k, 0, n}] // Flatten (* Robert G. Wilson v, Jul 06 2012 *)
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, k * (n - k) + 1)}; /* Michael Somos, Mar 20 2011 */

Formula

As a square array read by antidiagonals, a(n, k) = 1 + n*k. a(n, k) = a(n-1, k) + k. Row n has g.f. (1+(n-1)x)/(1-x)^2, n >= 0. - Paul Barry, Feb 22 2003
Still thinking of square arrays. Let f:N->Z and g:N->Z be given and I an integer, then define a(n, k) = I + f(n)*g(k). Then a(n, k)*a(n-1, k-1) = a(n-1, k)*a(n, k-1) + I*(f(n) - f(n-1))*(g(k) - g(k-1)) for suitable n and k. S= (E*W + 1)/N. arises with I = 1, and f = g = id. - Terry Lindgren, Apr 10 2011
Using the above: Having just read J. Fleron's nice article in Discovering the Art of Mathematics on the rascal triangle, it is neat to note and straightforward to show that when I = 1, a(n, k) + a(n-1, k-1) = a(n-1, k) + a(n, k-1) + (f(n) - f(n-1))*(g(k) - g(k-1)), so with I = 1, and f = g = id, we have S+N = E+W + 1, as his students discovered. - Terry Lindgren, Nov 28 2016
T(n, k) = A128139(n-1, k-1). - Gary W. Adamson, Jul 02 2012
O.g.f. (1 - x*(1 + t) + 2*t*x^2)/((1 - x)^2*(1 - t*x)^2) = 1 + (1 + t)*x + (1 + 2*t + t^2)*x^2 + .... Cf. A105851. - Peter Bala, Jul 26 2015
T(n, k) = 0 if n < k, T(n, 0) = 1, T(n,n) = 1, for n >= 0, and T(n, k) = (T(n-1, k-1)*T(n-1, k) + 1)/(T(n-2, k-1)) for 0 < k < n. See the first comment referring to the triangle with its apex in the middle. - Wolfdieter Lang, Dec 19 2017
E.g.f. as square array: exp(x+y)*(1 + x*y). - Stefano Spezia, Aug 10 2025

Extensions

Better definition based on Murthy's comment of Jul 15 2005 and the Anggoro et al. paper. - N. J. A. Sloane, Mar 05 2011

A060354 The n-th n-gonal number: a(n) = n*(n^2 - 3*n + 4)/2.

Original entry on oeis.org

0, 1, 2, 6, 16, 35, 66, 112, 176, 261, 370, 506, 672, 871, 1106, 1380, 1696, 2057, 2466, 2926, 3440, 4011, 4642, 5336, 6096, 6925, 7826, 8802, 9856, 10991, 12210, 13516, 14912, 16401, 17986, 19670, 21456, 23347, 25346, 27456, 29680, 32021
Offset: 0

Views

Author

Hareendra Yalamanchili (hyalaman(AT)mit.edu), Apr 01 2001

Keywords

Comments

Binomial transform of (0,1,0,3,0,0,0,...). - Paul Barry, Sep 14 2006
Also the number of permutations of length n which can be sorted by a single cut-and-paste move (in the sense of Cranston, Sudborough, and West). - Vincent Vatter, Aug 21 2013
Main diagonal of A317302. - Omar E. Pol, Aug 11 2018
a(n) is the number of ternary strings of length n that contain exactly one 1, zero or two 2's and have no restriction on the number of 0's. For example, a(5) = 35 since the strings are 12200 (30 of this type) and 10000 (5 of this type). - Enrique Navarrete, May 08 2025

Crossrefs

First differences of A004255.

Programs

  • Magma
    [(n*(n-2)^2+n^2)/2: n in [0..50]]; // Vincenzo Librandi, Feb 16 2015
  • Maple
    A060354 := proc(n)
        (n*(n-2)^2+n^2)/2 ;
    end proc: # R. J. Mathar, Jul 28 2016
  • Mathematica
    Table[(n (n-2)^2+n^2)/2,{n,0,50}] (* Harvey P. Dale, Aug 05 2011 *)
    CoefficientList[Series[x (1 - 2 x + 4 x^2) / (1 - x)^4, {x, 0, 50}], x] (* Vincenzo Librandi, Feb 16 2015 *)
    Table[PolygonalNumber[n,n],{n,0,50}] (* Harvey P. Dale, Mar 07 2016 *)
    LinearRecurrence[{4,-6,4,-1},{0,1,2,6},50] (* Harvey P. Dale, Mar 07 2016 *)
  • PARI
    a(n) = { (n*(n - 2)^2 + n^2)/2 } \\ Harry J. Smith, Jul 04 2009
    

Formula

a(n) = (n*(n-2)^2 + n^2)/2.
E.g.f.: exp(x)*x*(1+x^2/2). - Paul Barry, Sep 14 2006
G.f.: x*(1-2*x+4*x^2)/(1-x)^4. - R. J. Mathar, Sep 02 2008
a(n) = A057145(n,n). - R. J. Mathar, Jul 28 2016
a(n) = A000124(n-2) * n. - Bruce J. Nicholson, Jul 13 2018
a(n) = Sum_{i=0..n-1} (i*(n-2) + 1). - Ivan N. Ianakiev, Sep 25 2020
Showing 1-10 of 82 results. Next