A069099 -id:A069099 - cp's OEIS frontend

A003215 Hex (or centered hexagonal) numbers: 3n(n+1)+1 (crystal ball sequence for hexagonal lattice).

1, 7, 19, 37, 61, 91, 127, 169, 217, 271, 331, 397, 469, 547, 631, 721, 817, 919, 1027, 1141, 1261, 1387, 1519, 1657, 1801, 1951, 2107, 2269, 2437, 2611, 2791, 2977, 3169, 3367, 3571, 3781, 3997, 4219, 4447, 4681, 4921, 5167, 5419, 5677, 5941, 6211, 6487, 6769

Offset: 0

N. J. A. Sloane

The hexagonal lattice is the familiar 2-dimensional lattice in which each point has 6 neighbors. This is sometimes called the triangular lattice.
Crystal ball sequence for A_2 lattice. - Michael Somos, Jun 03 2012
Sixth spoke of hexagonal spiral (cf. A056105-A056109).
Number of ordered integer triples (a,b,c), -n <= a,b,c <= n, such that a+b+c=0. - Benoit Cloitre, Jun 14 2003
Also the number of partitions of 6n into at most 3 parts, A001399(6n). - R. K. Guy, Oct 20 2003
Also, a(n) is the number of partitions of 6(n+1) into exactly 3 distinct parts. - William J. Keith, Jul 01 2004
Number of dots in a centered hexagonal figure with n+1 dots on each side.
Values of second Bessel polynomial y_2(n) (see A001498).
First differences of cubes (A000578). - Cecilia Rossiter (cecilia(AT)noticingnumbers.net), Dec 15 2004
Final digits of Hex numbers (hex(n) mod 10) are periodic with palindromic period of length 5 {1, 7, 9, 7, 1}. Last two digits of Hex numbers (hex(n) mod 100) are periodic with palindromic period of length 100. - Alexander Adamchuk, Aug 11 2006
All divisors of a(n) are congruent to 1, modulo 6. Proof: If p is an odd prime different from 3 then 3n^2 + 3n + 1 = 0 (mod p) implies 9(2n + 1)^2 = -3 (mod p), whence p = 1 (mod 6). - Nick Hobson, Nov 13 2006
For n>=1, a(n) is the side of Outer Napoleon Triangle whose reference triangle is a right triangle with legs (3a(n))^(1/2) and 3n(a(n))^(1/2). - Tom Schicker (tschicke(AT)email.smith.edu), Apr 25 2007
Number of triples (a,b,c) where 0<=(a,b)<=n and c=n (at least once the term n). E.g., for n = 1: (0,0,1), (0,1,0), (1,0,0), (0,1,1), (1,0,1), (1,1,0), (1,1,1), so a(1)=7. - Philippe Lallouet (philip.lallouet(AT)wanadoo.fr), Aug 20 2007
Equals the triangular numbers convolved with [1, 4, 1, 0, 0, 0, ...]. - Gary W. Adamson and Alexander R. Povolotsky, May 29 2009
From Terry Stickels, Dec 07 2009: (Start)
Also the maximum number of viewable cubes from any one static point while viewing a cube stack of identical cubes of varying magnitude.
For example, viewing a 2 X 2 X 2 stack will yield 7 maximum viewable cubes.
If the stack is 3 X 3 X 3, the maximum number of viewable cubes from any one static position is 19, and so on.
The number of cubes in the stack must always be the same number for width, length, height (at true regular cubic stack) and the maximum number of visible cubes can always be found by taking any cubic number and subtracting the number of the cube that is one less.
Examples: 125 - 64 = 61, 64 - 27 = 37, 27 - 8 = 19. (End)
The sequence of digital roots of the a(n) is period 3: repeat [1,7,1]. - Ant King, Jun 17 2012
The average of the first n (n>0) centered hexagonal numbers is the n-th square. - Philippe Deléham, Feb 04 2013
A002024 is the following array A read along antidiagonals:
  1,  2,  3,  4,  5,  6, ...
  2,  3,  4,  5,  6,  7, ...
  3,  4,  5,  6,  7,  8, ...
  4,  5,  6,  7,  8,  9, ...
  5,  6,  7,  8,  9, 10, ...
  6,  7,  8,  9, 10, 11, ...
and a(n) is the hook sum Sum_{k=0..n} A(n,k) + Sum_{r=0..n-1} A(r,n). - R. J. Mathar, Jun 30 2013
a(n) is the sum of the terms in the n+1 X n+1 matrices minus those in n X n matrices in an array formed by considering A158405 an array (the beginning terms in each row are 1,3,5,7,9,11,...). - J. M. Bergot, Jul 05 2013
The formula also equals the product of the three distinct combinations of two consecutive numbers: n^2, (n+1)^2, and n*(n+1). - J. M. Bergot, Mar 28 2014
The sides of any triangle ABC are divided into 2n + 1 equal segments by 2n points: A_1, A_2, ..., A_2n in side a, and also on the sides b and c cyclically. If A'B'C' is the triangle delimited by AA_n, BB_n and CC_n cevians, we have (ABC)/(A'B'C') = a(n) (see Java applet link). - Ignacio Larrosa Cañestro, Jan 02 2015
a(n) is the maximal number of parts into which (n+1) triangles can intersect one another. - Ivan N. Ianakiev, Feb 18 2015
((2^m-1)n)^t mod a(n) = ((2^m-1)(n+1))^t mod a(n) = ((2^m-1)(2n+1))^t mod a(n), where m any positive integer, and t = 0(mod 6). - Alzhekeyev Ascar M, Oct 07 2016
((2^m-1)n)^t mod a(n) = ((2^m-1)(n+1))^t mod a(n) = a(n) - (((2^m-1)(2n+1))^t mod a(n)), where m any positive integer, and t = 3(mod 6). - Alzhekeyev Ascar M, Oct 07 2016
(3n+1)^(a(n)-1) mod a(n) = (3n+2)^(a(n)-1) mod a(n) = 1. If a(n) not prime, then always strong pseudoprime. - Alzhekeyev Ascar M, Oct 07 2016
Every positive integer is the sum of 8 hex numbers (zero included), at most 3 of which are greater than 1. - Mauro Fiorentini, Jan 01 2018
Area enclosed by the segment of Archimedean spiral between n*Pi/2 and (n+1)*Pi/2 in Pi^3/48 units. - Carmine Suriano, Apr 10 2018
This sequence contains all numbers k such that 12*k - 3 is a square. - Klaus Purath, Oct 19 2021
The continued fraction expansion of sqrt(3*a(n)) is [3n+1; {1, 1, 2n, 1, 1, 6n+2}]. For n = 0, this collapses to [1; {1, 2}]. - Magus K. Chu, Sep 12 2022

			G.f. = 1 + 7*x + 19*x^2 + 37*x^3 + 61*x^4 + 91*x^5 + 127*x^6 + 169*x^7 + 217*x^8 + ...
From _Omar E. Pol_, Aug 21 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 o o o      o o o o o o
.                   o o o        o o o o o
.                                 o o o o
.
.   1      7          19             37
.
(End)
From _Klaus Purath_, Dec 03 2021: (Start)
(1) a(19) is not a prime number, because besides a(19) = a(9) + P(29), a(19) = a(15) + P(20) = a(2) + P(33) is also true.
(2) a(25) is prime, because except for a(25) = a(12) + P(38) there is no other equation of this pattern. (End)

John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 81.
M. Gardner, Time Travel and Other Mathematical Bewilderments. Freeman, NY, 1988, p. 18.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 0..1000
G. L. Alexanderson and John E. Wetzel, Dissections of a tetrahedron, J. Combinatorial Theory Ser. B 11 (1971), 58--66. MR0303412 (46 #2549). See p. 58.
B. T. Bennett and R. B. Potts, Arrays and brooks, J. Austral. Math. Soc., 7 (1967), 23-31 (see p. 30).
B. T. Bennett and R. B. Potts, Arrays and brooks, J. Austral. Math. Soc., 7 (1967), 23-31. [Annotated scanned copy]
Aran Bingham, Commutative n-ary Arithmetic, University of New Orleans Theses and Dissertations, Paper 1959, 2015.
Henry Bottomley, Illustration of initial terms.
John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 41.
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
Dominique Désérable, Versatile Topology for Two-Dimensional Cellular Automata, Advances in Cellular Automata, Emergence, Complexity and Computation (ECC Vol 52) Springer, Cham (2025), Ch. 6, pp. 151-186.
M. Gardner & N. J. A. Sloane, Correspondence, 1973-74
R. K. Guy, Letter to N. J. A. Sloane, 1987
R. K. Guy, The strong law of small numbers. Amer. Math. Monthly 95 (1988), no. 8, 697-712.
R. K. Guy, The strong law of small numbers. Amer. Math. Monthly 95 (1988), no. 8, 697-712. [Annotated scanned copy]
Md. Towhidul Islam, Extending triangle
G. S. Kazandzidis, On a Conjecture of Moessner and a General Problem, Bull. Soc. Math. Grèce, Nouvelle Série - vol. 2, fasc. 1-2, pp. 23-30 (1961).
Ignacio Larrosa Cañestro, Hexágono y estrella determinados por tres pares de cevianas simétricas, (java applet).
G. Nebe and N. J. A. Sloane, Home page for hexagonal (or triangular) lattice A2
Vladimir Pletser, Global Generalized Mersenne Numbers: Definition, Decomposition, and Generalized Theorems, Preprints.org, 2024. See p. 20.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
B. K. Teo and N. J. A. Sloane, Magic numbers in polygonal and polyhedral clusters, Inorgan. Chem. 24 (1985), 4545-4558.
Eric Weisstein's World of Mathematics, Hex Number
Eric Weisstein's World of Mathematics, Nexus Number
Eric Weisstein's World of Mathematics, Outer Napoleon Triangle.
Index entries for sequences related to centered polygonal numbers
Index entries for crystal ball sequences
Index entries for sequences related to A2 = hexagonal = triangular lattice
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A000124, A000166, A000217, A000290, A000578 (the cubes, or partial sums), A001263, A001498, A002061, A002378, A002407 (primes), A003514, A005408, A005449, A005891, A028896, A048766, A056105, A056106, A056107, A056108, A056109, A063496, A056220, A130298, A132111 (second diagonal), A158405, A215630, A239449, A243201.
Column k=3 of A080853, and column k=2 of A047969.
See also A220083 for a list of numbers of the form n*P(s,n)-(n-1)*P(s,n-1), where P(s,n) is the n-th polygonal number with s sides.
Cf. A287326(A000124(n), 1).
Cf. A008292.
Cf. A005448, A140091, A001844, A000384, A060544, A045943.
Cf. A154105.

a003215 n = 3 * n * (n + 1) + 1  -- Reinhard Zumkeller, Oct 22 2011

[3*n*(n+1)+1: n in [0..50]]; // G. C. Greubel, Nov 04 2017

A003215:=n->3*n*(n+1)+1; seq(A003215(n), n=0..100); # Wesley Ivan Hurt, Mar 28 2014

FoldList[#1 + #2 &, 1, 6 Range@ 50] (* Robert G. Wilson v, Feb 02 2011 *)
LinearRecurrence[{3, -3, 1}, {1, 7, 19}, 47] (* Robert G. Wilson v, Jul 06 2013 *)

makelist(3*n*(n+1)+1, n, 0, 30); /* Martin Ettl, Nov 12 2012 */

PARI
```
{a(n) = 3*n*(n+1) + 1};
```

[3*n*(n+1)+1 for n in range(47)] # Michael S. Branicky, Jan 07 2021

a(n) = 3*n*(n+1) + 1, n >= 0 (see the name).
a(n) = (n+1)^3 - n^3 = a(-1-n).
G.f.: (1 + 4*x + x^2) / (1 - x)^3. - Simon Plouffe in his 1992 dissertation
a(n) = 6*A000217(n) + 1.
a(n) = a(n-1) + 6*n = 2a(n-1) - a(n-2) + 6 = 3*a(n-1) - 3*a(n-2) + a(n-3) = A056105(n) + 5n = A056106(n) + 4*n = A056107(n) + 3*n = A056108(n) + 2*n = A056108(n) + n.
n-th partial arithmetic mean is n^2. - Amarnath Murthy, May 27 2003
a(n) = 1 + Sum_{j=0..n} (6*j). E.g., a(2)=19 because 1+ 6*0 + 6*1 + 6*2 = 19. - Xavier Acloque, Oct 06 2003
The sum of the first n hexagonal numbers is n^3. That is, Sum_{n>=1} (3*n*(n-1) + 1) = n^3. - Edward Weed (eweed(AT)gdrs.com), Oct 23 2003
a(n) = right term in M^n * [1 1 1], where M = the 3 X 3 matrix [1 0 0 / 2 1 0 / 3 3 1]. M^n * [1 1 1] = [1 2n+1 a(n)]. E.g., a(4) = 61, right term in M^4 * [1 1 1], since M^4 * [1 1 1] = [1 9 61] = [1 2n+1 a(4)]. - Gary W. Adamson, Dec 22 2004
Row sums of triangle A130298. - Gary W. Adamson, Jun 07 2007
a(n) = 3*n^2 + 3*n + 1. Proof: 1) If n occurs once, it may be in 3 positions; for the two other ones, n terms are independently possible, then we have 3*n^2 different triples. 2) If the term n occurs twice, the third one may be placed in 3 positions and have n possible values, then we have 3*n more different triples. 3) The term n may occurs 3 times in one way only that gives the formula. - Philippe Lallouet (philip.lallouet(AT)wanadoo.fr), Aug 20 2007
Binomial transform of [1, 6, 6, 0, 0, 0, ...]; Narayana transform (A001263) of [1, 6, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
a(n) = (n-1)*A000166(n) + (n-2)*A000166(n-1) = (n-1)floor(n!*e^(-1)+1) + (n-2)*floor((n-1)!*e^(-1)+1) (with offset 0). - Gary Detlefs, Dec 06 2009
a(n) = A028896(n) + 1. - Omar E. Pol, Oct 03 2011
a(n) = integral( (sin((n+1/2)x)/sin(x/2))^3, x=0..Pi)/Pi. - Yalcin Aktar, Dec 03 2011
Sum_{n>=0} 1/a(n) = Pi/sqrt(3)*tanh(Pi/(2*sqrt(3))) = 1.305284153013581... - Ant King, Jun 17 2012
a(n) = A000290(n) + A000217(2n+1). - Ivan N. Ianakiev, Sep 24 2013
a(n) = A002378(n+1) + A056220(n) = A005408(n) + 2*A005449(n) = 6*A000217(n) + 1. - Ivan N. Ianakiev, Sep 26 2013
a(n) = 6*A000124(n) - 5. - Ivan N. Ianakiev, Oct 13 2013
a(n) = A239426(n+1) / A239449(n+1) = A215630(2*n+1,n+1). - Reinhard Zumkeller, Mar 19 2014
a(n) = A243201(n) / A002061(n + 1). - Mathew Englander, Jun 03 2014
a(n) = A101321(6,n). - R. J. Mathar, Jul 28 2016
E.g.f.: (1 + 6*x + 3*x^2)*exp(x). - Ilya Gutkovskiy, Jul 28 2016
a(n) = (A001844(n) + A016754(n))/2. - Bruce J. Nicholson, Aug 06 2017
a(n) = A045943(2n+1). - Miquel Cerda, Jan 22 2018
a(n) = 3*Integral_{x=n..n+1} x^2 dx. - Carmine Suriano, Apr 10 2018
a(n) = A287326(A000124(n), 1). - Kolosov Petro, Oct 22 2018
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=0} a(n)/n! = 10*e.
Sum_{n>=0} (-1)^(n+1)*a(n)/n! = 2/e. (End)
G.f.: polylog(-3, x)*(1-x)/x. See the Simon Plouffe formula above, and the g.f. of the rows of A008292 by Vladeta Jovovic, Sep 02 2002. - Wolfdieter Lang, May 08 2021
a(n) = T(n-1)^2 - 2*T(n)^2 + T(n+1)^2, n >= 1, T = triangular number A000217. - Klaus Purath, Oct 11 2021
a(n) = 1 + 2*Sum_{j=n..2n} j. - Klaus Purath, Oct 19 2021
a(n) = A069099(n+1) - A000217(n). - Klaus Purath, Nov 03 2021
From Leo Tavares, Dec 03 2021: (Start)
a(n) = A005448(n) + A140091(n);
a(n) = A001844(n) + A002378(n);
a(n) = A005891(n) + A000217(n);
a(n) = A000290(n) + A000384(n+1);
a(n) = A060544(n-1) + 3*A000217(n);
a(n) = A060544(n-1) + A045943(n).
a(2*n+1) = A154105(n).
(End)

Partially edited by Joerg Arndt, Mar 11 2010

A000566 Heptagonal numbers (or 7-gonal numbers): n(5n-3)/2.

Original entry on oeis.org

0, 1, 7, 18, 34, 55, 81, 112, 148, 189, 235, 286, 342, 403, 469, 540, 616, 697, 783, 874, 970, 1071, 1177, 1288, 1404, 1525, 1651, 1782, 1918, 2059, 2205, 2356, 2512, 2673, 2839, 3010, 3186, 3367, 3553, 3744, 3940, 4141, 4347, 4558, 4774, 4995, 5221, 5452, 5688

Offset: 0

N. J. A. Sloane

Binomial transform of (0, 1, 5, 0, 0, 0, ...). Binomial transform is A084899. - Paul Barry, Jun 10 2003
Row sums of triangle A131413. - Gary W. Adamson, Jul 08 2007
Sequence starting (1, 7, 18, 34, ...) = binomial transform of (1, 6, 5, 0, 0, 0, ...). Also row sums of triangle A131896. - Gary W. Adamson, Jul 24 2007
Also the partial sums of A016861, a zero added in front; therefore a(n) = n (mod 5). - R. J. Mathar, Mar 19 2008
Also sequence found by reading the line from 0, in the direction 0, 7, ..., and the line from 1, in the direction 1, 18, ..., in the square spiral whose edges have length A195013 and whose vertices are the numbers A195014. These parallel lines are the semi-axes perpendicular to the main axis A195015 in the same spiral. - Omar E. Pol, Oct 14 2011
Also sequence found by reading the line from 0, in the direction 0, 7, ... and the parallel line from 1, in the direction 1, 18, ..., in the square spiral whose vertices are the generalized heptagonal numbers A085787. - Omar E. Pol, Jul 18 2012
Partial sums give A002413. - Omar E. Pol, Jan 12 2013
The n-th heptagonal number equals the sum of the n consecutive integers starting at 2*n-1; for example, 1, 3+4, 5+6+7, 7+8+9+10, etc. In general, the n-th (2k+1)-gonal number is the sum of the n consecutive integers starting at (k-1)*n - (k-2). When k = 1 and 2, this result generates the triangular numbers, A000217, and the pentagonal numbers, A000326, respectively. - Charlie Marion, Mar 02 2022

			G.f. = x + 7*x^2 + 18*x^3 + 34*x^4 + 55*x^5 + 81*x^6 + 112*x^7 + ... - _Michael Somos_, Jan 25 2019

Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 38.
E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6.
Leonard 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).
David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 123.

Daniel Mondot, Table of n, a(n) for n = 0..10000 (first 1000 terms by T. D. Noe)
S. Barbero, U. Cerruti, and N. Murru, Transforming Recurrent Sequences by Using the Binomial and Invert Operators, J. Int. Seq. 13 (2010) # 10.7.7., section 4.4.
C. K. Cook and M. R. Bacon, Some polygonal number summation formulas, Fib. Q., 52 (2014), 336-343.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 341.
Bir Kafle, Florian Luca, and Alain Togbé, Pentagonal and heptagonal repdigits, Annales Mathematicae et Informaticae, pp. 137-145.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992.
Omar E. Pol, Illustration of initial terms of A000217, A000290, A000326, A000384, A000566, A000567.
B. Srinivasa Rao, Heptagonal Numbers in the Pell Sequence and Diophantine Equations 2x^2 = y^2(5y - 3)^2 ± 2, Fib. Quarterly, 43 (2005), 194-201.
B. Srinivasa Rao, Heptagonal numbers in the associated Pell sequence and Diophantine equations x^2(5x - 3)^2 = 8y^2 ± 4, Fib. Quarterly, 43 (2005), 302-306.
Leo Tavares, Illustration.
Eric Weisstein's World of Mathematics, Heptagonal Number.
Index to sequences related to polygonal numbers.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A014637, A014640, A014773, A014792, A069099, A131413, A131896, A134483, A000384.
a(n)= A093562(n+1, 2), (5, 1)-Pascal column.
Cf. A006564, A147875, A244639.
Cf. sequences listed in A254963.

a000566 n = n * (5 * (n - 1) + 2) `div` 2
a000566_list = scanl (+) 0 a016861_list  -- Reinhard Zumkeller, Jun 16 2013

a000566:=func< n | n*(5*n-3) div 2 >; [ a000566(n): n in [0..50] ];

A000566 := proc(n)
        n*(5*n-3)/2 ;
end proc:
seq(A000566(n),n=0..30); # R. J. Mathar, Oct 02 2020

Table[n (5n - 3)/2, {n, 0, 50}] (* or *) LinearRecurrence[{3, -3, 1}, {0, 1, 7}, 50] (* Harvey P. Dale, Oct 13 2011 *)
Join[{0},Accumulate[Range[1,315,5]]] (* Harvey P. Dale, Mar 26 2016 *)
(* For Mathematica 10.4+ *) Table[PolygonalNumber[RegularPolygon[7], n], {n, 0, 48}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
PolygonalNumber[7,Range[0,50]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Jan 23 2021 *)

makelist(n*(5*n-3)/2, n, 0, 20); /* Martin Ettl, Dec 11 2012 */

PARI
```
a(n) = n * (5*n - 3) / 2
```

# Intended to compute the initial segment of the sequence, not isolated terms.
def aList():
     x, y = 1, 1
     yield 0
     while True:
         yield x
         x, y = x + y + 5, y + 5
A000566 = aList()
print([next(A000566) for i in range(49)]) # Peter Luschny, Aug 04 2019

[n*(5*n-3)//2 for n in range(50)] # Gennady Eremin, Mar 24 2022

G.f.: x*(1 + 4*x)/(1 - x)^3. Simon Plouffe in his 1992 dissertation.
a(n) = C(n, 1) + 5*C(n, 2). - Paul Barry, Jun 10 2003
a(n) = Sum_{k = 1..n} (4*n - 3*k). - Paul Barry, Sep 06 2005
a(n) = n + 5*A000217(n-1) - Floor van Lamoen, Oct 14 2005
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for a(0) = 0, a(1) = 1, a(2) = 7. -  Jaume Oliver Lafont, Dec 02 2008
a(n+1) = A153126(n) + n mod 2; a(2*n+1) = A033571(n); a(2*(n+1)) = A153127(n) + 1. - Reinhard Zumkeller, Dec 20 2008
40*a(n)+ 9 = A017354(n-1). - Ken Rosenbaum, Dec 02 2009.
a(n) = 2*a(n-1) - a(n-2) + 5, with a(0) = 0 and a(1) = 1. - Mohamed Bouhamida, May 05 2010
a(n) = A000217(n) + 4*A000217(n-1). - Vincenzo Librandi, Nov 20 2010
a(n) = a(n-1) + 5*n - 4, with a(0) = 0. - Vincenzo Librandi, Nov 20 2010
a(n) = A130520(5*n). - Philippe Deléham, Mar 26 2013
a(5*a(n) + 11*n + 1) = a(5*a(n) + 11*n) + a(5*n + 1). - Vladimir Shevelev, Jan 24 2014
Sum_{n>=1} 1/a(n) = sqrt(1 - 2/sqrt(5))*Pi/3 + 5*log(5)/6 - sqrt(5)*log((1 + sqrt(5))/2)/3 = 1.32277925312238885674944226131... . See A244639. - Vaclav Kotesovec, Apr 27 2016
E.g.f.: x*(2 + 5*x)*exp(x)/2. - Ilya Gutkovskiy, Aug 27 2016
From Charlie Marion, May 02 2017: (Start)
a(n+m) = a(n) + 5*n*m + a(m);
a(n-m) = a(n) - 5*n*m + a(m) + 3*m;
a(n) - a(m) = (5*(n + m) - 3)*(n - m)/2.
In general, let P(k,n) be the n-th k-gonal number. Then
P(k,n+m) = P(k,n) + (k - 2)*n*m + P(k,m);
P(k,n-m) = P(k,n) - (k - 2)*n*m + P(k,m) + (k - 4)*m;
P(k,n) - P(k,m) = ((k-2)*(n + m) + 4 - k)*(n - m)/2.
(End)
a(n) = A147875(-n) for all n in Z. - Michael Somos, Jan 25 2019
a(n) = A000217(n-1) + A000217(2*n-1). - Charlie Marion, Dec 19 2019
Product_{n>=2} (1 - 1/a(n)) = 5/7. - Amiram Eldar, Jan 21 2021
a(n) + a(n+1) = (2*n+1)^2 + n^2 - 2*n. In general, if we let P(k,n) = the n-th k-gonal number, then P(k^2-k+1,n)+ P(k^2-k+1,n+1) = ((k-1)*n+1)^2 + (k-2)*(n^2-2*n) = ((k-1)*n+1)^2 + (k-2)*A005563(n-2). When k = 2, this formula reduces to the well-known triangular number formula: T(n) + T(n+1) = (n+1)^2. - Charlie Marion, Jul 01 2021

Partially edited by Joerg Arndt, Mar 11 2010

A001106 9-gonal (or enneagonal or nonagonal) numbers: a(n) = n(7n-5)/2.

Original entry on oeis.org

0, 1, 9, 24, 46, 75, 111, 154, 204, 261, 325, 396, 474, 559, 651, 750, 856, 969, 1089, 1216, 1350, 1491, 1639, 1794, 1956, 2125, 2301, 2484, 2674, 2871, 3075, 3286, 3504, 3729, 3961, 4200, 4446, 4699, 4959, 5226, 5500, 5781, 6069, 6364

Offset: 0

N. J. A. Sloane

Sequence found by reading the line from 0, in the direction 0, 9, ... and the parallel line from 1, in the direction 1, 24, ..., in the square spiral whose vertices are the generalized 9-gonal (enneagonal) numbers A118277. Also sequence found by reading the same lines in the square spiral whose edges have length A195019 and whose vertices are the numbers A195020. - Omar E. Pol, Sep 10 2011
Number of ordered pairs of integers (x,y) with abs(x) < n, abs(y) < n and x+y <= n. - Reinhard Zumkeller, Jan 23 2012
Partial sums give A007584. - Omar E. Pol, Jan 15 2013

Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe and William A. Tedeschi, Table of n, a(n) for n = 0..10000 (1000 terms were computed by T. D. Noe)
S. Barbero, U. Cerruti and N. Murru, Transforming Recurrent Sequences by Using the Binomial and Invert Operators, J. Int. Seq. 13 (2010) # 10.7.7, section 4.4.
C. K. Cook and M. R. Bacon, Some polygonal number summation formulas, Fib. Q., 52 (2014), 336-343.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 343
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Eric Weisstein's World of Mathematics, Nonagonal Number.
Index to sequences related to polygonal numbers
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A093564 ((7, 1) Pascal, column m=2). Partial sums of A016993.

Cf. A131875, A057655, A069099, A244646.

a001106 n = length [(x,y) | x <- [-n+1..n-1], y <- [-n+1..n-1], x + y <= n]
-- Reinhard Zumkeller, Jan 23 2012

a001106 n = n*(7*n-5) `div` 2 -- James Spahlinger, Oct 18 2012

Table[n(7n - 5)/2, {n, 0, 50}] (* or *) LinearRecurrence[{3, -3, 1}, {0, 1, 9}, 50] (* Harvey P. Dale, Nov 06 2011 *)
(* For Mathematica 10.4+ *) Table[PolygonalNumber[RegularPolygon[9], n], {n, 0, 43}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
PolygonalNumber[9,Range[0,50]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Nov 19 2019 *)

a(n)=n*(7*n-5)/2 \\ Charles R Greathouse IV, Jun 10 2011

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 + 7, y + 7
A001106 = aList()
print([next(A001106) for i in range(49)]) # Peter Luschny, Aug 04 2019

a(n) = (7*n - 5)*n/2.
G.f.: x*(1+6*x)/(1-x)^3. - Simon Plouffe in his 1992 dissertation.
a(n) = n + 7*A000217(n-1). - Floor van Lamoen, Oct 14 2005
Starting (1, 9, 24, 46, 75, ...) gives the binomial transform of (1, 8, 7, 0, 0, 0, ...). - Gary W. Adamson, Jul 22 2007
Row sums of triangle A131875 starting (1, 9, 24, 46, 75, 111, ...). A001106 = binomial transform of (1, 8, 7, 0, 0, 0, ...). - Gary W. Adamson, Jul 22 2007
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0) = 0, a(1) = 1, a(2) = 9. - Jaume Oliver Lafont, Dec 02 2008
a(n) = 2*a(n-1) - a(n-2) + 7. - Mohamed Bouhamida, May 05 2010
a(n) = a(n-1) + 7*n - 6 (with a(0) = 0). - Vincenzo Librandi, Nov 12 2010
a(n) = A174738(7n). - Philippe Deléham, Mar 26 2013
a(7*a(n) + 22*n + 1) =  a(7*a(n) + 22*n) + a(7*n+1). - Vladimir Shevelev, Jan 24 2014
E.g.f.: x*(2 + 7*x)*exp(x)/2. - Ilya Gutkovskiy, Jul 28 2016
a(n+2) + A000217(n) = (2*n+3)^2. - Ezhilarasu Velayutham, Mar 18 2020
Product_{n>=2} (1 - 1/a(n)) = 7/9. - Amiram Eldar, Jan 21 2021
Sum_{n>=1} 1/a(n) = A244646. - Amiram Eldar, Nov 12 2021
a(n) = A000217(3*n-2) - (n-1)^2. - Charlie Marion, Feb 27 2022
a(n) = 3*A000217(n) + 2*A005563(n-2). In general, if P(k,n) = the n-th k-gonal number, then P(m*k,n) = m*P(k,n) + (m-1)*A005563(n-2). - Charlie Marion, Feb 21 2023

A062786 Centered 10-gonal numbers.

Original entry on oeis.org

1, 11, 31, 61, 101, 151, 211, 281, 361, 451, 551, 661, 781, 911, 1051, 1201, 1361, 1531, 1711, 1901, 2101, 2311, 2531, 2761, 3001, 3251, 3511, 3781, 4061, 4351, 4651, 4961, 5281, 5611, 5951, 6301, 6661, 7031, 7411, 7801, 8201, 8611, 9031, 9461, 9901, 10351, 10811

Offset: 1

Jason Earls, Jul 19 2001

Deleting the least significant digit yields the (n-1)-st triangular number: a(n) = 10*A000217(n-1) + 1. - Amarnath Murthy, Dec 11 2003
All divisors of a(n) are congruent to 1 or -1, modulo 10; that is, they end in the decimal digit 1 or 9. Proof: If p is an odd prime different from 5 then 5n^2 - 5n + 1 == 0 (mod p) implies 25(2n - 1)^2 == 5 (mod p), whence p == 1 or -1 (mod 10). - Nick Hobson, Nov 13 2006
Centered decagonal numbers. - Omar E. Pol, Oct 03 2011
The partial sums of this sequence give A004466. - Leo Tavares, Oct 04 2021
The continued fraction expansion of sqrt(5*a(n)) is [5n-3; {2, 2n-2, 2, 10n-6}]. For n=1, this collapses to [2; {4}]. - Magus K. Chu, Sep 12 2022
Numbers m such that 20*m + 5 is a square. Also values of the Fibonacci polynomial y^2 - x*y - x^2 for x = n and y = 3*n - 1. This is a subsequence of A089270. - Klaus Purath, Oct 30 2022
All terms can be written as a difference of two consecutive squares a(n) = A005891(n-1)^2 - A028895(n-1)^2, and they can be represented by the forms (x^2 + 2mxy + (m^2-1)y^2) and (3x^2 + (6m-2)xy + (3m^2-2m)y^2), both of discriminant 4. - Klaus Purath, Oct 17 2023

T. D. Noe, Table of n, a(n) for n = 1..1000
Leo Tavares, Illustration: Pentagonal Stars.
Leo Tavares, Illustration: Mid-section Stars.
Leo Tavares, Illustration: Mid-section Star Pillars.
Leo Tavares, Illustration: Trapezoidal Rays.
R. Yin, J. Mu, and T. Komatsu, The p-Frobenius Number for the Triple of the Generalized Star Numbers, Preprints 2024, 2024072280. See p. 2.
Index entries for sequences related to centered polygonal numbers.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001263, A124080, A101321, A028387, A016861, A003154, A005891, A000217, A004466, A144390, A000326, A060544.

Cf. also A016754, A002378, A069099, A045943, A003215, A046092, A001844, A028896, A005448, A024966, A082970.

List([1..50], n-> 1+5*n*(n-1)); # G. C. Greubel, Mar 30 2019

[1+5*n*(n-1): n in [1..50]]; // G. C. Greubel, Mar 30 2019

FoldList[#1+#2 &, 1, 10Range@ 45] (* Robert G. Wilson v, Feb 02 2011 *)
1+5*Pochhammer[Range[50]-1, 2] (* G. C. Greubel, Mar 30 2019 *)

j=[]; for(n=1,75,j=concat(j,(5*n*(n-1)+1))); j

for (n=1, 1000, write("b062786.txt", n, " ", 5*n*(n - 1) + 1) ) \\ Harry J. Smith, Aug 11 2009

def a(n): return(5*n**2-5*n+1) # Torlach Rush, May 10 2024

[1+5*rising_factorial(n-1, 2) for n in (1..50)] # G. C. Greubel, Mar 30 2019

a(n) = 5*n*(n-1) + 1.
From  Gary W. Adamson, Dec 29 2007: (Start)
Binomial transform of [1, 10, 10, 0, 0, 0, ...];
Narayana transform (A001263) of [1, 10, 0, 0, 0, ...]. (End)
G.f.: x*(1+8*x+x^2) / (1-x)^3. - R. J. Mathar, Feb 04 2011
a(n) = A124080(n-1) + 1. - Omar E. Pol, Oct 03 2011
a(n) = A101321(10,n-1). - R. J. Mathar, Jul 28 2016
a(n) = A028387(A016861(n-1))/5 for n > 0. - Art Baker, Mar 28 2019
E.g.f.: (1+5*x^2)*exp(x) - 1. - G. C. Greubel, Mar 30 2019
Sum_{n>=1} 1/a(n) = Pi * tan(Pi/(2*sqrt(5))) / sqrt(5). - Vaclav Kotesovec, Jul 23 2019
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=1} a(n)/n! = 6*e - 1.
Sum_{n>=1} (-1)^n * a(n)/n! = 6/e - 1. (End)
a(n) = A005891(n-1) + 5*A000217(n-1). - Leo Tavares, Jul 14 2021
a(n) = A003154(n) - 2*A000217(n-1). See Mid-section Stars illustration. - Leo Tavares, Sep 06 2021
From Leo Tavares, Oct 06 2021: (Start)
a(n) = A144390(n-1) + 2*A028387(n-1). See Mid-section Star Pillars illustration.
a(n) = A000326(n) + A000217(n) + 3*A000217(n-1). See Trapezoidal Rays illustration.
a(n) = A060544(n) + A000217(n-1). (End)
From Leo Tavares, Oct 31 2021: (Start)
a(n) = A016754(n-1) + 2*A000217(n-1).
a(n) = A016754(n-1) + A002378(n-1).
a(n) = A069099(n) + 3*A000217(n-1).
a(n) = A069099(n) + A045943(n-1).
a(n) = A003215(n-1) + 4*A000217(n-1).
a(n) = A003215(n-1) + A046092(n-1).
a(n) = A001844(n-1) + 6*A000217(n-1).
a(n) = A001844(n-1) + A028896(n-1).
a(n) = A005448(n) + 7*A000217(n).
a(n) = A005448(n) + A024966(n). (End)
From Klaus Purath, Oct 30 2022: (Start)
a(n) = a(n-2) + 10*(2*n-3).
a(n) = 2*a(n-1) - a(n-2) + 10.
a(n) = A135705(n-1) + n.
a(n) = A190816(n) - n.
a(n) = 2*A005891(n-1) - 1. (End)

Better description from Terrel Trotter, Jr., Apr 06 2002

A069129 Centered 16-gonal numbers.

Original entry on oeis.org

1, 17, 49, 97, 161, 241, 337, 449, 577, 721, 881, 1057, 1249, 1457, 1681, 1921, 2177, 2449, 2737, 3041, 3361, 3697, 4049, 4417, 4801, 5201, 5617, 6049, 6497, 6961, 7441, 7937, 8449, 8977, 9521, 10081, 10657, 11249, 11857, 12481, 13121, 13777, 14449, 15137, 15841

Offset: 1

Terrel Trotter, Jr., Apr 07 2002

Also, sequence found by reading the line from 1, in the direction 1, 17, ..., in the square spiral whose vertices are the triangular numbers A000217. Opposite numbers to the members of A139098 in the same spiral. - Omar E. Pol, Apr 26 2008
The subsequence of primes begins: 17, 97, 241, 337, 449, 577, 881, 1249, 3041, 3361, 3697, 4049, 4801, 6961, 7937, 9521,  10657, 13121, 14449. See A184899: n such that the n-th centered 12-gonal number is prime. Indices of prime star numbers. - Jonathan Vos Post, Feb 27 2011
Binomial transform of [1, 16, 16, 0, 0, 0, ...] and Narayana transform (A001263) of [1, 16, 0, 0, 0, ...]. - Gary W. Adamson, Jul 28 2011
Centered hexadecagonal numbers or centered hexakaidecagonal numbers. - Omar E. Pol, Oct 03 2011
a(n) = m(n,n) for an array constructed by using the terms in A016813 as the antidiagonals; the first few antidiagonals are 1; 5,9; 13,17,21; 25,29,33,37. - J. M. Bergot, Jul 05 2013
[The first five rows begin: 1,9,21,37,57; 5,17,33,53,77; 13,29,49,73,101; 25,45,69,97,129; 41,65,93,125,161.]

			a(5) = 161 because 8*5^2 - 8*5 + 1 = 200 - 40 + 1 = 161.

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Omar E. Pol, Determinacion geometrica de los numeros primos y perfectos.
Leo Tavares, Illustration: Octagonal Stars.
Eric Weisstein's World of Mathematics, Centered Polygonal Numbers.
R. Yin, J. Mu, and T. Komatsu, The p-Frobenius Number for the Triple of the Generalized Star Numbers, Preprints 2024, 2024072280. See p. 2.
Index entries for sequences related to centered polygonal numbers.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A005448, A001844, A005891, A003215, A069099, A000217, A035008, A139098.

Bisection of A077221.

[8*n^2-8*n+1: n in [0..50]]; // Vincenzo Librandi, Feb 05 2013

FoldList[#1 + #2 &, 1, 16 Range@ 45] (* Robert G. Wilson v, Feb 02 2011 *)
Rest[CoefficientList[Series[-x(1+14x+x^2)/(x-1)^3,{x,0,50}],x]]  (* Harvey P. Dale, Apr 22 2011 *)

a(n)=8*n^2-8*n+1 \\ Charles R Greathouse IV, Sep 24 2015

a(n) = 8*n^2 - 8*n + 1.
a(n) = A035008(n-1) + 1. - Omar E. Pol, Apr 26 2008
a(n) = 16*n + a(n-1) - 16 with n > 1, a(1)=1. - Vincenzo Librandi, Aug 08 2010
G.f.: -x*(1+14*x+x^2) / (x-1)^3. - R. J. Mathar, Feb 04 2011
E.g.f.: (8*x^2 + 1)*exp(x). - G. C. Greubel, Jul 18 2017
a(n) = A056220(2n-1). - Bruce J. Nicholson, Aug 31 2017
Sum_{n>=1} 1/a(n) = Pi * tan(Pi/(2*sqrt(2))) / (4*sqrt(2)). - Vaclav Kotesovec, Jul 23 2019
From Amiram Eldar, Jun 21 2020: (Start)
Sum_{n>=1} a(n)/n! = 9*e - 1.
Sum_{n>=1} (-1)^n * a(n)/n! = 9/e - 1. (End)
Product_{n>=2} (a(n) - 1) / (a(n) + 1) = Pi/4. - Dimitris Valianatos, Jun 27 2020
a(n) = A016754(n-1) + 8*A000217(n-1). - Leo Tavares, Jul 19 2021

A101321 Table T(n,m) = 1 + nm(m+1)/2 read by antidiagonals: centered polygonal numbers.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 7, 7, 4, 1, 1, 11, 13, 10, 5, 1, 1, 16, 21, 19, 13, 6, 1, 1, 22, 31, 31, 25, 16, 7, 1, 1, 29, 43, 46, 41, 31, 19, 8, 1, 1, 37, 57, 64, 61, 51, 37, 22, 9, 1, 1, 46, 73, 85, 85, 76, 61, 43, 25, 10, 1, 1, 56, 91, 109, 113, 106, 91, 71, 49, 28, 11, 1, 1, 67

Offset: 0

Eugene McDonnell (eemcd(AT)mac.com), Dec 24 2004

Row n gives the centered figurate numbers of the n-gon.

Antidiagonal sums are in A101338.

			The upper left corner of the infinite array T is
|0| 1   1   1   1   1   1   1   1   1   1 ... A000012
|1| 1   2   4   7  11  16  22  29  37  46 ... A000124
|2| 1   3   7  13  21  31  43  57  73  91 ... A002061
|3| 1   4  10  19  31  46  64  85 109 136 ... A005448
|4| 1   5  13  25  41  61  85 113 145 181 ... A001844
|5| 1   6  16  31  51  76 106 141 181 226 ... A005891
|6| 1   7  19  37  61  91 127 169 217 271 ... A003215
|7| 1   8  22  43  71 106 148 197 253 316 ... A069099
|8| 1   9  25  49  81 121 169 225 289 361 ... A016754
|9| 1  10  28  55  91 136 190 253 325 406 ... A060544

Cf. A000217, A001263, A101338.

A101321 := proc(n,k)
        n*k*(k+1)/2+1 ;
end proc: # R. J. Mathar, Jul 28 2016

T[n_, m_] := 1 + n m (m + 1)/2;
Table[T[n - m, m], {n, 0, 12}, {m, n, 0, -1}] // Flatten (* Jean-François Alcover, Mar 23 2020 *)

T(n,m) = 1 + n*m*(m+1)/2 \\ Charles R Greathouse IV, Jul 28 2016

T(n,2) = A016777(n). T(n,3) = A016921(n). T(n,4) = A017281(n).
T(10,m) = A062786(m+1).
T(11,m) = A069125(m+1).
T(12,m) = A003154(m+1).
T(13,m) = A069126(m+1).
T(14,m) = A069127(m+1).
T(15,m) = A069128(m+1).
T(16,m) = A069129(m+1).
T(17,m) = A069130(m+1).
T(18,m) = A069131(m+1).
T(19,m) = A069132(m+1).
T(20,m) = A069133(m+1).
T(n+1,m) = T(n,m) + m*(m+1)/2. - Gary W. Adamson and Michel Marcus, Oct 13 2015

Edited by R. J. Mathar, Oct 21 2009

A004126 a(n) = n(7n^2 - 1)/6.

Original entry on oeis.org

0, 1, 9, 31, 74, 145, 251, 399, 596, 849, 1165, 1551, 2014, 2561, 3199, 3935, 4776, 5729, 6801, 7999, 9330, 10801, 12419, 14191, 16124, 18225, 20501, 22959, 25606, 28449, 31495, 34751, 38224, 41921, 45849, 50015, 54426, 59089, 64011

Offset: 0

Albert D. Rich (Albert_Rich(AT)msn.com)

3-dimensional analog of centered polygonal numbers.
Sum of n triangular numbers starting from T(n), where T = A000217. E.g., a(4) = T(4) + T(5) + T(6) + T(7) = 10 + 15 + 21 + 28 = 74. - Amarnath Murthy, Jul 16 2004
Also as a(n) = (1/6)*(7*n^3-n), n>0: structured heptagonal diamond numbers (vertex structure 8). Cf. A100179 = alternate vertex; A000447 = structured diamonds; A100145 for more on structured numbers. - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Partial sums of A069099, centered heptagonal numbers (A000566). - Jonathan Vos Post, Mar 16 2006
Binomial transform of (0, 1, 7, 7, 0, 0, 0, ...) and third partial sum of (0, 1, 6, 7, 7, 7, ...). - Gary W. Adamson, Oct 05 2015

E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 140.

Vincenzo Librandi, Table of n, a(n) for n = 0..5000
T. P. Martin, Shells of atoms, Phys. Reports, 273 (1996), 199-241, eq. (11).
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

1/12*t*(n^3-n)+n for t = 2, 4, 6, ... gives A004006, A006527, A006003, A005900, A004068, A000578, A004126, A000447, A004188, A004466, A004467, A007588, A062025, A063521, A063522, A063523.
Cf. A000217, A000566, A016993, A069099.
Cf. A000447, A000292.

[n*(7*n^2-1)/6: n in [0..50]]; // Vincenzo Librandi, May 15 2011

seq(binomial(2*n+1,3)-binomial(n+1,3), n=0..38); # Zerinvary Lajos, Jan 21 2007

Table[n (7 n^2 - 1)/6, {n, 0, 80}] (* Vladimir Joseph Stephan Orlovsky, Apr 18 2011 *)

makelist(n*(7*n^2-1)/6,n,0,30); /* Martin Ettl, Jan 08 2013 */

vector(100, n, n--; n*(7*n^2 - 1)/6) \\ Altug Alkan, Oct 06 2015

a(n) = C(2*n+1,3)-C(n+1,3), n>=0. - Zerinvary Lajos, Jan 21 2007
a(n) = A000447(n) - A000292(n). - Zerinvary Lajos, Jan 21 2007
G.f.: x*(1+5*x+x^2)/(1-x)^4. - Colin Barker, Mar 02 2012
E.g.f.: (x/6)*(7*x^2 + 21*x + 6)*exp(x). - G. C. Greubel, Oct 05 2015
a(n) = Sum_{i = n..2*n-1} A000217(i). - Bruno Berselli, Sep 06 2017
a(n) = n^3 + Sum_{k=0..n-1} k*(k+1)/2. Alternately, a(n) = A000578(n) + A000292(n-1) for n>0. - Bruno Berselli, May 23 2018

A024966 7 times triangular numbers: 7n(n+1)/2.

Original entry on oeis.org

0, 7, 21, 42, 70, 105, 147, 196, 252, 315, 385, 462, 546, 637, 735, 840, 952, 1071, 1197, 1330, 1470, 1617, 1771, 1932, 2100, 2275, 2457, 2646, 2842, 3045, 3255, 3472, 3696, 3927, 4165, 4410, 4662, 4921, 5187, 5460, 5740, 6027, 6321, 6622

Offset: 0

Joe Keane (jgk(AT)jgk.org), Dec 11 1999

Sequence found by reading the line from 0, in the direction 0, 7, ... and the same line from 0, in the direction 1, 21, ..., in the square spiral whose edges have length A195019 and whose vertices are the numbers A195020. This is the main diagonal in the spiral. - Omar E. Pol, Sep 09 2011
Also sequence found by reading the same line mentioned above in the square spiral whose vertices are the generalized enneagonal numbers A118277. Axis perpendicular to A195145 in the same spiral. - Omar E. Pol, Sep 18 2011
Sequence provides all integers m such that 56*m + 49 is a square. - Bruno Berselli, Oct 07 2015
Sum of the numbers from 3*n to 4*n. - Wesley Ivan Hurt, Dec 22 2015

Ivan Panchenko, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A000217, A001106, A002378, A022264, A022265, A028895, A028896, A033996.
Cf. A045943, A046092, A069099, A118277, A174738, A179986, A186029, A193053.
Cf. A195019, A195020, A195145, A218471.

[ (7*n^2 + 7*n)/2 : n in [0..50] ]; // Wesley Ivan Hurt, Jun 09 2014

[seq(7*binomial(n,2), n=1..44)]; # Zerinvary Lajos, Nov 24 2006

7 Table[n (n + 1)/2, {n, 0, 43}] (* or *)
Table[Sum[i, {i, 3 n, 4 n}], {n, 0, 43}] (* or *)
Table[SeriesCoefficient[7 x/(1 - x)^3, {x, 0, n}], {n, 0, 43}] (* Michael De Vlieger, Dec 22 2015 *)
7*Accumulate[Range[0,50]] (* or *) LinearRecurrence[{3,-3,1},{0,7,21},50] (* Harvey P. Dale, Jul 20 2025 *)

x='x+O('x^100); concat(0, Vec(7*x/(1-x)^3)) \\ Altug Alkan, Dec 23 2015

a(n) = (7/2)*n*(n+1).
G.f.: 7*x/(1-x)^3.
a(n) = (7*n^2 + 7*n)/2 = 7*A000217(n). - Omar E. Pol, Dec 12 2008
a(n) = a(n-1) + 7*n with n > 0, a(0)=0. - Vincenzo Librandi, Nov 19 2010
a(n) = A069099(n+1) - 1. - Omar E. Pol, Oct 03 2011
a(n) = a(-n-1), a(n+2) = A193053(n+2) + 2*A193053(n+1) + A193053(n). - Bruno Berselli, Oct 21 2011
From Philippe Deléham, Mar 26 2013: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) with a(0) = 0, a(1) = 7, a(2) = 21.
a(n) = A174738(7*n+6).
a(n) = A179986(n) + n = A186029(n) + 2*n = A022265(n) + 3*n = A022264(n) + 4*n = A218471(n) + 5*n = A001106(n) + 6*n. (End)
a(n) = Sum_{i=3*n..4*n} i. - Wesley Ivan Hurt, Dec 22 2015
E.g.f.: (7/2)*x*(x+2)*exp(x). - G. C. Greubel, Aug 19 2017
From Amiram Eldar, Feb 25 2022: (Start)
Sum_{n>=1} 1/a(n) = 2/7.
Sum_{n>=1} (-1)^(n+1)/a(n) = (2/7)*(2*log(2) - 1). (End)
From Amiram Eldar, Feb 21 2023: (Start)
Product_{n>=1} (1 - 1/a(n)) = -(7/(2*Pi))*cos(sqrt(15/7)*Pi/2).
Product_{n>=1} (1 + 1/a(n)) = (7/(2*Pi))*cosh(Pi/(2*sqrt(7))). (End)

cp's OEIS Frontend

A003215 Hex (or centered hexagonal) numbers: 3*n*(n+1)+1 (crystal ball sequence for hexagonal lattice).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

Maxima

PARI

Python

Formula

Extensions

A000566 Heptagonal numbers (or 7-gonal numbers): n*(5*n-3)/2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

Maxima

PARI

Python

Python

Formula

Extensions

A001106 9-gonal (or enneagonal or nonagonal) numbers: a(n) = n*(7*n-5)/2.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Haskell

Mathematica

PARI

Python

Formula

A062786 Centered 10-gonal numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Mathematica

PARI

PARI

Python

Sage

Formula

Extensions

A069129 Centered 16-gonal numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

A003215 Hex (or centered hexagonal) numbers: 3n(n+1)+1 (crystal ball sequence for hexagonal lattice).

A000566 Heptagonal numbers (or 7-gonal numbers): n(5n-3)/2.

A001106 9-gonal (or enneagonal or nonagonal) numbers: a(n) = n(7n-5)/2.

A101321 Table T(n,m) = 1 + nm(m+1)/2 read by antidiagonals: centered polygonal numbers.

A004126 a(n) = n(7n^2 - 1)/6.

A024966 7 times triangular numbers: 7n(n+1)/2.

A069125 a(n) = (11n^2 - 11n + 2)/2.