A051869 -id:A051869 - cp's OEIS frontend

A137878 Perfect squares among 17-gonal numbers A051869(k) = k(15k - 13)/2.

1, 17689, 378225, 4109707449, 87870152041, 954775454112481, 20414169462254569, 221815343046210267025, 4742660677722035990769, 51532584126226886201833161, 1101824413949324675985344641, 11972153009151467313136073526409, 255978051492792346696545201859225

Offset: 1

Alexander Adamchuk, Feb 19 2008

Corresponding square roots sqrt(a(n)) are listed in A137879.

Indices of perfect squares among the 17-gonal numbers A051869(k) = k*(15*k - 13)/2 are listed in A137880. Note that all such indices are also perfect squares, their square roots are listed in A137881(k) = sqrt(A137880(k)).

Colin Barker, Table of n, a(n) for n = 1..350
Index entries for linear recurrences with constant coefficients, signature (1,232322,-232322,-1,1).

Cf. A051869 (17-gonal numbers), A137879, A137880, A137881.

Vec(x*(1+17688*x+128214*x^2+17688*x^3+x^4)/((1-x)*(1-482*x+x^2)*(1+482*x+x^2)) + O(x^20)) \\ Colin Barker, Jun 19 2016

a(n) = A137879(n)^2 = A051869( A137880(n) ) = A051869( A137881(n)^2 ).
From Colin Barker, Jun 19 2016: (Start)
a(n) = a(n-1) + 232322*a(n-2) - 232322*a(n-3) - a(n-4) + a(n-5) for n > 5.
G.f.: x*(1 + 17688*x + 128214*x^2 + 17688*x^3 + x^4) / ((1-x)*(1 - 482*x + x^2)*(1 + 482*x + x^2)).
(End)

Edited and extended by Max Alekseyev, Oct 19 2008

A137880 Indices k of perfect squares among 17-gonal numbers A051869(k) = k(15k - 13)/2.

Original entry on oeis.org

1, 49, 225, 23409, 108241, 11282881, 52171729, 5438325025, 25146664929, 2621261378961, 12120640323841, 1263442546333969, 5842123489426225, 608976686071593889, 2815891401263116401, 293525499243961920321, 1357253813285332678849, 141478681658903574000625

Offset: 1

Alexander Adamchuk, Feb 19 2008

Corresponding perfect squares are listed in A137878.

Note that all a(n) are perfect squares themselves, their square roots are listed in A137881.

Matthew House, Table of n, a(n) for n = 1..746
Index entries for linear recurrences with constant coefficients, signature (1,482,-482,-1,1).

Cf. A051869 (17-gonal numbers), A137878 (17-gonal numbers that are perfect squares), A137879, A137881.

Rest@ CoefficientList[Series[x (1 + 48 x - 306 x^2 + 48 x^3 + x^4)/((1 - x) (1 - 22 x + x^2) (1 + 22 x + x^2)), {x, 0, 18}], x] (* Michael De Vlieger, Jun 18 2016 *)

Vec(x*(1+48*x-306*x^2+48*x^3+x^4)/((1-x)*(1-22*x+x^2)*(1+22*x+x^2)) + O(x^20)) \\ Colin Barker, Jun 18 2016

A051869( a(n) ) = A137878(n); a(n) = A137881(n)^2.
From Max Alekseyev, Oct 19 2008: (Start)
a(n) = 482*a(n-2) - a(n-4) - 208.
a(2n) = ( (15 - sqrt(30))/30 * (11 + 2*sqrt(30))^n + (15 + sqrt(30))/30 * (11 - 2*sqrt(30))^n )^2.
a(2n+1) = ( (15 + sqrt(30))/30 * (11 + 2*sqrt(30))^n + (15 - sqrt(30))/30 * (11 - 2*sqrt(30))^n )^2. (End)
a(n) = a(n-1) + 482*a(n-2) - 482*a(n-3) - a(n-4) + a(n-5). - Matthew House, Jun 18 2016
G.f.: x*(1 + 48*x - 306*x^2 + 48*x^3 + x^4) / ((1-x)*(1 - 22*x + x^2)*(1 + 22*x + x^2)). - Colin Barker, Jun 18 2016

Edited and extended by Max Alekseyev, Oct 19 2008

A139600 Square array T(n,k) = n(k-1)k/2+k, of nonnegative numbers together with polygonal numbers, read by antidiagonals upwards.

Original entry on oeis.org

0, 0, 1, 0, 1, 2, 0, 1, 3, 3, 0, 1, 4, 6, 4, 0, 1, 5, 9, 10, 5, 0, 1, 6, 12, 16, 15, 6, 0, 1, 7, 15, 22, 25, 21, 7, 0, 1, 8, 18, 28, 35, 36, 28, 8, 0, 1, 9, 21, 34, 45, 51, 49, 36, 9, 0, 1, 10, 24, 40, 55, 66, 70, 64, 45, 10, 0, 1, 11, 27, 46, 65, 81, 91, 92, 81, 55, 11

Offset: 0

Omar E. Pol, Apr 27 2008

A general formula for polygonal numbers is P(n,k) = (n-2)*(k-1)*k/2 + k, where P(n,k) is the k-th n-gonal number.

The triangle sums, see A180662 for their definitions, link this square array read by antidiagonals with twelve different sequences, see the crossrefs. Most triangle sums are linear sums of shifted combinations of a sequence, see e.g. A189374. - Johannes W. Meijer, Apr 29 2011

			The square array of nonnegatives together with polygonal numbers begins:
=========================================================
....................... A   A   .   .   A    A    A    A
....................... 0   0   .   .   0    0    1    1
....................... 0   0   .   .   1    1    3    3
....................... 0   0   .   .   6    7    9    9
....................... 0   0   .   .   9    3    6    6
....................... 0   1   .   .   5    2    0    0
....................... 4   2   .   .   7    9    6    7
=========================================================
Nonnegatives . A001477: 0,  1,  2,  3,  4,   5,   6,   7, ...
Triangulars .. A000217: 0,  1,  3,  6, 10,  15,  21,  28, ...
Squares ...... A000290: 0,  1,  4,  9, 16,  25,  36,  49, ...
Pentagonals .. A000326: 0,  1,  5, 12, 22,  35,  51,  70, ...
Hexagonals ... A000384: 0,  1,  6, 15, 28,  45,  66,  91, ...
Heptagonals .. A000566: 0,  1,  7, 18, 34,  55,  81, 112, ...
Octagonals ... A000567: 0,  1,  8, 21, 40,  65,  96, 133, ...
9-gonals ..... A001106: 0,  1,  9, 24, 46,  75, 111, 154, ...
10-gonals .... A001107: 0,  1, 10, 27, 52,  85, 126, 175, ...
11-gonals .... A051682: 0,  1, 11, 30, 58,  95, 141, 196, ...
12-gonals .... A051624: 0,  1, 12, 33, 64, 105, 156, 217, ...
...
=========================================================
The column with the numbers 2, 3, 4, 5, 6, ... is formed by the numbers > 1 of A000027. The column with the numbers 3, 6, 9, 12, 15, ... is formed by the positive members of A008585.

Alois P. Heinz, Rows n = 0..140, flattened
Peter Luschny, Figurate number — a very short introduction. With plots from Stefan Friedrich Birkner.
Omar E. Pol, Polygonal numbers, An alternative illustration of initial terms.
Index to sequences related to polygonal numbers

Cf. A001477, A057145, A139601, A139617, A139618, A139620, A180662.
A formal extension negative n is in A326728.
Triangle sums (see the comments): A055795 (Row1), A080956 (Row2; terms doubled), A096338 (Kn11, Kn12, Kn13, Fi1, Ze1), A002624 (Kn21, Kn22, Kn23, Fi2, Ze2), A000332 (Kn3, Ca3, Gi3), A134393 (Kn4), A189374 (Ca1, Ze3), A011779 (Ca2, Ze4), A101357 (Ca4), A189375 (Gi1), A189376 (Gi2), A006484 (Gi4). - Johannes W. Meijer, Apr 29 2011
Sequences of m-gonal numbers: A000217 (m=3), A000290 (m=4), A000326 (m=5), A000384 (m=6), A000566 (m=7), A000567 (m=8), A001106 (m=9), A001107 (m=10), A051682 (m=11), A051624 (m=12), A051865 (m=13), A051866 (m=14), A051867 (m=15), A051868 (m=16), A051869 (m=17), A051870 (m=18), A051871 (m=19), A051872 (m=20), A051873 (m=21), A051874 (m=22), A051875 (m=23), A051876 (m=24), A255184 (m=25), A255185 (m=26), A255186 (m=27), A161935 (m=28), A255187 (m=29), A254474 (m=30).

T:= func< n,k | k*(n*(k-1)+2)/2 >;
A139600:= func< n,k | T(n-k, k) >;
[A139600(n,k): k in  [0..n], n in [0..12]]; // G. C. Greubel, Jul 12 2024

T:= (n, k)-> n*(k-1)*k/2+k:
seq(seq(T(d-k, k), k=0..d), d=0..14);  # Alois P. Heinz, Oct 14 2018

T[n_, k_] := (n + 1)*(k - 1)*k/2 + k; Table[T[n - k - 1, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Robert G. Wilson v, Jul 12 2009 *)

def A139600Row(n):
    x, y = 1, 1
    yield 0
    while True:
        yield x
        x, y = x + y + n, y + n
for n in range(8):
    R = A139600Row(n)
    print([next(R) for  in range(11)]) # _Peter Luschny, Aug 04 2019

def T(n,k): return k*(n*(k-1)+2)/2
def A139600(n,k): return T(n-k, k)
flatten([[A139600(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Jul 12 2024

T(n,k) = n*(k-1)*k/2+k.
T(n,k) = A057145(n+2,k). - R. J. Mathar, Jul 28 2016
From Stefano Spezia, Apr 12 2024: (Start)
G.f.: y*(1 - x - y + 2*x*y)/((1 - x)^2*(1 - y)^3).
E.g.f.: exp(x+y)*y*(2 + x*y)/2. (End)

Edited by Omar E. Pol, Jan 05 2009

A303305 Generalized 17-gonal (or heptadecagonal) numbers: m(15m - 13)/2 with m = 0, +1, -1, +2, -2, +3, -3, ...

Original entry on oeis.org

0, 1, 14, 17, 43, 48, 87, 94, 146, 155, 220, 231, 309, 322, 413, 428, 532, 549, 666, 685, 815, 836, 979, 1002, 1158, 1183, 1352, 1379, 1561, 1590, 1785, 1816, 2024, 2057, 2278, 2313, 2547, 2584, 2831, 2870, 3130, 3171, 3444, 3487, 3773, 3818, 4117, 4164, 4476, 4525, 4850

Offset: 0

Omar E. Pol, Jun 06 2018

120*a(n) + 169 is a square. - Bruno Berselli, Jun 08 2018
Partial sums of A317313. - Omar E. Pol, Jul 28 2018
Generalized k-gonal numbers are second k-gonal numbers and positive terms of k-gonal numbers interleaved, k >= 5. They are also the partial sums of the sequence formed by the multiples of (k - 4) and the odd numbers (A005408) interleaved, k >= 5. In this case k = 17. - Omar E. Pol, Apr 25 2021

			From _Omar E. Pol_, Apr 24 2021: (Start)
Illustration of initial terms as vertices of a rectangular spiral:
        43_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _17
         |                                                   |
         |                         0                         |
         |                         |_ _ _ _ _ _ _ _ _ _ _ _ _|
         |                         1                         14
         |
        48
More generally, all generalized k-gonal numbers can be represented with this kind of spirals, k >= 5". (End)

Colin Barker, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (1,2,-2,-1,1).

Cf. A051869, A194715, A244636, A317313.

Sequences of generalized k-gonal numbers: A001318 (k=5), A000217 (k=6), A085787 (k=7), A001082 (k=8), A118277 (k=9), A074377 (k=10), A195160 (k=11), A195162 (k=12), A195313 (k=13), A195818 (k=14), A277082 (k=15), A274978 (k=16), this sequence (k=17), A274979 (k=18), A303813 (k=19), A218864 (k=20), A303298 (k=21), A303299 (k=22), A303303 (k=23), A303814 (k=24), A303304 (k=25), A316724 (k=26), A316725 (k=27), A303812 (k=28), A303815 (k=29), A316729 (k=30).

With[{pp = 17, nn = 55}, {0}~Join~Riffle[Array[PolygonalNumber[pp, #] &, Ceiling[nn/2]], Array[PolygonalNumber[pp, -#] &, Ceiling[nn/2]]]] (* Michael De Vlieger, Jun 06 2018 *)
Table[(30 n (n + 1) + 11 (2 n + 1) (-1)^n - 11)/16, {n, 0, 60}] (* Bruno Berselli, Jun 08 2018 *)
CoefficientList[ Series[-x (x^2 + 13x + 1)/((x - 1)^3 (x + 1)^2), {x, 0, 50}], x] (* or *)
LinearRecurrence[{1, 2, -2, -1, 1}, {0, 1, 14, 17, 43}, 51] (* Robert G. Wilson v, Jul 28 2018 *)

concat(0, Vec(x*(1 + 13*x + x^2)/((1 + x)^2*(1 - x)^3) + O(x^40))) \\ Colin Barker, Jun 12 2018

From Bruno Berselli, Jun 08 2018: (Start)
G.f.: x*(1 + 13*x + x^2)/((1 + x)^2*(1 - x)^3).
a(n) = a(-n-1) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5).
a(n) = (30*n*(n + 1) + 11*(2*n + 1)*(-1)^n - 11)/16. Therefore:
a(n) = n*(15*n + 26)/8, if n is even, or (n + 1)*(15*n - 11)/8 otherwise.
2*(2*n - 1)*a(n) + 2*(2*n + 1)*a(n-1) - n*(15*n^2 - 13) = 0. (End)

A237616 a(n) = n(n + 1)(5*n - 4)/2.

Original entry on oeis.org

0, 1, 18, 66, 160, 315, 546, 868, 1296, 1845, 2530, 3366, 4368, 5551, 6930, 8520, 10336, 12393, 14706, 17290, 20160, 23331, 26818, 30636, 34800, 39325, 44226, 49518, 55216, 61335, 67890, 74896, 82368, 90321, 98770, 107730, 117216, 127243, 137826, 148980, 160720

Offset: 0

Bruno Berselli, Feb 10 2014

Also 17-gonal (or heptadecagonal) pyramidal numbers.

This sequence is related to A226489 by 2*a(n) = n*A226489(n) - Sum_{i=0..n-1} A226489(i).

			After 0, the sequence is provided by the row sums of the triangle:
   1;
   2,  16;
   3,  32,  31;
   4,  48,  62,  46;
   5,  64,  93,  92,  61;
   6,  80, 124, 138, 122,  76;
   7,  96, 155, 184, 183, 152,  91;
   8, 112, 186, 230, 244, 228, 182, 106;
   9, 128, 217, 276, 305, 304, 273, 212, 121;
  10, 144, 248, 322, 366, 380, 364, 318, 242, 136; etc.,
where (r = row index, c = column index):
T(r,r) = T(c,c) = 15*r-14 and T(r,c) = T(r-1,c)+T(r,r) = (r-c+1)*T(r,r), with r>=c>0.

E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93 (fifteenth row of the table).

Bruno Berselli, Table of n, a(n) for n = 0..1000
Eric Weisstein's World of Mathematics, Pyramidal Number.
Index to sequences related to pyramidal numbers.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Cf. A051869, A104728.

Cf. sequences with formula n*(n+1)*(k*n-k+3)/6: A000217 (k=0), A000292 (k=1), A000330 (k=2), A002411 (k=3), A002412 (k=4), A002413 (k=5), A002414 (k=6), A007584 (k=7), A007585 (k=8), A007586 (k=9), A007587 (k=10), A050441 (k=11), A172073 (k=12), A177890 (k=13), A172076 (k=14), this sequence (k=15), A172078(k=16), A237617 (k=17), A172082 (k=18), A237618 (k=19), A172117(k=20), A256718 (k=21), A256716 (k=22), A256645 (k=23), A256646(k=24), A256647 (k=25), A256648 (k=26), A256649 (k=27), A256650(k=28).

List([0..40], n-> n*(n+1)*(5*n-4)/2); # G. C. Greubel, Aug 30 2019

Magma
```
[n*(n+1)*(5*n-4)/2: n in [0..40]];
```

I:=[0,1,18,66]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..50]]; // Vincenzo Librandi, Feb 12 2014

seq(n*(n+1)*(5*n-4)/2, n=0..40); # G. C. Greubel, Aug 30 2019

Table[n(n+1)(5n-4)/2, {n, 0, 40}]
CoefficientList[Series[x (1+14x)/(1-x)^4, {x, 0, 40}], x] (* Vincenzo Librandi, Feb 12 2014 *)
LinearRecurrence[{4,-6,4,-1},{0,1,18,66},50] (* Harvey P. Dale, Jan 11 2015 *)

a(n)=n*(n+1)*(5*n-4)/2 \\ Charles R Greathouse IV, Sep 24 2015

[n*(n+1)*(5*n-4)/2 for n in (0..40)] # G. C. Greubel, Aug 30 2019

G.f.: x*(1 + 14*x)/(1 - x)^4.
For n>0, a(n) = Sum_{i=0..n-1} (n-i)*(15*i+1). More generally, the sequence with the closed form n*(n+1)*(k*n-k+3)/6 is also given by Sum_{i=0..n-1} (n-i)*(k*i+1) for n>0.
a(n) = A104728(A001844(n-1)) for n>0.
Sum_{n>=1} 1/a(n) = (2*sqrt(5*(5 + 2*sqrt(5)))*Pi + 10*sqrt(5)*arccoth(sqrt(5)) + 25*log(5) - 16)/72 = 1.086617842136293176... . - Vaclav Kotesovec, Dec 07 2016
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) for n >= 4. - Wesley Ivan Hurt, Dec 18 2020
E.g.f.: exp(x)*x*(2 + 16*x + 5*x^2)/2. - Elmo R. Oliveira, Aug 04 2025

A139601 Square array of polygonal numbers read by ascending antidiagonals: T(n, k) = (n + 1)(k - 1)k/2 + k.

Original entry on oeis.org

0, 0, 1, 0, 1, 3, 0, 1, 4, 6, 0, 1, 5, 9, 10, 0, 1, 6, 12, 16, 15, 0, 1, 7, 15, 22, 25, 21, 0, 1, 8, 18, 28, 35, 36, 28, 0, 1, 9, 21, 34, 45, 51, 49, 36, 0, 1, 10, 24, 40, 55, 66, 70, 64, 45, 0, 1, 11, 27, 46, 65, 81, 91, 92, 81, 55, 0, 1, 12, 30, 52, 75, 96, 112, 120, 117, 100, 66

Offset: 0

Omar E. Pol, Apr 27 2008

A general formula for polygonal numbers is P(n,k) = (n-2)(k-1)k/2 + k, where P(n,k) is the k-th n-gonal number. - Omar E. Pol, Dec 21 2008

			The square array of polygonal numbers begins:
========================================================
Triangulars .. A000217: 0, 1,  3,  6, 10,  15,  21,  28,
Squares ...... A000290: 0, 1,  4,  9, 16,  25,  36,  49,
Pentagonals .. A000326: 0, 1,  5, 12, 22,  35,  51,  70,
Hexagonals ... A000384: 0, 1,  6, 15, 28,  45,  66,  91,
Heptagonals .. A000566: 0, 1,  7, 18, 34,  55,  81, 112,
Octagonals ... A000567: 0, 1,  8, 21, 40,  65,  96, 133,
9-gonals ..... A001106: 0, 1,  9, 24, 46,  75, 111, 154,
10-gonals .... A001107: 0, 1, 10, 27, 52,  85, 126, 175,
11-gonals .... A051682: 0, 1, 11, 30, 58,  95, 141, 196,
12-gonals .... A051624: 0, 1, 12, 33, 64, 105, 156, 217,
And so on ..............................................
========================================================

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
Peter Luschny, Figurate number — a very short introduction. With plots from Stefan Friedrich Birkner.
Omar E. Pol, Polygonal numbers, An alternative illustration of initial terms.

Sequences of m-gonal numbers: A000217 (m=3), A000290 (m=4), A000326 (m=5), A000384 (m=6), A000566 (m=7), A000567 (m=8), A001106 (m=9), A001107 (m=10), A051682 (m=11), A051624 (m=12), A051865 (m=13), A051866 (m=14), A051867 (m=15), A051868 (m=16), A051869 (m=17), A051870 (m=18), A051871 (m=19), A051872 (m=20), A051873 (m=21), A051874 (m=22), A051875 (m=23), A051876 (m=24), A255184 (m=25), A255185 (m=26), A255186 (m=27), A161935 (m=28), A255187 (m=29), A254474 (m=30).
Cf. A000007, A000012, A000027, A002411, A006000, A006003, A006522.
Cf. A008585, A016957, A017329, A086270, A117142, A139606, A139607.
Cf. A139608, A139609, A139610, A139611, A139612, A139613, A139614.
Cf. A139615, A139616, A057145, A086271, A139600, A139617, A139618.
Cf. A139619, A139620.

T:= func< n,k | k*((n+1)*(k-1) +2)/2 >;
A139601:= func< n,k | T(n-k, k) >;
[A139601(n,k): k in  [0..n], n in [0..12]]; // G. C. Greubel, Jul 12 2024

T[n_, k_] := (n + 1)*(k - 1)*k/2 + k; Table[ T[n - k, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Robert G. Wilson v, Jul 12 2009 *)

def T(n,k): return k*((n+1)*(k-1)+2)/2
def A139601(n,k): return T(n-k, k)
flatten([[A139601(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Jul 12 2024

T(n,k) = A086270(n,k), k>0. - R. J. Mathar, Aug 06 2008
T(n,k) = (n+1)*(k-1)*k/2 +k, n>=0, k>=0. - Omar E. Pol, Jan 07 2009
From G. C. Greubel, Jul 12 2024: (Start)
t(n, k) = (k/2)*( (k-1)*(n-k+1) + 2), where t(n,k) is this array read by rising antidiagonals.
t(2*n, n) = A006003(n).
t(2*n+1, n) = A002411(n).
t(2*n-1, n) = A006000(n-1).
Sum_{k=0..n} t(n, k) = A006522(n+2).
Sum_{k=0..n} (-1)^k*t(n, k) = (-1)^n * A117142(n).
Sum_{k=0..n} t(n-k, k) = (2*n^4 + 34*n^2 + 48*n - 15 + 3*(-1)^n*(2*n^2 + 16*n + 5))/384. (End)

A255184 25-gonal numbers: a(n) = n(23n-21)/2.

Original entry on oeis.org

0, 1, 25, 72, 142, 235, 351, 490, 652, 837, 1045, 1276, 1530, 1807, 2107, 2430, 2776, 3145, 3537, 3952, 4390, 4851, 5335, 5842, 6372, 6925, 7501, 8100, 8722, 9367, 10035, 10726, 11440, 12177, 12937, 13720, 14526, 15355, 16207, 17082, 17980

Offset: 0

Luciano Ancora, Apr 03 2015

If b(n,k) = n*((k-2)*n-(k-4))/2 is n-th k-gonal number, then b(n,k) = A000217(n) + (k-3)* A000217(n-1) (see Deza in References section, page 21, where the formula is attributed to Bachet de Méziriac).
Also, b(n,k) = b(n,k-1) + A000217(n-1) (see Deza and Picutti in References section, page 20 and 137 respectively, where the formula is attributed to Nicomachus). Some examples:
for k=4, A000290(n) = A000217(n) + A000217(n-1);
for k=5, A000326(n) = A000290(n) + A000217(n-1);
for k=6, A000384(n) = A000326(n) + A000217(n-1), etc.
This is the case k=25.

E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6 (23rd row of the table).
E. Picutti, Sul numero e la sua storia, Feltrinelli Economica (1977), pages 131-147.

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

Cf. k-gonal numbers: A000217 (k=3), A000290 (k=4), A000326 (k=5), A000384 (k=6), A000566 (k=7), A000567 (k=8), A001106 (k=9), A001107 (k=10), A051682 (k=11), A051624 (k=12), A051865 (k=13), A051866 (k=14), A051867 (k=15), A051868 (k=16), A051869 (k=17), A051870 (k=18), A051871 (k=19), A051872 (k=20), A051873 (k=21), A051874 (k=22), A051875 (k=23), A051876 (k=24), this sequence (k=25), A255185 (k=26), A255186 (k=27), A161935 (k=28), A255187 (k=29), A254474 (k=30).

k:=25; [n*((k-2)*n-(k-4))/2: n in [0..40]]; // Bruno Berselli, Apr 10 2015

Mathematica
```
Table[n (23 n - 21)/2, {n, 40}]
```

a(n)=n*(23*n-21)/2 \\ Charles R Greathouse IV, Oct 07 2015

G.f.: x*(-1 - 22*x)/(-1 + x)^3.
a(n) = A000217(n) + 22*A000217(n-1) = A051876(n) + A000217(n-1), see comments.
Product_{n>=2} (1 - 1/a(n)) = 23/25. - Amiram Eldar, Jan 22 2021
E.g.f.: exp(x)*(x + 23*x^2/2). - Nikolaos Pantelidis, Feb 05 2023

A316672 Numbers k for which 120*k + 169 is a square.

Original entry on oeis.org

-1, 0, 1, 3, 10, 14, 17, 22, 36, 43, 48, 56, 77, 87, 94, 105, 133, 146, 155, 169, 204, 220, 231, 248, 290, 309, 322, 342, 391, 413, 428, 451, 507, 532, 549, 575, 638, 666, 685, 714, 784, 815, 836, 868, 945, 979, 1002, 1037, 1121, 1158, 1183, 1221, 1312, 1352, 1379, 1420

Offset: 1

Bruno Berselli, Jul 10 2018

All terms of A303305 belong to this sequence.

Colin Barker, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (1,0,0,2,-2,0,0,-1,1).

Cf. A051869, A303305.
Subsequence of A047283.
Cf. Numbers k for which 8*(2*h+1)*k + (2*h-1)^2 is a square: A000217 (h=0), A001318 (h=1), A085787 (h=2), A118277 (h=3), A195160 (h=4), A195313 (h=5), A277082 (h=6), this sequence (h=7), A303813 (h=8), A303298 (h=9); A303815 (h=13).

[k: k in [0..1500] | IsSquare(120*k+169)];

select(k->issqr(120*k+169),[$-1..1500]); # Muniru A Asiru, Jul 10 2018

LinearRecurrence[{1, 0, 0, 2, -2, 0, 0, -1, 1}, {-1, 0, 1, 3, 10, 14, 17, 22, 36}, 60]

isok(n) = issquare(120*n+169); \\ Michel Marcus, Jul 11 2018

Vec(x*(-1 + x + x^2 + 2*x^3 + 9*x^4 + 2*x^5 + x^6 + x^7 - x^8)/((1 + x)^2*(1 - x)^3*(1 + x^2)^2) + O(x^40)) \\ Colin Barker, Jul 18 2018

print([k for k in (0..1500) if is_square(120*k+169)])

O.g.f.: x*(-1 + x + x^2 + 2*x^3 + 9*x^4 + 2*x^5 + x^6 + x^7 - x^8)/((1 + x)^2*(1 - x)^3*(1 + x^2)^2).
a(n) = a(1-n) = a(n-1) + 2*a(n-4) - 2*a(n-5) - a(n-8) + a(n-9).
a(n) = (30*n^2 - 2*(15 + 3*(-1)^n + 10*i^(n*(n+1)))*n + 2*(5 + (-1)^n)*i^(n*(n+1)) + 3*(-1)^n - 79)/64, with i = sqrt(-1). Therefore:
a(4*k+1) = (3*k + 2)*(5*k - 1)/2;
a(4*k+2) = k*(15*k + 13)/2, first bisection of A303305;
a(4*k+3) = (k + 1)*(15*k + 2)/2, second bisection of A303305 (see A051869);
a(4*k+4) = (3*k + 1)*(5*k + 6)/2.

A261696 a(n) is least number > 0 such that the concatenation of a(1) ... a(n) is 17-gonal: (15n^2 - 13n)/2.

Original entry on oeis.org

1, 7, 689, 6797, 67984832, 6798483348333332, 8455610150480042707742277762479, 707328322040172689545426423113211907561874137758547957769721082

Offset: 1

Anders Hellström, Nov 26 2015

From Chai Wah Wu, Mar 16 2018: (Start)
There are some interesting patterns observed in the terms. Terms a(5), a(6), a(9), a(10), a(11), a(12), ... share the same prefix of 6798483...
From terms a(n) for n > 5, there seems to a pattern of how they are constructed from previous terms. a(6) is formed by inserting 3483...3 between the penultimate digit and the last digit of a(5). Then a(7) and (8) do not follow this pattern.
The digits of a(9) and a(6) match until the last digit of a(6). Next, a(10), a(11) and (12) are formed from a(9), a(10) and a(11) resp. by inserting 3483...3. Then this pattern is interrupted by a(13) and a(14), and continue again for a(15) ..., etc.
(End)

			1, 17, 17689, 176896797 are 17-gonal.

Chai Wah Wu, Table of n, a(n) for n = 1..11
Wikipedia, Polygonal number
Chai Wah Wu, terms < 10^70000 to illustrate the pattern observed (see comments)

Cf. A051671, A051869 (17-gonal numbers), A061109, A061110, A264733, A264738, A264776, A264777, A264842, A264848, A264849, A264804.

heptadecagonal(n)=ispolygonal(n, 17)
first(m)=my(s=""); s="1"; print1(1, ", ");for(i=2, m, n=1; while(!heptadecagonal(eval(concat(s, Str(n)))), n++); print1(n, ", "); s=concat(s, Str(n)))

a(6)-a(8) from Chai Wah Wu, Mar 16 2018

A137881 a(n) = sqrt(A137880(n)).

Original entry on oeis.org

1, 7, 15, 153, 329, 3359, 7223, 73745, 158577, 1619031, 3481471, 35544937, 76433785, 780369583, 1678061799, 17132585889, 36840925793, 376136519975, 808822305647, 8257870853561, 17757249798441, 181297022258367, 389850673260055, 3980276618830513, 8558957561922769

Offset: 1

Alexander Adamchuk, Feb 19 2008

A137880 gives the indices m (= a(n)^2) of perfect squares in 17-gonal numbers A051869(m) = m(15m -13)/2. Corresponding 17-gonal numbers are listed in A137878(n) = A051869( a(n)^2 ).

Positive values of x (or y) satisfying x^2 - 22xy + y^2 + 104 = 0. - Colin Barker, Feb 19 2014

Vincenzo Librandi, Table of n, a(n) for n = 1..200
Index entries for linear recurrences with constant coefficients, signature (0,22,0,-1).

Cf. A051869 (17-gonal numbers), A137878 (17-gonal numbers that are perfect squares), A137879, A137880.

I:=[1,7,15,153]; [n le 4 select I[n] else 22*Self(n-2)-Self(n-4): n in [1..30]]; // Vincenzo Librandi, Feb 21 2014

CoefficientList[Series[(1 - x) (x^2 + 8 x + 1)/(x^4 - 22 x^2 + 1), {x, 0, 40}], x] (* Vincenzo Librandi, Feb 21 2014 *)

a(n) = sqrt(A137880(n)). A051869( a(n)^2 ) = A137878(n).
For n>=5, a(n) = 22*a(n-2) - a(n-4). [Alekseyev]
a(2n) = (15 - sqrt(30))/30 * (11 + 2*sqrt(30))^n + (15 + sqrt(30))/30 * (11 - 2*sqrt(30))^n. [Alekseyev]
a(2n+1) = (15 + sqrt(30))/30 * (11 + 2*sqrt(30))^n + (15 - sqrt(30))/30 * (11 - 2*sqrt(30))^n. [Alekseyev]
G.f.: -x*(x-1)*(x^2+8*x+1) / (x^4-22*x^2+1). - Colin Barker, Feb 19 2014

Edited and extended by Max Alekseyev, Oct 19 2008

cp's OEIS Frontend

A137878 Perfect squares among 17-gonal numbers A051869(k) = k*(15*k - 13)/2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

Extensions

A137880 Indices k of perfect squares among 17-gonal numbers A051869(k) = k*(15*k - 13)/2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A139600 Square array T(n,k) = n*(k-1)*k/2+k, of nonnegative numbers together with polygonal numbers, read by antidiagonals upwards.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Python

SageMath

Formula

Extensions

A303305 Generalized 17-gonal (or heptadecagonal) numbers: m*(15*m - 13)/2 with m = 0, +1, -1, +2, -2, +3, -3, ...

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A237616 a(n) = n*(n + 1)*(5*n - 4)/2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Magma

Magma

Maple

Mathematica

PARI

Sage

Formula

A139601 Square array of polygonal numbers read by ascending antidiagonals: T(n, k) = (n + 1)*(k - 1)*k/2 + k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

A137878 Perfect squares among 17-gonal numbers A051869(k) = k(15k - 13)/2.

A137880 Indices k of perfect squares among 17-gonal numbers A051869(k) = k(15k - 13)/2.

A139600 Square array T(n,k) = n(k-1)k/2+k, of nonnegative numbers together with polygonal numbers, read by antidiagonals upwards.

A303305 Generalized 17-gonal (or heptadecagonal) numbers: m(15m - 13)/2 with m = 0, +1, -1, +2, -2, +3, -3, ...

A237616 a(n) = n(n + 1)(5*n - 4)/2.

A139601 Square array of polygonal numbers read by ascending antidiagonals: T(n, k) = (n + 1)(k - 1)k/2 + k.

A255184 25-gonal numbers: a(n) = n(23n-21)/2.