A069128 -id:A069128 - cp's OEIS frontend

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

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

A069190 Centered 24-gonal numbers.

Original entry on oeis.org

1, 25, 73, 145, 241, 361, 505, 673, 865, 1081, 1321, 1585, 1873, 2185, 2521, 2881, 3265, 3673, 4105, 4561, 5041, 5545, 6073, 6625, 7201, 7801, 8425, 9073, 9745, 10441, 11161, 11905, 12673, 13465, 14281, 15121, 15985, 16873, 17785, 18721, 19681, 20665, 21673

Offset: 1

Terrel Trotter, Jr., Apr 10 2002

Sequence found by reading the line from 1, in the direction 1, 25, ..., in the square spiral whose vertices are the generalized octagonal numbers A001082. Semi-axis opposite to A135453 in the same spiral. - Omar E. Pol, Sep 16 2011

			a(5) = 241 because 12*5^2 - 12*5 + 1 = 300 - 60 + 1 = 241.

Ivan Panchenko, Table of n, a(n) for n = 1..1000
John Elias, Illustration: Odd Ordered Star Perimeters.
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. centered k-gonal numbers with k=3..25: A005448, A001844, A005891, A003215, A069099, A016754, A060544, A062786, A069125, A003154, A069126, A069127, A069128, A069129, A069130, A069131, A069132, A069133, A069178, A069173, A069174, A069190, A262221.

FoldList[#1 + #2 &, 1, 24 Range@ 45] (* Robert G. Wilson v *)
Table[12n^2-12n+1,{n,50}] (* or *) LinearRecurrence[{3,-3,1},{1,25,73},50] (* Harvey P. Dale, Jul 17 2011 *)

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

a(n) = 12*n^2 - 12*n + 1.
a(n) = 24*n + a(n-1) - 24 with a(1)=1. - Vincenzo Librandi, Aug 08 2010
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3); a(1)=1, a(2)=25, a(3)=73. - Harvey P. Dale, Jul 17 2011
G.f.: x*(1+22*x+x^2)/(1-x)^3. - Harvey P. Dale, Jul 17 2011
Binomial transform of [1, 24, 24, 0, 0, 0, ...] and Narayana transform (cf. A001263) of [1, 24, 0, 0, 0, ...]. - Gary W. Adamson, Jul 26 2011
From Amiram Eldar, Jun 21 2020: (Start)
Sum_{n>=1} 1/a(n) = Pi*tan(Pi/sqrt(6))/(4*sqrt(6)).
Sum_{n>=1} a(n)/n! = 13*e - 1.
Sum_{n>=1} (-1)^n * a(n)/n! = 13/e - 1. (End)
E.g.f.: exp(x)*(1 + 12*x^2) - 1. - Stefano Spezia, May 31 2022

More terms from Harvey P. Dale, Jul 17 2011

A262221 a(n) = 25n(n + 1)/2 + 1.

Original entry on oeis.org

1, 26, 76, 151, 251, 376, 526, 701, 901, 1126, 1376, 1651, 1951, 2276, 2626, 3001, 3401, 3826, 4276, 4751, 5251, 5776, 6326, 6901, 7501, 8126, 8776, 9451, 10151, 10876, 11626, 12401, 13201, 14026, 14876, 15751, 16651, 17576, 18526, 19501, 20501, 21526, 22576, 23651

Offset: 0

Bruno Berselli, Sep 15 2015

Also centered 25-gonal (or icosipentagonal) numbers.
This is the case k=25 of the formula (k*n*(n+1) - (-1)^k + 1)/2. See table in Links section for similar sequences.
For k=2*n, the formula shown above gives A011379.
Primes in sequence: 151, 251, 701, 1951, 3001, 4751, 10151, 12401, ...

E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 51 (23rd row of the table).

Bruno Berselli, Table of n, a(n) for n = 0..1000
Bruno Berselli, Table of numbers of the form (k*n*(n+1)-(-1)^k+1)/2.
Eric Weisstein's World of Mathematics, Centered Polygonal Numbers.
Index entries for sequences related to centered polygonal numbers.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A000217, A008607 (first differences), A011379, A034856, A054254, A080956, A123296, A186349, A255184.
Cf. centered polygonal numbers listed in A069190.
Similar sequences of the form (k*n*(n+1) - (-1)^k + 1)/2 with -1 <= k <= 26: A000004, A000124, A002378, A005448, A005891, A028896, A033996, A035008, A046092, A049598, A060544, A064200, A069099, A069125, A069126, A069128, A069130, A069132, A069174, A069178, A080956, A124080, A163756, A163758, A163761, A164136, A173307.

Magma
```
[25*n*(n+1)/2+1: n in [0..50]];
```

Table[25 n (n + 1)/2 + 1, {n, 0, 50}]
25*Accumulate[Range[0,50]]+1 (* or *) LinearRecurrence[{3,-3,1},{1,26,76},50] (* Harvey P. Dale, Jan 29 2023 *)

PARI
```
vector(50, n, n--; 25*n*(n+1)/2+1)
```
Sage
```
[25*n*(n+1)/2+1 for n in (0..50)]
```

G.f.: (1 + 23*x + x^2)/(1 - x)^3.
a(n) = a(-n-1) = 3*a(n-1) - 3*a(n-2) + a(n-3).
a(n) = A123296(n) + 1.
a(n) = A000217(5*n+2) - 2.
a(n) = A034856(5*n+1).
a(n) = A186349(10*n+1).
a(n) = A054254(5*n+2) with n>0, a(0)=1.
a(n) = A000217(n+1) + 23*A000217(n) + A000217(n-1) with A000217(-1)=0.
Sum_{i>=0} 1/a(i) = 1.078209111... = 2*Pi*tan(Pi*sqrt(17)/10)/(5*sqrt(17)).
From Amiram Eldar, Jun 21 2020: (Start)
Sum_{n>=0} a(n)/n! = 77*e/2.
Sum_{n>=0} (-1)^(n+1) * a(n)/n! = 23/(2*e). (End)
E.g.f.: exp(x)*(2 + 50*x + 25*x^2)/2. - Elmo R. Oliveira, Dec 24 2024

A194715 15 times triangular numbers.

Original entry on oeis.org

0, 15, 45, 90, 150, 225, 315, 420, 540, 675, 825, 990, 1170, 1365, 1575, 1800, 2040, 2295, 2565, 2850, 3150, 3465, 3795, 4140, 4500, 4875, 5265, 5670, 6090, 6525, 6975, 7440, 7920, 8415, 8925, 9450, 9990, 10545, 11115, 11700, 12300, 12915, 13545, 14190, 14850, 15525

Offset: 0

Omar E. Pol, Oct 03 2011

Sequence found by reading the line from 0, in the direction 0, 15, ... and the same line from 0, in the direction 0, 45, ..., in the square spiral whose vertices are the generalized 17-gonal numbers.
Sum of the numbers from 7*n to 8*n. - Wesley Ivan Hurt, Dec 23 2015
Also the number of 4-cycles in the (n+6)-triangular honeycomb obtuse knight graph. - Eric W. Weisstein, Jul 28 2017

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
M. Janjic and B. Petkovic, A Counting Function, arXiv 1301.4550 [math.CO], 2013.
Eric Weisstein's World of Mathematics, Graph Cycle.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A000217, A028895, A035008, A045943, A069128, A163756.

Cf. A001105 (3-cycles in the triangular honeycomb obtuse knight graph), A290391 (5-cycles), A290392 (6-cycles). - Eric W. Weisstein, Jul 29 2017

[15*n*(n+1)/2: n in [0..50]]; // Vincenzo Librandi, Oct 04 2011

A194715:=n->15*n*(n+1)/2: seq(A194715(n), n=0..60); # Wesley Ivan Hurt, Dec 23 2015

15*Accumulate[Range[0, 60]] (* Harvey P. Dale, Feb 12 2012 *)
Table[15 n (n + 1)/2, {n, 0, 60}] (* Wesley Ivan Hurt, Dec 23 2015 *)
15 Binomial[Range[20], 2] (* Eric W. Weisstein, Jul 28 2017 *)
15 PolygonalNumber[Range[0, 20]] (* Eric W. Weisstein, Jul 28 2017 *)

a(n)=15*n*(n+1)/2 \\ Charles R Greathouse IV, Jun 17 2017

a(n) = 15*n*(n+1)/2 = 15*A000217(n) = 5*A045943(n) = 3*A028895(n) = A069128(n+1) - 1.
From Wesley Ivan Hurt, Dec 23 2015: (Start)
G.f.: 15*x/(1-x)^3.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2.
a(n) = Sum_{i=7*n..8*n} i. (End)
From Amiram Eldar, Feb 21 2023: (Start)
Sum_{n>=1} 1/a(n) = 2/15.
Sum_{n>=1} (-1)^(n+1)/a(n) = (4*log(2) - 2)/15.
Product_{n>=1} (1 - 1/a(n)) = -(15/(2*Pi))*cos(sqrt(23/15)*Pi/2).
Product_{n>=1} (1 + 1/a(n)) = (15/(2*Pi))*cos(sqrt(7/15)*Pi/2). (End)
E.g.f.: 15*exp(x)*x*(2 + x)/2. - Elmo R. Oliveira, Dec 25 2024

A264822 Centered 15-gonal (or pentadecagonal) primes.

Original entry on oeis.org

151, 421, 541, 991, 1171, 1801, 2851, 6091, 11701, 12301, 14851, 16921, 19891, 30241, 34171, 42751, 43891, 52291, 53551, 58741, 62791, 64171, 80341, 81901, 93241, 107101, 121921, 131671, 156601, 163171, 165391, 183691, 193201, 210421, 231001, 233641, 241651, 244351

Offset: 1

Ilya Gutkovskiy, Nov 26 2015

nonn

Primes of the form (15*k^2 - 15*k + 2)/2.
All the terms in this sequence are congruent to 1 (mod 10). - K. D. Bajpai, Nov 29 2015
The associated k-values are 5, 8, 9, 12, 13, 16, 20, 29, 40, 41, 45, 48, 52, 64, 68, 76, 77, 84, 85, 89, ... - Danny Rorabaugh, Jan 18 2016

K. D. Bajpai, Table of n, a(n) for n = 1..10000
OEIS Wiki, Figurate numbers
Eric Weisstein's World of Mathematics, Centered Polygonal Number and Prime Number

Cf. A000040, A069128.

[k: n in [1..10000] | IsPrime(k) where k is (15*n^2-15*n+2) div 2]; // K. D. Bajpai, Nov 29 2015

select(isprime, [seq((15*k^2 - 15*k + 2) / 2, k=0..1000)]); # K. D. Bajpai, Nov 29 2015

Select[Table[(15n^2 - 15n + 2) / 2, {n, 500}], PrimeQ] (* K. D. Bajpai, Nov 29 2015 *)

for(n=1, 1e3, if(isprime(k=(15*n^2-15*n+2)/2), print1(k,", "))) \\ Altug Alkan, Nov 26 2015

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A069190 Centered 24-gonal numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A262221 a(n) = 25*n*(n + 1)/2 + 1.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Sage

Formula

A194715 15 times triangular numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A264822 Centered 15-gonal (or pentadecagonal) primes.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

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

A262221 a(n) = 25n(n + 1)/2 + 1.