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 31 results. Next

A050441 Partial sums of A051865.

Original entry on oeis.org

0, 1, 14, 50, 120, 235, 406, 644, 960, 1365, 1870, 2486, 3224, 4095, 5110, 6280, 7616, 9129, 10830, 12730, 14840, 17171, 19734, 22540, 25600, 28925, 32526, 36414, 40600, 45095, 49910, 55056, 60544, 66385, 72590, 79170, 86136, 93499, 101270
Offset: 0

Views

Author

Barry E. Williams, Dec 23 1999

Keywords

Comments

This sequence is related to A180223 by 2*a(n) = n*A180223(n) - Sum_{i=0..n-1} A180223(i). Also, 13-gonal (or tridecagonal) pyramidal numbers. - Bruno Berselli, Dec 14 2010

Examples

			After 0, the sequence is provided by the row sums of the triangle (see above, fourth formula):
  1;
  2, 12;
  3, 24, 23;
  4, 36, 46, 34;
  5, 48, 69, 68, 45; ... - _Vincenzo Librandi_, Feb 12 2014
		

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, pp. 189-196.
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93.

Crossrefs

Similar sequences are listed in A237616.

Programs

  • GAP
    List([0..40], n-> n*(n+1)*(11*n-8)/6); # G. C. Greubel, Aug 30 2019
  • Magma
    I:=[0,1,14,50]; [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
    
  • Maple
    seq(n*(n+1)*(11*n-8)/6, n=0..40); # G. C. Greubel, Aug 30 2019
  • Mathematica
    Accumulate[Table[n (11n-9)/2,{n,0,40}]] (* or *) LinearRecurrence[ {4,-6,4,-1},{0,1,14,50},40] (* Harvey P. Dale, Nov 14 2011 *)
    CoefficientList[Series[x (1 + 10 x)/(1 - x)^4, {x, 0, 40}], x] (* Vincenzo Librandi, Feb 12 2014 *)
  • PARI
    a(n)=n*(n+1)*(11*n-8)/6 \\ Charles R Greathouse IV, Oct 07 2015
    
  • Sage
    [n*(n+1)*(11*n-8)/6 for n in (0..40)] # G. C. Greubel, Aug 30 2019
    

Formula

a(n) = n*(n+1)*(11*n-8)/6.
G.f.: x*(1+10*x)/(1-x)^4. - Bruno Berselli, Aug 19 2010
a(n) = 4*a(n-1)-6*a(n-2)+4*a(n-3)-a(n-4). - Bruno Berselli, Aug 19 2010
a(n) = Sum_{i=0..n-1} (n-i)*(11*i+1), with a(0)=0. - Bruno Berselli, Feb 10 2014
E.g.f.: exp(x)*x*(6 + 36*x + 11*x^2)/6. - Stefano Spezia, May 04 2022

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

A195313 Generalized 13-gonal numbers: m*(11*m-9)/2 with m = 0, 1, -1, 2, -2, 3, -3, ...

Original entry on oeis.org

0, 1, 10, 13, 31, 36, 63, 70, 106, 115, 160, 171, 225, 238, 301, 316, 388, 405, 486, 505, 595, 616, 715, 738, 846, 871, 988, 1015, 1141, 1170, 1305, 1336, 1480, 1513, 1666, 1701, 1863, 1900, 2071, 2110, 2290, 2331, 2520, 2563, 2761, 2806, 3013, 3060, 3276
Offset: 0

Views

Author

Omar E. Pol, Sep 14 2011

Keywords

Comments

Also generalized tridecagonal numbers or generalized triskaidecagonal numbers.
Also A211013 and positive terms of A051865 interleaved. - Omar E. Pol, Aug 04 2012
Numbers k for which 88*k + 81 is a square. - Bruno Berselli, Jul 10 2018

Crossrefs

Partial sums of A195312.
Column 9 of A195152.
Cf. A316672.
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), this sequence (k=13), A195818 (k=14), A277082 (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

  • Magma
    [(22*n*(n+1)+7*(2*n+1)*(-1)^n-7)/16: n in [0..50]]; // Vincenzo Librandi, Sep 16 2011
    
  • Magma
    A195313:=func; [0] cat [A195313(n*m): m in [1,-1], n in [1..25]]; // Bruno Berselli, Nov 13 2012
    
  • Maple
    a:= n-> (m-> m*(11*m-9)/2)(-ceil(n/2)*(-1)^n):
    seq(a(n), n=0..60);  # Alois P. Heinz, Jul 10 2018
  • Mathematica
    lim = 50; Sort[Table[n*(11*n - 9)/2, {n, -lim, lim}]] (* T. D. Noe, Sep 15 2011 *)
    Accumulate[With[{nn=30},Riffle[9Range[0,nn],Range[1,2nn+1,2]]]] (* Harvey P. Dale, Sep 24 2011 *)
  • PARI
    a(n)=(22*n*(n+1)+7*(2*n+1)*(-1)^n-7)/16 \\ Charles R Greathouse IV, Sep 24 2015

Formula

From Bruno Berselli, Sep 15 2011: (Start)
G.f.: x*(1 + 9*x + x^2)/((1 + x)^2*(1 - x)^3).
a(n) = (22*n*(n + 1) + 7*(2*n + 1)*(-1)^n - 7)/16.
a(n) - a(n-2) = A175885(n). (End)
Sum_{n>=1} 1/a(n) = 22/81 + 2*Pi*cot(2*Pi/11)/9. - 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)

A254963 a(n) = n*(11*n + 3)/2.

Original entry on oeis.org

0, 7, 25, 54, 94, 145, 207, 280, 364, 459, 565, 682, 810, 949, 1099, 1260, 1432, 1615, 1809, 2014, 2230, 2457, 2695, 2944, 3204, 3475, 3757, 4050, 4354, 4669, 4995, 5332, 5680, 6039, 6409, 6790, 7182, 7585, 7999, 8424, 8860, 9307, 9765, 10234, 10714, 11205, 11707
Offset: 0

Views

Author

Bruno Berselli, Feb 11 2015

Keywords

Comments

This sequence provides the first differences of A254407 and the partial sums of A017473.
Also:
a(n) - n = A022269(n);
a(n) + n = n*(11*n+5)/2: 0, 8, 27, 57, 98, 150, 213, 287, ...;
a(n) - 2*n = A022268(n);
a(n) + 2*n = n*(11*n+7)/2: 0, 9, 29, 60, 102, 155, 219, 294, ...;
a(n) - 3*n = n*(11*n-3)/2: 0, 4, 19, 45, 82, 130, 189, 259, ...;
a(n) + 3*n = A211013(n);
a(n) - 4*n = A226492(n);
a(n) + 4*n = A152740(n);
a(n) - 5*n = A180223(n);
a(n) + 5*n = n*(11*n+13)/2: 0, 12, 35, 69, 114, 170, 237, 315, ...;
a(n) - 6*n = A051865(n);
a(n) + 6*n = n*(11*n+15)/2: 0, 13, 37, 72, 118, 175, 243, 322, ...;
a(n) - 7*n = A152740(n-1) with A152740(-1) = 0;
a(n) + 7*n = n*(11*n+17)/2: 0, 14, 39, 75, 122, 180, 249, 329, ...;
a(n) - n*(n-1)/2 = A168668(n);
a(n) + n*(n-1)/2 = A049453(n);
a(n) - n*(n+1)/2 = A202803(n);
a(n) + n*(n+1)/2 = A033580(n).

Crossrefs

Cf. A008729 and A218530 (seventh column); A017473, A254407.
Cf. similar sequences of the type 4*n^2 + k*n*(n+1)/2: A055999 (k=-7, n>6), A028552 (k=-6, n>2), A095794 (k=-5, n>1), A046092 (k=-4, n>0), A000566 (k=-3), A049450 (k=-2), A022264 (k=-1), A016742 (k=0), A022267 (k=1), A202803 (k=2), this sequence (k=3), A033580 (k=4).
Cf. A069125: (2*n+1)^2 + 3*n*(n+1)/2; A147875: n^2 + 3*n*(n+1)/2.

Programs

  • Magma
    [n*(11*n+3)/2: n in [0..50]];
    
  • Mathematica
    Table[n (11 n + 3)/2, {n, 0, 50}]
    LinearRecurrence[{3,-3,1},{0,7,25},50] (* Harvey P. Dale, Mar 25 2018 *)
  • Maxima
    makelist(n*(11*n+3)/2, n, 0, 50);
  • PARI
    vector(50, n, n--; n*(11*n+3)/2)
    
  • Sage
    [n*(11*n+3)/2 for n in (0..50)]
    

Formula

G.f.: x*(7 + 4*x)/(1 - x)^3.
From Elmo R. Oliveira, Dec 15 2024: (Start)
E.g.f.: exp(x)*x*(14 + 11*x)/2.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2. (End)

A152734 5 times pentagonal numbers: 5*n*(3*n-1)/2.

Original entry on oeis.org

0, 5, 25, 60, 110, 175, 255, 350, 460, 585, 725, 880, 1050, 1235, 1435, 1650, 1880, 2125, 2385, 2660, 2950, 3255, 3575, 3910, 4260, 4625, 5005, 5400, 5810, 6235, 6675, 7130, 7600, 8085, 8585, 9100, 9630, 10175, 10735, 11310, 11900, 12505, 13125, 13760, 14410
Offset: 0

Views

Author

Omar E. Pol, Dec 11 2008

Keywords

Comments

a(n) can be represented as a figurate number using n concentric pentagons (see example). - Omar E. Pol, Aug 21 2011

Examples

			From _Omar E. Pol_, Aug 22 2011 (Start):
Illustration of initial terms as concentric pentagons (in a precise representation the pentagons should be strictly concentric):
.
.                                          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
.
.  5             25                        60
(End)
		

Crossrefs

Cf. sequences of the form n*(d*n+10-d)/2 indexed in A140090.

Programs

Formula

a(n) = 5*A000326(n).
a(n) = a(n-1)+15*n-10 (with a(0)=0). - Vincenzo Librandi, Nov 26 2010
G.f.: 5*x*(1+2*x)/(1-x)^3. a(n) = 4*A000217(n)+A051865(n). - Bruno Berselli, Feb 11 2011
E.g.f.: (5/2)*(3*x^2 + 2*x)*exp(x). - G. C. Greubel, Jul 17 2017
From Amiram Eldar, Feb 26 2022: (Start)
Sum_{n>=1} 1/a(n) = (9*log(3) - sqrt(3)*Pi)/15.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*(sqrt(3)*Pi- 6*log(2))/15. (End)

A226492 a(n) = n*(11*n-5)/2.

Original entry on oeis.org

0, 3, 17, 42, 78, 125, 183, 252, 332, 423, 525, 638, 762, 897, 1043, 1200, 1368, 1547, 1737, 1938, 2150, 2373, 2607, 2852, 3108, 3375, 3653, 3942, 4242, 4553, 4875, 5208, 5552, 5907, 6273, 6650, 7038, 7437, 7847, 8268, 8700, 9143, 9597, 10062, 10538, 11025, 11523
Offset: 0

Views

Author

Bruno Berselli, Jun 11 2013

Keywords

Comments

Sequences of numbers of the form n*(n*k - k + 6)/2:
. k from 0 to 10, respectively: A008585, A055998, A005563, A045943, A014105, A005475, A033428, A022264, A033991, A062741, A147874;
. k=11: a(n);
. k=12: A094159;
. k=13: 0, 3, 19, 48, 90, 145, 213, 294, 388, 495, 615, 748, 894, ...;
. k=14: 0, 3, 20, 51, 96, 155, 228, 315, 416, 531, 660, 803, 960, ...;
. k=15: A152773;
. k=16: A139272;
. k=17: 0, 3, 23, 60, 114, 185, 273, 378, 500, 639, 795, 968, ...;
. k=18: A152751;
. k=19: 0, 3, 25, 66, 126, 205, 303, 420, 556, 711, 885, 1078, ...;
. k=20: 0, 3, 26, 69, 132, 215, 318, 441, 584, 747, 930, 1133, ...;
. k=21: A152759;
. k=22: 0, 3, 28, 75, 144, 235, 348, 483, 640, 819, 1020, 1243, ...;
. k=23: 0, 3, 29, 78, 150, 245, 363, 504, 668, 855, 1065, 1298, ...;
. k=24: A152767;
. k=25: 0, 3, 31, 84, 162, 265, 393, 546, 724, 927, 1155, 1408, ...;
. k=26: 0, 3, 32, 87, 168, 275, 408, 567, 752, 963, 1200, 1463, ...;
. k=27: A153783;
. k=28: A195021;
. k=29: 0, 3, 35, 96, 186, 305, 453, 630, 836, 1071, 1335, 1628, ...;
. k=30: A153448;
. k=31: 0, 3, 37, 102, 198, 325, 483, 672, 892, 1143, 1425, 1738, ...;
. k=32: 0, 3, 38, 105, 204, 335, 498, 693, 920, 1179, 1470, 1793, ...;
. k=33: A153875.
Also:
a(n) - n = A180223(n);
a(n) + n = n*(11*n-3)/2 = 0, 4, 19, 45, 82, 130, 189, 259, ...;
a(n) - 2*n = A051865(n);
a(n) + 2*n = A022268(n);
a(n) - 3*n = A152740(n-1);
a(n) + 3*n = A022269(n);
a(n) - 4*n = n*(11*n-13)/2 = 0, -1, 9, 30, 62, 105, 159, 224, ...;
a(n) + 4*n = A254963(n);
a(n) - n*(n-1)/2 = A147874(n+1);
a(n) + n*(n-1)/2 = A094159(n) (case k=12);
a(n) - n*(n-1) = A062741(n) (see above, this is the case k=9);
a(n) + n*(n-1) = n*(13*n-7)/2 (case k=13);
a(n) - n*(n+1)/2 = A135706(n);
a(n) + n*(n+1)/2 = A033579(n);
a(n) - n*(n+1) = A051682(n);
a(n) + n*(n+1) = A186030(n);
a(n) - n^2 = A062708(n);
a(n) + n^2 = n*(13*n-5)/2 = 0, 4, 21, 51, 94, 150, 219, ..., etc.
Sum of reciprocals of a(n), for n > 0: 0.47118857003113149692081665034891...

Crossrefs

Cf. sequences in Comments lines.
First differences are in A017425.

Programs

  • Magma
    [n*(11*n-5)/2: n in [0..50]];
    
  • Magma
    I:=[0,3,17]; [n le 3 select I[n] else 3*Self(n-1)-3*Self(n-2)+Self(n-3): n in [1..46]]; // Vincenzo Librandi, Aug 18 2013
    
  • Mathematica
    Table[n (11 n - 5)/2, {n, 0, 50}]
    CoefficientList[Series[x (3 + 8 x) / (1 - x)^3, {x, 0, 45}], x] (* Vincenzo Librandi, Aug 18 2013 *)
    LinearRecurrence[{3,-3,1},{0,3,17},50] (* Harvey P. Dale, Jan 14 2019 *)
  • PARI
    a(n)=n*(11*n-5)/2 \\ Charles R Greathouse IV, Sep 24 2015

Formula

G.f.: x*(3+8*x)/(1-x)^3.
a(n) + a(-n) = A033584(n).
From Elmo R. Oliveira, Dec 27 2024: (Start)
E.g.f.: exp(x)*x*(6 + 11*x)/2.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2.
a(n) = n + A180223(n). (End)

A220212 Convolution of natural numbers (A000027) with tetradecagonal numbers (A051866).

Original entry on oeis.org

0, 1, 16, 70, 200, 455, 896, 1596, 2640, 4125, 6160, 8866, 12376, 16835, 22400, 29240, 37536, 47481, 59280, 73150, 89320, 108031, 129536, 154100, 182000, 213525, 248976, 288666, 332920, 382075, 436480, 496496, 562496, 634865, 714000, 800310, 894216, 996151
Offset: 0

Views

Author

Bruno Berselli, Dec 08 2012

Keywords

Comments

Partial sums of A172073.
Apart from 0, all terms are in A135021: a(n) = A135021(A034856(n+1)) with n>0.

Crossrefs

Cf. convolution of the natural numbers (A000027) with the k-gonal numbers (* means "except 0"):
k= 2 (A000027 ): A000292;
k= 3 (A000217 ): A000332 (after the third term);
k= 4 (A000290 ): A002415 (after the first term);
k= 5 (A000326 ): A001296;
k= 6 (A000384*): A002417;
k= 7 (A000566 ): A002418;
k= 8 (A000567*): A002419;
k= 9 (A001106*): A051740;
k=10 (A001107*): A051797;
k=11 (A051682*): A051798;
k=12 (A051624*): A051799;
k=13 (A051865*): A055268.
Cf. similar sequences with formula n*(n+1)*(n+2)*(k*n-k+2)/12 listed in A264850.

Programs

  • Magma
    A051866:=func; [&+[(n-k+1)*A051866(k): k in [0..n]]: n in [0..37]];
    
  • Magma
    I:=[0,1,16,70,200]; [n le 5 select I[n] else 5*Self(n-1)-10*Self(n-2)+10*Self(n-3)-5*Self(n-4)+Self(n-5): n in [1..50]]; // Vincenzo Librandi, Aug 18 2013
  • Mathematica
    A051866[k_] := k (6 k - 5); Table[Sum[(n - k + 1) A051866[k], {k, 0, n}], {n, 0, 37}]
    CoefficientList[Series[x (1 + 11 x) / (1 - x)^5, {x, 0, 40}], x] (* Vincenzo Librandi, Aug 18 2013 *)

Formula

G.f.: x*(1+11*x)/(1-x)^5.
a(n) = n*(n+1)*(n+2)*(3*n-2)/6.
From Amiram Eldar, Feb 15 2022: (Start)
Sum_{n>=1} 1/a(n) = 3*(3*sqrt(3)*Pi + 27*log(3) - 17)/80.
Sum_{n>=1} (-1)^(n+1)/a(n) = 3*(6*sqrt(3)*Pi - 64*log(2) + 37)/80. (End)

A153875 3 times 13-gonal (or tridecagonal) numbers: a(n) = 3*n*(11*n - 9)/2.

Original entry on oeis.org

0, 3, 39, 108, 210, 345, 513, 714, 948, 1215, 1515, 1848, 2214, 2613, 3045, 3510, 4008, 4539, 5103, 5700, 6330, 6993, 7689, 8418, 9180, 9975, 10803, 11664, 12558, 13485, 14445, 15438, 16464, 17523, 18615, 19740, 20898, 22089
Offset: 0

Views

Author

Omar E. Pol, Jan 03 2009

Keywords

Crossrefs

Cf. numbers of the form n*(n*k-k+6)/2, this sequence is the case k=33: see Comments lines of A226492.

Programs

Formula

a(n) = (33*n^2 - 27*n)/2 = A051865(n)*3.
a(n) = a(n-1) + 33*n - 30, with n>0, a(0)=0. - Vincenzo Librandi, Dec 14 2010
G.f.: 3*x*(1 + 10*x)/(1-x)^3. - Bruno Berselli, Jan 21 2011
From G. C. Greubel, Aug 31 2016: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
E.g.f.: (3/2)*x*(2 + 11*x)*exp(x). (End)

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
Showing 1-10 of 31 results. Next