cp's OEIS Frontend

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

Previous Showing 31-40 of 286 results. Next

A016921 a(n) = 6*n + 1.

Original entry on oeis.org

1, 7, 13, 19, 25, 31, 37, 43, 49, 55, 61, 67, 73, 79, 85, 91, 97, 103, 109, 115, 121, 127, 133, 139, 145, 151, 157, 163, 169, 175, 181, 187, 193, 199, 205, 211, 217, 223, 229, 235, 241, 247, 253, 259, 265, 271, 277, 283, 289, 295, 301, 307, 313, 319, 325, 331
Offset: 0

Views

Author

Keywords

Comments

Apart from initial term(s), dimension of the space of weight 2n cusp forms for Gamma_0( 22 ).
Also solutions to 2^x + 3^x == 5 (mod 7). - Cino Hilliard, May 10 2003
Except for 1, exponents n > 1 such that x^n - x^2 - 1 is reducible. - N. J. A. Sloane, Jul 19 2005
Let M(n) be the n X n matrix m(i,j) = min(i,j); then the trace of M(n)^(-2) is a(n-1) = 6*n - 5. - Benoit Cloitre, Feb 09 2006
If Y is a 3-subset of an (2n+1)-set X then, for n >= 3, a(n-1) is the number of 3-subsets of X having at least two elements in common with Y. - Milan Janjic, Dec 16 2007
All composite terms belong to A269345 as shown in there. - Waldemar Puszkarz, Apr 13 2016
First differences of the number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 773", based on the 5-celled von Neumann neighborhood. - Robert Price, May 23 2016
For b(n) = A103221(n) one has b(a(n)-1) = b(a(n)+1) = b(a(n)+2) = b(a(n)+3) = b(a(n)+4) = n+1 but b(a(n)) = n. So-called "dips" in A103221. See the Avner and Gross remark on p. 178. - Wolfdieter Lang, Sep 16 2016
A (n+1,n) pebbling move involves removing n + 1 pebbles from a vertex in a simple graph and placing n pebbles on an adjacent vertex. A two-player impartial (n+1,n) pebbling game involves two players alternating (n+1,n) pebbling moves. The first player unable to make a move loses. The sequence a(n) is also the minimum number of pebbles such that any assignment of those pebbles on a complete graph with 3 vertices is a next-player winning game in the two player impartial (k+1,k) pebbling game. These games are represented by A347637(3,n). - Joe Miller, Oct 18 2021
Interleaving of A017533 and A017605. - Leo Tavares, Nov 16 2021

Examples

			From _Ilya Gutkovskiy_, Apr 15 2016: (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
n=0  n=1     n=2     n=3
(End)
		

References

  • Avner Ash and Robert Gross, Summing it up, Princeton University Press, 2016, p. 178.

Crossrefs

Cf. A093563 ((6, 1) Pascal, column m=1).
a(n) = A007310(2*(n+1)); complement of A016969 with respect to A007310.
Cf. A287326 (second column).

Programs

Formula

a(n) = 6*n + 1, n >= 0 (see the name).
G.f.: (1+5*x)/(1-x)^2.
A008615(a(n)) = n. - Reinhard Zumkeller, Feb 27 2008
A157176(a(n)) = A013730(n). - Reinhard Zumkeller, Feb 24 2009
a(n) = 4*(3*n-1) - a(n-1) (with a(0)=1). - Vincenzo Librandi, Nov 20 2010
E.g.f.: (1 + 6*x)*exp(x). - G. C. Greubel, Sep 18 2019
a(n) = A003215(n) - 6*A000217(n-1). See Hexagonal Lines illustration. - Leo Tavares, Sep 10 2021
From Leo Tavares, Oct 27 2021: (Start)
a(n) = 6*A001477(n-1) + 7
a(n) = A016813(n) + 2*A001477(n)
a(n) = A017605(n-1) + A008588(n-1)
a(n) = A016933(n) - 1
a(n) = A008588(n) + 1. (End)
Sum_{n>=0} (-1)^n/a(n) = Pi/6 + sqrt(3)*arccoth(sqrt(3))/3. - Amiram Eldar, Dec 10 2021

A001107 10-gonal (or decagonal) numbers: a(n) = n*(4*n-3).

Original entry on oeis.org

0, 1, 10, 27, 52, 85, 126, 175, 232, 297, 370, 451, 540, 637, 742, 855, 976, 1105, 1242, 1387, 1540, 1701, 1870, 2047, 2232, 2425, 2626, 2835, 3052, 3277, 3510, 3751, 4000, 4257, 4522, 4795, 5076, 5365, 5662, 5967, 6280, 6601, 6930, 7267, 7612, 7965, 8326
Offset: 0

Views

Author

Keywords

Comments

Write 0, 1, 2, ... in a square spiral, with 0 at the origin and 1 immediately below it; sequence gives numbers on the negative y-axis (see Example section).
Number of divisors of 48^(n-1) for n > 0. - J. Lowell, Aug 30 2008
a(n) is the Wiener index of the graph obtained by connecting two copies of the complete graph K_n by an edge (for n = 3, approximately: |>-<|). The Wiener index of a connected graph is the sum of the distances between all unordered pairs of vertices in the graph. - Emeric Deutsch, Sep 20 2010
This sequence does not contain any squares other than 0 and 1. See A188896. - T. D. Noe, Apr 13 2011
For n > 0: right edge of the triangle A033293. - Reinhard Zumkeller, Jan 18 2012
Sequence found by reading the line from 0, in the direction 0, 10, ... and the parallel line from 1, in the direction 1, 27, ..., in the square spiral whose vertices are the generalized decagonal numbers A074377. - Omar E. Pol, Jul 18 2012
Partial sums give A007585. - Omar E. Pol, Jan 15 2013
This is also a star pentagonal number: a(n) = A000326(n) + 5*A000217(n-1). - Luciano Ancora, Mar 28 2015
Also the number of undirected paths in the n-sunlet graph. - Eric W. Weisstein, Sep 07 2017
After 0, a(n) is the sum of 2*n consecutive integers starting from n-1. - Bruno Berselli, Jan 16 2018
Number of corona of an H0 hexagon with a T(n) triangle. - Craig Knecht, Dec 13 2024

Examples

			On a square lattice, place the nonnegative integers at lattice points forming a spiral as follows: place "0" at the origin; then move one step downward (i.e., in the negative y direction) and place "1" at the lattice point reached; then turn 90 degrees in either direction and place a "2" at the next lattice point; then make another 90-degree turn in the same direction and place a "3" at the lattice point; etc. The terms of the sequence will lie along the negative y-axis, as seen in the example below:
  99  64--65--66--67--68--69--70--71--72
   |   |                               |
  98  63  36--37--38--39--40--41--42  73
   |   |   |                       |   |
  97  62  35  16--17--18--19--20  43  74
   |   |   |   |               |   |   |
  96  61  34  15   4---5---6  21  44  75
   |   |   |   |   |       |   |   |   |
  95  60  33  14   3  *0*  7  22  45  76
   |   |   |   |   |   |   |   |   |   |
  94  59  32  13   2--*1*  8  23  46  77
   |   |   |   |           |   |   |   |
  93  58  31  12--11-*10*--9  24  47  78
   |   |   |                   |   |   |
  92  57  30--29--28-*27*-26--25  48  79
   |   |                           |   |
  91  56--55--54--53-*52*-51--50--49  80
   |                                   |
  90--89--88--87--86-*85*-84--83--82--81
[Edited by _Jon E. Schoenfield_, Jan 02 2017]
		

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
  • Bruce C. Berndt, Ramanujan's Notebooks, Part II, Springer; see p. 23.
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6.
  • S. M. Ellerstein, The square spiral, J. Recreational Mathematics 29 (#3, 1998) 188; 30 (#4, 1999-2000), 246-250.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 2nd ed., 1994, p. 99.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A093565 ((8, 1) Pascal, column m = 2). Partial sums of A017077.
Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.
Cf. A003215.

Programs

  • Magma
    [4*n^2-3*n : n in [0..50] ]; // Wesley Ivan Hurt, Jun 05 2014
    
  • Maple
    A001107:=-(1+7*z)/(z-1)**3; # Simon Plouffe in his 1992 dissertation
  • Mathematica
    LinearRecurrence[{3, -3, 1}, {0, 1, 10}, 60] (* Harvey P. Dale, May 08 2012 *)
    Table[PolygonalNumber[RegularPolygon[10], n], {n, 0, 46}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
    Table[4 n^2 - 3 n, {n, 0, 49}] (* Alonso del Arte, Jan 24 2017 *)
    PolygonalNumber[10, Range[0, 20]] (* Eric W. Weisstein, Sep 07 2017 *)
    LinearRecurrence[{3, -3, 1}, {1, 10, 27}, {0, 20}] (* Eric W. Weisstein, Sep 07 2017 *)
  • PARI
    a(n)=4*n^2-3*n
    
  • Python
    a=lambda n: 4*n**2-3*n # Indranil Ghosh, Jan 01 2017
    def aList(): # Intended to compute the initial segment of the sequence, not isolated terms.
         x, y = 1, 1
         yield 0
         while True:
             yield x
             x, y = x + y + 8, y + 8
    A001107 = aList()
    print([next(A001107) for i in range(49)]) # Peter Luschny, Aug 04 2019

Formula

a(n) = A033954(-n) = A074377(2*n-1).
a(n) = n + 8*A000217(n-1). - Floor van Lamoen, Oct 14 2005
G.f.: x*(1 + 7*x)/(1 - x)^3.
Partial sums of odd numbers 1 mod 8, i.e., 1, 1 + 9, 1 + 9 + 17, ... . - Jon Perry, Dec 18 2004
1^3 + 3^3*(n-1)/(n+1) + 5^3*((n-1)*(n-2))/((n+1)*(n+2)) + 7^3*((n-1)*(n-2)*(n-3))/((n+1)*(n+2)*(n+3)) + ... = n*(4*n-3) [Ramanujan]. - Neven Juric, Apr 15 2008
Starting (1, 10, 27, 52, ...), this is the binomial transform of [1, 9, 8, 0, 0, 0, ...]. - Gary W. Adamson, Apr 30 2008
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n>2, a(0)=0, a(1)=1, a(2)=10. - Jaume Oliver Lafont, Dec 02 2008
a(n) = 8*n + a(n-1) - 7 for n>0, a(0)=0. - Vincenzo Librandi, Jul 10 2010
a(n) = 8 + 2*a(n-1) - a(n-2). - Ant King, Sep 04 2011
a(n) = A118729(8*n). - Philippe Deléham, Mar 26 2013
a(8*a(n) + 29*n+1) = a(8*a(n) + 29*n) + a(8*n + 1). - Vladimir Shevelev, Jan 24 2014
Sum_{n >= 1} 1/a(n) = Pi/6 + log(2) = 1.216745956158244182494339352... = A244647. - Vaclav Kotesovec, Apr 27 2016
From Ilya Gutkovskiy, Aug 28 2016: (Start)
E.g.f.: x*(1 + 4*x)*exp(x).
Sum_{n >= 1} (-1)^(n+1)/a(n) = (sqrt(2)*Pi - 2*log(2) + 2*sqrt(2)*log(1 + sqrt(2)))/6 = 0.92491492293323294695... (End)
a(n) = A000217(3*n-2) - A000217(n-2). In general, if P(k,n) be the n-th k-gonal number and T(n) be the n-th triangular number, A000217(n), then P(T(k),n) = T((k-1)*n - (k-2)) - T(k-3)*T(n-2). - Charlie Marion, Sep 01 2020
Product_{n>=2} (1 - 1/a(n)) = 4/5. - Amiram Eldar, Jan 21 2021
a(n) = A003215(n-1) + A000290(n) - 1. - Leo Tavares, Jul 23 2022

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

Original entry on oeis.org

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

Views

Author

Keywords

Comments

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

Examples

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

References

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

Crossrefs

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

Programs

Formula

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

A005898 Centered cube numbers: n^3 + (n+1)^3.

Original entry on oeis.org

1, 9, 35, 91, 189, 341, 559, 855, 1241, 1729, 2331, 3059, 3925, 4941, 6119, 7471, 9009, 10745, 12691, 14859, 17261, 19909, 22815, 25991, 29449, 33201, 37259, 41635, 46341, 51389, 56791, 62559, 68705, 75241, 82179, 89531, 97309, 105525, 114191, 123319, 132921
Offset: 0

Views

Author

Keywords

Comments

Write the natural numbers in groups: 1; 2,3,4; 5,6,7,8,9; 10,11,12,13,14,15,16; ..... and add the groups, i.e., a(n) = Sum_{j=n^2-2(n-1)..n^2} j. - Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Sep 05 2001
The numbers 1, 9, 35, 91, etc. are divisible by 1, 3, 5, 7, etc. Therefore there are no prime numbers in this list. 9 is divisible by 3 and every third number after 9 is also divisible by 3. 35 is divisible by 5 and 7 and every fifth number after 35 is also divisible by 5 and every seventh number after 35 is also divisible by 7. This pattern continues indefinitely. - Howard Berman (howard_berman(AT)hotmail.com), Nov 07 2008
n^3 + (n+1)^3 = (2n+1)*(n^2+n+1), hence all terms are composite. - Zak Seidov, Feb 08 2011
This is the order of an n-ball centered at a node in the Kronecker product (or direct product) of three cycles, each of whose lengths is at least 2n+2. - Pranava K. Jha, Oct 10 2011
Positive y values of 4*x^3 - 3*x^2 = y^2. - Bruno Berselli, Apr 28 2018

References

  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 52.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

(1/12)*t*(2*n^3 - 3*n^2 + n) + 2*n - 1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A005917, A063493, A063494, A063495, A063496.
Partial sums of A005897.
The 28 uniform 3D tilings: cab: A299266, A299267; crs: A299268, A299269; fcu: A005901, A005902; fee: A299259, A299265; flu-e: A299272, A299273; fst: A299258, A299264; hal: A299274, A299275; hcp: A007899, A007202; hex: A005897, A005898; kag: A299256, A299262; lta: A008137, A299276; pcu: A005899, A001845; pcu-i: A299277, A299278; reo: A299279, A299280; reo-e: A299281, A299282; rho: A008137, A299276; sod: A005893, A005894; sve: A299255, A299261; svh: A299283, A299284; svj: A299254, A299260; svk: A010001, A063489; tca: A299285, A299286; tcd: A299287, A299288; tfs: A005899, A001845; tsi: A299289, A299290; ttw: A299257, A299263; ubt: A299291, A299292; bnn: A007899, A007202. See the Proserpio link in A299266 for overview.

Programs

Formula

a(n) = Sum_{i=0..n} A005897(i), partial sums. - Jonathan Vos Post, Feb 06 2011
G.f.: (x^2+4*x+1)*(1+x)/(1-x)^3. - Simon Plouffe (see MAPLE section) and Colin Barker, Jan 02 2012; edited by N. J. A. Sloane, Feb 07 2018
a(n) = A037270(n+1) - A037270(n). - Ivan N. Ianakiev, May 13 2012
a(n) = A000217(n+1)^2 - A000217(n-1)^2. - Bob Selcoe, Mar 25 2016
a(n) = A005408(n) * A002061(n+1). - Miquel Cerda, Oct 05 2016
From Ilya Gutkovskiy, Oct 06 2016: (Start)
E.g.f.: (1 + 8*x + 9*x^2 + 2*x^3)*exp(x).
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). (End)
a(n) = (A081435(n))^2 - (A081435(n) - 1)^2. - Sergey Pavlov, Mar 01 2017

A033428 a(n) = 3*n^2.

Original entry on oeis.org

0, 3, 12, 27, 48, 75, 108, 147, 192, 243, 300, 363, 432, 507, 588, 675, 768, 867, 972, 1083, 1200, 1323, 1452, 1587, 1728, 1875, 2028, 2187, 2352, 2523, 2700, 2883, 3072, 3267, 3468, 3675, 3888, 4107, 4332, 4563, 4800, 5043, 5292, 5547, 5808, 6075, 6348
Offset: 0

Views

Author

Keywords

Comments

The number of edges of a complete tripartite graph of order 3n, K_n,n,n. - Roberto E. Martinez II, Oct 18 2001
From Floor van Lamoen, Jul 21 2001: (Start)
Write 1,2,3,4,... in a hexagonal spiral around 0; then a(n) is the sequence found by reading the line from 0 in the direction 0,3,.... The spiral begins:
.
33--32--31--30
/ \
34 16--15--14 29
/ / \ \
35 17 5---4 13 28
/ / / \ \ \
36 18 6 0---3--12--27--48-->
/ / / / / / / /
37 19 7 1---2 11 26 47
\ \ \ / / /
38 20 8---9--10 25 46
\ \ / /
39 21--22--23--24 45
\ /
40--41--42--43--44
(End)
Number of edges of the complete bipartite graph of order 4n, K_n,3n. - Roberto E. Martinez II, Jan 07 2002
Also the number of partitions of 6n + 3 into at most 3 parts. - R. K. Guy, Oct 23 2003
Also the number of partitions of 6n into exactly 3 parts. - Colin Barker, Mar 23 2015
Numbers n such that the imaginary quadratic field Q[sqrt(-n)] has six units. - Marc LeBrun, Apr 12 2006
The denominators of Hoehn's sequence (recalled by G. L. Honaker, Jr.) and the numerators of that sequence reversed. The sequence is 1/3, (1+3)/(5+7), (1+3+5)/(7+9+11), (1+3+5+7)/(9+11+13+15), ...; reduced to 1/3, 4/12, 9/27, 16/48, ... . For the reversal, the reduction is 3/1, 12/4, 27/9, 48/16, ... . - Enoch Haga, Oct 05 2007
Right edge of tables in A200737 and A200741: A200737(n, A000292(n)) = A200741(n, A100440(n)) = a(n). - Reinhard Zumkeller, Nov 21 2011
The Wiener index of the crown graph G(n) (n>=3). The crown graph G(n) is the graph with vertex set {x(1), x(2), ..., x(n), y(1), y(2), ..., y(n)} and edge set {(x(i), y(j)): 1<=i, j<=n, i/=j} (= the complete bipartite graph K(n,n) with horizontal edges removed). Example: a(3)=27 because G(3) is the cycle C(6) and 6*1 + 6*2 + 3*3 = 27. The Hosoya-Wiener polynomial of G(n) is n(n-1)(t+t^2)+nt^3. - Emeric Deutsch, Aug 29 2013
From Michel Lagneau, May 04 2015: (Start)
Integer area A of equilateral triangles whose side lengths are in the commutative ring Z[3^(1/4)] = {a + b*3^(1/4) + c*3^(1/2) + d*3^(3/4), a,b,c and d in Z}.
The area of an equilateral triangle of side length k is given by A = k^2*sqrt(3)/4. In the ring Z[3^(1/4)], if k = q*3^(1/4), then A = 3*q^2/4 is an integer if q is even. Example: 27 is in the sequence because the area of the triangle (6*3^(1/4), 6*3^(1/4), 6*3^(1/4)) is 27. (End)
a(n) is 2*sqrt(3) times the area of a 30-60-90 triangle with short side n. Also, 3 times the area of an n X n square. - Wesley Ivan Hurt, Apr 06 2016
Consider the hexagonal tiling of the plane. Extract any four hexagons adjacent by edge. This can be done in three ways. Fold the four hexagons so that all opposite faces occupy parallel planes. For all parallel projections of the resulting object, at least two correspond to area a(n) for side length of n of the original hexagons. - Torlach Rush, Aug 17 2022
The sequence terms are the exponents in the expansion of Product_{n >= 1} (1 - q^(3*n))/(1 + q^(3*n)) = ( Sum_{n in Z} q^(n*(3*n+1)/2) ) / ( Product_{n >= 1} 1 + q^n ) = 1 - 2*q^3 + 2*q^12 - 2*q^27 + 2*q^48 - 2*q^75 + - .... - Peter Bala, Dec 30 2024

Examples

			From _Ilya Gutkovskiy_, Apr 13 2016: (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   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
. n=1      n=2            n=3                 n=4
(End)
		

Crossrefs

Programs

  • Haskell
    a033428 = (* 3) . (^ 2)
    a033428_list = 0 : 3 : 12 : zipWith (+) a033428_list
       (map (* 3) $ tail $ zipWith (-) (tail a033428_list) a033428_list)
    -- Reinhard Zumkeller, Jul 11 2013
    
  • Magma
    [3*n^2: n in [0..50]]; // Vincenzo Librandi, May 18 2015
    
  • Maple
    seq(3*n^2, n=0..46); # Nathaniel Johnston, Jun 26 2011
  • Mathematica
    3 Range[0, 50]^2
    LinearRecurrence[{3, -3, 1}, {0, 3, 12}, 50] (* Harvey P. Dale, Feb 16 2013 *)
  • Maxima
    makelist(3*n^2,n,0,30); /* Martin Ettl, Nov 12 2012 */
    
  • PARI
    a(n)=3*n^2
    
  • Python
    def a(n): return 3 * (n**2) # Torlach Rush, Aug 25 2022

Formula

a(n) = 3*a(n-1)-3*a(n-2)+a(n-3) for n>2.
G.f.: 3*x*(1+x)/(1-x)^3. - R. J. Mathar, Sep 09 2008
Main diagonal of triangle in A132111: a(n) = A132111(n,n). - Reinhard Zumkeller, Aug 10 2007
A214295(a(n)) = -1. - Reinhard Zumkeller, Jul 12 2012
a(n) = A215631(n,n) for n > 0. - Reinhard Zumkeller, Nov 11 2012
a(n) = A174709(6n+2). - Philippe Deléham, Mar 26 2013
a(n) = a(n-1) + 6*n - 3, with a(0)=0. - Jean-Bernard François, Oct 04 2013
E.g.f.: 3*x*(1 + x)*exp(x). - Ilya Gutkovskiy, Apr 13 2016
a(n) = t(3*n) - 3*t(n), where t(i) = i*(i+k)/2 for any k. Special case (k=1): A000217(3*n) - 3*A000217(n). - Bruno Berselli, Aug 31 2017
a(n) = A000326(n) + A005449(n). - Bruce J. Nicholson, Jan 10 2020
From Amiram Eldar, Jul 03 2020: (Start)
Sum_{n>=1} 1/a(n) = Pi^2/18 (A086463).
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/36. (End)
From Amiram Eldar, Feb 03 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = sqrt(3)*sinh(Pi/sqrt(3))/Pi.
Product_{n>=1} (1 - 1/a(n)) = sqrt(3)*sin(Pi/sqrt(3))/Pi. (End)
a(n) = A003215(n) - A016777(n). - Leo Tavares, Apr 29 2023

Extensions

Better description from N. J. A. Sloane, May 15 1998

A211422 Number of ordered triples (w,x,y) with all terms in {-n,...,0,...,n} and w^2 + x*y = 0.

Original entry on oeis.org

1, 9, 17, 25, 41, 49, 57, 65, 81, 105, 113, 121, 137, 145, 153, 161, 193, 201, 225, 233, 249, 257, 265, 273, 289, 329, 337, 361, 377, 385, 393, 401, 433, 441, 449, 457, 505, 513, 521, 529, 545, 553, 561, 569, 585, 609, 617, 625, 657, 713, 753, 761
Offset: 0

Views

Author

Clark Kimberling, Apr 10 2012

Keywords

Comments

Suppose that S={-n,...,0,...,n} and that f(w,x,y,n) is a function, where w,x,y are in S. The number of ordered triples (w,x,y) satisfying f(w,x,y,n)=0, regarded as a function of n, is a sequence t of nonnegative integers. Sequences such as t/4 may also be integer sequences for all except certain initial values of n. In the following guide, such sequences are indicated in the related sequences column and may be included in the corresponding Mathematica programs.
...
sequence... f(w,x,y,n) ..... related sequences
A211415 ... w^2+x*y-1 ...... t+2, t/4, (t/4-1)/4
A211422 ... w^2+x*y ........ (t-1)/8, A120486
A211423 ... w^2+2x*y ....... (t-1)/4
A211424 ... w^2+3x*y ....... (t-1)/4
A211425 ... w^2+4x*y ....... (t-1)/4
A211426 ... 2w^2+x*y ....... (t-1)/4
A211427 ... 3w^2+x*y ....... (t-1)/4
A211428 ... 2w^2+3x*y ...... (t-1)/4
A211429 ... w^3+x*y ........ (t-1)/4
A211430 ... w^2+x+y ........ (t-1)/2
A211431 ... w^3+(x+y)^2 .... (t-1)/2
A211432 ... w^2-x^2-y^2 .... (t-1)/8
A003215 ... w+x+y .......... (t-1)/2, A045943
A202253 ... w+2x+3y ........ (t-1)/2, A143978
A211433 ... w+2x+4y ........ (t-1)/2
A211434 ... w+2x+5y ........ (t-1)/4
A211435 ... w+4x+5y ........ (t-1)/2
A211436 ... 2w+3x+4y ....... (t-1)/2
A211435 ... 2w+3x+5y ....... (t-1)/2
A211438 ... w+2x+2y ....... (t-1)/2, A118277
A001844 ... w+x+2y ......... (t-1)/4, A000217
A211439 ... w+3x+3y ........ (t-1)/2
A211440 ... 2w+3x+3y ....... (t-1)/2
A028896 ... w+x+y-1 ........ t/6, A000217
A211441 ... w+x+y-2 ........ t/3, A028387
A182074 ... w^2+x*y-n ...... t/4, A028387
A000384 ... w+x+y-n
A000217 ... w+x+y-2n
A211437 ... w*x*y-n ........ t/4, A007425
A211480 ... w+2x+3y-1
A211481 ... w+2x+3y-n
A211482 ... w*x+w*y+x*y-w*x*y
A211483 ... (n+w)^2-x-y
A182112 ... (n+w)^2-x-y-w
...
For the following sequences, S={1,...,n}, rather than
{-n,...,0,...n}. If f(w,x,y,n) is linear in w,x,y,n, then the sequence is a linear recurrence sequence.
A132188 ... w^2-x*y
A211506 ... w^2-x*y-n
A211507 ... w^2-x*y+n
A211508 ... w^2+x*y-n
A211509 ... w^2+x*y-2n
A211510 ... w^2-x*y+2n
A211511 ... w^2-2x*y ....... t/2
A211512 ... w^2-3x*y ....... t/2
A211513 ... 2w^2-x*y ....... t/2
A211514 ... 3w^2-x*y ....... t/2
A211515 ... w^3-x*y
A211516 ... w^2-x-y
A211517 ... w^3-(x+y)^2
A063468 ... w^2-x^2-y^2 .... t/2
A000217 ... w+x-y
A001399 ... w-2x-3y
A211519 ... w-2x+3y
A008810 ... w+2x-3y
A001399 ... w-2x-3y
A008642 ... w-2x-4y
A211520 ... w-2x+4y
A211521 ... w+2x-4y
A000115 ... w-2x-5y
A211522 ... w-2x+5y
A211523 ... w+2x-5y
A211524 ... w-3x-5y
A211533 ... w-3x+5y
A211523 ... w+3x-5y
A211535 ... w-4x-5y
A211536 ... w-4x+5y
A008812 ... w+4x-5y
A055998 ... w+x+y-2n
A074148 ... 2w+x+y-2n
A211538 ... 2w+2x+y-2n
A211539 ... 2w+2x-y-2n
A211540 ... 2w-3x-4y
A211541 ... 2w-3x+4y
A211542 ... 2w+3x-4y
A211543 ... 2w-3x-5y
A211544 ... 2w-3x+5y
A008812 ... 2w+3x-5y
A008805 ... w-2x-2y (repeated triangular numbers)
A001318 ... w-2x+2y
A000982 ... w+x-2y
A211534 ... w-3x-3y
A211546 ... w-3x+3y (triply repeated triangular numbers)
A211547 ... 2w-3x-3y (triply repeated squares)
A082667 ... 2w-3x+3y
A055998 ... w-x-y+2
A001399 ... w-2x-3y+1
A108579 ... w-2x-3y+n
...
Next, S={-n,...-1,1,...,n}, and the sequence counts the cases (w,x,y) satisfying the indicated inequality. If f(w,x,y,n) is linear in w,x,y,n, then the sequence is a linear recurrence sequence.
A211545 ... w+x+y>0; recurrence degree: 4
A211612 ... w+x+y>=0
A211613 ... w+x+y>1
A211614 ... w+x+y>2
A211615 ... |w+x+y|<=1
A211616 ... |w+x+y|<=2
A211617 ... 2w+x+y>0; recurrence degree: 5
A211618 ... 2w+x+y>1
A211619 ... 2w+x+y>2
A211620 ... |2w+x+y|<=1
A211621 ... w+2x+3y>0
A211622 ... w+2x+3y>1
A211623 ... |w+2x+3y|<=1
A211624 ... w+2x+2y>0; recurrence degree: 6
A211625 ... w+3x+3y>0; recurrence degree: 8
A211626 ... w+4x+4y>0; recurrence degree: 10
A211627 ... w+5x+5y>0; recurrence degree: 12
A211628 ... 3w+x+y>0; recurrence degree: 6
A211629 ... 4w+x+y>0; recurrence degree: 7
A211630 ... 5w+x+y>0; recurrence degree: 8
A211631 ... w^2>x^2+y^2; all terms divisible by 8
A211632 ... 2w^2>x^2+y^2; all terms divisible by 8
A211633 ... w^2>2x^2+2y^2; all terms divisible by 8
...
Next, S={1,...,n}, and the sequence counts the cases (w,x,y) satisfying the indicated relation.
A211634 ... w^2<=x^2+y^2
A211635 ... w^2A211790
A211636 ... w^2>=x^2+y^2
A211637 ... w^2>x^2+y^2
A211638 ... w^2+x^2+y^2
A211639 ... w^2+x^2+y^2<=n
A211640 ... w^2+x^2+y^2>n
A211641 ... w^2+x^2+y^2>=n
A211642 ... w^2+x^2+y^2<2n
A211643 ... w^2+x^2+y^2<=2n
A211644 ... w^2+x^2+y^2>2n
A211645 ... w^2+x^2+y^2>=2n
A211646 ... w^2+x^2+y^2<3n
A211647 ... w^2+x^2+y^2<=3n
A063691 ... w^2+x^2+y^2=n
A211649 ... w^2+x^2+y^2=2n
A211648 ... w^2+x^2+y^2=3n
A211650 ... w^3A211790
A211651 ... w^3>x^3+y^3; see Comments at A211790
A211652 ... w^4A211790
A211653 ... w^4>x^4+y^4; see Comments at A211790

Examples

			a(1) counts these 9 triples: (-1,-1,1), (-1, 1,-1), (0, -1, 0), (0, 0, -1), (0,0,0), (0,0,1), (0,1,0), (1,-1,1), (1,1,-1).
		

Crossrefs

Cf. A120486.

Programs

  • Mathematica
    t[n_] := t[n] = Flatten[Table[w^2 + x*y, {w, -n, n}, {x, -n, n}, {y, -n, n}]]
    c[n_] := Count[t[n], 0]
    t = Table[c[n], {n, 0, 70}] (* A211422 *)
    (t - 1)/8                   (* A120486 *)

A005891 Centered pentagonal numbers: (5n^2+5n+2)/2; crystal ball sequence for 3.3.3.4.4. planar net.

Original entry on oeis.org

1, 6, 16, 31, 51, 76, 106, 141, 181, 226, 276, 331, 391, 456, 526, 601, 681, 766, 856, 951, 1051, 1156, 1266, 1381, 1501, 1626, 1756, 1891, 2031, 2176, 2326, 2481, 2641, 2806, 2976, 3151, 3331, 3516, 3706, 3901, 4101, 4306, 4516, 4731, 4951, 5176, 5406
Offset: 0

Keywords

Comments

Equals the triangular numbers convolved with [1, 3, 1, 0, 0, 0, ...]. - Gary W. Adamson and Alexander R. Povolotsky, May 29 2009
From Ant King, Jun 15 2012: (Start)
a(n) == 1 (mod 5) for all n.
The digital roots of the a(n) form a purely periodic palindromic 9-cycle 1, 6, 7, 4, 6, 4, 7, 6, 1.
The units' digits of the a(n) form a purely periodic palindromic 4-cycle 1, 6, 6, 1.
(End)
Binomial transform of (1, 5, 5, 0, 0, 0, ...) and second partial sum of (1, 4, 5, 5, 5, ...). - Gary W. Adamson, Sep 09 2015
a(n) = a(-1-n) for all n in Z. - Michael Somos, Jan 25 2019
On the plane start with a single regular pentagon, and repeat the following procedure, "For each edge of any pentagon not already connected to an existing pentagon create a mirror image such that the mirror image does not overlap with an existing pentagon." a(n) is the number of pentagons occupying the plane after n repetitions. - Torlach Rush, Sep 14 2022

Examples

			a(2)= 5*T(2) + 1 = 5*3 + 1 = 16, a(4) = 5*T(4) + 1 = 5*10 + 1 = 51. - _Thomas M. Green_, Nov 16 2009
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • B. K. Teo and N. J. A. Sloane, Magic numbers in polygonal and polyhedral clusters, Inorgan. Chem. 24 (1985), 4545-4558.

Crossrefs

Cf. A028895, A001844, A003215, A004068 (partial sums), A006322, A001263.
Partial sums of A008706.
Equals second row of A167546 divided by 2.

Programs

Formula

G.f.: (1 + 3*x + x^2)/(1 - x)^3. Simon Plouffe in his 1992 dissertation
Narayana transform (A001263) of [1, 5, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0)=1, a(1)=6, a(2)=16. - Jaume Oliver Lafont, Dec 02 2008
a(n) = 5*A000217(n) + 1 = 5*T(n) + 1, for n = 0, 1, 2, 3, ... and where T(n) = n*(n+1)/2 = n-th triangular number. - Thomas M. Green, Nov 25 2009
a(n) = a(n-1) + 5*n, with a(0)=1. - Vincenzo Librandi, Nov 18 2010
a(n) = A028895(n) + 1. - Omar E. Pol, Oct 03 2011
a(n) = 2*a(n-1) - a(n-2) + 5. - Ant King, Jun 12 2012
Sum_{n>=0} 1/a(n) = 2*Pi /sqrt(15) *tanh(Pi/2*sqrt(3/5)) = 1.360613169863... - Ant King, Jun 15 2012
a(n) = A101321(5,n). - R. J. Mathar, Jul 28 2016
E.g.f.: (2 + 10*x + 5*x^2)*exp(x)/2. - Ilya Gutkovskiy, Jul 28 2016
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=0} a(n)/n! = 17*e/2.
Sum_{n>=0} (-1)^(n+1)*a(n)/n! = 3/(2*e). (End)

Extensions

Formula corrected and relocated by Johannes W. Meijer, Nov 07 2009

A056220 a(n) = 2*n^2 - 1.

Original entry on oeis.org

-1, 1, 7, 17, 31, 49, 71, 97, 127, 161, 199, 241, 287, 337, 391, 449, 511, 577, 647, 721, 799, 881, 967, 1057, 1151, 1249, 1351, 1457, 1567, 1681, 1799, 1921, 2047, 2177, 2311, 2449, 2591, 2737, 2887, 3041, 3199, 3361, 3527, 3697, 3871, 4049, 4231, 4417, 4607, 4801
Offset: 0

Author

N. J. A. Sloane, Aug 06 2000

Keywords

Comments

Image of squares (A000290) under "little Hankel" transform that sends [c_0, c_1, ...] to [d_0, d_1, ...] where d_n = c_n^2 - c_{n+1}*c_{n-1}. - Henry Bottomley, Dec 12 2000
Surround numbers of an n X n square. - Jason Earls, Apr 16 2001
Numbers n such that 2*n + 2 is a perfect square. - Cino Hilliard, Dec 18 2003, Juri-Stepan Gerasimov, Apr 09 2016
The sums of the consecutive integer sequences 2n^2 to 2(n+1)^2-1 are cubes, as 2n^2 + ... + 2(n+1)^2-1 = (1/2)(2(n+1)^2 - 1 - 2n^2 + 1)(2(n+1)^2 - 1 + 2n^2) = (2n+1)^3. E.g., 2+3+4+5+6+7 = 27 = 3^3, then 8+9+10+...+17 = 125 = 5^3. - Andras Erszegi (erszegi.andras(AT)chello.hu), Apr 29 2005
X values (other than 0) of solutions to the equation 2*X^3 + 2*X^2 = Y^2. To find Y values: b(n) = 2n*(2*n^2 - 1). - Mohamed Bouhamida, Nov 06 2007
Average of the squares of two consecutive terms is also a square. In fact: (2*n^2 - 1)^2 + (2*(n+1)^2 - 1)^2 = 2*(2*n^2 + 2*n + 1)^2. - Matias Saucedo (solomatias(AT)yahoo.com.ar), Aug 18 2008
Equals row sums of triangle A143593 and binomial transform of [1, 6, 4, 0, 0, 0, ...] with n > 1. - Gary W. Adamson, Aug 26 2008
Start a spiral of square tiles. Trivially the first tile fits in a 1 X 1 square. 7 tiles fit in a 3 X 3 square, 17 tiles fit in a 5 X 5 square and so on. - Juhani Heino, Dec 13 2009
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=-2, A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n >= 1, a(n) = coeff(charpoly(A,x),x^(n-2)). - Milan Janjic, Jan 26 2010
For each n > 0, the recursive series, formula S(b) = 6*S(b-1) - S(b-2) - 2*a(n) with S(0) = 4n^2-4n+1 and S(1) = 2n^2, has the property that every even term is a perfect square and every odd term is twice a perfect square. - Kenneth J Ramsey, Jul 18 2010
Fourth diagonal of A154685 for n > 2. - Vincenzo Librandi, Aug 07 2010
First integer of (2*n)^2 consecutive integers, where the last integer is 3 times the first + 1. As example, n = 2: term = 7; (2*n)^2 = 16; 7, 8, 9, ..., 20, 21, 22: 7*3 + 1 = 22. - Denis Borris, Nov 18 2012
Chebyshev polynomial of the first kind T(2,n). - Vincenzo Librandi, May 30 2014
For n > 0, number of possible positions of a 1 X 2 rectangle in a (n+1) X (n+2) rectangular integer lattice. - Andres Cicuttin, Apr 07 2016
This sequence also represents the best solution for Ripà's n_1 X n_2 X n_3 dots problem, for any 0 < n_1 = n_2 < n_3 = floor((3/2)*(n_1 - 1)) + 1. - Marco Ripà, Jul 23 2018

Examples

			a(0) = 0^2-1*1 = -1, a(1) = 1^2 - 4*0 = 1, a(2) = 2^2 - 9*1 = 7, etc.
a(4) = 31 = (1, 3, 3, 1) dot (1, 6, 4, 0) = (1 + 18 + 12 + 0). - _Gary W. Adamson_, Aug 29 2008
		

Crossrefs

Cf. A066049 (indices of prime terms)
Column 2 of array A188644 (starting at offset 1).

Programs

Formula

G.f.: (-1 + 4*x + x^2)/(1-x)^3. - Henry Bottomley, Dec 12 2000
a(n) = A119258(n+1,2) for n > 0. - Reinhard Zumkeller, May 11 2006
From Doug Bell, Mar 08 2009: (Start)
a(0) = -1,
a(n) = sqrt(A001844(n)^2 - A069074(n-1)),
a(n+1) = sqrt(A001844(n)^2 + A069074(n-1)) = sqrt(a(n)^2 + A069074(n-1)*2). (End)
a(n) + a(n+1) + 1 = (2n+1)^2. - Doug Bell, Mar 09 2009
a(n) = a(n-1) + 4*n - 2 (with a(0)=-1). - Vincenzo Librandi, Dec 25 2010
a(n) = A188653(2*n) for n > 0. - Reinhard Zumkeller, Apr 13 2011
a(n) = A162610(2*n-1,n) for n > 0. - Reinhard Zumkeller, Jan 19 2013
a(n) = ( Sum_{k=0..2} (C(n+k,3)-C(n+k-1,3))*(C(n+k,3)+C(n+k+1,3)) ) - (C(n+2,3)-C(n-1,3))*(C(n,3)+C(n+3,3)), for n > 3. - J. M. Bergot, Jun 16 2014
a(n) = j^2 + k^2 - 2 or 2*j*k if n >= 2 and j = n + sqrt(2)/2 and k = n - sqrt(2)/2. - Avi Friedlich, Mar 30 2015
a(n) = A002593(n)/n^2. - Bruce J. Nicholson, Apr 03 2017
a(n) = A000384(n) + n - 1. - Bruce J. Nicholson, Nov 12 2017
a(n)*a(n+k) + 2k^2 = m^2 (a perfect square), m = a(n) + (2n*k), for n>=1. - Ezhilarasu Velayutham, May 13 2019
From Amiram Eldar, Aug 10 2020: (Start)
Sum_{n>=1} 1/a(n) = 1/2 - sqrt(2)*Pi*cot(Pi/sqrt(2))/4.
Sum_{n>=1} (-1)^(n+1)/a(n) = sqrt(2)*Pi*csc(Pi/sqrt(2))/4 - 1/2. (End)
From Amiram Eldar, Feb 04 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = (Pi/sqrt(2))*csc(Pi/sqrt(2)).
Product_{n>=2} (1 - 1/a(n)) = (Pi/(4*sqrt(2)))*csc(Pi/sqrt(2)). (End)
a(n) = A003215(n) - A000217(n-2)*2. - Leo Tavares, Jun 29 2021
Let T(n) = n*(n+1)/2. Then a(n)^2 = T(2n-2)*T(2n+1) + n^2. - Charlie Marion, Feb 12 2023
E.g.f.: exp(x)*(2*x^2 + 2*x - 1). - Stefano Spezia, Jul 08 2023

A010057 a(n) = 1 if n is a cube, else 0.

Original entry on oeis.org

1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Author

N. J. A. Sloane, Mar 15 1996

Keywords

Comments

Multiplicative with a(p^e) = 1 if 3 divides e, 0 otherwise. - Mitch Harris, Jun 09 2005
a(A000578(n)) = 1; a(A007412(n)) = 0. - Reinhard Zumkeller, Oct 22 2011
a(n) = A000007(sum(A010872(A124010(n,k))): k = 1..A001221(n)) for n > 0. - Reinhard Zumkeller, Jun 21 2013
If n has 4 divisors, a(n) = bigomega(n) - 2. - Wesley Ivan Hurt, Jun 06 2014

References

  • E. Landau, Elementary Number Theory, translation by Jacob E. Goodman of Elementare Zahlentheorie (Vol. I_1 (1927) of Vorlesungen ueber Zahlentheorie), by Edmund Landau, with added exercises by Paul T. Bateman and E. E. Kohlbecker, Chelsea Publishing Co., New York, 1958, pp. 31-32.

Crossrefs

Cf. A000578.
Cf. A003215. - Reinhard Zumkeller, Sep 27 2008

Programs

  • Haskell
    a010057 0 = 1
    a010057 n = fromEnum $ all ((== 0) . (`mod` 3)) $ a124010_row n
    a010057_list = concatMap (\x -> 1 : replicate (a003215 x - 1) 0) [0..]
    -- Reinhard Zumkeller, Jun 21 2013, Oct 22 2011
    
  • Maple
    A010057 := proc(n)
        if n = 0 then
            1;
        else
            for pe in ifactors(n)[2] do
                if modp(op(2,pe),3) <> 0 then
                    return 0 ;
                end if;
            end do:
        end if;
        1 ;
    end proc: # R. J. Mathar, Feb 07 2023
  • Mathematica
    Table[ Boole[ IntegerQ[n^(1/3)]], {n, 0, 80}] (* Jean-François Alcover, Jun 10 2013 *)
  • PARI
    a(n) = ispower(n, 3); \\ Michel Marcus, Feb 24 2015
    
  • Python
    from sympy import integer_nthroot
    def A010057(n): return int(integer_nthroot(n,3)[1]) # Chai Wah Wu, Apr 02 2021

Formula

Dirichlet generating function: zeta(3s). - Franklin T. Adams-Watters, Sep 11 2005
a(n) = f(n,0) with f(x,y) = if x>0 then f(x-3*y*(y+1),y+1) else 0^(-x). - Reinhard Zumkeller, Sep 27 2008
a(n) = 1 + floor(n^(1/3)) - ceiling(n^(1/3)). - Wesley Ivan Hurt, Jun 06 2014
a(n) = floor(n^(1/3)) - floor((n-1)^(1/3)). - Mikael Aaltonen, Feb 24 2015

A003154 Centered 12-gonal numbers, or centered dodecagonal numbers: numbers of the form 6*k*(k-1) + 1.

Original entry on oeis.org

1, 13, 37, 73, 121, 181, 253, 337, 433, 541, 661, 793, 937, 1093, 1261, 1441, 1633, 1837, 2053, 2281, 2521, 2773, 3037, 3313, 3601, 3901, 4213, 4537, 4873, 5221, 5581, 5953, 6337, 6733, 7141, 7561, 7993, 8437, 8893, 9361, 9841, 10333, 10837, 11353, 11881, 12421
Offset: 1

Keywords

Comments

Binomial transform of [1, 12, 12, 0, 0, 0, ...]. Narayana transform (A001263) of [1, 12, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
Numbers k such that 6*k+3 is a square, these squares are given in A016946. - Gary Detlefs and Vincenzo Librandi, Aug 08 2010
Odd numbers of the form floor(n^2/6). - Juri-Stepan Gerasimov, Jul 27 2011
Bisection of A032528. - Omar E. Pol, Aug 20 2011
Sequence found by reading the line from 1, in the direction 1, 13, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. Opposite numbers to the members of A033581 in the same spiral. - Omar E. Pol, Sep 08 2011
The digital root has period 3 (1, 4, 1) (A146325), the same digital root as the centered triangular numbers A005448(n). - Peter M. Chema, Dec 20 2023

Examples

			From _Omar E. Pol_, Aug 21 2011: (Start)
1. Classic illustration of initial terms of the star numbers:
.
.                                     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            13                 37
.
2. Alternative illustration of initial terms using n-1 concentric hexagons around a central element:
.
.                                 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
(End)
		

References

  • Martin Gardner, Time Travel and Other Mathematical Bewilderments. Freeman, NY, 1988, p. 20.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Programs

Formula

G.f.: x*(1+10*x+x^2)/(1-x)^3. Simon Plouffe in his 1992 dissertation
a(n) = 1 + Sum_{j=0..n} (12*j). E.g., a(2)=37 because 1 + 12*0 + 12*1 + 12*2 = 37. - Xavier Acloque, Oct 06 2003
a(n) = numerator in B_2(x) = (1/2)x^2 - (1/2)x + 1/12 = Bernoulli polynomial of degree 2. - Gary W. Adamson, May 30 2005
a(n) = 12*(n-1) + a(n-1), with n>1, a(1)=1. - Vincenzo Librandi, Aug 08 2010
a(n) = A049598(n-1) + 1. - Omar E. Pol, Oct 03 2011
Sum_{n>=1} 1/a(n) = A306980 = Pi * tan(Pi/(2*sqrt(3))) / (2*sqrt(3)). - Vaclav Kotesovec, Jul 23 2019
From Amiram Eldar, Jun 21 2020: (Start)
Sum_{n>=1} a(n)/n! = 7*e - 1.
Sum_{n>=1} (-1)^n * a(n)/n! = 7/e - 1. (End)
a(n) = 2*A003215(n-1) - 1. - Leo Tavares, Jul 30 2021
E.g.f.: exp(x)*(1 + 6*x^2) - 1. - Stefano Spezia, Aug 19 2022

Extensions

More terms from Michael Somos
Previous Showing 31-40 of 286 results. Next