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 41-50 of 124 results. Next

A003518 a(n) = 8*binomial(2*n+1,n-3)/(n+5).

Original entry on oeis.org

1, 8, 44, 208, 910, 3808, 15504, 62016, 245157, 961400, 3749460, 14567280, 56448210, 218349120, 843621600, 3257112960, 12570420330, 48507033744, 187187399448, 722477682080, 2789279908316, 10772391370048, 41620603020640, 160878516023680, 622147386185325
Offset: 3

Views

Author

Keywords

Comments

a(n-6) is the number of n-th generation nodes in the tree of sequences with unit increase labeled by 7 (cf. Zoran Sunic reference). - Benoit Cloitre, Oct 07 2003
Number of standard tableaux of shape (n+4,n-3). - Emeric Deutsch, May 30 2004

Examples

			G.f. = x^3 + 8*x^4 + 44*x^5 + 208*x^6 + 910*x^7 + 3808*x^8 + 15504*x^9 + ...
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A002057.
First differences are in A026018.
A diagonal of any of the essentially equivalent arrays A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • Magma
    [8*Binomial(2*n+1,n-3)/(n+5): n in [3..30]]; // Vincenzo Librandi, Jan 23 2017
  • Mathematica
    Table[8 Binomial[2 n + 1, n - 3]/(n + 5), {n, 3, 25}] (* Michael De Vlieger, Oct 26 2016 *)
    CoefficientList[Series[((1 - Sqrt[1 - 4 x])/(2 x))^8, {x, 0, 30}], x] (* Vincenzo Librandi, Jan 23 2017 *)
  • PARI
    {a(n) = if( n<3, 0, 8 * binomial(2*n + 1, n-3) / (n + 5))}; /* Michael Somos, Mar 14 2011 */
    
  • PARI
    my(x='x+O('x^50)); Vec(x^3*((1-(1-4*x)^(1/2))/(2*x))^8) \\ Altug Alkan, Nov 01 2015
    

Formula

G.f.: x^3*C(x)^8, where C(x)=(1-sqrt(1-4*x))/(2*x) is g.f. for the Catalan numbers (A000108). - Emeric Deutsch, May 30 2004
The convolution of A002057 with itself. - Gerald McGarvey, Nov 08 2007
Let A be the Toeplitz matrix of order n defined by: A[i,i-1]=-1, A[i,j]=Catalan(j-i), (i<=j), and A[i,j]=0, otherwise. Then, for n>=7, a(n-4)=(-1)^(n-7)*coeff(charpoly(A,x),x^7). - Milan Janjic, Jul 08 2010
a(n) = A214292(2*n,n-4) for n > 3. - Reinhard Zumkeller, Jul 12 2012
Integral representation as the n-th moment of the signed weight function W(x) on (0,4), i.e.: a(n+3) = Integral_{x=0..4} x^n*W(x) dx, n >= 0, with W(x) = (1/2)*x^(7/2)*(x-2)*(x^2-4*x+2)*sqrt(4-x)/Pi. - Karol A. Penson, Oct 26 2016
From Ilya Gutkovskiy, Jan 22 2017: (Start)
E.g.f.: 4*BesselI(4,2*x)*exp(2*x)/x.
a(n) ~ 4^(n+2)/(sqrt(Pi)*n^(3/2)). (End)
D-finite with recurrence: -(n+5)*(n-3)*a(n) +2*n*(2*n+1)*a(n-1)=0. - R. J. Mathar, Feb 20 2020
From Amiram Eldar, Jan 02 2022: (Start)
Sum_{n>=3} 1/a(n) = 43*Pi/(36*sqrt(3)) - 81/80.
Sum_{n>=3} (-1)^(n+1)/a(n) = 6213*log(phi)/(50*sqrt(5)) - 10339/400, where phi is the golden ratio (A001622). (End)

Extensions

More terms from Jon E. Schoenfield, May 06 2010

A030237 Catalan's triangle with right border removed (n > 0, 0 <= k < n).

Original entry on oeis.org

1, 1, 2, 1, 3, 5, 1, 4, 9, 14, 1, 5, 14, 28, 42, 1, 6, 20, 48, 90, 132, 1, 7, 27, 75, 165, 297, 429, 1, 8, 35, 110, 275, 572, 1001, 1430, 1, 9, 44, 154, 429, 1001, 2002, 3432, 4862, 1, 10, 54, 208, 637, 1638, 3640, 7072, 11934, 16796, 1, 11, 65, 273, 910, 2548, 6188, 13260, 25194, 41990, 58786
Offset: 1

Views

Author

Keywords

Comments

This triangle appears in the totally asymmetric exclusion process as Y(alpha=1,beta=1,n,m), written in the Derrida et al. reference as Y_n(m) for alpha=1, beta=1. - Wolfdieter Lang, Jan 13 2006

Examples

			Triangle begins as:
  1;
  1, 2;
  1, 3,  5;
  1, 4,  9,  14;
  1, 5, 14,  28,  42;
  1, 6, 20,  48,  90,  132;
  1, 7, 27,  75, 165,  297,  429;
  1, 8, 35, 110, 275,  572, 1001, 1430;
  1, 9, 44, 154, 429, 1001, 2002, 3432, 4862;
		

Crossrefs

Alternate versions of (essentially) the same Catalan triangle: A009766, A033184, A047072, A059365, A099039, A106566, A130020.
Row sums give A071724.

Programs

  • Haskell
    a030237 n k = a030237_tabl !! n !! k
    a030237_row n = a030237_tabl !! n
    a030237_tabl = map init $ tail a009766_tabl
    -- Reinhard Zumkeller, Jul 12 2012
    
  • Magma
    [(n-k+1)*Binomial(n+k, k)/(n+1): k in [0..n-1], n in [1..12]]; // G. C. Greubel, Mar 17 2021
  • Maple
    A030237 := proc(n,m)
        (n-m+1)*binomial(n+m,m)/(n+1) ;
    end proc: # R. J. Mathar, May 31 2016
    # Compare the analogue algorithm for the Bell numbers in A011971.
    CatalanTriangle := proc(len) local P, T, n; P := [1]; T := [[1]];
    for n from 1 to len-1 do P := ListTools:-PartialSums([op(P), P[-1]]);
    T := [op(T), P] od; T end: CatalanTriangle(6):
    ListTools:-Flatten(%); # Peter Luschny, Mar 26 2022
    # Alternative:
    ogf := n -> (1 - 2*x)/(1 - x)^(n + 2):
    ser := n -> series(ogf(n), x, n):
    row := n -> seq(coeff(ser(n), x, k), k = 0..n-1):
    seq(row(n), n = 1..11); # Peter Luschny, Mar 27 2022
  • Mathematica
    T[n_, k_]:= T[n, k] = Which[k==0, 1, k>n, 0, True, T[n-1, k] + T[n, k-1]];
    Table[T[n, k], {n,1,12}, {k,0,n-1}] // Flatten (* Jean-François Alcover, Nov 14 2017 *)
  • PARI
    T(n,k) = (n-k+1)*binomial(n+k, k)/(n+1) \\ Andrew Howroyd, Feb 23 2018
    
  • Sage
    flatten([[(n-k+1)*binomial(n+k, k)/(n+1) for k in (0..n-1)] for n in (1..12)]) # G. C. Greubel, Mar 17 2021
    

Formula

T(n, k) = (n-k+1)*binomial(n+k, k)/(n+1).
Sum_{k=0..n-1} T(n,k) = A000245(n). - G. C. Greubel, Mar 17 2021
T(n, k) = [x^k] ((1 - 2*x)/(1 - x)^(n + 2)). - Peter Luschny, Mar 27 2022

Extensions

Missing a(8) = T(7,0) = 1 inserted by Reinhard Zumkeller, Jul 12 2012

A115728 Number of subpartitions of partitions in Abramowitz and Stegun order.

Original entry on oeis.org

1, 2, 3, 3, 4, 5, 4, 5, 7, 6, 7, 5, 6, 9, 9, 10, 9, 9, 6, 7, 11, 12, 13, 10, 14, 13, 10, 12, 11, 7, 8, 13, 15, 16, 14, 19, 17, 16, 16, 19, 16, 14, 15, 13, 8, 9, 15, 18, 19, 18, 24, 21, 15, 23, 22, 26, 21, 19, 22, 23, 24, 19, 15, 18, 18, 15, 9, 10, 17, 21, 22, 22, 29
Offset: 0

Views

Author

Keywords

Comments

subpart([n^k]) = C(n+k,k); subpart([1,2,3,...,n]) = C_n = A000108(n). The b(i,j) defined in the formula for sequences [1,2,3,...] form A009766.
Row sums are A297388. Row lengths are A000041. - Geoffrey Critzer, Jan 10 2021

Examples

			Partition 5 in A&S order is [2,1]; it has 5 subpartitions: [], [1], [2], [1^2] and [2,1] itself.
1
2
3, 3
4, 5, 4
5, 7, 6, 7,  5
6, 9, 9, 10, 9, 9, 6
		

Crossrefs

Programs

  • PARI
    /* Expects input as vector in increasing order - e.g. [1,1,2,3] */
    subpart(p)=local(i,j,v,n);n=matsize(p)[2];if(n==0,1,v=vector(p[n]+1);v[1] =1;for(i=1,n,for(j=1,p[i],v[j+1]+=v[j]));for(j=1,p[n],v[j+1]+=v[j]);v[p[n ]+1])
    
  • PARI
    /* Given Partition p(), Find Subpartitions s(): */ {s(n)=polcoeff(x^n-sum(k=0, n-1, s(k)*x^k*(1-x+x*O(x^n))^p(k)),n)} \\ Paul D. Hanna, Jul 03 2006

Formula

For a partition P = [p_1,...,p_n] with the p_i in increasing order, define b(i,j) to be the number of subpartitions of [p_1,...,p_i] with the i-th part = j (b(i,0) is subpartitions with less than i parts). Then b(1,j)=1 for j<=p_1, b(i+1,j) = Sum_{k=0..j} b(i,k) for 0<=j<=p_{i+1}; and the total number of subpartitions is sum_{k=1..p_n} b(n,k).
For a partition P = {p(n)}, the number of subpartitions {s(n)} of P can be determined by the g.f.: 1/(1-x) = Sum_{n>=0} s(n)*x^n*(1-x)^p(n). - Paul D. Hanna, Jul 03 2006

A193091 Augmentation of the triangular array A158405. See Comments.

Original entry on oeis.org

1, 1, 3, 1, 6, 14, 1, 9, 37, 79, 1, 12, 69, 242, 494, 1, 15, 110, 516, 1658, 3294, 1, 18, 160, 928, 3870, 11764, 22952, 1, 21, 219, 1505, 7589, 29307, 85741, 165127, 1, 24, 287, 2274, 13355, 61332, 224357, 638250, 1217270, 1, 27, 364, 3262, 21789, 115003
Offset: 0

Views

Author

Clark Kimberling, Jul 30 2011

Keywords

Comments

Suppose that P is an infinite triangular array of numbers:
p(0,0)
p(1,0)...p(1,1)
p(2,0)...p(2,1)...p(2,2)
p(3,0)...p(3,1)...p(3,2)...p(3,3)...
...
Let w(0,0)=1, w(1,0)=p(1,0), w(1,1)=p(1,1), and define
W(n)=(w(n,0), w(n,1), w(n,2),...w(n,n-1), w(n,n)) recursively by W(n)=W(n-1)*PP(n), where PP(n) is the n X (n+1) matrix given by
...
row 0 ... p(n,0) ... p(n,1) ...... p(n,n-1) ... p(n,n)
row 1 ... 0 ..... p(n-1,0) ..... p(n-1,n-2) .. p(n-1,n-1)
row 2 ... 0 ..... 0 ............ p(n-2,n-3) .. p(n-2,n-2)
...
row n-1 . 0 ..... 0 ............. p(2,1) ..... p(2,2)
row n ... 0 ..... 0 ............. p(1,0) ..... p(1,1)
...
The augmentation of P is here introduced as the triangular array whose n-th row is W(n), for n>=0. The array P may be represented as a sequence of polynomials; viz., row n is then the vector of coefficients: p(n,0), p(n,1),...,p(n,n), from p(n,0)*x^n+p(n,1)*x^(n-1)+...+p(n,n). For example, (C(n,k)) is represented by ((x+1)^n); using this choice of P (that is, Pascal's triangle), the augmentation of P is calculated one row at a time, either by the above matrix products or by polynomial substitutions in the following manner:
...
row 0 of W: 1, by decree
row 1 of W: 1 augments to 1,1
...polynomial version: 1 -> x+1
row 2 of W: 1,1 augments to 1,3,2
...polynomial version: x+1 -> (x^2+2x+1)+(x+1)=x^2+3x+2
row 3 to W: 1,3,2 augments to 1,6,11,6
...polynomial version:
x^2+3x+2 -> (x+1)^3+3(x+1)^2+2(x+1)=(x+1)(x+2)(x+3)
...
Examples of augmented triangular arrays:
(p(n,k)=1) augments to A009766, Catalan triangle.
Catalan triangle augments to A193560.
Pascal triangle augments to A094638, Stirling triangle.
A002260=((k+1)) augments to A023531.
A154325 augments to A033878.
A158405 augments to A193091.
((k!)) augments to A193092.
A094727 augments to A193093.
A130296 augments to A193094.
A004736 augments to A193561.
...
Regarding the specific augmentation W=A193091: w(n,n)=A003169.
From Peter Bala, Aug 02 2012: (Start)
This is the table of g(n,k) in the notation of Carlitz (p. 124). The triangle enumerates two-line arrays of positive integers
............a_1 a_2 ... a_n..........
............b_1 b_2 ... b_n..........
such that
1) max(a_i, b_i) <= min(a_(i+1), b_(i+1)) for 1 <= i <= n-1
2) max(a_i, b_i) <= i for 1 <= i <= n
3) max(a_n, b_n) = k.
See A071948 and A211788 for other two-line array enumerations.
(End)

Examples

			The triangle P, at A158405, is given by rows
  1
  1...3
  1...3...5
  1...3...5...7
  1...3...5...7...9...
The augmentation of P is the array W starts with w(0,0)=1, by definition of W.
Successive polynomials (rows of W) arise from P as shown here:
  ...
  1->x+3, so that W has (row 1)=(1,3);
  ...
  x+3->(x^2+3x+5)+3*(x+3), so that W has (row 2)=(1,6,14);
  ...
  x^2+6x+14->(x^3+3x^2+5x+7)+6(x^2+3x+5)+14(x+3), so that (row 3)=(1,9,37,79).
  ...
First 7 rows of W:
  1
  1    3
  1    6    14
  1    9    37    79
  1   12    69   242    494
  1   15   110   516   1658    3294
  1   18   160   928   3870   11764   22952
		

Crossrefs

Programs

  • Mathematica
    p[n_, k_] := 2 k + 1
    Table[p[n, k], {n, 0, 5}, {k, 0, n}] (* A158405 *)
    m[n_] := Table[If[i <= j, p[n + 1 - i, j - i], 0], {i, n}, {j, n + 1}]
    TableForm[m[4]]
    w[0, 0] = 1; w[1, 0] = p[1, 0]; w[1, 1] = p[1, 1];
    v[0] = w[0, 0]; v[1] = {w[1, 0], w[1, 1]};
    v[n_] := v[n - 1].m[n]
    TableForm[Table[v[n], {n, 0, 6}]] (* A193091 *)
    Flatten[Table[v[n], {n, 0, 9}]]

Formula

From Peter Bala, Aug 02 2012: (Start)
T(n,k) = (n-k+1)/n*Sum_{i=0..k} C(n+1,n-k+i+1)*C(2*n+i+1,i) for 0 <= k <= n.
Recurrence equation: T(n,k) = Sum_{i=0..k} (2*k-2*i+1)*T(n-1,i).
(End)

A003519 a(n) = 10*C(2n+1, n-4)/(n+6).

Original entry on oeis.org

1, 10, 65, 350, 1700, 7752, 33915, 144210, 600875, 2466750, 10015005, 40320150, 161280600, 641886000, 2544619500, 10056336264, 39645171810, 155989499540, 612815891050, 2404551645100, 9425842448792, 36921502679600, 144539291740025, 565588532895750, 2212449261033375
Offset: 4

Views

Author

Keywords

Comments

Number of standard tableaux of shape (n+5,n-4). - Emeric Deutsch, May 30 2004
a(n) is the number of North-East paths from (0,0) to (n,n) that cross the diagonal y = x horizontally exactly twice. By symmetry, it is also the number of North-East paths from (0,0) to (n,n) that cross the diagonal y = x vertically exactly twice. Details can be found in Section 3.3 in Pan and Remmel's link. - Ran Pan, Feb 02 2016

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

A diagonal of any of the essentially equivalent arrays A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • Magma
    [10*Binomial(2*n+1, n-4)/(n+6): n in [4..35]]; // Vincenzo Librandi, Feb 03 2016
  • Maple
    seq(10*binomial(2*n+1,n-4)/(n+6), n=4..50); # Robert Israel, Feb 02 2016
  • Mathematica
    Table[10 Binomial[2 n + 1, n - 4]/(n + 6), {n, 4, 28}] (* Michael De Vlieger, Feb 03 2016 *)
  • PARI
    a(n) = 10*binomial(2*n+1, n-4)/(n+6); \\ Michel Marcus, Feb 02 2016
    

Formula

G.f.: x^4*C(x)^10, where C(x)=[1-sqrt(1-4x)]/(2x) is g.f. for the Catalan numbers (A000108). - Emeric Deutsch, May 30 2004
Let A be the Toeplitz matrix of order n defined by: A[i,i-1]=-1, A[i,j]=Catalan(j-i), (i<=j), and A[i,j]=0, otherwise. Then, for n>=9, a(n-5)=(-1)^(n-9)*coeff(charpoly(A,x),x^9). [Milan Janjic, Jul 08 2010]
a(n) = A214292(2*n,n-5) for n > 4. - Reinhard Zumkeller, Jul 12 2012
From Robert Israel, Feb 02 2016: (Start)
D-finite with recurrence a(n+1) = 2*(n+1)*(2n+3)/((n+7)*(n-3)) * a(n).
a(n) ~ 20 * 4^n/sqrt(Pi*n^3). (End)
E.g.f.: 5*BesselI(5,2*x)*exp(2*x)/x. - Ilya Gutkovskiy, Jan 23 2017
From Amiram Eldar, Jan 02 2022: (Start)
Sum_{n>=4} 1/a(n) = 34*Pi/(45*sqrt(3)) - 44/175.
Sum_{n>=4} (-1)^n/a(n) = 53004*log(phi)/(125*sqrt(5)) - 79048/875, where phi is the golden ratio (A001622). (End)

A122848 Exponential Riordan array (1, x(1+x/2)).

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 3, 1, 0, 0, 3, 6, 1, 0, 0, 0, 15, 10, 1, 0, 0, 0, 15, 45, 15, 1, 0, 0, 0, 0, 105, 105, 21, 1, 0, 0, 0, 0, 105, 420, 210, 28, 1, 0, 0, 0, 0, 0, 945, 1260, 378, 36, 1, 0, 0, 0, 0, 0, 945, 4725, 3150, 630, 45, 1, 0, 0, 0, 0, 0, 0, 10395, 17325, 6930, 990, 55, 1, 0, 0
Offset: 0

Views

Author

Paul Barry, Sep 14 2006

Keywords

Comments

Entries are Bessel polynomial coefficients. Row sums are A000085. Diagonal sums are A122849. Inverse is A122850. Product of A007318 and A122848 gives A100862.
T(n,k) is the number of self-inverse permutations of {1,2,...,n} having exactly k cycles. - Geoffrey Critzer, May 08 2012
Bessel numbers of the second kind. For relations to the Hermite polynomials and the Catalan (A033184 and A009766) and Fibonacci (A011973, A098925, and A092865) matrices, see Yang and Qiao. - Tom Copeland, Dec 18 2013.
Also the inverse Bell transform of the double factorial of odd numbers Product_{k= 0..n-1} (2*k+1) (A001147). For the definition of the Bell transform see A264428 and for cross-references A265604. - Peter Luschny, Dec 31 2015

Examples

			Triangle begins:
    1
    0    1
    0    1    1
    0    0    3    1
    0    0    3    6    1
    0    0    0   15   10    1
    0    0    0   15   45   15    1
    0    0    0    0  105  105   21    1
    0    0    0    0  105  420  210   28    1
    0    0    0    0    0  945 1260  378   36    1
From _Gus Wiseman_, Jan 12 2021: (Start)
As noted above, a(n) is the number of set partitions of {1..n} into k singletons or pairs. This is also the number of set partitions of subsets of {1..n} into n - k pairs. In the first case, row n = 5 counts the following set partitions:
  {{1},{2,3},{4,5}}  {{1},{2},{3},{4,5}}  {{1},{2},{3},{4},{5}}
  {{1,2},{3},{4,5}}  {{1},{2},{3,4},{5}}
  {{1,2},{3,4},{5}}  {{1},{2,3},{4},{5}}
  {{1,2},{3,5},{4}}  {{1,2},{3},{4},{5}}
  {{1},{2,4},{3,5}}  {{1},{2},{3,5},{4}}
  {{1},{2,5},{3,4}}  {{1},{2,4},{3},{5}}
  {{1,3},{2},{4,5}}  {{1},{2,5},{3},{4}}
  {{1,3},{2,4},{5}}  {{1,3},{2},{4},{5}}
  {{1,3},{2,5},{4}}  {{1,4},{2},{3},{5}}
  {{1,4},{2},{3,5}}  {{1,5},{2},{3},{4}}
  {{1,4},{2,3},{5}}
  {{1,4},{2,5},{3}}
  {{1,5},{2},{3,4}}
  {{1,5},{2,3},{4}}
  {{1,5},{2,4},{3}}
In the second case, we have:
  {{1,2},{3,4}}  {{1,2}}  {}
  {{1,2},{3,5}}  {{1,3}}
  {{1,2},{4,5}}  {{1,4}}
  {{1,3},{2,4}}  {{1,5}}
  {{1,3},{2,5}}  {{2,3}}
  {{1,3},{4,5}}  {{2,4}}
  {{1,4},{2,3}}  {{2,5}}
  {{1,4},{2,5}}  {{3,4}}
  {{1,4},{3,5}}  {{3,5}}
  {{1,5},{2,3}}  {{4,5}}
  {{1,5},{2,4}}
  {{1,5},{3,4}}
  {{2,3},{4,5}}
  {{2,4},{3,5}}
  {{2,5},{3,4}}
(End)
		

Crossrefs

Row sums are A000085.
Column sums are A001515.
Same as A049403 but with a first column k = 0.
The same set partitions counted by number of pairs are A100861.
Reversing rows gives A111924 (without column k = 0).
A047884 counts standard Young tableaux by size and greatest row length.
A238123 counts standard Young tableaux by size and least row length.
A320663/A339888 count unlabeled multiset partitions into singletons/pairs.
A322661 counts labeled covering half-loop-graphs.
A339742 counts factorizations into distinct primes or squarefree semiprimes.

Programs

  • Maple
    # The function BellMatrix is defined in A264428.
    BellMatrix(n -> `if`(n<2,1,0), 9); # Peter Luschny, Jan 27 2016
  • Mathematica
    t[n_, k_] := k!*Binomial[n, k]/((2 k - n)!*2^(n - k)); Table[ t[n, k], {n, 0, 11}, {k, 0, n}] // Flatten
    (* Second program: *)
    rows = 12;
    t = Join[{1, 1}, Table[0, rows]];
    T[n_, k_] := BellY[n, k, t];
    Table[T[n, k], {n, 0, rows}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 23 2018,after Peter Luschny *)
    sbs[{}]:={{}};sbs[set:{i_,_}]:=Join@@Function[s,(Prepend[#1,s]&)/@sbs[Complement[set,s]]]/@Cases[Subsets[set],{i}|{i,_}];
    Table[Length[Select[sbs[Range[n]],Length[#]==k&]],{n,0,6},{k,0,n}] (* Gus Wiseman, Jan 12 2021 *)
  • PARI
    {T(n,k)=if(2*kn, 0, n!/(2*k-n)!/(n-k)!*2^(k-n))} /* Michael Somos, Oct 03 2006 */
    
  • Sage
    # uses[inverse_bell_transform from A265605]
    multifact_2_1 = lambda n: prod(2*k + 1 for k in (0..n-1))
    inverse_bell_matrix(multifact_2_1, 9) # Peter Luschny, Dec 31 2015

Formula

Number triangle T(n,k) = k!*C(n,k)/((2k-n)!*2^(n-k)).
T(n,k) = A001498(k,n-k). - Michael Somos, Oct 03 2006
E.g.f.: exp(y(x+x^2/2)). - Geoffrey Critzer, May 08 2012
Triangle equals the matrix product A008275*A039755. Equivalently, the n-th row polynomial R(n,x) is given by the Type B Dobinski formula R(n,x) = exp(-x/2)*Sum_{k>=0} P(n,2*k+1)*(x/2)^k/k!, where P(n,x) = x*(x-1)*...*(x-n+1) denotes the falling factorial polynomial. Cf. A113278. - Peter Bala, Jun 23 2014
From Daniel Checa, Aug 28 2022: (Start)
E.g.f. for the m-th column: (x^2/2+x)^m/m!.
T(n,k) = T(n-1,k-1) + (n-1)*T(n-2,k-1) for n>1 and k=1..n, T(0,0) = 1. (End)

A008282 Triangle of Euler-Bernoulli or Entringer numbers read by rows: T(n,k) is the number of down-up permutations of n+1 starting with k+1.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 2, 4, 5, 5, 5, 10, 14, 16, 16, 16, 32, 46, 56, 61, 61, 61, 122, 178, 224, 256, 272, 272, 272, 544, 800, 1024, 1202, 1324, 1385, 1385, 1385, 2770, 4094, 5296, 6320, 7120, 7664, 7936, 7936, 7936, 15872, 23536, 30656, 36976, 42272, 46366, 49136, 50521, 50521
Offset: 1

Views

Author

Keywords

Examples

			Triangle T(n,k) (with rows n >= 1 and columns k = 1..n) begins
   1
   1  1
   1  2  2
   2  4  5  5
   5 10 14 16 16
  16 32 46 56 61 61
  ...
Each row is constructed by forming the partial sums of the previous row, reading from the right and repeating the final term.
T(4,3) = 5 because we have 41325, 41523, 42314, 42513 and 43512. All these permutations have length n+1 = 5, start with k+1 = 4, and they are down-up permutations.
		

References

  • R. C. Entringer, A combinatorial interpretation of the Euler and Bernoulli numbers, Nieuw Archief voor Wiskunde, 14 (1966), 241-246.

Crossrefs

Programs

  • Haskell
    a008282 n k = a008282_tabl !! (n-1) !! (k-1)
    a008282_row n = a008282_tabl !! (n-1)
    a008282_tabl = iterate f [1] where
       f xs = zs ++ [last zs] where zs = scanl1 (+) (reverse xs)
    -- Reinhard Zumkeller, Dec 28 2011
  • Maple
    f:=series(sec(x)+tan(x),x=0,25): E[0]:=1: for n from 1 to 20 do E[n]:=n!*coeff(f,x^n) od: T:=proc(n,k) if kPeter Luschny, Aug 03 2017
    # Third program:
    T := proc(n, k) local w: if 0 = n mod 2 then w := coeftayl(cos(x)/cos(x + y), [x, y] = [0, 0], [n - k, k]): end if: if 1 = n mod 2 then w := coeftayl(sin(x)/cos(x + y), [x, y] = [0, 0], [k, n - k]): end if: w*(n - k)!*k!: end proc:
    for n from 1 to 6 do seq(T(n,k), k=1..n) od; # Petros Hadjicostas, Feb 17 2021
  • Mathematica
    ro[1] = {1}; ro[n_] := ro[n] = (s = Accumulate[ Reverse[ ro[n-1]]]; Append[ s, Last[s]]); Flatten[ Table[ ro[n], {n, 1, 10}]] (* Jean-François Alcover, Oct 03 2011 *)
    nxt[lst_]:=Module[{lst2=Accumulate[Reverse[lst]]},Flatten[Join[ {lst2,Last[ lst2]}]]]; Flatten[NestList[nxt,{1},10]] (* Harvey P. Dale, Aug 17 2014 *)

Formula

From Emeric Deutsch, May 15 2004: (Start)
Let E[j] = A000111(j) = j! * [x^j](sec(x) + tan(x)) be the up/down or Euler numbers. For 1 <= k < n,
T(n, k) = Sum_{i=0..floor((k-1)/2)} (-1)^i * binomial(k, 2*i+1) * E[n-2*i-1];
T(n,k) = Sum_{i=0..floor((n-k)/2)} (-1)^i * binomial(n-k, 2*i) * E[n-2*i];
T(n, k) = Sum_{i=0..floor((n-k)/2)} (-1)^i * binomial(n-k, 2*i) * E[n-2*i]; and
T(n, n) = E[n] for n >= 1. (End)
From Petros Hadjicostas, Feb 17 2021: (Start)
If n is even, then T(n,k) = k!*(n-k)!*[x^(n-k),y^k] cos(x)/cos(x + y).
If n is odd, then T(n,k) = k!*(n-k)!*[x^k,y^(n-k)] sin(x)/cos(x + y).
(These were adapted and corrected from the formulas in Corollary 1.3 in Foata and Guo-Niu Han (2014).) (End)
Comment from Masanobu Kaneko: (Start)
A generating function that applies for all n, both even and odd:
Sum_{n=0..oo} Sum_{k=0..n} T(n,k) x^(n-k)/(n-k)! * y^k/k! = {cos x + sin y}/cos(x + y).
(End) - N. J. A. Sloane, Feb 06 2022

Extensions

Example and Formula sections edited by Petros Hadjicostas, Feb 17 2021

A100100 Triangle T(n,k) = binomial(2*n-k-1, n-k) read by rows.

Original entry on oeis.org

1, 1, 1, 3, 2, 1, 10, 6, 3, 1, 35, 20, 10, 4, 1, 126, 70, 35, 15, 5, 1, 462, 252, 126, 56, 21, 6, 1, 1716, 924, 462, 210, 84, 28, 7, 1, 6435, 3432, 1716, 792, 330, 120, 36, 8, 1, 24310, 12870, 6435, 3003, 1287, 495, 165, 45, 9, 1, 92378, 48620, 24310, 11440, 5005, 2002
Offset: 0

Views

Author

Paul Barry, Nov 08 2004

Keywords

Comments

Sometimes called a Catalan triangle, although there are many other triangles that carry that name - see A009766, A008315, A028364, A033184, A053121, A059365, A062103.
Number of nodes of outdegree k in all ordered trees with n edges. Equivalently, number of ascents of length k in all Dyck paths of semilength n. Example: T(3,2) = 3 because the Dyck paths of semilength 3 are UDUDUD, UD(UU)DD, (UU)DDUD, (UU)DUDD and UUUDDD, where U = (1,1), D = (1,-1), the ascents of length 2 being shown between parentheses. - Emeric Deutsch, Nov 19 2006
Riordan array (f(x), x*g(x)) where f(x) is the g.f. of A088218 and g(x) is the g.f. of A000108. - Philippe Deléham, Jan 23 2010
T(n,k) is the number of nonnegative paths of upsteps U = (1,1) and downsteps D = (1,-1) of length 2*n with k returns to ground level, the horizontal line through the initial vertex. Example: T(2,1) = 2 counts UDUU, UUDD. Also, T(n,k) = number of these paths whose last descent has length k, that is, k downsteps follow the last upstep. Example: T(2,1) = 2 counts UUUD, UDUD. - David Callan, Nov 21 2011
Belongs to the hitting-time subgroup of the Riordan group. Multiplying this triangle by the square Pascal matrix gives A092392 read as a square array. See the example below. - Peter Bala, Nov 03 2015

Examples

			From _Paul Barry_, Mar 15 2010: (Start)
Triangle begins in row n=0 with columns 0<=k<=n as:
    1;
    1,   1;
    3,   2,   1;
   10,   6,   3,  1;
   35,  20,  10,  4,  1;
  126,  70,  35, 15,  5, 1;
  462, 252, 126, 56, 21, 6, 1;
Production matrix begins
  1, 1;
  2, 1, 1;
  3, 1, 1, 1;
  4, 1, 1, 1, 1;
  5, 1, 1, 1, 1, 1;
  6, 1, 1, 1, 1, 1, 1;
  7, 1, 1, 1, 1, 1, 1, 1;
(End)
A092392 as a square array = A100100 * square Pascal matrix:
/1   1  1  1 ...\   / 1          \/1 1  1  1 ...\
|2   3  4  5 ...|   | 1 1        ||1 2  3  4 ...|
|6  10 15 21 ...| = | 3 2 1      ||1 3  6 10 ...|
|20 35 56 84 ...|   |10 6 3 1    ||1 4 10 20 ...|
|70 ...         |   |35 ...      ||1 ...        |
- _Peter Bala_, Nov 03 2015
		

Crossrefs

Row sums are A000984. Equivalent to A092392, to which A088218 has been added as a first column. Columns include A088218, A000984, A001700, A001791, A002054, A002694. Diagonal sums are A100217. Matrix inverse is A100218.
Cf. A059481 (mirrored). Cf. A033184, A094527, A113955.

Programs

  • Haskell
    a100100 n k = a100100_tabl !! n !! n
    a100100_row n = a100100_tabl !! n
    a100100_tabl = [1] : f a092392_tabl where
       f (us : wss'@(vs : wss)) = (vs !! 1 : us) : f wss'
    -- Reinhard Zumkeller, Jan 15 2014
    
  • Magma
    /* As triangle */ [[Binomial(2*n - k - 1, n - k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Nov 21 2018
  • Maple
    A100100 := proc(n,k)
        binomial(2*n-k-1,n-1) ;
    end proc:
    seq(seq(A100100(n,k),k=0..n),n=0..10) ; # R. J. Mathar, Feb 06 2015
  • Mathematica
    Flatten[Table[Binomial[2 n - k - 1, n - k], {n, 0, 11}, {k, 0, n}]] (* Vincenzo Librandi, Nov 21 2018 *)
  • PARI
    T(n,k)=binomial(2*n-k-1,n-k) \\ Charles R Greathouse IV, Jan 16 2012
    

Formula

From Peter Bala, Sep 06 2015: (Start)
Matrix product A094527 * P^(-1) = A113955 * P^(-2), where P denotes Pascal's triangle A007318.
Essentially, the logarithmic derivative of A033184. (End)
Let column(k) = [T(n, k), n >= k], then the generating function for column(k) is (2/(sqrt(1-4*x)+1))^(k-1)/sqrt(1-4*x). - Peter Luschny, Mar 19 2021
O.g.f. row polynomials R(n, x) := Sum_{k=0..n} T(n, k)*x^k, i.e. o.g.f. of the triangle, G(z,x) = 1/((2 - c(z))*(1 - x*z*c(z))), with c the o.g.f. of A000108 (Catalan). See the Riordan coomment by Philippe Deléham above. - Wolfdieter Lang, Apr 06 2021

A047072 Array A read by diagonals: A(h,k)=number of paths consisting of steps from (0,0) to (h,k) such that each step has length 1 directed up or right and no step touches the line y=x unless x=0 or x=h.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 3, 2, 2, 3, 1, 1, 4, 5, 4, 5, 4, 1, 1, 5, 9, 5, 5, 9, 5, 1, 1, 6, 14, 14, 10, 14, 14, 6, 1, 1, 7, 20, 28, 14, 14, 28, 20, 7, 1, 1, 8, 27, 48, 42, 28, 42, 48, 27, 8, 1, 1, 9, 35, 75, 90, 42, 42, 90, 75, 35, 9, 1
Offset: 0

Views

Author

Keywords

Examples

			Array, A(n, k), begins as:
  1, 1,  1,  1,  1,   1,   1,   1, ...;
  1, 2,  1,  2,  3,   4,   5,   6, ...;
  1, 1,  2,  2,  5,   9,  14,  20, ...;
  1, 2,  2,  4,  5,  14,  28,  48, ...;
  1, 3,  5,  5, 10,  14,  42,  90, ...;
  1, 4,  9, 14, 14,  28,  42, 132, ...;
  1, 5, 14, 28, 42,  42,  84, 132, ...;
  1, 6, 20, 48, 90, 132, 132, 264, ...;
Antidiagonals, T(n, k), begins as:
  1;
  1,  1;
  1,  2,  1;
  1,  1,  1,  1;
  1,  2,  2,  2,  1;
  1,  3,  2,  2,  3,  1;
  1,  4,  5,  4,  5,  4,  1;
  1,  5,  9,  5,  5,  9,  5,  1;
  1,  6, 14, 14, 10, 14, 14,  6,  1;
		

Crossrefs

The following are all versions of (essentially) the same Catalan triangle: A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • Magma
    b:= func< n | n eq 0 select 1 else 2*Catalan(n-1) >;
    function A(n,k)
      if k eq n then return b(n);
      elif k gt n then return Binomial(n+k-1, n) - Binomial(n+k-1, n-1);
      else return Binomial(n+k-1, k) - Binomial(n+k-1, k-1);
      end if; return A;
    end function;
    // [[A(n,k): k in [0..12]]: n in [0..12]];
    T:= func< n,k | A(n-k, k) >;
    [T(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Oct 13 2022
    
  • Mathematica
    A[, 0]= 1; A[0, ]= 1; A[h_, k_]:= A[h, k]= If[(k-1>h || k-1Jean-François Alcover, Mar 06 2019 *)
  • SageMath
    def A(n,k):
        if (k==n): return 2*catalan_number(n-1) + 2*int(n==0)
        elif (k>n): return binomial(n+k-1, n) - binomial(n+k-1, n-1)
        else: return binomial(n+k-1, k) - binomial(n+k-1, k-1)
    def T(n,k): return A(n-k, k)
    # [[A(n,k) for k in range(12)] for n in range(12)]
    flatten([[T(n,k) for k in range(n+1)] for n in range(12)]) # G. C. Greubel, Oct 13 2022

Formula

A(n, n) = 2*[n=0] - A002420(n),
A(n, n+1) = 2*A000108(n-1), n >= 1.
From G. C. Greubel, Oct 13 2022: (Start)
T(n, n-1) = A000027(n-2) + 2*[n<3], n >= 1.
T(n, n-2) = A000096(n-4) + 2*[n<5], n >= 2.
T(n, n-3) = A005586(n-6) + 4*[n<7] - 2*[n=3], n >= 3.
T(2*n, n) = 2*A000108(n-1) + 3*[n=0].
T(2*n-1, n-1) = T(2*n+1, n+1) = A000180(n).
T(3*n, n) = A025174(n) + [n=0]
Sum_{k=0..n} T(n, k) = 2*A063886(n-2) + [n=0] - 2*[n=1]
Sum_{k=0..n} (-1)^k * T(n, k) = A000007(n).
Sum_{k=0..floor(n/2)} T(n, k) = A047079(n). (End)

A005586 a(n) = n*(n+4)*(n+5)/6.

Original entry on oeis.org

0, 5, 14, 28, 48, 75, 110, 154, 208, 273, 350, 440, 544, 663, 798, 950, 1120, 1309, 1518, 1748, 2000, 2275, 2574, 2898, 3248, 3625, 4030, 4464, 4928, 5423, 5950, 6510, 7104, 7733, 8398, 9100, 9840, 10619, 11438, 12298, 13200, 14145, 15134, 16168, 17248
Offset: 0

Views

Author

Keywords

Comments

Number of walks on square lattice.
Number of standard tableaux of shape (n+2,3) (n >= 1). - Emeric Deutsch, May 20 2004
Number of left factors of Dyck paths from (0,0) to (n+5,n-1). E.g. a(1)=5 because we have UDUDUD, UDUUDD, UUDDUD, UUDUDD and UUUDDD, where U=(1,1) and D=(1,-1). - Emeric Deutsch, Jan 25 2005
Column 4 of Catalan triangle A009766. - Zerinvary Lajos, Nov 25 2006
Sum of first n triangular numbers minus next triangular number. - Vladimir Joseph Stephan Orlovsky, Oct 13 2009
Number of packed increasing tableaux of shape 3 X (n+1) with alphabet [n+4]. - Oliver Pechenik, Jan 03 2022

Examples

			G.f. = 5*x + 14*x^2 + 28*x^3 + 48*x^4 + 75*x^5 + 110*x^6 + 154*x^7 + ...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 796.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n)=A053121(n+5,n-1).

Programs

  • Magma
    [n*(n+4)*(n+5)/6: n in [0..50]]; // Vincenzo Librandi, Jun 09 2013
  • Maple
    [seq(binomial(n,3 )-binomial(n,1),n=4..48)]; # Zerinvary Lajos, Nov 25 2006
    a:=n->sum ((j-3)*j/2,j=0..n): seq(a(n),n=4..48); # Zerinvary Lajos, Dec 17 2006
    A005586:=z*(5-6*z+2*z**2)/(z-1)**4; # conjectured by Simon Plouffe in his 1992 dissertation
    seq(sum(binomial(n,m), m=1..3)-n^2,n=5..49); # Zerinvary Lajos, Jun 19 2008
  • Mathematica
    Clear[lst,n,a,f]; f[n_]:=n*(n+1)/2; a=0;lst={};Do[a+=f[n];AppendTo[lst,a-f[n+1]],{n,5!}];lst (* Vladimir Joseph Stephan Orlovsky, Oct 13 2009 *)
    CoefficientList[Series[x (5 - 6 x + 2 x^2) / (1 - x)^4, {x, 0, 50}], x] (* Vincenzo Librandi, Jun 09 2013 *)
    Table[(n(n+4)(n+5))/6,{n,0,50}] (* or *) LinearRecurrence[{4,-6,4,-1},{0,5,14,28},50] (* Harvey P. Dale, Jul 14 2018 *)
  • PARI
    {a(n) = n * (n+4) * (n+5) / 6}; /* Michael Somos, Apr 13 2007 */
    

Formula

G.f.: x * (5 - 6*x + 2*x^2) / (1 - x)^4.
E.g.f.: (5*x + 2*x^2 + x^3/6) * exp(x). - Michael Somos, Apr 13 2007
Let t(n) = n*(n+1)/2, te(n) = (n+1)*(n+2)*(n+3)/6. Then a(n-4) = -2*t(n) + te(n-1), e.g., a(2) = -2*t(6) + te(5) = -2*21 + 56 = 14, where te(n) are the tetrahedral numbers A000292 and t(n) are the triangular numbers A000217. - Jon Perry, Jul 23 2003
a(n) = C(5+n, 3)-C(5+n, 2). - Zerinvary Lajos, Jan 09 2006
a(n) = C(n,3) - C(n,1), n>=4. - Zerinvary Lajos, Nov 25 2006
a(n) = - A005581(-4-n) for all n in Z. - Michael Somos, Apr 13 2007
a(n) = A214292(n+4,2). - Reinhard Zumkeller, Jul 12 2012
From Amiram Eldar, Feb 20 2022: (Start)
Sum_{n>=1} 1/a(n) = 77/200.
Sum_{n>=1} (-1)^(n+1)/a(n) = 363/200 - 12*log(2)/5. (End)
a(n) = A005581(n+2)-2. - R. J. Mathar, Nov 22 2024

Extensions

M3842=A005555 in the 1995 EIS was the same sequence as this.
More terms from Zerinvary Lajos, Jan 09 2006
Previous Showing 41-50 of 124 results. Next