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

A300656 Triangle read by rows: T(n,k) = 30*k^2*(n-k)^2 + 1; n >= 0, 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 31, 1, 1, 121, 121, 1, 1, 271, 481, 271, 1, 1, 481, 1081, 1081, 481, 1, 1, 751, 1921, 2431, 1921, 751, 1, 1, 1081, 3001, 4321, 4321, 3001, 1081, 1, 1, 1471, 4321, 6751, 7681, 6751, 4321, 1471, 1, 1, 1921, 5881, 9721, 12001, 12001, 9721, 5881, 1921, 1
Offset: 0

Views

Author

Kolosov Petro, Mar 10 2018

Keywords

Comments

From Kolosov Petro, Apr 12 2020: (Start)
Let A(m, r) = A302971(m, r) / A304042(m, r).
Let L(m, n, k) = Sum_{r=0..m} A(m, r) * k^r * (n - k)^r.
Then T(n, k) = L(2, n, k).
Fifth power can be expressed as row sum of triangle T(n, k).
T(n, k) is symmetric: T(n, k) = T(n, n-k). (End)

Examples

			Triangle begins:
--------------------------------------------------------------------------
k=    0     1     2      3      4      5      6      7     8     9    10
--------------------------------------------------------------------------
n=0:  1;
n=1:  1,    1;
n=2:  1,   31,    1;
n=3:  1,  121,  121,     1;
n=4:  1,  271,  481,   271,     1;
n=5:  1,  481, 1081,  1081,   481,     1;
n=6:  1,  751, 1921,  2431,  1921,   751,     1;
n=7:  1, 1081, 3001,  4321,  4321,  3001,  1081,     1;
n=8:  1, 1471, 4321,  6751,  7681,  6751,  4321,  1471,    1;
n=9:  1, 1921, 5881,  9721, 12001, 12001,  9721,  5881, 1921,    1;
n=10: 1, 2431, 7681, 13231, 17281, 18751, 17281, 13231, 7681, 2431,   1;
		

Crossrefs

Various cases of L(m, n, k): A287326(m=1), This sequence (m=2), A300785(m=3). See comments for L(m, n, k).
Row sums give the nonzero terms of A002561.

Programs

  • GAP
    T:=Flat(List([0..9],n->List([0..n],k->30*k^2*(n-k)^2+1))); # Muniru A Asiru, Oct 24 2018
    
  • Magma
    [[30*k^2*(n-k)^2+1: k in [0..n]]: n in [0..12]]; // G. C. Greubel, Dec 14 2018
    
  • Maple
    a:=(n,k)->30*k^2*(n-k)^2+1: seq(seq(a(n,k),k=0..n),n=0..9); # Muniru A Asiru, Oct 24 2018
  • Mathematica
    T[n_, k_] := 30 k^2 (n - k)^2 + 1; Column[
    Table[T[n, k], {n, 0, 10}, {k, 0, n}], Center] (* Kolosov Petro, Apr 12 2020 *)
    f[n_]:=Table[SeriesCoefficient[(1 + 26 y + 336 y^2 + 326 y^3 + 31 y^4 + x^2 (1 + 116 y + 486 y^2 + 116 y^3 + y^4) + x (-2 - 82 y - 882 y^2 - 502 y^3 + 28 y^4))/((-1 + x)^3 (-1 + y)^5), {x, 0, i}, {y, 0, j}], {i, n, n}, {j, 0, n}]; Flatten[Array[f, 11, 0]] (* Stefano Spezia, Oct 30 2018 *)
  • PARI
    t(n, k) = 30*k^2*(n-k)^2+1
    trianglerows(n) = for(x=0, n-1, for(y=0, x, print1(t(x, y), ", ")); print(""))
    /* Print initial 9 rows of triangle as follows */ trianglerows(9)
    
  • Sage
    [[30*k^2*(n-k)^2+1 for k in range(n+1)] for n in range(12)] # G. C. Greubel, Dec 14 2018

Formula

From Kolosov Petro, Apr 12 2020: (Start)
T(n, k) = 30 * k^2 * (n-k)^2 + 1.
T(n, k) = 30 * A094053(n,k)^2 + 1.
T(n, k) = A158558((n-k) * k).
T(n+2, k) = 3*T(n+1, k) - 3*T(n, k) + T(n-1, k), for n >= k.
Sum_{k=1..n} T(n, k) = A000584(n).
Sum_{k=0..n-1} T(n, k) = A000584(n).
Sum_{k=0..n} T(n, k) = A002561(n).
Sum_{k=1..n-1} T(n, k) = A258807(n).
Sum_{k=1..n-1} T(n, k) = -A024003(n), n > 1.
Sum_{k=1..r} T(n, k) = A316349(2,r,0)*n^0 - A316349(2,r,1)*n^1 + A316349(2,r,2)*n^2. (End)
G.f.: (1 + 26*y + 336*y^2 + 326*y^3 + 31*y^4 + x^2*(1 + 116*y + 486*y^2 + 116*y^3 + y^4) + x*(-2 - 82*y - 882*y^2 - 502*y^3 + 28*y^4))/((-1 + x)^3*(-1 + y)^5). - Stefano Spezia, Oct 30 2018

A300785 Triangle read by rows: T(n,k) = 140*k^3*(n-k)^3 - 14*k*(n-k) + 1; n >= 0, 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 127, 1, 1, 1093, 1093, 1, 1, 3739, 8905, 3739, 1, 1, 8905, 30157, 30157, 8905, 1, 1, 17431, 71569, 101935, 71569, 17431, 1, 1, 30157, 139861, 241753, 241753, 139861, 30157, 1, 1, 47923, 241753, 472291, 573217, 472291, 241753, 47923, 1, 1, 71569, 383965, 816229, 1119721, 1119721, 816229, 383965, 71569, 1
Offset: 0

Views

Author

Kolosov Petro, Mar 12 2018

Keywords

Comments

From Kolosov Petro, Apr 12 2020: (Start)
Let A(m, r) = A302971(m, r) / A304042(m, r).
Let L(m, n, k) = Sum_{r=0..m} A(m, r) * k^r * (n - k)^r.
Then T(n, k) = L(3, n, k).
T(n, k) is symmetric: T(n, k) = T(n, n-k). (End)

Examples

			Triangle begins:
--------------------------------------------------------------------
k=   0      1       2       3       4       5       6      7     8
--------------------------------------------------------------------
n=0: 1;
n=1: 1,     1;
n=2: 1,   127,      1;
n=3: 1,  1093,   1093,      1;
n=4: 1,  3739,   8905,   3739,      1;
n=5: 1,  8905,  30157,  30157,   8905,      1;
n=6: 1, 17431,  71569, 101935,  71569,  17431,      1;
n=7: 1, 30157, 139861, 241753, 241753, 139861,  30157,     1;
n=8: 1, 47923, 241753, 472291, 573217, 472291, 241753, 47923,    1;
		

Crossrefs

Various cases of L(m, n, k): A287326 (m=1), A300656 (m=2), This sequence (m=3). See comments for L(m, n, k).
Row sums give A258806.

Programs

  • GAP
    T:=Flat(List([0..9], n->List([0..n], k->140*k^3*(n-k)^3 - 14*k*(n-k)+1))); # G. C. Greubel, Dec 14 2018
  • Magma
    /* As triangle */ [[140*k^3*(n-k)^3-14*k*(n-k)+1: k in [0..n]]: n in [0..10]]; // Bruno Berselli, Mar 21 2018
    
  • Maple
    T:=(n,k)->140*k^3*(n-k)^3-14*k*(n-k)+1: seq(seq(T(n,k),k=0..n),n=0..9); # Muniru A Asiru, Dec 14 2018
  • Mathematica
    T[n_, k_] := 140*k^3*(n - k)^3 - 14*k*(n - k) + 1; Column[
    Table[T[n, k], {n, 0, 10}, {k, 0, n}], Center] (* From Kolosov Petro, Apr 12 2020 *)
  • PARI
    t(n, k) = 140*k^3*(n-k)^3-14*k*(n-k)+1
    trianglerows(n) = for(x=0, n-1, for(y=0, x, print1(t(x, y), ", ")); print(""))
    /* Print initial 9 rows of triangle as follows */ trianglerows(9)
    
  • Sage
    [[140*k^3*(n-k)^3 - 14*k*(n-k)+1 for k in range(n+1)] for n in range(12)] # G. C. Greubel, Dec 14 2018
    

Formula

From Kolosov Petro, Apr 12 2020: (Start)
T(n, k) = 140*k^3*(n-k)^3 - 14*k*(n-k) + 1.
T(n, k) = 140*A094053(n, k)^3 + 0*A094053(n, k)^2 - 14*A094053(n, k)^1 + 1.
T(n+3, k) = 4*T(n+2, k) - 6*T(n+1, k) + 4*T(n, k) - T(n-1, k), for n >= k.
Sum_{k=1..n} T(n, k) = A001015(n).
Sum_{k=0..n} T(n, k) = A258806(n).
Sum_{k=0..n-1} T(n, k) = A001015(n).
Sum_{k=1..n-1} T(n, k) = A258808(n).
Sum_{k=1..n-1} T(n, k) = -A024005(n).
Sum_{k=1..r} T(n, k) = -A316387(3, r, 0)*n^0 + A316387(3, r, 1)*n^1 - A316387(3, r, 2)*n^2 + A316387(3, r, 3)*n^3. (End)
G.f.: (1 + 127*x^6*y^3 - 3*x*(1 + y) + 585*x^5*y^2*(1 + y) + 129*x^4*y*(1 + 17*y + y^2) + 3*x^2*(1 + 45*y + y^2) - x^3*(1 - 579*y - 579*y^2 + y^3))/((1 - x)^4*(1 - x*y)^4). - Stefano Spezia, Sep 14 2024

A302971 Triangle read by rows: T(n,k) is the numerator of R(n,k) defined implicitly by the identity Sum_{i=0..l-1} Sum_{j=0..m} R(m,j)*(l-i)^j*i^j = l^(2*m+1) holding for all l,m >= 0.

Original entry on oeis.org

1, 1, 6, 1, 0, 30, 1, -14, 0, 140, 1, -120, 0, 0, 630, 1, -1386, 660, 0, 0, 2772, 1, -21840, 18018, 0, 0, 0, 12012, 1, -450054, 491400, -60060, 0, 0, 0, 51480, 1, -11880960, 15506040, -3712800, 0, 0, 0, 0, 218790, 1, -394788954, 581981400, -196409840, 8817900, 0, 0, 0, 0, 923780, 1, -16172552880, 26003271294, -10863652800, 1031151660, 0, 0, 0, 0, 0, 3879876
Offset: 0

Views

Author

Kolosov Petro, Apr 16 2018

Keywords

Examples

			Triangle begins:
------------------------------------------------------------------------
k=   0          1         2         3    4     5      6      7       8
------------------------------------------------------------------------
n=0: 1;
n=1: 1,         6;
n=2: 1,         0,       30;
n=3: 1,       -14,        0,      140;
n=4: 1,      -120,        0,        0, 630;
n=5: 1,     -1386,      660,        0,   0, 2772;
n=6: 1,    -21840,    18018,        0,   0,    0, 12012;
n=7: 1,   -450054,   491400,   -60060,   0,    0,     0, 51480;
n=8: 1, -11880960, 15506040, -3712800,   0,    0,     0,     0, 218790;
		

Crossrefs

Items of second row are the coefficients in the definition of A287326.
Items of third row are the coefficients in the definition of A300656.
Items of fourth row are the coefficients in the definition of A300785.
T(n,n) gives A002457(n).
Denominators of R(n,k) are shown in A304042.
Row sums return A000079(2n+1) - 1.

Programs

  • Maple
    R := proc(n, k) if k < 0 or k > n then return 0 fi; (2*k+1)*binomial(2*k, k);
    if n = k then % else -%*add((-1)^j*R(n, j)*binomial(j, 2*k+1)*
    bernoulli(2*j-2*k)/(j-k), j=2*k+1..n) fi end: T := (n, k) -> numer(R(n, k)):
    seq(print(seq(T(n, k), k=0..n)), n=0..12);
    # Numerical check that S(m, n) = n^(2*m+1):
    S := (m, n) -> add(add(R(m, j)*(n-k)^j*k^j, j=0..m), k=0..n-1):
    seq(seq(S(m, n) - n^(2*m+1), n=0..12), m=0..12); # Peter Luschny, Apr 30 2018
  • Mathematica
    R[n_, k_] := 0
    R[n_, k_] := (2 k + 1)*Binomial[2 k, k]*
       Sum[R[n, j]*Binomial[j, 2 k + 1]*(-1)^(j - 1)/(j - k)*
       BernoulliB[2 j - 2 k], {j, 2 k + 1, n}] /; 2 k + 1 <= n
    R[n_, k_] := (2 n + 1)*Binomial[2 n, n] /; k == n;
    T[n_, k_] := Numerator[R[n, k]];
    (* Print Fifteen Initial rows of Triangle A302971 *)
    Column[ Table[ T[n, k], {n, 0, 15}, {k, 0, n}], Center]
  • PARI
    T(n, k) = if ((n>k) || (n<0), 0, if (k==n, (2*n+1)*binomial(2*n, n), if (2*n+1>k, 0, if (n==0, 1, (2*n+1)*binomial(2*n, n)*sum(j=2*n+1, k+1, T(j, k)*binomial(j, 2*n+1)*(-1)^(j-1)/(j-n)*bernfrac(2*j-2*n))))));
    tabl(nn) = for (n=0, nn, for (k=0, n, print1(numerator(T(k,n)), ", ")); print); \\ Michel Marcus, Apr 27 2018

Formula

Recurrence given by Max Alekseyev (see the MathOverflow link):
R(n, k) = 0 if k < 0 or k > n.
R(n, k) = (2k+1)*binomial(2k, k) if k = n.
R(n, k) = (2k+1)*binomial(2k, k)*Sum_{j=2k+1..n} R(n, j)*binomial(j, 2k+1)*(-1)^(j-1)/(j-k)*Bernoulli(2j-2k), otherwise.
T(n, k) = numerator(R(n, k)).

A304042 Triangle read by rows: T(n,k) is the denominator of R(n,k) defined implicitly by the identity Sum_{i=0..l-1} Sum_{j=0..m} R(m,j)*(l-i)^j*i^j = l^(2*m+1) holding for all l,m >= 0.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 3, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 0

Views

Author

Kolosov Petro, May 05 2018

Keywords

Examples

			Triangle begins:
-----------------------------------------------------
k=    0  1  2  3  4  5  6  7  8  9  10 11 12 13 14 15
-----------------------------------------------------
n=0:  1;
n=1:  1, 1;
n=2:  1, 1, 1;
n=3:  1, 1, 1, 1;
n=4:  1, 1, 1, 1, 1;
n=5:  1, 1, 1, 1, 1, 1;
n=6:  1, 1, 1, 1, 1, 1, 1;
n=7:  1, 1, 1, 1, 1, 1, 1, 1;
n=8:  1, 1, 1, 1, 1, 1, 1, 1, 1;
n=9:  1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n=10: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n=11: 1, 5, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1;
n=12: 1, 3, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n=13: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n=14: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n=15: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
		

Crossrefs

Programs

  • Mathematica
    R[n_, k_] := 0
    R[n_, k_] := (2 k + 1)*Binomial[2 k, k]*
       Sum[R[n, j]*Binomial[j, 2 k + 1]*(-1)^(j - 1)/(j - k)*
       BernoulliB[2 j - 2 k], {j, 2 k + 1, n}] /; 2 k + 1 <= n
    R[n_, k_] := (2 n + 1)*Binomial[2 n, n] /; k == n;
    T[n_, k_] := Denominator[R[n, k]];
    (* Print Fifteen Initial rows of Triangle A304042 *)
    Column[ Table[ T[n, k], {n, 0, 15}, {k, 0, n}], Center]
  • PARI
    up_to = 1274; \\ = binomial(50+1,2)-1
    A304042aux(n, k) = if((k<0)||(k>n),0,(k+k+1)*binomial(2*k, k)*if(k==n,1,sum(j=k+k+1,n, A304042aux(n, j)*binomial(j, k+k+1)*((-1)^(j-1))/(j-k)*bernfrac(2*(j-k)))));
    A304042tr(n, k) = denominator(A304042aux(n, k));
    A304042list(up_to) = { my(v = vector(up_to), i=0); for(n=0,oo, for(k=0,n, if(i++ > up_to, return(v)); v[i] = A304042tr(n,k))); (v); };
    v304042 = A304042list(1+up_to);
    A304042(n) = v304042[1+n]; \\ Antti Karttunen, Nov 07 2018

Formula

Recurrence given by Max Alekseyev (see the MathOverflow link):
R(n, k) = 0 if k < 0 or k > n.
R(n, k) = (2k+1)*binomial(2k, k) if k = n.
R(n, k) = (2k+1)*binomial(2k, k)*Sum_{j=2k+1..n} R(n, j)*binomial(j, 2k+1)*(-1)^(j-1)/(j-k)*Bernoulli(2j-2k), otherwise.
T(n, k) = denominator(R(n, k)).

A130154 Triangle read by rows: T(n, k) = 1 + 2*(n-k)*(k-1) (1 <= k <= n).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 5, 1, 1, 7, 9, 7, 1, 1, 9, 13, 13, 9, 1, 1, 11, 17, 19, 17, 11, 1, 1, 13, 21, 25, 25, 21, 13, 1, 1, 15, 25, 31, 33, 31, 25, 15, 1, 1, 17, 29, 37, 41, 41, 37, 29, 17, 1, 1, 19, 33, 43, 49, 51, 49, 43, 33, 19, 1, 1, 21, 37, 49, 57, 61, 61, 57, 49, 37, 21, 1
Offset: 1

Views

Author

Emeric Deutsch, May 22 2007

Keywords

Comments

Column k, except for the initial k-1 0's, is an arithmetic progression with first term 1 and common difference 2(k-1). Row sums yield A116731. First column of the inverse matrix is A129779.
Studied by Paul Curtz circa 1993.
From Rogério Serôdio, Dec 19 2017: (Start)
T(n, k) gives the number of distinct sums of 2(k-1) elements in {1,1,2,2,...,n-1,n-1}. For example, T(6, 2) = the number of distinct sums of 2 elements in {1,1,2,2,3,3,4,4,5,5}, and because each sum from the smallest 1 + 1 = 2 to the largest 5 + 5 = 10 appears, T(6, 2) = 10 - 1 = 9. [In general: 2*(Sum_{j=1..(k-1)} n-j) - (2*(Sum_{j=1..k-1} j) - 1) = 2*n*(k-1) - 4*(k-1)*k/2 + 1 = 2*(k-1)*(n-k) + 1 = T(n, k). - Wolfdieter Lang, Dec 20 2017]
T(n, k) is the number of lattice points with abscissa x = 2*(k-1) and even ordinate in the closed region bounded by the parabola y = x*(2*(n-1) - x) and the x axis. [That is, (1/2)*y(2*(k-1)) + 1 = T(n, k). - Wolfdieter Lang, Dec 20 2017]
Pascal's triangle (A007318, but with apex in the middle) is formed using the rule South = West + East; the rascal triangle A077028 uses the rule South = (West*East + 1)/North; the present triangle uses a similar rule: South = (West*East + 2)/North. See the formula section for this recurrence. (End)

Examples

			The triangle T(n, k) starts:
  n\k  1  2  3  4  5  6  7  8  9 10 ...
  1:   1
  2:   1  1
  3:   1  3  1
  4:   1  5  5  1
  5:   1  7  9  7  1
  6:   1  9 13 13  9  1
  7:   1 11 17 19 17 11  1
  8:   1 13 21 25 25 21 13  1
  9:   1 15 25 31 33 31 25 15  1
 10:   1 17 29 37 41 41 37 29 17  1
 ... reformatted. - _Wolfdieter Lang_, Dec 19 2017
		

Crossrefs

Column sequences (no leading zeros): A000012, A016813, A016921, A017077, A017281, A017533, A131877, A158057, A161705, A215145.

Programs

  • GAP
    Flat(List([1..12], n-> List([1..n], k-> 1 + 2*(n-k)*(k-1) ))); # G. C. Greubel, Nov 25 2019
  • Magma
    [1 + 2*(n-k)*(k-1): k in [1..n], n in [1..12]]; // G. C. Greubel, Nov 25 2019
    
  • Maple
    T:=proc(n,k) if k<=n then 2*(n-k)*(k-1)+1 else 0 fi end: for n from 1 to 14 do seq(T(n,k),k=1..n) od; # yields sequence in triangular form
  • Mathematica
    Flatten[Table[1+2(n-k)(k-1),{n,0,20},{k,n}]] (* Harvey P. Dale, Jul 13 2013 *)
  • PARI
    T(n, k) = 1 + 2*(n-k)*(k-1) \\ Iain Fox, Dec 19 2017
    
  • PARI
    first(n) = my(res = vector(binomial(n+1,2)), i = 1); for(r=1, n, for(k=1, r, res[i] = 1 + 2*(r-k)*(k-1); i++)); res \\ Iain Fox, Dec 19 2017
    
  • Sage
    [[1 + 2*(n-k)*(k-1) for k in (1..n)] for n in (1..12)] # G. C. Greubel, Nov 25 2019
    

Formula

T(n, k) = 1 + 2*(n-k)*(k-1) (1 <= k <= n).
G.f.: G(t,z) = t*z*(3*t*z^2 - z - t*z + 1)/((1-t*z)*(1-z))^2.
Equals = 2 * A077028 - A000012 as infinite lower triangular matrices. - Gary W. Adamson, Oct 23 2007
T(n, 1) = 1 and T(n, n) = 1 for n >= 1; T(n, k) = (T(n-1, k-1)*T(n-1, k) + 2)/T(n-2, k-1), for n > 2 and 1 < k < n. See a comment above. - Rogério Serôdio, Dec 19 2017
G.f. column k (with leading zeros): (x^k/(1-x)^2)*(1 + (2*k-3)*x), k >= 1. See the g.f. of the triangle G(t,z) above: (d/dt)^k G(t,x)/k!|{t=0}. - _Wolfdieter Lang, Dec 20 2017

Extensions

Edited by Wolfdieter Lang, Dec 19 2017

A051744 a(n) = n*(n+1)*(n^2+5*n+18)/24.

Original entry on oeis.org

2, 8, 21, 45, 85, 147, 238, 366, 540, 770, 1067, 1443, 1911, 2485, 3180, 4012, 4998, 6156, 7505, 9065, 10857, 12903, 15226, 17850, 20800, 24102, 27783, 31871, 36395, 41385, 46872, 52888, 59466, 66640, 74445, 82917, 92093, 102011, 112710, 124230, 136612
Offset: 1

Views

Author

Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Dec 07 1999

Keywords

Comments

a(n) is the number of binary words with length <= n+1 which contain at least one 0 and one 1 and have at most one ascent. - Amelia Gibbs, May 21 2024

Programs

Formula

a(n) = binomial(n+3, n-1) + binomial(n+1, n-1).
G.f.: x*(2-2*x+x^2)/(1-x)^5. - Colin Barker, Mar 19 2012
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). - Vincenzo Librandi, Apr 27 2012
a(n) = sum_{k=1..n} sum{j=1..k} sum{i=1..j} (i + binomial(j,k)). - Wesley Ivan Hurt, Nov 01 2014
E.g.f.: (1/24)*x*(x^3+12*x^2+48*x+48)*exp(x). - Robert Israel, Nov 02 2014
a(n) = Sum_{i=1..n+1} Sum_{j=1...i-1} A077028(i,j). - Amelia Gibbs, May 21 2024

A300401 Array T(n,k) = n*(binomial(k, 2) + 1) + k*(binomial(n, 2) + 1) read by antidiagonals.

Original entry on oeis.org

0, 1, 1, 2, 2, 2, 3, 4, 4, 3, 4, 7, 8, 7, 4, 5, 11, 14, 14, 11, 5, 6, 16, 22, 24, 22, 16, 6, 7, 22, 32, 37, 37, 32, 22, 7, 8, 29, 44, 53, 56, 53, 44, 29, 8, 9, 37, 58, 72, 79, 79, 72, 58, 37, 9, 10, 46, 74, 94, 106, 110, 106, 94, 74, 46, 10, 11, 56, 92, 119
Offset: 0

Views

Author

Keywords

Comments

Antidiagonal sums are given by 2*A055795.
Rows/columns n are binomial transform of {n, A152947(n+1), n, 0, 0, 0, ...}.
Some primes in the array are
n = 1: {2, 7, 11, 29, 37, 67, 79, 137, 191, 211, 277, 379, ...} = A055469, primes of the form k*(k + 1)/2 + 1;
n = 3: {3, 7, 37, 53, 479, 653, 1249, 1619, 2503, 3727, 4349, 5737, 7109, 8179, 9803, 11839, 12107, ...};
n = 4: {11, 37, 79, 137, 211, 821, 991, 1597, 1831, 2081, 2347, ...} = A188382, primes of the form 8*(2*k - 1)^2 + 2*(2*k - 1) + 1.

Examples

			The array T(n,k) begins
0     1    2    3    4     5     6     7     8     9    10    11  ...
1     2    4    7   11    16    22    29    37    46    56    67  ...
2     4    8   14   22    32    44    58    74    92   112   134  ...
3     7   14   24   37    53    72    94   119   147   178   212  ...
4    11   22   37   56    79   106   137   172   211   254   301  ...
5    16   32   53   79   110   146   187   233   284   340   401  ...
6    22   44   72  106   146   192   244   302   366   436   512  ...
7    29   58   94  137   187   244   308   379   457   542   634  ...
8    37   74  119  172   233   302   379   464   557   658   767  ...
9    46   92  147  211   284   366   457   557   666   784   911  ...
10   56  112  178  254   340   436   542   658   784   920  1066  ...
11   67  134  212  301   401   512   634   767   911  1066  1232  ...
12   79  158  249  352   467   594   733   884  1047  1222  1409  ...
13   92  184  289  407   538   682   839  1009  1192  1388  1597  ...
14  106  212  332  466   614   776   952  1142  1346  1564  1796  ...
15  121  242  378  529   695   876  1072  1283  1509  1750  2006  ...
16  137  274  427  596   781   982  1199  1432  1681  1946  2227  ...
17  154  308  479  667   872  1094  1333  1589  1862  2152  2459  ...
18  172  344  534  742   968  1212  1474  1754  2052  2368  2702  ...
19  191  382  592  821  1069  1336  1622  1927  2251  2594  2956  ...
20  211  422  653  904  1175  1466  1777  2108  2459  2830  3221  ...
...
The inverse binomial transforms of the columns are
0     1    2    3    4     5     6     7     8     9    10    11  ...  A001477
1     1    2    4    7    11    22    29    37    45    56    67  ...  A152947
0     1    2    3    4     5     6     7     8     9    10    11  ...  A001477
0     0    0    0    0     0     0     0     0     0     0     0  ...
0     0    0    0    0     0     0     0     0     0     0     0  ...
0     0    0    0    0     0     0     0     0     0     0     0  ...
...
		

References

  • Miklós Bóna, Introduction to Enumerative Combinatorics, McGraw-Hill, 2007.
  • L. Comtet, Advanced Combinatorics: The Art of Finite and Infinite Expansions, Reidel Publishing Company, 1974.
  • R. P. Stanley, Enumerative Combinatorics, second edition, Cambridge University Press, 2011.

Crossrefs

Programs

  • Maple
    T := (n, k) -> n*(binomial(k, 2) + 1) + k*(binomial(n, 2) + 1);
    for n from 0 to 20 do seq(T(n, k), k = 0 .. 20) od;
  • Mathematica
    T[n_, k_] := n (Binomial[k, 2] + 1) + k (Binomial[n, 2] + 1);
    Table[T[n - k, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Jean-François Alcover, Dec 07 2018 *)
  • Maxima
    T(n, k) := n*(binomial(k, 2) + 1) + k*(binomial(n, 2) + 1)$
    for n:0 thru 20 do
      print(makelist(T(n, k), k, 0, 20));
    
  • PARI
    T(n, k) = n*(binomial(k,2) + 1) + k*(binomial(n,2) + 1);
    tabl(nn) = for (n=0, nn, for (k=0, nn, print1(T(n, k), ", ")); print); \\ Michel Marcus, Mar 12 2018

Formula

T(n,k) = T(k,n) = n*A152947(k+1) + k*A152947(n+1).
T(n,0) = A001477(n).
T(n,1) = A000124(n).
T(n,2) = A014206(n).
T(n,3) = A273465(3*n+2).
T(n,4) = A084849(n+1).
T(n,n) = A179000(n-1,n), n >= 1.
T(2*n,2*n) = 8*A081436(n-1), n >= 1.
T(2*n+1,2*n+1) = 2*A006000(2*n+1).
T(n,n+1) = A188377(n+3).
T(n,n+2) = A188377(n+2), n >= 1.
Sum_{k=0..n} T(k,n-k) = 2*(binomial(n, 4) + binomial(n, 2)).
G.f.: -((2*x*y - y - x)*(2*x*y - y - x + 1))/(((x - 1)*(y - 1))^3).
E.g.f.: (1/2)*(x + y)*(x*y + 2)*exp(x + y).

A309555 Triangle read by rows: T(n,k) = 3 + k*(n-k) for n >= 0, 0 <= k <= n.

Original entry on oeis.org

3, 3, 3, 3, 4, 3, 3, 5, 5, 3, 3, 6, 7, 6, 3, 3, 7, 9, 9, 7, 3, 3, 8, 11, 12, 11, 8, 3, 3, 9, 13, 15, 15, 13, 9, 3, 3, 10, 15, 18, 19, 18, 15, 10, 3, 3, 11, 17, 21, 23, 23, 21, 17, 11, 3, 3, 12, 19, 24, 27, 28, 27, 24, 19, 12, 3, 3, 13, 21, 27, 31, 33, 33, 31, 27, 21, 13, 3, 3, 14, 23, 30, 35, 38, 39, 38, 35, 30, 23, 14, 3
Offset: 0

Views

Author

Philip K Hotchkiss, Aug 07 2019

Keywords

Comments

The rascal triangle (A077028) can be generated by either of the rules South = (East*West+1)/North or South = East+West+1-North; this number triangle can be generated by either of the rules South = (East*West+3)/North or South = East+West+1-North.
It is more suggestive to observe that N*S-E*W = 1 or 3 in the two cases, and (N+S)-(E+W) = 1 in both cases. In fact "3" in the present definition can be replaced by any integer c, and we get a triangle of integers with N*S-E*W = c and (N+S)-(E+W) = 1. I say "suggestive", because these rules also arise in frieze patterns. - N. J. A. Sloane, Aug 28 2019

Examples

			For the row n=3: a(3,0)=3, a(3,1)=5, a(3,2)=5, a(3,3)=3, ...
For the antidiagonal r=2: T(2,0)=3, T(2,1)=5, T(2,3)=7, T(2,4)=9, ...
The triangle begins:
..............3..
............3..3..
..........3..4..3..
........3..5...5..3..
......3..6...7...6..3..
....3..7...9...9..7..3..
..3..8..11..12..11..8..3..
3..9..13..15..15..13..9..3.
...
		

Crossrefs

Programs

  • Maple
    T:= proc(n, k)
        if n<0 or k<0 or k>n then
           0;
        else
           k*(n-k)+3 ;
        end if;
        end:
    seq(seq(T(n,k), k=0..n), n=0..12);
  • Mathematica
    T[n,k]:=k(n-k)+3;T[0,0] = 3; Table[T[n,k],{n,0,12},{k,0,n}]//Flatten

Formula

By rows: a(n,k) = 3 + k(n-k), n >= 0, 0 <= k <= n.
By antidiagonals: T(r,k) = 3 + r*k, r,k >= 0.

Extensions

Missing a(50)=23 inserted by Georg Fischer, Nov 08 2021

A114219 Number triangle T(n,k) = (k-(k-1)*0^(n-k))*[k<=n].

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 2, 3, 1, 0, 1, 2, 3, 4, 1, 0, 1, 2, 3, 4, 5, 1, 0, 1, 2, 3, 4, 5, 6, 1, 0, 1, 2, 3, 4, 5, 6, 7, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1
Offset: 0

Views

Author

Paul Barry, Nov 18 2005

Keywords

Comments

Row sums are n*(n-1)/2+1 (essentially A000124). Diagonal sums are A114220. First difference triangle of A077028, when this is viewed as a number triangle.
From R. J. Mathar, Mar 22 2013: (Start)
The matrix inverse is
1;
0, 1;
0, -1, 1;
0, 1, -2, 1;
0, -2, 4, -3, 1;
0, 6, -12, 9, -4, 1;
0, -24, 48, -36, 16, -5, 1;
0, 120, -240, 180, -80, 25, -6, 1;
0, -720, 1440, -1080, 480, -150, 36, -7, 1;
... apparently related to A208058. (End)
Number of permutations of length n avoiding simultaneously the patterns 132 and 321 with k left-to-right maxima (resp., right-to-left minima). A left-to-right maximum (resp., right-to-left minimum) in a permutation p(1)p(2)...p(n) is a position i such that p(j) < p(i) for all j < i (resp., p(j) > p(i) for all j > i). - Sergey Kitaev, Nov 18 2023

Examples

			Triangle begins
  1;
  0, 1;
  0, 1, 1;
  0, 1, 2, 1;
  0, 1, 2, 3, 1;
  0, 1, 2, 3, 4, 1;
  0, 1, 2, 3, 4, 5, 1;
  0, 1, 2, 3, 4, 5, 6, 1;
  0, 1, 2, 3, 4, 5, 6, 7, 1;
  ...
		

Crossrefs

Programs

  • Maple
    A114219 := proc(n,k)
        if k < 0 or k > n then
            0;
        elif n = k then
            1;
        else
            k ;
        end if;
    end proc: # R. J. Mathar, Mar 22 2013

Formula

G.f.: (1-x-u*x + 2u*x^2)/((1-x)(1-u*x)^2), where x records length and u records left-to-right maxima (or right-to-left minima). - Sergey Kitaev, Nov 18 2023

A309557 Number triangle with T(n,k) = 2 + 3*n - 2*k + 2*k*(n-k) for n >= 0, 0 <= k <= n.

Original entry on oeis.org

2, 5, 3, 8, 8, 4, 11, 13, 11, 5, 14, 18, 18, 14, 6, 17, 23, 25, 23, 17, 7, 20, 28, 32, 32, 28, 20, 8, 23, 33, 39, 41, 39, 33, 23, 9, 26, 38, 46, 50, 50, 46, 38, 26, 10, 29, 43, 53, 59, 61, 59, 53, 43, 29, 11, 32, 48, 60, 68, 72, 72, 68, 60, 48, 32, 12, 35, 53, 67, 77, 83, 85, 83, 77, 67, 53, 35, 13
Offset: 0

Views

Author

Philip K Hotchkiss, Aug 07 2019

Keywords

Comments

The rascal triangle (A077028) can be generated by South = (East*West+1)/North or South = East+West+1-North; this triangle can be generated by South = (East*West+1)/North, South = East+West+2-North.

Examples

			For row n=3: a(3,0)=11, a(3,1)=13, a(3,2)=11, a(3,3)=5, ...
For antidiagonal r=2: T(2,0)=4, T(2,1)=11, T(2,2)=18, ...
Triangle T begins:
              2
            5   3
          8   8   4
        11  13  11  5
      14  18  18  14  6
    17  23  25  23  17  7
  20  28  32  32  28  20  8
23  33  39  41  39  33  23  9
             ...
		

Crossrefs

Programs

  • Maple
    :=proc(n,k)
       if n<0 or k<0 or k>n then
           0;
       else
           2+3*n -2*k +2*k*(n-k);
       end if;
  • Mathematica
    T[n_,k_]:=2+3*n-2*k+2*k*(n-k); Table[T[n,k], {n,0,11}, {k,0,n}] // Flatten
    f[n_] := Table[SeriesCoefficient[(x*(-1-4*y+y^2)-2*(1-4*y+y^2))/((-1+x)^2*(-1+y)^3), {x, 0, i}, {y, 0, j}], {i, n, n}, {j, 0, n}]; Flatten[Array[f, 12,0]] (* Stefano Spezia, Sep 08 2019 *)

Formula

By rows: a(n,k) = 2 + 3*n - 2*k + 2*k*(n-k) n >= 0, 0 <= k <= n.
By antidiagonals: T(r,k) = 2 + 3*k + r + 2*r*k, r,k >= 0.
G.f.: (x*(-1-4*y+y^2)-2*(1-4*y+y^2))/((-1+x)^2*(-1+y)^3). - Stefano Spezia, Sep 08 2019
Previous Showing 11-20 of 42 results. Next