cp's OEIS Frontend

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

Showing 1-10 of 15 results. Next

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

Views

Author

Omar E. Pol, Apr 27 2008

Keywords

Comments

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

Examples

			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.
		

Crossrefs

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).

Programs

  • Magma
    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
    
  • Maple
    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
  • Mathematica
    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 *)
  • Python
    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
    
  • SageMath
    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

Formula

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)

Extensions

Edited by Omar E. Pol, Jan 05 2009

A277082 Generalized 15-gonal (or pentadecagonal) numbers: n*(13*n - 11)/2, n = 0,+1,-1,+2,-2,+3,-3, ...

Original entry on oeis.org

0, 1, 12, 15, 37, 42, 75, 82, 126, 135, 190, 201, 267, 280, 357, 372, 460, 477, 576, 595, 705, 726, 847, 870, 1002, 1027, 1170, 1197, 1351, 1380, 1545, 1576, 1752, 1785, 1972, 2007, 2205, 2242, 2451, 2490, 2710, 2751, 2982, 3025, 3267, 3312, 3565, 3612, 3876, 3925, 4200, 4251, 4537, 4590, 4887, 4942
Offset: 0

Views

Author

Ilya Gutkovskiy, Sep 29 2016

Keywords

Comments

More generally, the ordinary generating function for the generalized k-gonal numbers is x*(1 + (k - 4)*x + x^2)/((1 - x)^3*(1 + x)^2). A general formula for the generalized k-gonal numbers is given by (k*(2*n^2 + 2*((-1)^n + 1)*n + (-1)^n - 1) - 2*(2*n^2 + 2*(3*(-1)^n + 1)*n + 3*((-1)^n - 1)))/16.
For k>4, Sum_{n>=1} 1/a(k,n) = 2*(k-2)/(k-4)^2 + 2*Pi*cot(2*Pi/(k-2))/(k-4). - Vaclav Kotesovec, Oct 05 2016
Numbers k for which 104*k + 121 is a square. - Bruno Berselli, Jul 10 2018
Partial sums of A317311. - Omar E. Pol, Jul 28 2018

Crossrefs

Cf. A051867 (15-gonal numbers), A316672, A317311.
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), this sequence (k=15), A274978 (k=16), A303305 (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).

Programs

  • GAP
    a:=[0,1,12,15,37];;  for n in [6..60] do a[n]:=a[n-1]+2*a[n-2]-2*a[n-3]-a[n-4]+a[n-5]; od; a; # Muniru A Asiru, Jul 10 2018
  • Mathematica
    LinearRecurrence[{1, 2, -2, -1, 1}, {0, 1, 12, 15, 37}, 56]
    Table[(26 n^2 + 26 n + 9 (-1)^n (2 n + 1) - 9)/16, {n, 0, 55}]
  • PARI
    concat(0, Vec(x*(1+11*x+x^2)/((1-x)^3*(1+x)^2) + O(x^99))) \\ Altug Alkan, Oct 01 2016
    

Formula

G.f.: x*(1 + 11*x + x^2)/((1 - x)^3*(1 + x)^2).
a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5).
a(n) = (26*n^2 + 26*n + 9*(-1)^n*(2*n+1) - 9)/16.
Sum_{n>=1} 1/a(n) = 26/121 + 2*Pi*cot(2*Pi/13)/11 = 1.3032041594895857... . - Vaclav Kotesovec, Oct 05 2016

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

Views

Author

Omar E. Pol, Apr 27 2008

Keywords

Comments

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

Examples

			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 ..............................................
========================================================
		

Crossrefs

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).

Programs

  • Magma
    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
    
  • Mathematica
    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 *)
  • SageMath
    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

Formula

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)

A254474 30-gonal numbers: a(n) = n*(14*n-13).

Original entry on oeis.org

0, 1, 30, 87, 172, 285, 426, 595, 792, 1017, 1270, 1551, 1860, 2197, 2562, 2955, 3376, 3825, 4302, 4807, 5340, 5901, 6490, 7107, 7752, 8425, 9126, 9855, 10612, 11397, 12210, 13051, 13920, 14817, 15742, 16695, 17676, 18685, 19722, 20787, 21880
Offset: 0

Views

Author

Luciano Ancora, Apr 04 2015

Keywords

Comments

See comments in A255184.
Also star 15-gonal numbers.

References

  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6 (28th row of the table).

Crossrefs

Cf. similar sequences listed in A255184.

Programs

Formula

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

A255184 25-gonal numbers: a(n) = n*(23*n-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

Views

Author

Luciano Ancora, Apr 03 2015

Keywords

Comments

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.

References

  • 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.

Crossrefs

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).

Programs

  • Magma
    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}]
  • PARI
    a(n)=n*(23*n-21)/2 \\ Charles R Greathouse IV, Oct 07 2015

Formula

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

A190991 a(n) = 13*n + 1.

Original entry on oeis.org

1, 14, 27, 40, 53, 66, 79, 92, 105, 118, 131, 144, 157, 170, 183, 196, 209, 222, 235, 248, 261, 274, 287, 300, 313, 326, 339, 352, 365, 378, 391, 404, 417, 430, 443, 456, 469, 482, 495, 508, 521, 534, 547, 560, 573, 586, 599, 612, 625, 638, 651, 664, 677
Offset: 0

Views

Author

Keywords

Comments

Partial sums give A051867. - Leo Tavares, Mar 19 2023

Crossrefs

Programs

Formula

a(n) = 2*a(n-1) - a(n-2). - Vincenzo Librandi, Jun 11 2011
From Alejandro J. Becerra Jr., Jun 04 2020: (Start)
a(n) = 13*A001477(n) + A000012(n).
G.f.: (1 + 12*x)/(1 - x)^2. (End)
E.g.f.: (1+13*x)*exp(x). - G. C. Greubel, Sep 16 2022

A317302 Square array T(n,k) = (n - 2)*(k - 1)*k/2 + k, with n >= 0, k >= 0, read by antidiagonals upwards.

Original entry on oeis.org

0, 0, 1, 0, 1, 0, 0, 1, 1, -3, 0, 1, 2, 0, -8, 0, 1, 3, 3, -2, -15, 0, 1, 4, 6, 4, -5, -24, 0, 1, 5, 9, 10, 5, -9, -35, 0, 1, 6, 12, 16, 15, 6, -14, -48, 0, 1, 7, 15, 22, 25, 21, 7, -20, -63, 0, 1, 8, 18, 28, 35, 36, 28, 8, -27, -80, 0, 1, 9, 21, 34, 45, 51, 49, 36, 9, -35, -99, 0, 1, 10, 24, 40, 55, 66
Offset: 0

Views

Author

Omar E. Pol, Aug 09 2018

Keywords

Comments

Note that the formula gives several kinds of numbers, for example:
Row 0 gives 0 together with A258837.
Row 1 gives 0 together with A080956.
Row 2 gives A001477, the nonnegative numbers.
For n >= 3, row n gives the n-gonal numbers (see Crossrefs section).

Examples

			Array begins:
------------------------------------------------------------------------
n\k  Numbers       Seq. No.   0   1   2   3   4    5    6    7    8
------------------------------------------------------------------------
0    ............ (A258837):  0,  1,  0, -3, -8, -15, -24, -35, -48, ...
1    ............ (A080956):  0,  1,  1,  0, -2,  -5,  -9, -14, -20, ...
2    Nonnegatives  A001477:   0,  1,  2,  3,  4,   5,   6,   7,   8, ...
3    Triangulars   A000217:   0,  1,  3,  6, 10,  15,  21,  28,  36, ...
4    Squares       A000290:   0,  1,  4,  9, 16,  25,  36,  49,  64, ...
5    Pentagonals   A000326:   0,  1,  5, 12, 22,  35,  51,  70,  92, ...
6    Hexagonals    A000384:   0,  1,  6, 15, 28,  45,  66,  91, 120, ...
7    Heptagonals   A000566:   0,  1,  7, 18, 34,  55,  81, 112, 148, ...
8    Octagonals    A000567:   0,  1,  8, 21, 40,  65,  96, 133, 176, ...
9    9-gonals      A001106:   0,  1,  9, 24, 46,  75, 111, 154, 204, ...
10   10-gonals     A001107:   0,  1, 10, 27, 52,  85, 126, 175, 232, ...
11   11-gonals     A051682:   0,  1, 11, 30, 58,  95, 141, 196, 260, ...
12   12-gonals     A051624:   0,  1, 12, 33, 64, 105, 156, 217, 288, ...
13   13-gonals     A051865:   0,  1, 13, 36, 70, 115, 171, 238, 316, ...
14   14-gonals     A051866:   0,  1, 14, 39, 76, 125, 186, 259, 344, ...
15   15-gonals     A051867:   0,  1, 15, 42, 82, 135, 201, 280, 372, ...
...
		

Crossrefs

Column 0 gives A000004.
Column 1 gives A000012.
Column 2 gives A001477, which coincides with the row numbers.
Main diagonal gives A060354.
Row 0 gives 0 together with A258837.
Row 1 gives 0 together with A080956.
Row 2 gives A001477, the same as column 2.
For n >= 3, row n gives the n-gonal numbers: A000217 (n=3), A000290 (n=4), A000326 (n=5), A000384 (n=6), A000566 (n=7), A000567 (n=8), A001106 (n=9), A001107 (n=10), A051682 (n=11), A051624 (n=12), A051865 (n=13), A051866 (n=14), A051867 (n=15), A051868 (n=16), A051869 (n=17), A051870 (n=18), A051871 (n=19), A051872 (n=20), A051873 (n=21), A051874 (n=22), A051875 (n=23), A051876 (n=24), A255184 (n=25), A255185 (n=26), A255186 (n=27), A161935 (n=28), A255187 (n=29), A254474 (n=30).
Cf. A303301 (similar table but with generalized polygonal numbers).

Formula

T(n,k) = A139600(n-2,k) if n >= 2.
T(n,k) = A139601(n-3,k) if n >= 3.

A098923 33-gonal numbers: n(31n-29)/2.

Original entry on oeis.org

0, 1, 33, 96, 190, 315, 471, 658, 876, 1125, 1405, 1716, 2058, 2431, 2835, 3270, 3736, 4233, 4761, 5320, 5910, 6531, 7183, 7866, 8580, 9325, 10101, 10908, 11746, 12615, 13515, 14446, 15408, 16401, 17425, 18480, 19566, 20683, 21831, 23010
Offset: 0

Views

Author

Parthasarathy Nambi, Oct 18 2004

Keywords

Comments

Similar to 21-gonal and 15-gonal numbers (A051873, A051867).

Crossrefs

Programs

Formula

a(n) = n*(31*n-29)/2.
G.f.: x*(1+30*x)/(1-x)^3. - Bruno Berselli, Feb 04 2011
a(n) = 31*n + a(n-1) - 30 (with a(0)=0). - Vincenzo Librandi, Nov 16 2010
E.g.f.: exp(x)*(x + 31*x^2/2). - Nikolaos Pantelidis, Feb 10 2023

Extensions

More terms from Stefan Steinerberger, Feb 15 2006

A177890 15-gonal (or pentadecagonal) pyramidal numbers: a(n) = n*(n+1)*(13*n-10)/6.

Original entry on oeis.org

0, 1, 16, 58, 140, 275, 476, 756, 1128, 1605, 2200, 2926, 3796, 4823, 6020, 7400, 8976, 10761, 12768, 15010, 17500, 20251, 23276, 26588, 30200, 34125, 38376, 42966, 47908, 53215, 58900, 64976, 71456, 78353, 85680, 93450, 101676, 110371, 119548, 129220
Offset: 0

Views

Author

Bruno Berselli, Dec 14 2010

Keywords

Comments

Also a(n) = (15-m)*A000292(n-1) + n*(n+1)*((m-2)*n - (m-5))/6 being n*(n+1)*((m-2)*n - (m-5))/6 a m-gonal pyramidal number (1 < m < 15). For m=6, a(n) = 9*A000292(n-1) + A002412(n).
Inverse binomial transform of this sequence: 0, 1, 14, 13, 0, 0 (0 continued).

References

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

Crossrefs

Cf. similar sequences listed in A237616.

Programs

  • GAP
    List([0..40], n-> n*(n+1)*(13*n-10)/6); # G. C. Greubel, Aug 30 2019
  • Magma
    I:=[0,1,16,58]; [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..40]]; // Vincenzo Librandi, Jul 04 2012
    
  • Magma
    [n*(n+1)*(13*n-10)/6: n in [0..40]]; // G. C. Greubel, Aug 30 2019
    
  • Maple
    seq(n*(n+1)*(13*n-10)/6, n=0..40); # G. C. Greubel, Aug 30 2019
  • Mathematica
    CoefficientList[Series[x*(1+12*x)/(1-x)^4,{x,0,40}],x] (* Vincenzo Librandi, Jul 04 2012 *)
    Table[n*(n-1)*(13*n-23)/6, {n,40}] (* G. C. Greubel, Aug 30 2019 *)
    LinearRecurrence[{4,-6,4,-1},{0,1,16,58},40] (* Harvey P. Dale, Dec 21 2022 *)
  • PARI
    vector(40, n, n*(n-1)*(13*n-23)/6) \\ G. C. Greubel, Aug 30 2019
    
  • Sage
    [n*(n+1)*(13*n-10)/6 for n in (0..40)] # G. C. Greubel, Aug 30 2019
    

Formula

G.f.: x*(1+12*x)/(1-x)^4.
a(n) = Sum_{i=0..n} A051867(i).
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Vincenzo Librandi, Jul 04 2012
a(n) = Sum_{i=0..n-1} (n-i)*(13*i+1), with a(0)=0. - Bruno Berselli, Feb 10 2014
E.g.f.: x*(6 + 42*x + 13*x^2)*exp(x)/6. - G. C. Greubel, Aug 30 2019

A256648 28-gonal pyramidal numbers: a(n) = n*(n+1)*(26*n-23)/6.

Original entry on oeis.org

0, 1, 29, 110, 270, 535, 931, 1484, 2220, 3165, 4345, 5786, 7514, 9555, 11935, 14680, 17816, 21369, 25365, 29830, 34790, 40271, 46299, 52900, 60100, 67925, 76401, 85554, 95410, 105995, 117335, 129456, 142384, 156145, 170765, 186270, 202686, 220039, 238355
Offset: 0

Views

Author

Luciano Ancora, Apr 07 2015

Keywords

Comments

See comments in A256645.
This sequence is related to A051867 by a(n) = n*A051867(n) - Sum_{i=0..n-1} A051867(i). - Bruno Berselli, Apr 09 2015

References

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

Crossrefs

Partial sums of A161935.
Cf. similar sequences listed in A237616.

Programs

  • Magma
    [n*(n+1)*(26*n-23)/6: n in [0..50]]; // Vincenzo Librandi, Apr 08 2015
  • Mathematica
    Table[n (n + 1)(26 n - 23)/6, {n, 0, 40}]
    LinearRecurrence[{4, -6, 4, -1}, {0, 1, 29, 110}, 40] (* Vincenzo Librandi, Apr 08 2015 *)

Formula

G.f.: x*(1 + 25*x)/(1 - x)^4.
a(n) = A000292(n) + 25*A000292(n-1).
From Elmo R. Oliveira, Aug 04 2025: (Start)
E.g.f.: exp(x)*x*(6 + 81*x + 26*x^2)/6.
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). (End)
Showing 1-10 of 15 results. Next