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.

Previous Showing 11-20 of 69 results. Next

A280072 Indices of centered 11-gonal numbers (A060544) that are also 11-gonal numbers (A051682).

Original entry on oeis.org

1, 11, 210, 4180, 83381, 1663431, 33185230, 662041160, 13207637961, 263490718051, 5256606723050, 104868643742940, 2092116268135741, 41737456718971871, 832657018111301670, 16611402905507061520, 331395401092029928721, 6611296618935091512891
Offset: 1

Views

Author

Colin Barker, Dec 25 2016

Keywords

Comments

Also positive integers y in the solutions to 9*x^2 - 11*y^2 - 7*x + 11*y - 2 = 0, the corresponding values of x being A280071.

Examples

			11 is in the sequence because the 11th centered 11-gonal number is 606, which is also the 12th 11-gonal number.
		

Crossrefs

Programs

  • Mathematica
    LinearRecurrence[{21,-21,1},{1,11,210},20] (* Harvey P. Dale, Aug 19 2020 *)
  • PARI
    Vec(x*(1 - 10*x) / ((1 - x)*(1 - 20*x + x^2)) + O(x^30)) \\ Colin Barker, Dec 25 2016

Formula

a(n) = (6 - (3+sqrt(11))*(10+3*sqrt(11))^(-n) + (-3+sqrt(11))*(10+3*sqrt(11))^n)/12.
a(n) = 21*a(n-1) - 21*a(n-2) + a(n-3) for n>3.
G.f.: x*(1 - 10*x) / ((1 - x)*(1 - 20*x + x^2)).

A344377 Numbers that are both 11-gonal numbers (A051682) and 11-gonal pyramidal numbers (A007586).

Original entry on oeis.org

0, 1, 23725, 1519937678700, 7248070597636
Offset: 1

Views

Author

Seiichi Manyama, May 17 2021

Keywords

Comments

Intersection of A051682 and A007586.

Crossrefs

Programs

  • PARI
    for(k=0, 1e5, if(ispolygonal(m=k*(k+1)*(3*k-2)/2, 11), print1(m", ")))

A060544 Centered 9-gonal (also known as nonagonal or enneagonal) numbers. Every third triangular number, starting with a(1)=1.

Original entry on oeis.org

1, 10, 28, 55, 91, 136, 190, 253, 325, 406, 496, 595, 703, 820, 946, 1081, 1225, 1378, 1540, 1711, 1891, 2080, 2278, 2485, 2701, 2926, 3160, 3403, 3655, 3916, 4186, 4465, 4753, 5050, 5356, 5671, 5995, 6328, 6670, 7021, 7381, 7750, 8128, 8515, 8911, 9316
Offset: 1

Views

Author

Henry Bottomley, Apr 02 2001

Keywords

Comments

Triangular numbers not == 0 (mod 3). - Amarnath Murthy, Nov 13 2005
Shallow diagonal of triangular spiral in A051682. - Paul Barry, Mar 15 2003
Equals the triangular numbers convolved with [1, 7, 1, 0, 0, 0, ...]. - Gary W. Adamson & Alexander R. Povolotsky, May 29 2009
a(n) is congruent to 1 (mod 9) for all n. The sequence of digital roots of the a(n) is A000012(n). The sequence of units' digits of the a(n) is period 20: repeat [1, 0, 8, 5, 1, 6, 0, 3, 5, 6, 6, 5, 3, 0, 6, 1, 5, 8, 0, 1]. - Ant King, Jun 18 2012
Divide each side of any triangle ABC with area (ABC) into 2n + 1 equal segments by 2n points: A_1, A_2, ..., A_(2n) on side a, and similarly for sides b and c. If the hexagon with area (Hex(n)) delimited by AA_n, AA_(n+1), BB_n, BB_(n+1), CC_n and CC_(n+1) cevians, we have a(n+1) = (ABC)/(Hex(n)) for n >= 1, (see link with java applet). - Ignacio Larrosa Cañestro, Jan 02 2015; edited by Wolfdieter Lang, Jan 30 2015
For the case n = 1 see the link for Marion's Theorem (actually Marion Walter's Theorem, see the Cugo et al, reference). Also, the generalization considered here has been called there (Ryan) Morgan's Theorem. - Wolfdieter Lang, Jan 30 2015
Pollock states that every number is the sum of at most 11 terms of this sequence, but note that "1, 10, 28, 35, &c." has a typo (35 should be 55). - Michel Marcus, Nov 04 2017
a(n) is also the number of (nontrivial) paths as well as the Wiener sum index of the (n-1)-alkane graph. - Eric W. Weisstein, Jul 15 2021

Crossrefs

Programs

  • GAP
    List([1..50],n->(2*n-1)^2+(n-1)*n/2); # Muniru A Asiru, Mar 01 2019
    
  • Magma
    [(2*n-1)^2+(n-1)*n/2: n in [1..50]]; // Vincenzo Librandi, Nov 18 2015
    
  • Maple
    H := n -> simplify(1/hypergeom([-3*n,3*n+3,1],[3/2,2],3/4)); A060544 := n -> H(n-1); seq(A060544(i),i=1..19); # Peter Luschny, Jan 09 2012
  • Mathematica
    Take[Accumulate[Range[150]], {1, -1, 3}] (* Harvey P. Dale, Mar 11 2013 *)
    LinearRecurrence[{3, -3, 1}, {1, 10, 28}, 50] (* Harvey P. Dale, Mar 11 2013 *)
    FoldList[#1 + #2 &, 1, 9 Range @ 50] (* Robert G. Wilson v, Feb 02 2011 *)
    Table[(3 n - 1) (3 n - 2)/2, {n, 20}] (* Eric W. Weisstein, Jul 15 2021 *)
    Table[Binomial[3 n - 1, 2], {n, 20}] (* Eric W. Weisstein, Jul 15 2021 *)
    Table[PolygonalNumber[3 n - 2], {n, 20}] (* Eric W. Weisstein, Jul 15 2021 *)
  • PARI
    a(n)=(3*n-1)*(3*n-2)/2
    
  • Sage
    [(3*n-1)*(3*n-2)/2 for n in (1..50)] # G. C. Greubel, Mar 02 2019

Formula

a(n) = C(3*n, 3)/n = (3*n-1)*(3*n-2)/2 = A001504(n-1)/2.
a(n) = a(n-1) + 9*(n-1) = A060543(n, 3) = A006566(n)/n.
a(n) = A025035(n)/A025035(n-1) = A027468(n-1) + 1 = A000217(3*n-2).
a(1-n) = a(n).
From Paul Barry, Mar 15 2003: (Start)
a(n) = C(n-1, 0) + 9*C(n-1, 1) + 9*C(n-1, 2); binomial transform of (1, 9, 9, 0, 0, 0, ...).
a(n) = 9*A000217(n-1) + 1.
G.f.: x*(1 + 7*x + x^2)/(1-x)^3. (End)
Narayana transform (A001263) of [1, 9, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
a(n-1) = Pochhammer(4,3*n)/(Pochhammer(2,n)*Pochhammer(n+1,2*n)).
a(n-1) = 1/Hypergeometric([-3*n,3*n+3,1],[3/2,2],3/4). - Peter Luschny, Jan 09 2012
From Ant King, Jun 18 2012: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
a(n) = 2*a(n-1) - a(n-2) + 9.
a(n) = A000217(n) + 7*A000217(n-1) + A000217(n-2).
Sum_{n>=1} 1/a(n) = 2*Pi/(3*sqrt(3)) = A248897.
(End)
a(n) = (2*n-1)^2 + (n-1)*n/2. - Ivan N. Ianakiev, Nov 18 2015
a(n) = A101321(9,n-1). - R. J. Mathar, Jul 28 2016
E.g.f.: (2 + 9*x^2)*exp(x)/2 - 1. - G. C. Greubel, Mar 02 2019
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=1} a(n)/n! = 11*e/2 - 1.
Sum_{n>=1} (-1)^n * a(n)/n! = 11/(2*e) - 1. (End)
a(n) = A000567(n) + A005449(n-1) (see illustration in links). - John Elias, Nov 10 2020
a(n) = P(2*n,4)*P(3*n,3)/24 for n>=2, where P(s,k) = ((s - 2)*k^2 - (s - 4)*k)/2 is the k-th s-gonal number. - Lechoslaw Ratajczak, Jul 18 2021

Extensions

Additional description from Terrel Trotter, Jr., Apr 06 2002
Formulas by Paul Berry corrected for offset 1 by Wolfdieter Lang, Jan 30 2015

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

A195160 Generalized 11-gonal (or hendecagonal) numbers: m*(9*m - 7)/2 with m = 0, 1, -1, 2, -2, 3, -3, ...

Original entry on oeis.org

0, 1, 8, 11, 25, 30, 51, 58, 86, 95, 130, 141, 183, 196, 245, 260, 316, 333, 396, 415, 485, 506, 583, 606, 690, 715, 806, 833, 931, 960, 1065, 1096, 1208, 1241, 1360, 1395, 1521, 1558, 1691, 1730, 1870, 1911, 2058, 2101, 2255, 2300, 2461, 2508, 2676
Offset: 0

Views

Author

Omar E. Pol, Sep 10 2011

Keywords

Comments

Exponents of q in the expansion of Product_{n >= 1} (1 - q^(9*n))*(1 + q^(9*n-1))*(1 + q^(9*n-8)) = 1 + q + q^8 + q^11 + q^25 + q^30 + .... - Peter Bala, Nov 21 2024

Crossrefs

Partial sums of A195159.
Column 7 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), this sequence (k=11), A195162 (k=12), A195313 (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
    I:=[0, 1, 8, 11, 25]; [n le 5 select I[n] else Self(n-1)+2*Self(n-2)-2*Self(n-3)-Self(n-4)+Self(n-5): n in [1..50]]; // Vincenzo Librandi, Apr 09 2013
    
  • Mathematica
    CoefficientList[Series[x (1 + 7 x + x^2)/((1 + x)^2 (1 - x)^3), {x, 0, 60}], x] (* Vincenzo Librandi, Apr 09 2013 *)
  • PARI
    a(n)=(18*n*(n+1)+5*(2*n+1)*(-1)^n-5)/16 \\ Charles R Greathouse IV, Sep 24 2015

Formula

From Bruno Berselli, Sep 14 2011: (Start)
G.f.: x*(1+7*x+x^2)/((1+x)^2*(1-x)^3).
a(n) = (18*n*(n+1)+5*(2*n+1)*(-1)^n-5)/16.
a(2n) = A062728(n), a(2n-1) = A051682(n). (End)
Sum_{n>=1} 1/a(n) = 18/49 + 2*Pi*cot(2*Pi/9)/7. - Vaclav Kotesovec, Oct 05 2016

A027468 9 times the triangular numbers A000217.

Original entry on oeis.org

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

Views

Author

Keywords

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

Examples

			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
		

Crossrefs

Programs

  • Magma
    [9*n*(n+1)/2: n in [0..50]]; // Vincenzo Librandi, Dec 29 2012
    
  • 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 *)
  • PARI
    a(n)=9*n*(n+1)/2
    
  • Sage
    [9*binomial(n+1, 2) for n in (0..50)] # G. C. Greubel, May 20 2021

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)

Extensions

More terms from Patrick De Geest, Oct 15 1999

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

Original entry on oeis.org

0, 3, 15, 36, 66, 105, 153, 210, 276, 351, 435, 528, 630, 741, 861, 990, 1128, 1275, 1431, 1596, 1770, 1953, 2145, 2346, 2556, 2775, 3003, 3240, 3486, 3741, 4005, 4278, 4560, 4851, 5151, 5460, 5778, 6105, 6441, 6786, 7140, 7503, 7875, 8256, 8646, 9045
Offset: 0

Views

Author

Floor van Lamoen, Jul 21 2001

Keywords

Comments

Write 0,1,2,3,4,... in a triangular spiral; then a(n) is the sequence found by reading from 0 in the vertical upward direction.
Number of edges in the join of two complete graphs of order 2n and n, K_2n * K_n - Roberto E. Martinez II, Jan 07 2002

Examples

			The spiral begins:
            15
          16  14
        17   3  13
      18   4   2  12
    19   5   0   1  11
  20   6   7   8   9  10
		

Crossrefs

Programs

  • Magma
    [Binomial(3*n,2): n in [0..50]]; // G. C. Greubel, Dec 26 2023
    
  • Maple
    [seq(binomial(3*n,2),n=0..45)]; # Zerinvary Lajos, Jan 02 2007
  • Mathematica
    3*PolygonalNumber[5,Range[0,50]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Mar 06 2019 *)
  • PARI
    a(n)=3*n*(3*n-1)/2 \\ Charles R Greathouse IV, Sep 24 2015
    
  • SageMath
    [binomial(3*n,2) for n in range(51)] # G. C. Greubel, Dec 26 2023

Formula

a(n) = binomial(3*n, 2). - Zerinvary Lajos, Jan 02 2007
a(n) = (9*n^2 - 3*n)/2 = 3*n(3*n-1)/2 = A000326(n)*3. - Omar E. Pol, Dec 11 2008
a(n) = a(n-1) + 9*n - 6, with n > 0, a(0)=0. - Vincenzo Librandi, Aug 07 2010
G.f.: 3*x*(1+2*x)/(1-x)^3. - Bruno Berselli, Jan 21 2011
a(n) = A218470(9n+2). - Philippe Deléham, Mar 27 2013
a(n) = n*A008585(n) + Sum_{i=0..n-1} A008585(i) for n > 0. - Bruno Berselli, Dec 19 2013
From Amiram Eldar, Jan 10 2022: (Start)
Sum_{n>=1} 1/a(n) = log(3) - Pi/(3*sqrt(3)).
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*Pi/(3*sqrt(3)) - 4*log(2)/3. (End)
E.g.f.: (3/2)*x*(2 + 3*x)*exp(x). - G. C. Greubel, Dec 26 2023

Extensions

Better definition and edited by Omar E. Pol, Dec 11 2008

A004188 a(n) = n*(3*n^2 - 1)/2.

Original entry on oeis.org

0, 1, 11, 39, 94, 185, 321, 511, 764, 1089, 1495, 1991, 2586, 3289, 4109, 5055, 6136, 7361, 8739, 10279, 11990, 13881, 15961, 18239, 20724, 23425, 26351, 29511, 32914, 36569, 40485, 44671, 49136, 53889, 58939, 64295, 69966, 75961
Offset: 0

Views

Author

Albert D. Rich (Albert_Rich(AT)msn.com)

Keywords

Comments

3-dimensional analog of centered polygonal numbers.
(1), (4+7), (10+13+16), (19+22+25+28), ... - Jon Perry, Sep 10 2004

References

  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 140.
  • T. P. Martin, Shells of atoms, Phys. Reports, 273 (1996), 199-241, eq. (11).

Crossrefs

1/12*t*(n^3-n)+n for t = 2, 4, 6, ... gives A004006, A006527, A006003, A005900, A004068, A000578, A004126, A000447, A004188, A004466, A004467, A007588, A062025, A063521, A063522, A063523.
Cf. A236770 (partial sums).

Programs

Formula

Partial sums of n-1 3-spaced triangular numbers, e.g., a(4) = t(1) + t(4) + t(7) = 1 + 10 + 28 = 39. - Jon Perry, Jul 23 2003
a(n) = C(2*n+1,3) + C(n+1,3), n >= 0. - Zerinvary Lajos, Jan 21 2007
a(n) = A000447(n) + A000292(n). - Zerinvary Lajos, Jan 21 2007
G.f.: x*(1+7*x+x^2) / (x-1)^4. - R. J. Mathar, Oct 08 2011
From Miquel Cerda, Dec 25 2016: (Start)
a(n) = A000578(n) + A135503(n).
a(n) = A007588(n) - A135503(n). (End)
E.g.f.: (x/2)*(2 + 9*x + 3*x^2)*exp(x). - G. C. Greubel, Sep 01 2017

A022267 a(n) = n*(9*n + 1)/2.

Original entry on oeis.org

0, 5, 19, 42, 74, 115, 165, 224, 292, 369, 455, 550, 654, 767, 889, 1020, 1160, 1309, 1467, 1634, 1810, 1995, 2189, 2392, 2604, 2825, 3055, 3294, 3542, 3799, 4065, 4340, 4624, 4917, 5219, 5530, 5850, 6179
Offset: 0

Views

Author

Keywords

Comments

From Floor van Lamoen, Jul 21 2001: (Start)
Write 0, 1, 2, 3, 4, ... in a triangular spiral; then a(n) is the sequence found by reading the line from 0 in the direction 0, 5, ... . The spiral begins:
.
15
/ \
16 14
/ \
17 3 13
/ / \ \
18 4 2 12
/ / \ \
19 5 0---1 11
/ / \
20 6---7---8---9--10
.
(End)
a(n) is the sum of n consecutive integers starting from 4*n+1: (5), (9+10), (13+14+15), ... - Klaus Purath, Jul 07 2020
a(n) with n>0 are the numbers with the periodic length 3 in the Bulgarian and Mancala solitaire. - Paul Weisenhorn, Jan 29 2022

Crossrefs

Cf. similar sequences listed in A254963.
Cf. similar sequences listed in A022289.

Programs

  • Maple
    seq(binomial(9*n+1,2)/9, n=0..37); # Zerinvary Lajos, Jan 21 2007
  • Mathematica
    Table[ n (9 n + 1)/2, {n, 0, 40}] (* or *) LinearRecurrence[{3, -3, 1}, {0, 5, 19}, 40] (* Harvey P. Dale, Jul 01 2013 *)
  • PARI
    vector(100,n,(n-1)*(9*n-8)/2) \\ Derek Orr, Feb 06 2015

Formula

a(n) = A110449(n, 4) for n>3.
From Bruno Berselli, Feb 11 2011: (Start)
G.f.: x*(5 + 4*x)/(1 - x)^3.
a(n) = 4*A000217(n) + A000566(n). (End)
a(n) = 9*n + a(n-1) - 4 with n>0, a(0)=0. - Vincenzo Librandi, Aug 04 2010
a(n) = A218470(9*n+4). - Philippe Deléham, Mar 27 2013
a(n) = A000217(5*n) - A000217(4*n). - Bruno Berselli, Oct 13 2016
E.g.f.: (1/2)*(9*x^2 + 10*x)*exp(x). - G. C. Greubel, Jul 17 2017
a(n) = A060544(n+1) - A016813(n). - Leo Tavares, Mar 20 2022

A093644 (9,1) Pascal triangle.

Original entry on oeis.org

1, 9, 1, 9, 10, 1, 9, 19, 11, 1, 9, 28, 30, 12, 1, 9, 37, 58, 42, 13, 1, 9, 46, 95, 100, 55, 14, 1, 9, 55, 141, 195, 155, 69, 15, 1, 9, 64, 196, 336, 350, 224, 84, 16, 1, 9, 73, 260, 532, 686, 574, 308, 100, 17, 1, 9, 82, 333, 792, 1218, 1260, 882, 408, 117, 18, 1, 9, 91, 415
Offset: 0

Views

Author

Wolfdieter Lang, Apr 22 2004

Keywords

Comments

The array F(9;n,m) gives in the columns m>=1 the figurate numbers based on A017173, including the 11-gonal numbers A051682 (see the W. Lang link).
This is the ninth member, d=9, in the family of triangles of figurate numbers, called (d,1) Pascal triangles: A007318 (Pascal), A029653, A093560-5, for d=1..8.
This is an example of a Riordan triangle (see A093560 for a comment and A053121 for a comment and the 1991 Shapiro et al. reference on the Riordan group). Therefore the o.g.f. for the row polynomials p(n,x) := Sum_{m=0..n} a(n,m)*x^m is G(z,x) = (1+8*z)/(1-(1+x)*z).
The SW-NE diagonals give A022099(n-1) = Sum_{k=0..ceiling((n-1)/2)} a(n-1-k,k), n >= 1, with n=0 value 8. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs.
Triangle T(n,k), read by rows, given by (9,-8,0,0,0,0,0,0,0,...) DELTA (1,0,0,0,0,0,0,0,...) where DELTA is the operator defined in A084938. - Philippe Deléham, Oct 10 2011

Examples

			Triangle begins
  [1];
  [9,  1];
  [9, 10,  1];
  [9, 19, 11,  1];
  ...
		

References

  • Kurt Hawlitschek, Johann Faulhaber 1580-1635, Veroeffentlichung der Stadtbibliothek Ulm, Band 18, Ulm, Germany, 1995, Ch. 2.1.4. Figurierte Zahlen.
  • Ivo Schneider: Johannes Faulhaber 1580-1635, Birkhäuser, Basel, Boston, Berlin, 1993, ch.5, pp. 109-122.

Crossrefs

Row sums: A020714(n-1), n >= 1, 1 for n=0, alternating row sums are 1 for n=0, 8 for n=2 and 0 otherwise.
The column sequences give for m=1..9: A017173, A051682 (11-gonal), A007586, A051798, A051879, A050405, A052206, A056117, A056003.
Cf. A093645 (d=10).

Programs

  • Haskell
    a093644 n k = a093644_tabl !! n !! k
    a093644_row n = a093644_tabl !! n
    a093644_tabl = [1] : iterate
                   (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [9, 1]
    -- Reinhard Zumkeller, Aug 31 2014
  • Mathematica
    Join[{1},Table[Binomial[n,k]+8Binomial[n-1,k],{n,20},{k,0,n}]//Flatten] (* Harvey P. Dale, Aug 17 2024 *)

Formula

a(n, m) = F(9;n-m, m) for 0 <= m <= n, otherwise 0, with F(9;0, 0)=1, F(9;n, 0)=9 if n >= 1 and F(9;n, m):=(9*n+m)*binomial(n+m-1, m-1)/m if m >= 1.
Recursion: a(n, m)=0 if m > n, a(0, 0)= 1; a(n, 0)=9 if n >= 1; a(n, m) = a(n-1, m) + a(n-1, m-1).
G.f. column m (without leading zeros): (1+8*x)/(1-x)^(m+1), m >= 0.
T(n, k) = C(n, k) + 8*C(n-1, k). - Philippe Deléham, Aug 28 2005
Row n: Expansion of (9+x)*(1+x)^(n-1), n > 0. - Philippe Deléham, Oct 10 2011
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(9 + 19*x + 11*x^2/2! + x^3/3!) = 9 + 28*x + 58*x^2/2! + 100*x^3/3! + 155*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ). - Peter Bala, Dec 22 2014
G.f.: (-1-8*x)/(-1+x+x*y). - R. J. Mathar, Aug 11 2015
Previous Showing 11-20 of 69 results. Next