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

A250120 Coordination sequence for planar net 3.3.3.3.6 (also called the fsz net).

Original entry on oeis.org

1, 5, 9, 15, 19, 24, 29, 33, 39, 43, 48, 53, 57, 63, 67, 72, 77, 81, 87, 91, 96, 101, 105, 111, 115, 120, 125, 129, 135, 139, 144, 149, 153, 159, 163, 168, 173, 177, 183, 187, 192, 197, 201, 207, 211, 216, 221, 225, 231, 235
Offset: 0

Views

Author

N. J. A. Sloane, Nov 23 2014

Keywords

Comments

There are eleven uniform (or Archimedean) tilings (or planar nets), with vertex symbols 3^6, 3^4.6, 3^3.4^2, 3^2.4.3.4, 4^4, 3.4.6.4, 3.6.3.6, 6^3, 3.12^2, 4.6.12, and 4.8^2. Grünbaum and Shephard (1987) is the best reference.
a(n) is the number of vertices at graph distance n from any fixed vertex.
The Mathematica notebook can compute 30 or 40 iterations, and colors them with period 5. You could also change out images if you want to. These graphs are better for analyzing 5-iteration chunks of the pattern. You can see that under iteration all fragments of the circumferences are preserved in shape and translated outwards a distance approximately sqrt(21) (relative to small triangle edge), the length of a long diagonal of larger rhombus unit cell. The conjectured recurrence should follow from an analysis of how new pieces occur in between the translated pieces. - Bradley Klee, Nov 26 2014

References

  • Branko Grünbaum and G. C. Shephard, Tilings and Patterns. W. H. Freeman, New York, 1987, Fig. 2.1.5, p. 63.
  • Marjorie Senechal, Quasicrystals and geometry, Cambridge University Press, Cambridge, 1995, Fig. 1.10, Section 1.3, pp. 13-16.

Crossrefs

List of coordination sequences for uniform planar nets: A008458 (the planar net 3.3.3.3.3.3), A008486 (6^3), A008574 (4.4.4.4 and 3.4.6.4), A008576 (4.8.8), A008579 (3.6.3.6), A008706 (3.3.3.4.4), A072154 (4.6.12), A219529 (3.3.4.3.4), A250120 (3.3.3.3.6), A250122 (3.12.12).
For partial sums of the present sequence, see A250121.

Programs

  • C
    /* Comments on the C program (see link) from Maurizio Paolini, Nov 23 2014: Basically what I do is deform the net onto the integral lattice, connect nodes aligned either horizontally, vertically or diagonally from northeast to southwest, marking as UNREACHABLE the nodes with coordinates (i, j) satisfying i + 2*j = 0 mod 7. Then the code computes the distance from each node to the central node of the grid. */
  • Mathematica
    CoefficientList[Series[(x^2+x+1)(x^4+3x^3+3x+1)/((x^4+x^3+x^2+x+1)(x-1)^2), {x, 0, 80}], x] (* or *) LinearRecurrence[{1, 0, 0, 0, 1, -1}, {1, 5, 9, 15, 19, 24, 29}, 60] (* Harvey P. Dale, May 05 2018 *)

Formula

Based on the computations of Darrah Chavey, Bradley Klee, and Maurizio Paolini, there is a strong conjecture that the first differences of this sequence are 4, 4, 6, 4, 5, 5, 4, 6, 4, 5, 5, 4, 6, 4, 5, 5, ..., that is, 4 followed by (4,6,4,5,5) repeated.
This would imply that the sequence satisfies the recurrence:
for n > 2, a(n) = a(n-1) + { n == 0,3 (mod 5), 4; n == 4 (mod 5), 6; n == 1,2 (mod 5), 5 }
(from Darrah Chavey)
and has generating function
(x^2+x+1)*(x^4+3*x^3+3*x+1)/((x^4+x^3+x^2+x+1)*(x-1)^2)
All the above conjectures are true - for proof see link to my article with Chaim Goodman-Strauss. - N. J. A. Sloane, Jan 14 2018; link added Mar 26 2018
a(n) ~ 24*n/5. - Stefano Spezia, May 08 2022
For n>0, a(n) = 2*(12*n + sqrt(1+2/sqrt(5))*sin(4*Pi*n/5) - sqrt(1-2/sqrt(5))*sin(2*Pi*n/5))/5. - Natalia L. Skirrow, Apr 13 2025

Extensions

a(6)-a(10) from Bradley Klee, Nov 23 2014
a(11)-a(49) from Maurizio Paolini, Nov 23 2014

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

A005448 Centered triangular numbers: a(n) = 3*n*(n-1)/2 + 1.

Original entry on oeis.org

1, 4, 10, 19, 31, 46, 64, 85, 109, 136, 166, 199, 235, 274, 316, 361, 409, 460, 514, 571, 631, 694, 760, 829, 901, 976, 1054, 1135, 1219, 1306, 1396, 1489, 1585, 1684, 1786, 1891, 1999, 2110, 2224, 2341, 2461, 2584, 2710, 2839, 2971, 3106, 3244, 3385, 3529
Offset: 1

Views

Author

N. J. A. Sloane, R. K. Guy, Dec 12 1974

Keywords

Comments

These are Hogben's central polygonal numbers
2
.P
3 n
Also the sum of three consecutive triangular numbers (A000217); i.e., a(4) = 19 = T4 + T3 + T2 = 10 + 6 + 3. - Robert G. Wilson v, Apr 27 2001
For k>2, Sum_{n=1..k} a(n) gives the sum pertaining to the magic square of order k. E.g., Sum_{n=1..5} a(n) = 1 + 4 + 10 + 19 + 31 = 65. In general, Sum_{n=1..k} a(n) = k*(k^2 + 1)/2. - Amarnath Murthy, Dec 22 2001
Binomial transform of (1,3,3,0,0,0,...). - Paul Barry, Jul 01 2003
a(n) is the difference of two tetrahedral (or pyramidal) numbers: C(n+3,3) = (n+1)(n+2)(n+3)/6. a(n) = A000292(n) - A000292(n-3) = (n+1)(n+2)(n+3)/6 - (n-2)(n-1)(n)/6. - Alexander Adamchuk, May 20 2006
Partial sums are A006003(n) = n(n^2+1)/2. Finite differences are a(n+1) - a(n) = A008585(n) = 3n. - Alexander Adamchuk, Jun 03 2006
If X is an n-set and Y a fixed 3-subset of X then a(n-2) is equal to the number of 3-subsets of X intersecting Y. - Milan Janjic, Jul 30 2007
Equals (1, 2, 3, ...) convolved with (1, 2, 3, 3, 3, ...). a(4) = 19 = (1, 2, 3, 4) dot (3, 3, 2, 1) = (3 + 6 + 6 + 4). - Gary W. Adamson, May 01 2009
Equals the triangular numbers convolved with [1, 1, 1, 0, 0, 0, ...]. - Gary W. Adamson and Alexander R. Povolotsky, May 29 2009
a(n) is the number of triples (w,x,y) having all terms in {0,...,n} and min(w+x,x+y,y+w) = max(w,x,y). - Clark Kimberling, Jun 14 2012
a(n) = number of atoms at graph distance <= n from an atom in the graphite or graphene network (cf. A008486). - N. J. A. Sloane, Jan 06 2013
In 1826, Shiraishi gave a solution to the Diophantine equation a^3 + b^3 + c^3 = d^3 with b = a(n) for n > 1; see A226903. - Jonathan Sondow, Jun 22 2013
For n > 1, a(n) is the remainder of n^2 * (n-1)^2 mod (n^2 + (n-1)^2). - J. M. Bergot, Jun 27 2013
The equation A000578(x) - A000578(x-1) = A000217(y) - A000217(y-2) is satisfied by y=a(x). - Bruno Berselli, Feb 19 2014
A242357(a(n)) = n. - Reinhard Zumkeller, May 11 2014
A255437(a(n)) = 1. - Reinhard Zumkeller, Mar 23 2015
The first differences give A008486. a(n) seems to give the total number of triangles in the n-th generation of the six patterns of triangle expansion shown in the link. - Kival Ngaokrajang, Sep 12 2015
Number of binary shuffle squares of length 2n which contains exactly two 1's. - Bartlomiej Pawlik, Sep 07 2023
The digital root has period 3 (1, 4, 1) (A146325), the same digital root as the centered 12-gonal numbers, or centered dodecagonal numbers A003154(n). - Peter M. Chema, Dec 20 2023

Examples

			From _Seiichi Manyama_, Aug 12 2017: (Start)
a(1) = 1:
      *
     / \
    /   \
   /     \
  *-------*
.................................................
a(2) = 4:
            *
           / \
          /   \
         /     \
        *---*---*
           / \
      *   /   \   *
     / \ /     \ / \
    /   *-------*   \
   /     \     /     \
  *-------*   *-------*
.................................................
a(3) = 10:
                  *
                 / \
                /   \
               /     \
              *---*---*
                 / \
            *   /   \   *
           / \ /     \ / \
          /   *---*---*   \
         /     \ / \ /     \
        *---*---*   *---*---*
           / \ /     \ / \
      *   /   *---*---*   \   *
     / \ /     \ / \ /     \ / \
    /   *-------*   *-------*   \
   /     \     /     \     /     \
  *-------*   *-------*   *-------*
.................................................
a(4) = 19:
                        *
                       / \
                      /   \
                     /     \
                    *---*---*
                       / \
                  *   /   \   *
                 / \ /     \ / \
                /   *---*---*   \
               /     \ / \ /     \
              *---*---*   *---*---*
                 / \ /     \ / \
            *   /   \---*---*   \   *
           / \ /     \ / \ /     \ / \
          /   *---*---*   *---*---*   \
         /     \ / \ /     \ / \ /     \
        *---*---*   *---*---*   *---*---*
           / \ /     \ / \ /     \ / \
      *   /   *---*---*   *---*---*   \   *
     / \ /     \ / \ /     \ / \ /     \ / \
    /   *-------*   *-------*   *-------*   \
   /     \     /     \     /     \     /     \
  *-------*   *-------*   *-------*   *-------*
(End)
		

References

  • R. Reed, The Lemming Simulation Problem, Mathematics in School, 3 (#6, Nov. 1974), front cover and pp. 5-6.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Haskell
    a005448 n = 3 * n * (n - 1) `div` 2 + 1
    a005448_list = 1 : zipWith (+) a005448_list [3, 6 ..]
    -- Reinhard Zumkeller, Jun 20 2013
    
  • Magma
    I:=[1,4,10]; [n le 3 select I[n] else 3*Self(n-1)-3*Self(n-2)+Self(n-3): n in [1..60]]; // Vincenzo Librandi, Sep 13 2015
  • Maple
    A005448 := n->(3*(n-1)^2+3*(n-1)+2)/2: seq(A005448(n), n=1..100);
    A005448 := -(1+z+z**2)/(z-1)^3; # Simon Plouffe in his 1992 dissertation for offset 0
  • Mathematica
    FoldList[#1 + #2 &, 1, 3 Range@ 50] (* Robert G. Wilson v, Feb 02 2011 *)
    Join[{1,4},Total/@Partition[Accumulate[Range[50]],3,1]] (* Harvey P. Dale, Aug 17 2012 *)
    LinearRecurrence[{3, -3, 1}, {1, 4, 10}, 50] (* Vincenzo Librandi, Sep 13 2015 *)
    Table[ j! Coefficient[Series[Exp[x]*(1 + 3 x^2/2)-1, {x, 0, 20}], x, j], {j, 0, 20}] (* Nikolaos Pantelidis, Feb 07 2023 *)
    3#+1&/@Accumulate[Range[0,50]] (* Harvey P. Dale, Nov 20 2024 *)
  • PARI
    {a(n)=3*(n^2-n)/2+1} /* Michael Somos, Sep 23 2006 */
    
  • PARI
    isok(n) = my(k=(2*n-2)/3, m); (n==1) || ((denominator(k)==1) && (m=sqrtint(k)) && (m*(m+1)==k)); \\ Michel Marcus, May 20 2020
    

Formula

Expansion of x*(1-x^3)/(1-x)^4.
a(n) = C(n+3, 3)-C(n, 3) = C(n, 0)+3*C(n, 1)+3*C(n, 2). - Paul Barry, Jul 01 2003
a(n) = 1 + Sum_{j=0..n-1} (3*j). - Xavier Acloque, Oct 25 2003
a(n) = A000217(n) + A000290(n-1) = (3*A016754(n) + 5)/8. - Lekraj Beedassy, Nov 05 2005
Euler transform of length 3 sequence [4, 0, -1]. - Michael Somos, Sep 23 2006
a(1-n) = a(n). - Michael Somos, Sep 23 2006
a(n) = binomial(n+1,n-1) + binomial(n,n-2) + binomial(n-1,n-3). - Zerinvary Lajos, Sep 03 2006
Row sums of triangle A134482. - Gary W. Adamson, Oct 27 2007
Narayana transform (A001263) * [1, 3, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(1)=1, a(2)=4, a(3)=10. - Jaume Oliver Lafont, Dec 02 2008
a(n) = A000217(n-1)*3 + 1 = A045943(n-1) + 1. - Omar E. Pol, Dec 27 2008
a(n) = a(n-1) + 3*n-3. - Vincenzo Librandi, Nov 18 2010
Sum_{n>=1} 1/a(n) = A306324. - Ant King, Jun 12 2012
a(n) = 2*a(n-1) - a(n-2) + 3. - Ant King, Jun 12 2012
a(n) = A101321(3,n-1). - R. J. Mathar, Jul 28 2016
E.g.f.: -1 + (2 + 3*x^2)*exp(x)/2. - Ilya Gutkovskiy, Jul 28 2016
a(n) = A002061(n) + A000217(n-1). - Bruce J. Nicholson, Apr 20 2017
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=1} a(n)/n! = 5*e/2 - 1.
Sum_{n>=1} (-1)^n * a(n)/n! = 5/(2*e) - 1. (End)
a(n) = A000326(n) - n + 1. - Charlie Marion, Nov 21 2020

A008574 a(0) = 1, thereafter a(n) = 4n.

Original entry on oeis.org

1, 4, 8, 12, 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 56, 60, 64, 68, 72, 76, 80, 84, 88, 92, 96, 100, 104, 108, 112, 116, 120, 124, 128, 132, 136, 140, 144, 148, 152, 156, 160, 164, 168, 172, 176, 180, 184, 188, 192, 196, 200, 204, 208, 212, 216, 220, 224, 228, 232
Offset: 0

Views

Author

N. J. A. Sloane; entry revised Aug 24 2014

Keywords

Comments

Number of squares on the perimeter of an (n+1) X (n+1) board. - Jon Perry, Jul 27 2003
Coordination sequence for square lattice (or equivalently the planar net 4.4.4.4).
Apparently also the coordination sequence for the planar net 3.4.6.4. - Darrah Chavey, Nov 23 2014
From N. J. A. Sloane, Nov 26 2014: (Start)
I confirm that this is indeed the coordination sequence for the planar net 3.4.6.4. The points at graph distance n from a fixed point in this net essentially lie on a hexagon (see illustration in link).
If n = 3k, k >= 1, there are 2k + 1 nodes on each edge of the hexagon. This counts the corners of the hexagon twice, so the number of points in the shell is 6(2k + 1) - 6 = 4n. If n = 3k + 1, the numbers of points on the six edges of the hexagon are 2k + 2 (4 times) and 2k + 1 (twice), for a total of 12k + 10 - 6 = 4n. If n = 3k + 2 the numbers are 2k + 2 (4 times) and 2k + 3 twice, and again we get 4n points.
The illustration shows shells 0 through 12, as well as the hexagons formed by shells 9 (green, 36 points), 10 (black, 40 points), 11 (red, 44 points), and 12 (blue, 48 points).
It is clear from the net that this period-3 structure continues forever, and establishes the theorem.
In contrast, for the 4.4.4.4 planar net, the successive shells are diamonds instead of hexagons, and again the n-th shell (n > 0) contains 4n points.
Of course the two nets are very different, since 4.4.4.4 has the symmetry of the square, while 3.4.6.4 has only mirror symmetry (with respect to a point), and has the symmetry of a regular hexagon with respect to the center of any of the 12-gons. (End)
Also the coordination sequence for a 6.6.6.6 point in the 3-transitive tiling {4.6.6, 6.6.6, 6.6.6.6}, see A265045, A265046. - N. J. A. Sloane, Dec 27 2015
Also the coordination sequence for 2-dimensional cyclotomic lattice Z[zeta_4].
Susceptibility series H_1 for 2-dimensional Ising model (divided by 2).
Also the Engel expansion of exp^(1/4); cf. A006784 for the Engel expansion definition. - Benoit Cloitre, Mar 03 2002
This sequence differs from A008586, multiples of 4, only in its initial term. - Alonso del Arte, Apr 14 2011
Number of 2 X n binary matrices avoiding simultaneously the right angled numbered polyomino patterns (ranpp) (00,0), (00;1) and (10;1). 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 and j1 < j2 and these elements are in same relative order as those in the triple (x,y,z). - Sergey Kitaev, Nov 11 2004
Central terms of the triangle in A118013. - Reinhard Zumkeller, Apr 10 2006
Also the coordination sequence for the htb net. - N. J. A. Sloane, Mar 31 2018
This is almost certainly also the coordination sequence for Dual(3.3.4.3.4) with respect to a tetravalent node. - Tom Karzes, Apr 01 2020
Minimal number of segments (equivalently, corners) in a rook circuit of a 2n X 2n board (maximal number is A085622). - Ruediger Jehn, Jan 02 2021

Examples

			From _Omar E. Pol_, Aug 20 2011 (Start):
Illustration of initial terms as perimeters of squares (cf. Perry's comment above):
.                                         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    4      8        12         16           20
(End)
		

Crossrefs

Cf. A001844 (partial sums), A008586, A054275, A054410, A054389, A054764.
Convolution square of A040000.
Row sums of A130323 and A131032.
List of coordination sequences for uniform planar nets: A008458 (the planar net 3.3.3.3.3.3), A008486 (6^3), A008574 (4.4.4.4 and 3.4.6.4), A008576 (4.8.8), A008579(3.6.3.6), A008706 (3.3.3.4.4), A072154 (4.6.12), A219529(3.3.4.3.4), A250120 (3.3.3.3.6), A250122 (3.12.12).
List of coordination sequences for Laves tilings (or duals of uniform planar nets): [3,3,3,3,3.3] = A008486; [3.3.3.3.6] = A298014, A298015, A298016; [3.3.3.4.4] = A298022, A298024; [3.3.4.3.4] = A008574, A296368; [3.6.3.6] = A298026, A298028; [3.4.6.4] = A298029, A298031, A298033; [3.12.12] = A019557, A298035; [4.4.4.4] = A008574; [4.6.12] = A298036, A298038, A298040; [4.8.8] = A022144, A234275; [6.6.6] = A008458.
Coordination sequences for the 20 2-uniform tilings in the order in which they appear in the Galebach catalog, together with their names in the RCSR database (two sequences per tiling): #1 krt A265035, A265036; #2 cph A301287, A301289; #3 krm A301291, A301293; #4 krl A301298, A298024; #5 krq A301299, A301301; #6 krs A301674, A301676; #7 krr A301670, A301672; #8 krk A301291, A301293; #9 krn A301678, A301680; #10 krg A301682, A301684; #11 bew A008574, A296910; #12 krh A301686, A301688; #13 krf A301690, A301692; #14 krd A301694, A219529; #15 krc A301708, A301710; #16 usm A301712, A301714; #17 krj A219529, A301697; #18 kre A301716, A301718; #19 krb A301720, A301722; #20 kra A301724, A301726.
See also A265045, A265046.

Programs

  • Haskell
    a008574 0 = 1; a008574 n = 4 * n
    a008574_list = 1 : [4, 8 ..]  -- Reinhard Zumkeller, Apr 16 2015
  • Mathematica
    f[0] = 1; f[n_] := 4 n; Array[f, 59, 0] (* or *)
    CoefficientList[ Series[(1 + x)^2/(1 - x)^2, {x, 0, 58}], x] (* Robert G. Wilson v, Jan 02 2011 *)
    Join[{1},Range[4,232,4]] (* Harvey P. Dale, Aug 19 2011 *)
    a[ n_] := 4 n + Boole[n == 0]; (* Michael Somos, Jan 07 2019 *)
  • PARI
    {a(n) = 4*n + !n}; /* Michael Somos, Apr 16 2007 */
    

Formula

Binomial transform is A000337 (dropping the 0 there). - Paul Barry, Jul 21 2003
Euler transform of length 2 sequence [4, -2]. - Michael Somos, Apr 16 2007
G.f.: ((1 + x) / (1 - x))^2. E.g.f.: 1 + 4*x*exp(x). - Michael Somos, Apr 16 2007
a(-n) = -a(n) unless n = 0. - Michael Somos, Apr 16 2007
G.f.: exp(4*atanh(x)). - Jaume Oliver Lafont, Oct 20 2009
a(n) = a(n-1) + 4, n > 1. - Vincenzo Librandi, Dec 31 2010
a(n) = A005408(n-1) + A005408(n), n > 1. - Ivan N. Ianakiev, Jul 16 2012
a(n) = 4*n = A008586(n), n >= 1. - Tom Karzes, Apr 01 2020

A008458 Coordination sequence for hexagonal lattice.

Original entry on oeis.org

1, 6, 12, 18, 24, 30, 36, 42, 48, 54, 60, 66, 72, 78, 84, 90, 96, 102, 108, 114, 120, 126, 132, 138, 144, 150, 156, 162, 168, 174, 180, 186, 192, 198, 204, 210, 216, 222, 228, 234, 240, 246, 252, 258, 264, 270, 276, 282, 288, 294, 300, 306, 312, 318, 324, 330, 336, 342, 348
Offset: 0

Views

Author

Keywords

Comments

The hexagonal lattice is the familiar 2-dimensional lattice in which each point has 6 neighbors. This is sometimes called the triangular lattice. It is also the planar net 3.3.3.3.3.3.
Coordination sequence for 2-dimensional cyclotomic lattice Z[zeta_6].
Apart from initial term(s), dimension of the space of weight 2n cusp forms for Gamma_0( 20 ).
Also the Engel expansion of exp^(1/6); cf. A006784 for the Engel expansion definition. - Benoit Cloitre, Mar 03 2002
Numbers k such that k+floor(k/2) | k*floor(k/2). - Wesley Ivan Hurt, Dec 01 2020

Examples

			From _Omar E. Pol_, Aug 20 2011: (Start)
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
. 1             o o o       o       o       o           o
.       6                    o o o o         o         o
.                 12                          o o o o o
.                               18
.                                                 24
(End)
G.f. = 1 + 6*x + 12*x^2 + 18*x^3 + 24*x^4 + 30*x^5 + 36*x^6 + 42*x^7 + 48*x^8 + 54*x^9 + ...
		

Crossrefs

Essentially the same as A008588.
List of coordination sequences for uniform planar nets: A008458 (the planar net 3.3.3.3.3.3), A008486 (6^3), A008574(4.4.4.4 and 3.4.6.4), A008576 (4.8.8), A008579(3.6.3.6), A008706 (3.3.3.4.4), A072154 (4.6.12), A219529(3.3.4.3.4), A250120 (3.3.3.3.6), A250122 (3.12.12).
List of coordination sequences for Laves tilings (or duals of uniform planar nets): [3,3,3,3,3.3] = A008486; [3.3.3.3.6] = A298014, A298015, A298016; [3.3.3.4.4] = A298022, A298024; [3.3.4.3.4] = A008574, A296368; [3.6.3.6] = A298026, A298028; [3.4.6.4] = A298029, A298031, A298033; [3.12.12] = A019557, A298035; [4.4.4.4] = A008574; [4.6.12] = A298036, A298038, A298040; [4.8.8] = A022144, A234275; [6.6.6] = A008458.
Cf. A032528. - Omar E. Pol, Aug 20 2011
Cf. A048477 (binomial Transf.)

Programs

  • Magma
    [0^n+6*n: n in [0..60] ]; // Vincenzo Librandi, Aug 21 2011
    
  • Maple
    1, seq(6*n, n=1..65);
  • Mathematica
    Join[{1},6*Range[60]] (* Harvey P. Dale, Jul 21 2013 *)
    a[ n_] := Boole[n == 0] + 6 n; (* Michael Somos, May 21 2015 *)
  • Maxima
    makelist(if n=0 then 1 else 6*n,n,0,65); /* Martin Ettl, Nov 12 2012 */
    
  • PARI
    {a(n) = 6*n + (!n)};
    
  • SageMath
    [6*n+int(n==0) for n in range(66)] # G. C. Greubel, May 25 2023

Formula

G.f.: (1 + 4*x + x^2)/(1 - x)^2.
a(n) = A003215(n) - A003215(n-1), n > 0.
Equals binomial transform of [1, 5, 1, -1, 1, -1, 1, ...]. - Gary W. Adamson, Jul 08 2008
G.f.: Hypergeometric2F1([3,-2], [1], -x/(1-x)). - Paul Barry, Sep 18 2008
a(n) = 0^n + 6*n. - Vincenzo Librandi, Aug 21 2011
n*a(1) + (n-1)*a(2) + (n-2)*a(3) + ... + 2*a(n-1) + a(n) = n^3. - Warren Breslow, Oct 28 2013
E.g.f.: 1 + 6*x*exp(x). - Stefano Spezia, Jun 26 2022

A298024 Expansion of (x^4+3*x^3+6*x^2+3*x+1)/((1-x)*(1-x^3)).

Original entry on oeis.org

1, 4, 10, 14, 18, 24, 28, 32, 38, 42, 46, 52, 56, 60, 66, 70, 74, 80, 84, 88, 94, 98, 102, 108, 112, 116, 122, 126, 130, 136, 140, 144, 150, 154, 158, 164, 168, 172, 178, 182, 186, 192, 196, 200, 206, 210, 214, 220, 224, 228, 234, 238, 242, 248, 252, 256, 262
Offset: 0

Views

Author

N. J. A. Sloane, Jan 21 2018

Keywords

Comments

Coordination sequence for Dual(3^3.4^2) tiling with respect to a tetravalent node. This tiling is also called the prismatic pentagonal tiling, or the cem-d net. It is one of the 11 Laves tilings. (The identification of this coordination sequence with the g.f. in the definition was first conjectured by Colin Barker, Jan 22 2018.)
Also, coordination sequence for a tetravalent node in the "krl" 2-D tiling (or net).
Both of these identifications are easily established using the "coloring book" method - see the Goodman-Strauss & Sloane link.
For n>0, this is twice A047386 (numbers congruent to 0 or +-2 mod 7).
Linear recurrence and g.f. confirmed by Shutov/Maleev link. - Ray Chandler, Aug 31 2023

References

  • Branko Grünbaum and G. C. Shephard, Tilings and Patterns. W. H. Freeman, New York, 1987. See Table 2.2.1, page 66, 3rd row, second tiling. (For the krl tiling.)
  • B. Gruenbaum and G. C. Shephard, Tilings and Patterns, W. H. Freeman, New York, 1987. See p. 96. (For the Dual(3^3.4^2) tiling.)

Crossrefs

Cf. A301298.
See A298025 for partial sums, A298022 for a trivalent node.
See also A047486.
List of coordination sequences for Laves tilings (or duals of uniform planar nets): [3,3,3,3,3.3] = A008486; [3.3.3.3.6] = A298014, A298015, A298016; [3.3.3.4.4] = A298022, A298024; [3.3.4.3.4] = A008574, A296368; [3.6.3.6] = A298026, A298028; [3.4.6.4] = A298029, A298031, A298033; [3.12.12] = A019557, A298035; [4.4.4.4] = A008574; [4.6.12] = A298036, A298038, A298040; [4.8.8] = A022144, A234275; [6.6.6] = A008458.
Coordination sequences for the 20 2-uniform tilings in the order in which they appear in the Galebach catalog, together with their names in the RCSR database (two sequences per tiling): #1 krt A265035, A265036; #2 cph A301287, A301289; #3 krm A301291, A301293; #4 krl A301298, A298024; #5 krq A301299, A301301; #6 krs A301674, A301676; #7 krr A301670, A301672; #8 krk A301291, A301293; #9 krn A301678, A301680; #10 krg A301682, A301684; #11 bew A008574, A296910; #12 krh A301686, A301688; #13 krf A301690, A301692; #14 krd A301694, A219529; #15 krc A301708, A301710; #16 usm A301712, A301714; #17 krj A219529, A301697; #18 kre A301716, A301718; #19 krb A301720, A301722; #20 kra A301724, A301726.

Programs

  • Mathematica
    CoefficientList[Series[(x^4+3x^3+6x^2+3x+1)/((1-x)(1-x^3)),{x,0,60}],x] (* or *) LinearRecurrence[{1,0,1,-1},{1,4,10,14,18},80] (* Harvey P. Dale, Oct 03 2018 *)
  • PARI
    See Links section.

Formula

a(n) = a(n-1) + a(n-3) - a(n-4) for n>4. (Conjectured, correctly, by Colin Barker, Jan 22 2018.)

Extensions

More terms from Rémy Sigrist, Jan 21 2018
Entry revised by N. J. A. Sloane, Mar 25 2018

A219529 Coordination sequence for 3.3.4.3.4 Archimedean tiling.

Original entry on oeis.org

1, 5, 11, 16, 21, 27, 32, 37, 43, 48, 53, 59, 64, 69, 75, 80, 85, 91, 96, 101, 107, 112, 117, 123, 128, 133, 139, 144, 149, 155, 160, 165, 171, 176, 181, 187, 192, 197, 203, 208, 213, 219, 224, 229, 235, 240, 245, 251, 256, 261, 267, 272, 277, 283, 288, 293, 299
Offset: 0

Views

Author

Allan C. Wechsler, Nov 21 2012

Keywords

Comments

a(n) is the number of vertices of the 3.3.4.3.4 tiling (which has three triangles and two squares, in the given cyclic order, meeting at each vertex) whose shortest path connecting them to a given origin vertex contains n edges.
This is the dual tiling to the Cairo tiling (cf. A296368). - N. J. A. Sloane, Nov 02 2018
First few terms provided by Allan C. Wechsler; Fred Lunnon and Fred Helenius gave the next few; Fred Lunnon suggested that the recurrence was a(n+3) = a(n) + 16 for n > 1. [This conjecture is true - see the CGS-NJAS link for a proof. - N. J. A. Sloane, Dec 31 2017]
Appears also to be coordination sequence for node of type V2 in "krd" 2-D tiling (or net). This should be easy to prove by the coloring book method (see link). - N. J. A. Sloane, Mar 25 2018
Appears also to be coordination sequence for node of type V1 in "krj" 2-D tiling (or net). This also should be easy to prove by the coloring book method (see link). - N. J. A. Sloane, Mar 26 2018
First differences of A301696. - Klaus Purath, May 23 2020

References

  • Branko Grünbaum and G. C. Shephard, Tilings and Patterns. W. H. Freeman, New York, 1987. See Table 2.2.1, page 67, 1st row, 2nd tiling, also 2nd row, third tiling.

Crossrefs

List of coordination sequences for uniform planar nets: A008458 (the planar net 3.3.3.3.3.3), A008486 (6^3), A008574 (4.4.4.4 and 3.4.6.4), A008576 (4.8.8), A008579 (3.6.3.6), A008706 (3.3.3.4.4), A072154 (4.6.12), A219529 (3.3.4.3.4), A250120 (3.3.3.3.6), A250122 (3.12.12).
Coordination sequences for the 20 2-uniform tilings in the order in which they appear in the Galebach catalog, together with their names in the RCSR database (two sequences per tiling): #1 krt A265035, A265036; #2 cph A301287, A301289; #3 krm A301291, A301293; #4 krl A301298, A298024; #5 krq A301299, A301301; #6 krs A301674, A301676; #7 krr A301670, A301672; #8 krk A301291, A301293; #9 krn A301678, A301680; #10 krg A301682, A301684; #11 bew A008574, A296910; #12 krh A301686, A301688; #13 krf A301690, A301692; #14 krd A301694, A219529; #15 krc A301708, A301710; #16 usm A301712, A301714; #17 krj A219529, A301697; #18 kre A301716, A301718; #19 krb A301720, A301722; #20 kra A301724, A301726.

Programs

  • Haskell
    -- Very slow, could certainly be accelerated.  SST stands for Snub Square Tiling.
    setUnion [] l2 = l2
    setUnion (a:rst) l2 = if (elem a l2) then doRest else (a:doRest)
      where doRest = setUnion rst l2
    setDifference [] l2 = []
    setDifference (a:rst) l2 = if (elem a l2) then doRest else (a:doRest)
      where doRest = setDifference rst l2
    adjust k = (if (even k) then 1 else -1)
    weirdAdjacent (x,y) = (x+(adjust y),y+(adjust x))
    sstAdjacents (x,y) = [(x+1,y),(x-1,y),(x,y+1),(x,y-1),(weirdAdjacent (x,y))]
    sstNeighbors core = foldl setUnion core (map sstAdjacents core)
    sstGlob n core = if (n == 0) then core else (sstGlob (n-1) (sstNeighbors core))
    sstHalo core = setDifference (sstNeighbors core) core
    origin = [(0,0)]
    a219529 n = length (sstHalo (sstGlob (n-1) origin))
    -- Allan C. Wechsler, Nov 30 2012
    
  • Maple
    A219529:= n -> `if`(n=0, 1, (16*n +1 - `mod`(n+1,3))/3);
    seq(A219529(n), n = 0..60); # G. C. Greubel, May 27 2020
  • Mathematica
    Join[{1}, LinearRecurrence[{1,0,1,-1}, {5,11,16,21}, 60]] (* Jean-François Alcover, Dec 13 2018 *)
    Table[If[n==0, 1, (16*n +1 - Mod[n+1, 3])/3], {n, 0, 60}] (* G. C. Greubel, May 27 2020 *)
    CoefficientList[Series[(x+1)^4/((x^2+x+1)(x-1)^2),{x,0,70}],x] (* Harvey P. Dale, Jul 03 2021 *)
  • Sage
    [1]+[(16*n+1 -(n+1)%3)/3 for n in (1..60)] # G. C. Greubel, May 27 2020

Formula

Conjectured to be a(n) = floor((16n+1)/3) for n>0; a(0) = 1; this is a consequence of the suggested recurrence due to Lunnon (see comments). [This conjecture is true - see the CGS-NJAS link in A296368 for a proof. - N. J. A. Sloane, Dec 31 2017]
G.f.: (x+1)^4/((x^2+x+1)*(x-1)^2). - N. J. A. Sloane, Feb 07 2018
From G. C. Greubel, May 27 2020: (Start)
a(n) = (16*n - ChebyshevU(n-1, -1/2))/3 for n>0 with a(0)=1.
a(n) = (A008598(n) - A049347(n-1))/3 for n >0 with a(0)=1. (End)

Extensions

Corrected attributions and epistemological status in Comments; provided slow Haskell code - Allan C. Wechsler, Nov 30 2012
Extended by Joseph Myers, Dec 04 2014

A008576 Coordination sequence for planar net 4.8.8.

Original entry on oeis.org

1, 3, 5, 8, 11, 13, 16, 19, 21, 24, 27, 29, 32, 35, 37, 40, 43, 45, 48, 51, 53, 56, 59, 61, 64, 67, 69, 72, 75, 77, 80, 83, 85, 88, 91, 93, 96, 99, 101, 104, 107, 109, 112, 115, 117, 120, 123, 125, 128, 131, 133
Offset: 0

Views

Author

Keywords

Comments

Also, growth series for the affine Coxeter (or Weyl) groups B_2. - N. J. A. Sloane, Jan 11 2016

References

  • N. Bourbaki, Groupes et Algèbres de Lie, Chap. 4, 5 and 6, Hermann, Paris, 1968. See Chap. VI, Section 4, Problem 10b, page 231, W_a(t).
  • A. V. Shutov, On the number of words of a given length in plane crystallographic groups (Russian), Zap. Nauchn. Sem. S.-Peterburg. Otdel. Mat. Inst. Steklov. (POMI) 302 (2003), Anal. Teor. Chisel i Teor. Funkts. 19, 188--197, 203; translation in J. Math. Sci. (N.Y.) 129 (2005), no. 3, 3922-3926 [MR2023041]. See Table 1.

Crossrefs

List of coordination sequences for uniform planar nets: A008458 (the planar net 3.3.3.3.3.3), A008486 (6^3), A008574 (4.4.4.4 and 3.4.6.4), A008576 (4.8.8), A008579(3.6.3.6), A008706 (3.3.3.4.4), A072154 (4.6.12), A219529(3.3.4.3.4), A250120 (3.3.3.3.6), A250122 (3.12.12).
For partial sums see A008577.
The growth series for the finite Coxeter (or Weyl) groups B_3 through B_12 are A161696-A161699, A161716, A161717, A161733, A161755, A161776, A161858. These are all rows of A128084. The growth series for the affine Coxeter (or Weyl) groups B_2 through B_12 are A008576, A008137, A267167-A267175.

Programs

  • Maple
    if n mod 3 = 0 then 8*n/3 elif n mod 3 = 1 then 8*(n-1)/3+3 else 8*(n-2)/3+5 fi;
  • Mathematica
    cspn[n_]:=Module[{c=Mod[n,3]},Which[c==0,(8n)/3,c==1,(8(n-1))/3+3,True,(8(n-2))/3+5]]; Join[{1},Array[cspn,50]] (* or *) Join[{1}, LinearRecurrence[ {1,0,1,-1},{3,5,8,11},50]] (* Harvey P. Dale, Nov 24 2011 *)
  • PARI
    a(n)=([0,1,0,0; 0,0,1,0; 0,0,0,1; -1,1,0,1]^n*[1;3;5;8])[1,1] \\ Charles R Greathouse IV, Apr 08 2016

Formula

G.f.: ((1+x)^2*(1+x^2))/((1-x)^2*(1+x+x^2)). - Ralf Stephan, Apr 24 2004
a(0)=1, a(1)=3, a(2)=5, a(3)=8, a(4)=11, a(n) = a(n-1) + a(n-3) - a(n-4). - Harvey P. Dale, Nov 24 2011
a(0)=1; thereafter a(3k)=8k, a(3k+1)=8k+3, a(3k+2)=8k+5. - N. J. A. Sloane, Dec 22 2015
The above g.f. and recurrence were originally empirical observations, but I now have a proof (details will be added later). This also justifies the Maple and Mma programs and the b-file. - N. J. A. Sloane, Dec 22 2015
Sum of alternate terms of A042965 (numbers not congruent to 2 mod 4), such that A042965(n) = A042965(n+1) + A042965(n-1). - Gary W. Adamson, Sep 12 2007
a(n) = (2/9)*(12*n + (3/2)*A102283(n)) for n > 0. - Stefano Spezia, Aug 07 2022

A004277 1 together with positive even numbers.

Original entry on oeis.org

1, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42, 44, 46, 48, 50, 52, 54, 56, 58, 60, 62, 64, 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96, 98, 100, 102, 104, 106, 108, 110, 112, 114, 116, 118, 120, 122, 124, 126, 128, 130, 132
Offset: 0

Views

Author

Keywords

Comments

Also number of non-attacking bishops on n X n board. - Koksal Karakus (karakusk(AT)hotmail.com), May 27 2002
Engel expansion of e^(1/2) (see A006784 for definition) [when offset by 1]. - Henry Bottomley, Dec 18 2000
Numbers n such that a 2n-group (i.e., a group of order 2n) has subgroup C_2. - Lekraj Beedassy, Oct 14 2004
Image of 1/(1-2x) under the mapping g(x)->g(x/(1+x^2)). - Paul Barry, Jan 16 2005
Position of n in A113322: A113322(a(n-1)) = n for n>0. - Reinhard Zumkeller, Oct 26 2005
Incrementally largest terms in the continued fraction for e. - Nick Hobson, Jan 11 2007
Conjecturally, the differences of two consecutive primes (without repetition). - Juri-Stepan Gerasimov, Nov 09 2009
Equals (1, 2, 2, 2, ...) convolved with (1, 0, 2, 0, 2, 0, 2, ...). - Gary W. Adamson, Mar 03 2010
a(n) is the number of 0-dimensional elements (vertices) in an n-cross polytope. - Patrick J. McNab, Jul 06 2015
Numbers k such that in the symmetric representation of sigma(k) there is no pair bars as its ends (Cf. A237593). - Omar E. Pol, Sep 28 2018
Also, the coordination sequence of the L-lattice (see A332419). - Sean A. Irvine, Jul 29 2020

Crossrefs

INVERT transformation yields A098182 without A098182(0). - R. J. Mathar, Sep 11 2008

Programs

Formula

G.f.: (1+x^2)/(1-x)^2. - Paul Barry, Feb 28 2003
Inverse binomial transform of Cullen numbers A002064. a(n)=2n+0^n. - Paul Barry, Jun 12 2003
a(n) = Sum_{k=0..floor(n/2)} binomial(n-k-1)*(-1)^k*2^(n-2k). - Paul Barry, Jan 16 2005
Equals binomial transform of [1, 1, 1, -1, 1, -1, 1, ...]. - Gary W. Adamson, Jul 15 2008
E.g.f.: 1+x*sinh(x) (aerated sequence). - Paul Barry, Oct 11 2009
a(n) = 0^n + 2*n = A000007(n) + A005843(n). - Reinhard Zumkeller, Jan 11 2012

Extensions

Corrected by Charles R Greathouse IV, Mar 18 2010
Showing 1-10 of 85 results. Next