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 10 results.

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

A303299 Generalized 22-gonal (or icosidigonal) numbers: m*(10*m - 9) with m = 0, +1, -1, +2, -2, +3, -3, ...

Original entry on oeis.org

0, 1, 19, 22, 58, 63, 117, 124, 196, 205, 295, 306, 414, 427, 553, 568, 712, 729, 891, 910, 1090, 1111, 1309, 1332, 1548, 1573, 1807, 1834, 2086, 2115, 2385, 2416, 2704, 2737, 3043, 3078, 3402, 3439, 3781, 3820, 4180, 4221, 4599, 4642, 5038, 5083, 5497, 5544, 5976, 6025, 6475, 6526, 6994, 7047, 7533, 7588
Offset: 0

Views

Author

Omar E. Pol, Jun 23 2018

Keywords

Comments

Partial sums of A317318. - Omar E. Pol, Jul 28 2018
Exponents in expansion of Product_{n >= 1} (1 + x^(20*n-19))*(1 + x^(20*n-1))*(1 - x^(20*n)) = 1 + x + x^19 + x^22 + x^58 + .... - Peter Bala, Dec 10 2020

Crossrefs

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), A303305 (k=17), A274979 (k=18), A303813 (k=19), A218864 (k=20), A303298 (k=21), this sequence (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

  • Maple
    a:= n-> (m-> m*(10*m-9))(-ceil(n/2)*(-1)^n):
    seq(a(n), n=0..60);  # Alois P. Heinz, Jun 23 2018
  • Mathematica
    CoefficientList[ Series[-x (x^2 + 18x + 1)/((x - 1)^3 (x + 1)^2), {x, 0, 50}], x] (* or *)LinearRecurrence[{1, 2, -2, -1, 1}, {0, 1, 19, 22, 58}, 51] (* Robert G. Wilson v, Jul 28 2018 *)
    nn=30; Sort[Table[n (10 n - 9), {n, -nn, nn}]] (* Vincenzo Librandi, Jul 29 2018 *)
  • PARI
    a(n) = n++; my(m = (-1) ^ n * (n >> 1)); m * (10 * m - 9) \\ David A. Corneth, Jun 23 2018
    
  • PARI
    concat(0, Vec(x*(1 + 18*x + x^2) / ((1 - x)^3*(1 + x)^2) + O(x^60))) \\ Colin Barker, Jun 23 2018

Formula

From Colin Barker, Jun 23 2018: (Start)
G.f.: x*(1 + 18*x + x^2) / ((1 - x)^3*(1 + x)^2).
a(n) = (5*n^2 + 9*n)/2 for n even.
a(n) = (5*n^2 + n - 4)/2 for n odd.
a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5) for n>4.
(End)
Sum_{n>=1} 1/a(n) = (10 + 9*sqrt(5+2*sqrt(5))*Pi)/81. - Amiram Eldar, Mar 01 2022

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)

A051868 16-gonal (or hexadecagonal) numbers: a(n) = n*(7*n-6).

Original entry on oeis.org

0, 1, 16, 45, 88, 145, 216, 301, 400, 513, 640, 781, 936, 1105, 1288, 1485, 1696, 1921, 2160, 2413, 2680, 2961, 3256, 3565, 3888, 4225, 4576, 4941, 5320, 5713, 6120, 6541, 6976, 7425, 7888, 8365, 8856, 9361, 9880, 10413, 10960, 11521
Offset: 0

Views

Author

N. J. A. Sloane, Dec 15 1999

Keywords

Comments

Sequence found by reading the line from 0, in the direction 0, 16, ... and the parallel line from 1, in the direction 1, 45, ..., in the square spiral whose vertices are the generalized 16-gonal numbers. - Omar E. Pol, Jul 18 2012
This is also a star octagonal number: a(n) = A000567(n) + 8*A000217(n-1). - Luciano Ancora, Mar 29 2015
Let T(n) = A000217(n), the n-th triangular number. Then a(n) = T(n-1) + T(4n-3) - T(2n-4) + T(n-3). In general, let P(k,n) be the n-th k-gonal number. Then for k>1, P(T(k)+1,n) = T(n-1) + T((k-1)n-(k-2)) - T((k-3)n-2(k-3)) + T((k-4)n-3(k-4)) - ... + (-1)^(k+1)*T(n-(k-2)). - Charlie Marion, Dec 23 2019

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, p. 189.
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6.

Crossrefs

Programs

Formula

a(n) = 14*n + a(n-1) - 13, with n>0, a(0)=0. - Vincenzo Librandi, Aug 06 2010
G.f.: x*(1+13*x)/(1-x)^3. - Bruno Berselli, Feb 04 2011
a(0)=0, a(1)=1, a(2)=16; for n>2, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, May 07 2011
a(14*a(n) + 92*n + 1) = a(14*a(n) + 92*n) + a(14*n+1). - Vladimir Shevelev, Jan 24 2014
E.g.f.: exp(x)*x*(1 + 7*x). - Stefano Spezia, Dec 27 2019
a(n) = (4*n-3)^2 - (3*n-3)^2. In general, if we let P(k,n) be the n-th k-gonal number, then P(4k,n) = (k*n-k+1)^2 - ((k-1)*n-k+1)^2. In addition, {P(4k,n)} are the only polygonal number sequences each of whose terms can be written as the difference of two squares. - Charlie Marion, Feb 16 2020
Product_{n>=2} (1 - 1/a(n)) = 7/8. - Amiram Eldar, Jan 22 2021

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

A152965 Twice 12-gonal numbers: a(n) = 2*n*(5*n-4).

Original entry on oeis.org

0, 2, 24, 66, 128, 210, 312, 434, 576, 738, 920, 1122, 1344, 1586, 1848, 2130, 2432, 2754, 3096, 3458, 3840, 4242, 4664, 5106, 5568, 6050, 6552, 7074, 7616, 8178, 8760, 9362, 9984, 10626, 11288, 11970, 12672, 13394, 14136, 14898, 15680, 16482, 17304, 18146, 19008
Offset: 0

Views

Author

Omar E. Pol, Dec 21 2008

Keywords

Crossrefs

Cf. A051624 (12-gonal numbers), A051874.
Cf. numbers of the form n*(n*k - k + 4)/2 listed in A226488 (this sequence is the case k=20). - Bruno Berselli, Jun 10 2013

Programs

Formula

a(n) = 2*A051624(n).
From Vincenzo Librandi, Jul 10 2012: (Start)
G.f.: 2*x*(1+9*x)/(1-x)^3.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). (End)
From Elmo R. Oliveira, Dec 27 2024: (Start)
E.g.f.: 2*exp(x)*x*(1 + 5*x).
a(n) = n + A051874(n). (End)

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.

A131103 Rectangular array read by antidiagonals: a(n, k) is the number of ways to put k labeled objects into n labeled boxes so that there are no boxes with exactly one object (n, k >= 1).

Original entry on oeis.org

0, 0, 1, 0, 2, 1, 0, 3, 2, 1, 0, 4, 3, 8, 1, 0, 5, 4, 21, 22, 1, 0, 6, 5, 40, 63, 52, 1, 0, 7, 6, 65, 124, 243, 114, 1, 0, 8, 7, 96, 205, 664, 969, 240, 1, 0, 9, 8, 133, 306, 1405, 3196, 3657, 494, 1, 0, 10, 9, 176, 427, 2556, 7425, 15712, 12987, 1004, 1, 0, 11, 10, 225, 568, 4207
Offset: 1

Views

Author

David Wasserman, Jun 14 2007, Jun 15 2007

Keywords

Comments

Problem suggested by Brandon Zeidler. Columns four and five are A000567 and A051874. Second row is A130102.

Examples

			Array begins:
0 1 1 1 1 1 1
0 2 2 8 22 52 114
0 3 3 21 63 243 969
		

Crossrefs

Formula

a(n, k) = sum_{j=1..min(floor(k/2), n)} A008299(k, j)*n!/(n-j)!.

A326474 A(n, k) = (m*k)! [x^k] MittagLefflerE(m, x)^n, for m = 3, n >= 0, k >= 0; square array read by descending antidiagonals.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 22, 3, 1, 0, 1, 170, 63, 4, 1, 0, 1, 1366, 2187, 124, 5, 1, 0, 1, 10922, 59535, 7732, 205, 6, 1, 0, 1, 87382, 1594323, 599548, 18485, 306, 7, 1, 0, 1, 699050, 43033599, 39945364, 2416045, 36126, 427, 8, 1
Offset: 0

Views

Author

Peter Luschny, Jul 08 2019

Keywords

Examples

			Array starts:
[0] 1, 0,   0,     0,       0,          0,            0, ... A000007
[1] 1, 1,   1,     1,       1,          1,            1, ... A000012
[2] 1, 2,  22,   170,    1366,      10922,        87382, ... A007613
[3] 1, 3,  63,  2187,   59535,    1594323,     43033599, ...
[4] 1, 4, 124,  7732,  599548,   39945364,   2556712828, ...
[5] 1, 5, 205, 18485, 2416045,  352060805,  46660373965, ...
[6] 1, 6, 306, 36126, 6673266, 1544907006, 379696000626, ...
      A051874,
		

Crossrefs

Rows include: A000007, A000012, A007613.
Columns include: A051874.
Cf. A326476 (m=2, p>=0), A326327 (m=2, p<=0), this sequence (m=3, p>=0), A326475 (m=3, p<=0).

Programs

  • Mathematica
    (* The function MLPower is defined in A326327. *)
    For[n = 0, n < 8, n++, Print[MLPower[3, n, 8]]]
  • Sage
    # uses[MLPower from A326327]
    for n in (0..6): print(MLPower(3, n, 9))

A330892 Square array of polygonal numbers read by descending antidiagonals (the transpose of A317302).

Original entry on oeis.org

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

Views

Author

Robert G. Wilson v, Apr 27 2020

Keywords

Comments

\c 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 ...
r\
_0 0 1 0 -3 -8 -15 -24 -35 -48 -63 -80 -99 -120 -143 -168 -195 A067998
_1 0 1 1 0 -2 -5 -9 -14 -20 -27 -35 -44 -54 -65 -77 -90 A080956
_2 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 A001477
_3 0 1 3 6 10 15 21 28 36 45 55 66 78 91 105 120 A000217
_4 0 1 4 9 16 25 36 49 64 81 100 121 144 169 196 225 A000290
_5 0 1 5 12 22 35 51 70 92 117 145 176 210 247 287 330 A000326
_6 0 1 6 15 28 45 66 91 120 153 190 231 276 325 378 435 A000384
_7 0 1 7 18 34 55 81 112 148 189 235 286 342 403 469 540 A000566
_8 0 1 8 21 40 65 96 133 176 225 280 341 408 481 560 645 A000567
_9 0 1 9 24 46 75 111 154 204 261 325 396 474 559 651 750 A001106
10 0 1 10 27 52 85 126 175 232 297 370 451 540 637 742 855 A001107
11 0 1 11 30 58 95 141 196 260 333 415 506 606 715 833 960 A051682
12 0 1 12 33 64 105 156 217 288 369 460 561 672 793 924 1065 A051624
13 0 1 13 36 70 115 171 238 316 405 505 616 738 871 1015 1170 A051865
14 0 1 14 39 76 125 186 259 344 441 550 671 804 949 1106 1275 A051866
15 0 1 15 42 82 135 201 280 372 477 595 726 870 1027 1197 1380 A051867
...
Each row has a second forward difference of (r-2) and each column has a forward difference of c(c-1)/2.

Crossrefs

Cf. A317302 (the same array) but read by ascending antidiagonals.
Sub-arrays: A089000, A139600, A206735;
Number of times k>1 appears: A129654, First occurrence of k: A063778.

Programs

  • Mathematica
    Table[ PolygonalNumber[r - c, c], {r, 0, 11}, {c, r, 0, -1}] // Flatten

Formula

P(r, c) = (r - 2)(c(c-1)/2) + c.
Showing 1-10 of 10 results.