login
Search: a024966 -id:a024966
     Sort: relevance | references | number | modified | created      Format: long | short | data
8 times triangular numbers: a(n) = 4*n*(n+1).
+10
82
0, 8, 24, 48, 80, 120, 168, 224, 288, 360, 440, 528, 624, 728, 840, 960, 1088, 1224, 1368, 1520, 1680, 1848, 2024, 2208, 2400, 2600, 2808, 3024, 3248, 3480, 3720, 3968, 4224, 4488, 4760, 5040, 5328, 5624, 5928, 6240, 6560, 6888, 7224, 7568, 7920, 8280
OFFSET
0,2
COMMENTS
Write 0, 1, 2, ... in a clockwise spiral; sequence gives numbers on one of 4 diagonals.
Also, least m > n such that T(m)*T(n) is a square and more precisely that of A055112(n). {T(n) = A000217(n)}. - Lekraj Beedassy, May 14 2004
Also sequence found by reading the line from 0, in the direction 0, 8, ... and the same line from 0, in the direction 0, 24, ..., in the square spiral whose vertices are the generalized decagonal numbers A074377. Axis perpendicular to A195146 in the same spiral. - Omar E. Pol, Sep 18 2011
Number of diagonals with length sqrt(5) in an (n+1) X (n+1) square grid. Every 1 X 2 rectangle has two such diagonals. - Wesley Ivan Hurt, Mar 25 2015
Imagine a board made of squares (like a chessboard), one of whose squares is completely surrounded by square-shaped layers made of adjacent squares. a(n) is the total number of squares in the first to n-th layer. a(1) = 8 because there are 8 neighbors to the unit square; adding them gives a 3 X 3 square. a(2) = 24 = 8 + 16 because we need 16 more squares in the next layer to get a 5 X 5 square: a(n) = (2*n+1)^2 - 1 counting the (2n+1) X (2n+1) square minus the central square. - R. J. Cano, Sep 26 2015
The three platonic solids (the simplex, hypercube, and cross-polytope) with unit side length in n dimensions all have rational volume if and only if n appears in this sequence, after 0. - Brian T Kuhns, Feb 26 2016
The number of active (ON, black) cells in the n-th stage of growth of the two-dimensional cellular automaton defined by "Rule 645", based on the 5-celled von Neumann neighborhood. - Robert Price, May 19 2016
The square root of a(n), n>0, has continued fraction [2n; {1,4n}] with whole number part 2n and periodic part {1,4n}. - Ron Knott, May 11 2017
Numbers k such that k+1 is a square and k is a multiple of 4. - Bruno Berselli, Sep 28 2017
a(n) is the number of vertices of the octagonal network O(n,n); O(m,n) is defined by Fig. 1 of the Siddiqui et al. reference. - Emeric Deutsch, May 13 2018
a(n) is the number of vertices in conjoined n X n octagons which are arranged into a square array, a.k.a. truncated square tiling. - Donghwi Park, Dec 20 2020
a(n-2) is the number of ways to place 3 adjacent marks in a diagonal, horizontal, or vertical row on an n X n tic-tac-toe grid. - Matej Veselovac, May 28 2021
REFERENCES
Stuart M. Ellerstein, J. Recreational Math. 29 (3) 188, 1998.
R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 2nd ed., 1994, p. 99.
Stephen Wolfram, A New Kind of Science, Wolfram Media, 2002; p. 170.
LINKS
M. K. Siddiqui, M. Naeem, N. A. Rahman, and M. Imran, Computing topological indices of certain networks, J. of Optoelectronics and Advanced Materials, 18, No. 9-10, 2016, 884-892.
N. J. A. Sloane, On the Number of ON Cells in Cellular Automata, arXiv:1503.01168 [math.CO], 2015.
Eric Weisstein's World of Mathematics, Elementary Cellular Automaton.
Eric Weisstein's World of Mathematics, Hamiltonian Path.
Eric Weisstein's World of Mathematics, Knight Graph.
Stephen Wolfram, A New Kind of Science
FORMULA
a(n) = 4*n^2 + 4*n = (2*n+1)^2 - 1.
G.f.: 8*x/(1-x)^3.
a(n) = A016754(n) - 1 = 2*A046092(n) = 4*A002378(n). - Lekraj Beedassy, May 25 2004
a(n) = A049598(n) - A046092(n); a(n) = A124080(n) - A002378(n). - Zerinvary Lajos, Mar 06 2007
a(n) = 8*A000217(n). - Omar E. Pol, Dec 12 2008
a(n) = A005843(n) * A163300(n). - Juri-Stepan Gerasimov, Jul 26 2009
a(n) = a(n-1) + 8*n (with a(0)=0). - Vincenzo Librandi, Nov 17 2010
For n > 0, a(n) = A058031(n+1) - A062938(n-1). - Charlie Marion, Apr 11 2013
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Wesley Ivan Hurt, Mar 25 2015
a(n) = A000578(n+1) - A152618(n). - Bui Quang Tuan, Apr 01 2015
a(n) - a(n-1) = A008590(n), n > 0. - Altug Alkan, Sep 26 2015
From Ilya Gutkovskiy, May 19 2016: (Start)
E.g.f.: 4*x*(2 + x)*exp(x).
Sum_{n>=1} 1/a(n) = 1/4. (End)
Product_{n>=1} a(n)/A016754(n) = Pi/4. - Daniel Suteu, Dec 25 2016
a(n) = A056220(n) + A056220(n+1). - Bruce J. Nicholson, May 29 2017
sqrt(a(n)+1) - sqrt(a(n)) = (sqrt(n+1) - sqrt(n))^2. - Seiichi Manyama, Dec 23 2018
a(n)*a(n+k) + 4*k^2 = m^2 where m = (a(n) + a(n+k))/2 - 2*k^2; for k=1, m = 4*n^2 + 8*n + 2 = A060626(n). - Ezhilarasu Velayutham, May 22 2019
Sum_{n>=1} (-1)^n/a(n) = 1/4 - log(2)/2. - Vaclav Kotesovec, Dec 21 2020
From Amiram Eldar, Feb 21 2023: (Start)
Product_{n>=1} (1 - 1/a(n)) = -(4/Pi)*cos(Pi/sqrt(2)).
Product_{n>=1} (1 + 1/a(n)) = 4/Pi (A088538). (End)
EXAMPLE
Spiral with 0, 8, 24, 48, ... along lower right diagonal:
.
36--37--38--39--40--41--42
| |
35 16--17--18--19--20 43
| | | |
34 15 4---5---6 21 44
| | | | | |
33 14 3 0 7 22 45
| | | | \ | | |
32 13 2---1 8 23 46
| | | \ | |
31 12--11--10---9 24 47
| | \ |
30--29--28--27--26--25 48
\
[Reformatted by Jon E. Schoenfield, Dec 25 2016]
MAPLE
seq(8*binomial(n+1, 2), n=0..46); # Zerinvary Lajos, Nov 24 2006
[seq((2*n+1)^2-1, n=0..46)];
MATHEMATICA
Table[(2n - 1)^2 - 1, {n, 50}] (* Alonso del Arte, Mar 31 2013 *)
PROG
(PARI) nsqm1(n) = { forstep(x=1, n, 2, y = x*x-1; print1(y, ", ") ) }
(Magma) [ 4*n*(n+1) : n in [0..50] ]; // Wesley Ivan Hurt, Jun 09 2014
CROSSREFS
Cf. A000217, A016754, A002378, A024966, A027468, A028895, A028896, A045943, A046092, A049598, A088538, A124080, A008590 (first differences), A130809 (partial sums).
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.
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Dec 11 1999
STATUS
approved
Centered heptagonal numbers.
+10
63
1, 8, 22, 43, 71, 106, 148, 197, 253, 316, 386, 463, 547, 638, 736, 841, 953, 1072, 1198, 1331, 1471, 1618, 1772, 1933, 2101, 2276, 2458, 2647, 2843, 3046, 3256, 3473, 3697, 3928, 4166, 4411, 4663, 4922, 5188, 5461, 5741, 6028, 6322, 6623, 6931, 7246
OFFSET
1,2
COMMENTS
Equals the triangular numbers convolved with [ 1, 5, 1, 0, 0, 0, ...]. - Gary W. Adamson and Alexander R. Povolotsky, May 29 2009
Number of ordered pairs of integers (x,y) with abs(x) < n, abs(y) < n and abs(x + y) < n, counting twice pairs of equal numbers. - Reinhard Zumkeller, Jan 23 2012; corrected and extended by Mauro Fiorentini, Jan 01 2018
The number of pairs without repetitions is a(n) - 2n + 3 for n > 1. For example, there are 19 such pairs for n = 3: (-2, 0), (-2, 1), (-2, 2), (-1, -1), (-1, 0), (-1, 1), (-1, 2), (0, -2), (0, -1), (0, 0), (0, 1), (0, 2), (1, -2), (1, -1), (1, 0), (1, 1), (2, -2), (2, -1), (2, 0). - Mauro Fiorentini, Jan 01 2018
FORMULA
a(n) = (7*n^2 - 7*n + 2)/2.
a(n) = 1 + Sum_{k=1..n} 7*k. - Xavier Acloque, Oct 26 2003
Binomial transform of [1, 7, 7, 0, 0, 0, ...]; Narayana transform (A001263) of [1, 7, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
a(n) = 7*n + a(n-1) - 7 (with a(1)=1). - Vincenzo Librandi, Aug 08 2010
G.f.: x*(1+5*x+x^2) / (1-x)^3. - R. J. Mathar, Feb 04 2011
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3); a(0)=1, a(1)=8, a(2)=22. - Harvey P. Dale, Jun 04 2011
a(n) = A024966(n-1) + 1. - Omar E. Pol, Oct 03 2011
a(n) = 2*a(n-1) - a(n-2) + 7. - Ant King, Jun 17 2012
From Ant King, Jun 17 2012: (Start)
Sum_{n>=1} 1/a(n) = 2*Pi/sqrt(7)*tanh(Pi/(2*sqrt(7))) = 1.264723171685652...
a(n) == 1 (mod 7) for all n.
The sequence of digital roots of the a(n) is period 9: repeat [1, 8, 4, 7, 8, 7, 4, 8, 1] (the period is a palindrome).
The sequence of a(n) mod 10 is period 20: repeat [1, 8, 2, 3, 1, 6, 8, 7, 3, 6, 6, 3, 7, 8, 6, 1, 3, 2, 8, 1] (the period is a palindrome).
(End)
E.g.f.: -1 + (2 + 7*x^2)*exp(x)/2. - Ilya Gutkovskiy, Jun 30 2016
a(n) = A101321(7,n-1). - R. J. Mathar, Jul 28 2016
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=1} a(n)/n! = 9*e/2 - 1.
Sum_{n>=1} (-1)^n * a(n)/n! = 9/(2*e) - 1. (End)
a(n) = A003215(n-1) + A000217(n-1). - Leo Tavares, Jul 19 2022
EXAMPLE
a(5) = 71 because 71 = (7*5^2 - 7*5 + 2)/2 = (175 - 35 + 2)/2 = 142/2.
From Bruno Berselli, Oct 27 2017: (Start)
1 = -(0) + (1).
8 = -(0+1) + (2+3+4).
22 = -(0+1+2) + (3+4+5+6+7).
43 = -(0+1+2+3) + (4+5+6+7+8+9+10).
71 = -(0+1+2+3+4) + (5+6+7+8+9+10+11+12+13). (End)
MATHEMATICA
FoldList[#1 + #2 &, 1, 7 Range@ 50] (* Robert G. Wilson v, Feb 02 2011 *)
LinearRecurrence[{3, -3, 1}, {1, 8, 22}, 50] (* Harvey P. Dale, Jun 04 2011 *)
PROG
(Haskell)
a069099 n = length
[(x, y) | x <- [-n+1..n-1], y <- [-n+1..n-1], x + y <= n - 1]
-- Reinhard Zumkeller, Jan 23 2012
(PARI) a(n)=(7*n^2-7*n+2)/2 \\ Charles R Greathouse IV, Sep 24 2015
CROSSREFS
Cf. A000566 (heptagonal numbers).
KEYWORD
nonn,easy,nice
AUTHOR
Terrel Trotter, Jr., Apr 05 2002
STATUS
approved
Centered 10-gonal numbers.
+10
41
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,2
COMMENTS
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
FORMULA
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)
a(n) = 10*n + a(n-1) - 10; a(1)=1. - Vincenzo Librandi, Aug 07 2010
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)
MATHEMATICA
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 *)
PROG
(PARI) j=[]; for(n=1, 75, j=concat(j, (5*n*(n-1)+1))); j
(PARI) for (n=1, 1000, write("b062786.txt", n, " ", 5*n*(n - 1) + 1) ) \\ Harry J. Smith, Aug 11 2009
(Magma) [1+5*n*(n-1): n in [1..50]]; // G. C. Greubel, Mar 30 2019
(Sage) [1+5*rising_factorial(n-1, 2) for n in (1..50)] # G. C. Greubel, Mar 30 2019
(GAP) List([1..50], n-> 1+5*n*(n-1)); # G. C. Greubel, Mar 30 2019
(Python) def a(n): return(5*n**2-5*n+1) # Torlach Rush, May 10 2024
KEYWORD
easy,nonn
AUTHOR
Jason Earls, Jul 19 2001
EXTENSIONS
Better description from Terrel Trotter, Jr., Apr 06 2002
STATUS
approved
6 times triangular numbers: a(n) = 3*n*(n+1).
+10
39
0, 6, 18, 36, 60, 90, 126, 168, 216, 270, 330, 396, 468, 546, 630, 720, 816, 918, 1026, 1140, 1260, 1386, 1518, 1656, 1800, 1950, 2106, 2268, 2436, 2610, 2790, 2976, 3168, 3366, 3570, 3780, 3996, 4218, 4446, 4680, 4920, 5166, 5418, 5676
OFFSET
0,2
COMMENTS
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, 6, ...
The spiral begins:
85--84--83--82--81--80
/ \
86 56--55--54--53--52 79
/ / \ \
87 57 33--32--31--30 51 78
/ / / \ \ \
88 58 34 16--15--14 29 50 77
/ / / / \ \ \ \
89 59 35 17 5---4 13 28 49 76
/ / / / / \ \ \ \ \
<==90==60==36==18===6===0 3 12 27 48 75
/ / / / / / / / / /
61 37 19 7 1---2 11 26 47 74
\ \ \ \ / / / /
62 38 20 8---9--10 25 46 73
\ \ \ / / /
63 39 21--22--23--24 45 72
\ \ / /
64 40--41--42--43--44 71
\ /
65--66--67--68--69--70
(End)
If Y is a 4-subset of an n-set X then, for n >= 5, a(n-5) is the number of (n-4)-subsets of X having exactly two elements in common with Y. - Milan Janjic, Dec 28 2007
a(n) is the maximal number of points of intersection of n+1 distinct triangles drawn in the plane. For example, two triangles can intersect in at most a(1) = 6 points (as illustrated in the Star of David configuration). - Terry Stickels (Terrystickels(AT)aol.com), Jul 12 2008
Also sequence found by reading the line from 0, in the direction 0, 6, ... and the same line from 0, in the direction 0, 18, ..., in the square spiral whose vertices are the generalized octagonal numbers A001082. Axis perpendicular to A195143 in the same spiral. - Omar E. Pol, Sep 18 2011
Partial sums of A008588. - R. J. Mathar, Aug 28 2014
Also the number of 5-cycles in the (n+5)-triangular honeycomb acute knight graph. - Eric W. Weisstein, Jul 27 2017
a(n-4) is the maximum irregularity over all maximal 3-degenerate graphs with n vertices. The extremal graphs are 3-stars (K_3 joined to n-3 independent vertices). (The irregularity of a graph is the sum of the differences between the degrees over all edges of the graph.) - Allan Bickle, May 29 2023
LINKS
Allan Bickle and Zhongyuan Che, Irregularities of Maximal k-degenerate Graphs, Discrete Applied Math. 331 (2023) 70-87.
Allan Bickle, A Survey of Maximal k-degenerate Graphs and k-Trees, Theory and Applications of Graphs 0 1 (2024) Article 5.
Enrique Navarrete and Daniel Orellana, Finding Prime Numbers as Fixed Points of Sequences, arXiv:1907.10023 [math.NT], 2019.
Luis Manuel Rivera, Integer sequences and k-commuting permutations, arXiv preprint arXiv:1406.3081 [math.CO], 2014.
Eric Weisstein's World of Mathematics, Graph Cycle.
FORMULA
O.g.f.: 6*x/(1 - x)^3.
E.g.f.: 3*x*(x + 2)*exp(x). - G. C. Greubel, Aug 19 2017
a(n) = 6*A000217(n).
a(n) = polygorial(3, n+1). - Daniel Dockery (peritus(AT)gmail.com), Jun 16 2003
From Zerinvary Lajos, Mar 06 2007: (Start)
a(n) = A049598(n)/2.
a(n) = A124080(n) - A046092(n).
a(n) = A033996(n) - A002378(n). (End)
a(n) = A002378(n)*3 = A045943(n)*2. - Omar E. Pol, Dec 12 2008
a(n) = a(n-1) + 6*n for n>0, a(0)=0. - Vincenzo Librandi, Aug 05 2010
a(n) = A003215(n) - 1. - Omar E. Pol, Oct 03 2011
From Philippe Deléham, Mar 26 2013: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n>2, a(0)=0, a(1)=6, a(2)=18.
a(n) = A174709(6*n + 5). (End)
a(n) = A049450(n) + 4*n. - Lear Young, Apr 24 2014
a(n) = Sum_{i = n..2*n} 2*i. - Bruno Berselli, Feb 14 2018
a(n) = A320047(1, n, 1). - Kolosov Petro, Oct 04 2018
a(n) = T(3*n) - T(2*n-2) + T(n-2), 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 04 2020
From Amiram Eldar, Feb 15 2022: (Start)
Sum_{n>=1} 1/a(n) = 1/3.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/3 - 1/3. (End)
From Amiram Eldar, Feb 21 2023: (Start)
Product_{n>=1} (1 - 1/a(n)) = -(3/Pi)*cos(sqrt(7/3)*Pi/2).
Product_{n>=1} (1 + 1/a(n)) = (3/Pi)*cosh(Pi/(2*sqrt(3))). (End)
MAPLE
[seq(6*binomial(n, 2), n=1..44)]; # Zerinvary Lajos, Nov 24 2006
MATHEMATICA
6 Accumulate[Range[0, 50]] (* Harvey P. Dale, Mar 05 2012 *)
6 PolygonalNumber[Range[0, 20]] (* Eric W. Weisstein, Jul 27 2017 *)
LinearRecurrence[{3, -3, 1}, {0, 6, 18}, 20] (* Eric W. Weisstein, Jul 27 2017 *)
PROG
(Magma) [3*n*(n+1): n in [0..50]]; // Wesley Ivan Hurt, Jun 09 2014
(PARI) a(n)=3*n*(n+1) \\ Charles R Greathouse IV, Sep 24 2015
(PARI) first(n) = Vec(6*x/(1 - x)^3 + O(x^n), -n) \\ Iain Fox, Feb 14 2018
(GAP) List([0..44], n->3*n*(n+1)); # Muniru A Asiru, Mar 15 2019
CROSSREFS
Cf. A002378 (3-cycles in triangular honeycomb acute knight graph), A045943 (4-cycles), A152773 (6-cycles).
Cf. A007531.
The partial sums give A007531. - Leo Tavares, Jan 22 2022
Cf. A002378, A046092, A028896 (irregularities of maximal k-degenerate graphs).
KEYWORD
nonn,easy
AUTHOR
Joe Keane (jgk(AT)jgk.org), Dec 11 1999
STATUS
approved
9 times the triangular numbers A000217.
+10
28
0, 9, 27, 54, 90, 135, 189, 252, 324, 405, 495, 594, 702, 819, 945, 1080, 1224, 1377, 1539, 1710, 1890, 2079, 2277, 2484, 2700, 2925, 3159, 3402, 3654, 3915, 4185, 4464, 4752, 5049, 5355, 5670, 5994, 6327, 6669, 7020, 7380, 7749, 8127, 8514, 8910, 9315
OFFSET
0,2
COMMENTS
Staggered diagonal of triangular spiral in A051682, between (0,1,11) spoke and (0,8,25) spoke. - Paul Barry, Mar 15 2003
Number of permutations of n distinct letters (ABCD...) each of which appears thrice with n-2 fixed points. - Zerinvary Lajos, Oct 15 2006
Number of n permutations (n>=2) of 4 objects u, v, z, x with repetition allowed, containing n-2=0 u's. Example: if n=2 then n-2 =zero (0) u, a(1)=9 because we have vv, zz, xx, vx, xv, zx, xz, vz, zv. A027465 formatted as a triangular array: diagonal: 9, 27, 54, 90, 135, 189, 252, 324, ... . - Zerinvary Lajos, Aug 06 2008
a(n) is also the least weight of self-conjugate partitions having n different parts such that each part is a multiple of 3. - Augustine O. Munagi, Dec 18 2008
Also sequence found by reading the line from 0, in the direction 0, 9, ..., and the same line from 0, in the direction 0, 27, ..., in the square spiral whose vertices are the generalized hendecagonal numbers A195160. Axis perpendicular to A195147 in the same spiral. - Omar E. Pol, Sep 18 2011
Sum of the numbers from 4*n to 5*n. - Wesley Ivan Hurt, Nov 01 2014
LINKS
Augustine O. Munagi, Pairing conjugate partitions by residue classes, Discrete Math., Vol. 308, No. 12 (2008), pp. 2492-2501.
Enrique Navarrete and Daniel Orellana, Finding Prime Numbers as Fixed Points of Sequences, arXiv:1907.10023 [math.NT], 2019.
D. Zvonkine, Home Page.
FORMULA
Numerators of sequence a[n, n-2] in (a[i, j])^2 where a[i, j] = binomial(i-1, j-1)/2^(i-1) if j<=i, 0 if j>i.
a(n) = (9/2)*n*(n+1).
a(n) = 9*C(n, 1) + 9*C(n, 2) (binomial transform of (0, 9, 9, 0, 0, ...)). - Paul Barry, Mar 15 2003
G.f.: 9*x/(1-x)^3.
a(-1-n) = a(n).
a(n) = 9*C(n+1,2), n>=0. - Zerinvary Lajos, Aug 06 2008
a(n) = a(n-1) + 9*n (with a(0)=0). - Vincenzo Librandi, Nov 19 2010
a(n) = A060544(n+1) - 1. - Omar E. Pol, Oct 03 2011
a(n) = A218470(9*n+8). - Philippe Deléham, Mar 27 2013
E.g.f.: (9/2)*x*(x+2)*exp(x). - G. C. Greubel, Aug 22 2017
a(n) = A060544(n+1) - 1. See Centroid Triangles illustration. - Leo Tavares, Dec 27 2021
From Amiram Eldar, Feb 15 2022: (Start)
Sum_{n>=1} 1/a(n) = 2/9.
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*log(2)/9 - 2/9. (End)
From Amiram Eldar, Feb 21 2023: (Start)
Product_{n>=1} (1 - 1/a(n)) = -(9/(2*Pi))*cos(sqrt(17)*Pi/6).
Product_{n>=1} (1 + 1/a(n)) = 9*sqrt(3)/(4*Pi). (End)
EXAMPLE
The first such self-conjugate partitions, corresponding to a(n)=1,2,3,4 are 3+3+3, 6+6+6+3+3+3, 9+9+9+6+6+6+3+3+3, 12+12+12+9+9+9+6+6+6+3+3+3. - Augustine O. Munagi, Dec 18 2008
MAPLE
[seq(9*binomial(n+1, 2), n=0..50)]; # Zerinvary Lajos, Nov 24 2006
MATHEMATICA
Table[(9/2)*n*(n+1), {n, 0, 50}] (* G. C. Greubel, Aug 22 2017 *)
PROG
(PARI) a(n)=9*n*(n+1)/2
(Magma) [9*n*(n+1)/2: n in [0..50]]; // Vincenzo Librandi, Dec 29 2012
(Sage) [9*binomial(n+1, 2) for n in (0..50)] # G. C. Greubel, May 20 2021
KEYWORD
nonn,easy
EXTENSIONS
More terms from Patrick De Geest, Oct 15 1999
STATUS
approved
Vertex number of a square spiral in which the length of the first two edges are the legs of the primitive Pythagorean triple [3, 4, 5]. The edges of the spiral have length A195019.
+10
28
0, 3, 7, 13, 21, 30, 42, 54, 70, 85, 105, 123, 147, 168, 196, 220, 252, 279, 315, 345, 385, 418, 462, 498, 546, 585, 637, 679, 735, 780, 840, 888, 952, 1003, 1071, 1125, 1197, 1254, 1330, 1390, 1470, 1533, 1617, 1683, 1771, 1840, 1932, 2004, 2100
OFFSET
0,2
COMMENTS
Zero together with the partial sums of A195019.
The spiral contains infinitely many Pythagorean triples in which the hypotenuses on the main diagonal are the positives A008587. The vertices on the main diagonal are the numbers A024966 = (3+4)*A000217 = 7*A000217, where both 3 and 4 are the first two edges in the spiral. The distance "a" between nearest edges that are perpendicular to the initial edge of the spiral is 3, while the distance "b" between nearest edges that are parallel to the initial edge is 4, so the distance "c" between nearest vertices on the same axis is 5 because from the Pythagorean theorem we can write c = (a^2+b^2)^(1/2) = sqrt(3^2+4^2) = sqrt(9+16) = sqrt(25) = 5.
Let an array have m(0,n)=m(n,0)=n*(n-1)/2 and m(n,n)=n*(n+1)/2. The first n+1 terms in row(n) are the numbers in the closed interval m(0,n) to m(n,n). The terms in column(n) are the same from m(n,0) to m(n,n). The first few antidiagonals are 0; 0,0; 1,1,1; 3,2,2,3; 6,4,3,4,6; 10,7,5,5,7,10. a(n) is the difference between the sum of the terms in the n+1 X n+1 matrices and those in the n X n matrices. - J. M. Bergot, Jul 05 2013 [The first five rows are: 0,0,1,3,6; 0,1,2,4,7; 1,2,3,5,8; 3,4,5,6,9; 6,7,8,9,10]
LINKS
Eric Weisstein's World of Mathematics, Pythagorean Triple
FORMULA
From Bruno Berselli, Oct 13 2011: (Start)
G.f.: x*(3+4*x)/((1+x)^2*(1-x)^3).
a(n) = (1/2)*A004526(n+2)*A047335(n+1) = (2*n*(7*n+13) + (2*n-5)*(-1)^n+5)/16.
a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5).
a(n) - a(n-2) = A047355(n+1). (End)
MATHEMATICA
With[{r = Range[50]}, Join[{0}, Accumulate[Riffle[3*r, 4*r]]]] (* or *)
LinearRecurrence[{1, 2, -2, -1, 1}, {0, 3, 7, 13, 21}, 100] (* Paolo Xausa, Feb 09 2024 *)
PROG
(Magma) [(2*n*(7*n+13)+(2*n-5)*(-1)^n+5)/16: n in [0..50]]; // Vincenzo Librandi, Oct 14 2011
KEYWORD
nonn,easy
AUTHOR
Omar E. Pol, Sep 07 2011 - Sep 12 2011
STATUS
approved
Second 9-gonal (or nonagonal) numbers: a(n) = n*(7*n+5)/2.
+10
18
0, 6, 19, 39, 66, 100, 141, 189, 244, 306, 375, 451, 534, 624, 721, 825, 936, 1054, 1179, 1311, 1450, 1596, 1749, 1909, 2076, 2250, 2431, 2619, 2814, 3016, 3225, 3441, 3664, 3894, 4131, 4375, 4626, 4884, 5149, 5421, 5700, 5986, 6279, 6579, 6886
OFFSET
0,2
COMMENTS
This sequence is a bisection of A118277 (even part).
Sequence found by reading the line from 0, in the direction 0, 19... and the line from 6, in the direction 6, 39,..., in the square spiral whose vertices are the generalized 9-gonal numbers A118277. - Omar E. Pol, Jul 24 2012
The early part of this sequence is a strikingly close approximation to the early part of A100752. - Peter Munn, Nov 14 2019
FORMULA
G.f.: x*(6 + x)/(1 - x)^3.
a(n) = Sum_{i=0..(n-1)} A017053(i) for n>0.
a(-n) = A001106(n).
Sum_{i=0..n} (a(n)+i)^2 = ( Sum_{i=(n+1)..2*n} (a(n)+i)^2 ) + 21*A000217(n)^2 for n>0.
a(n) = a(n-1)+7*n-1 for n>0, with a(0)=0. - Vincenzo Librandi, Feb 05 2011
a(0)=0, a(1)=6, a(2)=19; for n>2, a(n) = 3*a(n-1)-3*a(n-2)+a(n-3). - Harvey P. Dale, Aug 19 2011
a(n) = A174738(7n+5). - Philippe Deléham, Mar 26 2013
a(n) = A001477(n) + 2*A000290(n) + 3*A000217(n). - J. M. Bergot, Apr 25 2014
a(n) = A055998(4*n) - A055998(3*n). - Bruno Berselli, Sep 23 2016
E.g.f.: (x/2)*(12 + 7*x)*exp(x). - G. C. Greubel, Aug 19 2017
MATHEMATICA
f[n_] := n (7 n + 5)/2; f[Range[0, 60]] (* Vladimir Joseph Stephan Orlovsky, Feb 05 2011*)
LinearRecurrence[{3, -3, 1}, {0, 6, 19}, 60] (* or *) Array[(#(7# + 5))/2&, 60, 0] (* Harvey P. Dale, Aug 19 2011 *)
CoefficientList[Series[x (6 + x)/(1 - x)^3, {x, 0, 60}], x] (* Vincenzo Librandi, Oct 15 2012 *)
PROG
(Magma) [n*(7*n+5)/2: n in [0..50]]; // Bruno Berselli, Sep 23 2016
(Magma) I:=[0, 6, 19]; [n le 3 select I[n] else 3*Self(n-1) -3*Self(n-2) +Self(n-3): n in [1..60]]; // Vincenzo Librandi, Oct 15 2012
(PARI) a(n)=n*(7*n+5)/2 \\ Charles R Greathouse IV, Sep 24 2015
CROSSREFS
Cf. second k-gonal numbers: A005449 (k=5), A014105 (k=6), A147875 (k=7), A045944 (k=8), this sequence (k=9), A033954 (k=10), A062728 (k=11), A135705 (k=12).
KEYWORD
nonn,easy
AUTHOR
Bruno Berselli, Jan 13 2011
STATUS
approved
Partial sums of floor(n/7).
+10
17
0, 0, 0, 0, 0, 0, 0, 1, 2, 3, 4, 5, 6, 7, 9, 11, 13, 15, 17, 19, 21, 24, 27, 30, 33, 36, 39, 42, 46, 50, 54, 58, 62, 66, 70, 75, 80, 85, 90, 95, 100, 105, 111, 117, 123, 129, 135, 141, 147, 154, 161, 168, 175, 182, 189, 196, 204, 212, 220, 228, 236
OFFSET
0,9
COMMENTS
Apart from the initial zeros, the same as A011867.
LINKS
Mircea Merca, Inequalities and Identities Involving Sums of Integer Functions, J. Integer Sequences, Vol. 14 (2011), Article 11.9.1.
FORMULA
a(n) = round(n*(n-5)/14).
a(n) = floor((n-2)*(n-3)/14).
a(n) = ceiling((n+1)*(n-6)/14).
a(n) = a(n-7) + n - 6, n > 6.
a(n) = +2*a(n-1) - a(n-2) + a(n-7) - 2*a(n-8) + a(n-9). - R. J. Mathar, Nov 30 2010
G.f.: x^7/( (1 + x + x^2 + x^3 + x^4 + x^5 + x^6)*(1-x)^3 ). - R. J. Mathar, Nov 30 2010
a(7n) = A001106(n), a(7n+1) = A218471(n), a(7n+2) = A022264(n), a(7n+3) = A022265(n), a(7n+4) = A186029(n), a(7n+5) = A179986(n), a(7n+6) = A024966(n). - Philippe Deléham, Mar 26 2013
EXAMPLE
a(9) = floor(0/7) + floor(1/7) + floor(2/7) + floor(3/7) + floor(4/7) + floor(5/7) + floor(6/7) + floor(7/7) + floor(8/7) + floor(9/7) = 3.
MAPLE
A174738 := proc(n) round(n*(n-5)/14) ; end proc:
seq(A174738(n), n=0..30) ;
MATHEMATICA
Table[Floor[(n - 2)*(n - 3)/14], {n, 0, 60}] (* G. C. Greubel, Dec 13 2016 *)
PROG
(Magma) [Round(n*(n-5)/14): n in [0..60]]; // Vincenzo Librandi, Jun 22 2011
(PARI) a(n)=(n-2)*(n-3)\14 \\ Charles R Greathouse IV, Sep 24 2015
(Sage) [floor((n-2)*(n-3)/14) for n in (0..60)] # G. C. Greubel, Aug 31 2019
(GAP) List([0..60], n-> Int((n-2)*(n-3)/14)); # G. C. Greubel, Aug 31 2019
KEYWORD
nonn,easy
AUTHOR
Mircea Merca, Nov 30 2010
STATUS
approved
a(n) = n*(7*n+3)/2.
+10
16
0, 5, 17, 36, 62, 95, 135, 182, 236, 297, 365, 440, 522, 611, 707, 810, 920, 1037, 1161, 1292, 1430, 1575, 1727, 1886, 2052, 2225, 2405, 2592, 2786, 2987, 3195, 3410, 3632, 3861, 4097, 4340, 4590, 4847, 5111, 5382, 5660, 5945, 6237, 6536, 6842, 7155, 7475
OFFSET
0,2
COMMENTS
This sequence is related to A050409 by A050409(n) = n*a(n) - Sum_{i=0..n-1} a(i).
FORMULA
G.f.: x*(5+2*x)/(1-x)^3.
a(n) - a(-n) = A008585(n).
a(n) + a(-n) = A033582(n).
n*a(n+1) - (n+1)*a(n) = A024966(n). - Bruno Berselli, May 30 2012
n*a(n+2) - (n+2)*a(n) = A067727(n) for n>0. - Bruno Berselli, May 30 2012
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n>2, a(0)=0, a(1)=5, a(2)=17. - Philippe Deléham, Mar 26 2013
a(n) = A174738(7*n+4). - Philippe Deléham, Mar 26 2013
E.g.f.: (1/2)*(7*x^2 + 10*x)*exp(x). - G. C. Greubel, Jul 17 2017
EXAMPLE
From Ilya Gutkovskiy, Mar 31 2016: (Start)
. o o o o o o o o o o o o
. o o
. o o o o o o o o o o o o o o o o o o o o o o o o
. o o o o o o o o o o o o
. o o o o o o o o o o o o o o o o o o o o
. o o o o o o o o o o o o o o o o o o o o o o o o
. o o 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)
MATHEMATICA
Table[(n - 1) (7 n - 4)/2, {n, 100}] (* Vladimir Joseph Stephan Orlovsky, Jul 06 2011 *)
LinearRecurrence[{3, -3, 1}, {0, 5, 17}, 50] (* Harvey P. Dale, Sep 07 2022 *)
PROG
(Magma) [n*(7*n+3)/2: n in [0..44]];
(PARI) a(n)=n*(7*n+3)/2 \\ Charles R Greathouse IV, Sep 24 2015
CROSSREFS
Cf. numbers of the form n*(d*n+10-d)/2 indexed in A140090.
Cf. A017041 (first differences).
KEYWORD
nonn,easy
AUTHOR
Bruno Berselli, Feb 11 2011
STATUS
approved
a(n) = n*(7*n + 1)/2.
+10
15
0, 4, 15, 33, 58, 90, 129, 175, 228, 288, 355, 429, 510, 598, 693, 795, 904, 1020, 1143, 1273, 1410, 1554, 1705, 1863, 2028, 2200, 2379, 2565, 2758, 2958, 3165, 3379, 3600, 3828, 4063, 4305, 4554, 4810
OFFSET
0,2
COMMENTS
For n >= 4, a(n) is the sum of the numbers appearing in the 4th row of an n X n square array whose elements are the numbers from 1..n^2, listed in increasing order by rows. - Wesley Ivan Hurt, May 17 2021
FORMULA
a(n) = A110449(n, 3) for n>2.
a(n) = A049453(n) - A005475(n). - Zerinvary Lajos, Jan 21 2007
a(n) = 7*n + a(n-1) - 3 for n>0, a(0)=0. - Vincenzo Librandi, Aug 04 2010
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) with a(0)=0, a(1)=4, a(2)=15. - Philippe Deléham, Mar 26 2013
a(n) = A174738(7n+3). - Philippe Deléham, Mar 26 2013
a(n) = A000217(4*n) - A000217(3*n). - Bruno Berselli, Oct 13 2016
G.f.: x*(4 + 3*x)/(1 - x)^3. - Ilya Gutkovskiy, Oct 13 2016
E.g.f.: (x/2)*(7*x + 8)*exp(x). - G. C. Greubel, Aug 23 2017
EXAMPLE
From Bruno Berselli, Oct 27 2017: (Start)
After 0:
4 = -(1) + (2 + 3).
15 = -(1 + 2) + (3 + 4 + 5 + 6).
33 = -(1 + 2 + 3) + (4 + 5 + 6 + 7 + 8 + 9).
58 = -(1 + 2 + 3 + 4) + (5 + 6 + 7 + 8 + 9 + 10 + 11 + 12). (End)
MAPLE
seq(binomial(7*n+1, 2)/7, n=0..37); # Zerinvary Lajos, Jan 21 2007
seq(binomial(6*n+1, 2)/3-binomial(5*n+1, 2)/5, n=0..42); # Zerinvary Lajos, Jan 21 2007
MATHEMATICA
Table[n (7 n + 1)/2, {n, 0, 40}] (* Bruno Berselli, Oct 13 2016 *)
LinearRecurrence[{3, -3, 1}, {0, 4, 15}, 40] (* Harvey P. Dale, Oct 09 2018 *)
PROG
(PARI) a(n)=n*(7*n+1)/2 \\ Charles R Greathouse IV, Oct 07 2015
CROSSREFS
Cf. similar sequences listed in A022289.
KEYWORD
nonn,easy
STATUS
approved

Search completed in 0.086 seconds