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-5 of 5 results.

A123164 Row sums of A123160.

Original entry on oeis.org

1, 2, 8, 38, 192, 1002, 5336, 28814, 157184, 864146, 4780008, 26572086, 148321344, 830764794, 4666890936, 26283115038, 148348809216, 838944980514, 4752575891144, 26964373486406, 153196621856192, 871460014012682, 4962895187697048, 28292329581548718, 161439727075246592
Offset: 0

Views

Author

Roger L. Bagula, Oct 02 2006

Keywords

Comments

Coefficient of x^n in ((1 + x)/(1 - x))^n. - Paul Barry, Jan 20 2008
a(n) is also the number of order-preserving partial transformations (of an n-element chain). Equivalently, it is the order of the semigroup (monoid) of order-preserving partial transformations (of an n-element chain), PO sub n. - Abdullahi Umar, Aug 25 2008
Hankel transform is A180966. - Paul Barry, Sep 29 2010

Crossrefs

Essentially identical to A002003.

Programs

  • Magma
    [1] cat [n le 2 select 2*4^(n-1) else (4*(3*(n-1)^2-1)*Self(n-1) - (2*n-1)*(n-2)*Self(n-2))/((2*n-3)*(n)): n in [1..30]]; // G. C. Greubel, Jul 19 2023
    
  • Mathematica
    a[n_]:= a[n]= Sum[Binomial[n+k-1,k]*Binomial[n,k], {k,0,n}];
    Table[a[n], {n,0,30}]
  • SageMath
    def A123164(n): return sum(binomial(n,j)*binomial(n+j-1,j) for j in range(n+1))
    [A123164(n) for n in range(31)] # G. C. Greubel, Jul 19 2023

Formula

a(n) = A122542(2*n,n). - Philippe Deléham, May 28 2007
a(n) = Sum_{k=0..n} C(n, k)*C(n+k-1, k). - Paul Barry, Aug 22 2007
(2*n-1)*(n+1)*a(n+1) = 4*(3*n^2-1)*a(n) - (2*n+1)*(n-1)*a(n-1) for n >= 1 with a(0) = 1 and a(1) = 2. - Abdullahi Umar, Aug 25 2008
a(n) = Jacobi_P(n, 0, -1, 3). - Paul Barry, Sep 27 2009
G.f.: (1 + x + sqrt(1 - 6*x + x^2))/(2*sqrt(1 - 6*x + x^2)). - Paul Barry, Sep 29 2010
From Abdullahi Umar, Oct 11 2008: (Start)
a(n+1) - a(n) = (2*n + 1)*A006318 (n >= 0);
2*a(n) = (n + 1)*A006318(n) - (n - 1)*A006318(n-1) (n > 0). (End)
a(n) = Hypergeometric2F1([-n, n], [1], -1). - Peter Luschny, Aug 02 2014
a(n) ~ (1 + sqrt(2))^(2*n) / (2^(3/4) * sqrt(Pi*n)). - Vaclav Kotesovec, Feb 14 2021
From Peter Bala, Oct 07 2021: (Start)
a(n) = Sum_{k = 0..floor(n/2)} (-1)^k*C(n, k)*C(3*n-2*k-1, n-2*k).
a(p) == 2 (mod p^3) for prime p >= 5.
Conjecture: a(n*p^k) == a(n*p^(k-1)) mod( p^(3*k) ) for prime p >= 5 and all positive integers n and k. (End)

Extensions

Edited by N. J. A. Sloane, Oct 04 2006
Offset changed (a(0)=1) by Michael Somos, Feb 07 2011

A088617 Triangle read by rows: T(n,k) = C(n+k,n)*C(n,k)/(k+1), for n >= 0, k = 0..n.

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 6, 10, 5, 1, 10, 30, 35, 14, 1, 15, 70, 140, 126, 42, 1, 21, 140, 420, 630, 462, 132, 1, 28, 252, 1050, 2310, 2772, 1716, 429, 1, 36, 420, 2310, 6930, 12012, 12012, 6435, 1430, 1, 45, 660, 4620, 18018, 42042, 60060, 51480, 24310, 4862
Offset: 0

Views

Author

N. J. A. Sloane, Nov 23 2003

Keywords

Comments

Row sums: A006318 (Schroeder numbers). Essentially same as triangle A060693 transposed.
T(n,k) is number of Schroeder paths (i.e., consisting of steps U=(1,1), D=(1,-1), H=(2,0) and never going below the x-axis) from (0,0) to (2n,0), having k U's. E.g., T(2,1)=3 because we have UHD, UDH and HUD. - Emeric Deutsch, Dec 06 2003
Little Schroeder numbers A001003 have a(n) = Sum_{k=0..n} A088617(n,k)*(-1)^(n-k)*2^k. - Paul Barry, May 24 2005
Conjecture: The expected number of U's in a Schroeder n-path is asymptotically Sqrt[1/2]*n for large n. - David Callan, Jul 25 2008
T(n, k) is also the number of order-preserving and order-decreasing partial transformations (of an n-chain) of width k (width(alpha) = |Dom(alpha)|). - Abdullahi Umar, Oct 02 2008
The antidiagonals of this lower triangular matrix are the rows of A055151. - Tom Copeland, Jun 17 2015

Examples

			Triangle begins:
  [0] 1;
  [1] 1,  1;
  [2] 1,  3,   2;
  [3] 1,  6,  10,    5;
  [4] 1, 10,  30,   35,    14;
  [5] 1, 15,  70,  140,   126,    42;
  [6] 1, 21, 140,  420,   630,   462,   132;
  [7] 1, 28, 252, 1050,  2310,  2772,  1716,   429;
  [8] 1, 36, 420, 2310,  6930, 12012, 12012,  6435,  1430;
  [9] 1, 45, 660, 4620, 18018, 42042, 60060, 51480, 24310, 4862;
		

References

  • Charles Jordan, Calculus of Finite Differences, Chelsea 1965, p. 449.

Crossrefs

Programs

  • Magma
    [[Binomial(n+k,n)*Binomial(n,k)/(k+1): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Jun 18 2015
    
  • Maple
    R := n -> simplify(hypergeom([-n, n + 1], [2], -x)):
    Trow := n -> seq(coeff(R(n, x), x, k), k = 0..n):
    seq(print(Trow(n)), n = 0..9); # Peter Luschny, Apr 26 2022
  • Mathematica
    Table[Binomial[n+k, n] Binomial[n, k]/(k+1), {n,0,10}, {k,0,n}]//Flatten (* Michael De Vlieger, Aug 10 2017 *)
  • PARI
    {T(n, k)= if(k+1, binomial(n+k, n)*binomial(n, k)/(k+1))}
    
  • SageMath
    flatten([[binomial(n+k, 2*k)*catalan_number(k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, May 22 2022

Formula

Triangle T(n, k) read by rows; given by [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...] DELTA [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...] where DELTA is Deléham's operator defined in A084938.
T(n, k) = A085478(n, k)*A000108(k); A000108 = Catalan numbers. - Philippe Deléham, Dec 05 2003
Sum_{k=0..n} T(n, k)*x^k*(1-x)^(n-k) = A000108(n), A001003(n), A007564(n), A059231(n), A078009(n), A078018(n), A081178(n), A082147(n), A082181(n), A082148(n), A082173(n) for x = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11. - Philippe Deléham, Aug 18 2005
Sum_{k=0..n} T(n,k)*x^k = (-1)^n*A107841(n), A080243(n), A000007(n), A000012(n), A006318(n), A103210(n), A103211(n), A133305(n), A133306(n), A133307(n), A133308(n), A133309(n) for x = -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8 respectively. - Philippe Deléham, Oct 18 2007
O.g.f. (with initial 1 excluded) is the series reversion with respect to x of (1-t*x)*x/(1+x). Cf. A062991 and A089434. - Peter Bala, Jul 31 2012
G.f.: 1 + (1 - x - T(0))/y, where T(k) = 1 - x*(1+y)/( 1 - x*y/T(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Nov 03 2013
From Peter Bala, Jul 20 2015: (Start)
O.g.f. A(x,t) = ( 1 - x - sqrt((1 - x)^2 - 4*x*t) )/(2*x*t) = 1 + (1 + t)*x + (1 + 3*t + 2*t^2)*x^2 + ....
1 + x*(dA(x,t)/dx)/A(x,t) = 1 + (1 + t)*x + (1 + 4*t + 3*t^2)*x^2 + ... is the o.g.f. for A123160.
For n >= 1, the n-th row polynomial equals (1 + t)/(n+1)*Jacobi_P(n-1,1,1,2*t+1). Removing a factor of 1 + t from the row polynomials gives the row polynomials of A033282. (End)
From Tom Copeland, Jan 22 2016: (Start)
The o.g.f. G(x,t) = {1 - (2t+1) x - sqrt[1 - (2t+1) 2x + x^2]}/2x = (t + t^2) x + (t + 3t^2 + 2t^3) x^2 + (t + 6t^2 + 10t^3 + 5t^3) x^3 + ... generating shifted rows of this entry, excluding the first, was given in my 2008 formulas for A033282 with an o.g.f. f1(x,t) = G(x,t)/(1+t) for A033282. Simple transformations presented there of f1(x,t) are related to A060693 and A001263, the Narayana numbers. See also A086810.
The inverse of G(x,t) is essentially given in A033282 by x1, the inverse of f1(x,t): Ginv(x,t) = x [1/(t+x) - 1/(1+t+x)] = [((1+t) - t) / (t(1+t))] x - [((1+t)^2 - t^2) / (t(1+t))^2] x^2 + [((1+t)^3 - t^3) / (t(1+t))^3] x^3 - ... . The coefficients in t of Ginv(xt,t) are the o.g.f.s of the diagonals of the Pascal triangle A007318 with signed rows and an extra initial column of ones. The numerators give the row o.g.f.s of signed A074909.
Rows of A088617 are shifted columns of A107131, whose reversed rows are the Motzkin polynomials of A055151, related to A011973. The diagonals of A055151 give the rows of A088671, and the antidiagonals (top to bottom) of A088617 give the rows of A107131 and reversed rows of A055151. The diagonals of A107131 give the columns of A055151. The antidiagonals of A088617 (bottom to top) give the rows of A055151.
(End)
T(n, k) = [x^k] hypergeom([-n, 1 + n], [2], -x). - Peter Luschny, Apr 26 2022

A178301 Triangle T(n,k) = binomial(n,k)*binomial(n+k+1,n+1) read by rows, 0 <= k <= n.

Original entry on oeis.org

1, 1, 3, 1, 8, 10, 1, 15, 45, 35, 1, 24, 126, 224, 126, 1, 35, 280, 840, 1050, 462, 1, 48, 540, 2400, 4950, 4752, 1716, 1, 63, 945, 5775, 17325, 27027, 21021, 6435, 1, 80, 1540, 12320, 50050, 112112, 140140, 91520, 24310, 1, 99, 2376, 24024, 126126, 378378, 672672, 700128, 393822, 92378
Offset: 0

Views

Author

Alford Arnold, May 30 2010

Keywords

Comments

Antidiagonal sums are given by A113682. - Johannes W. Meijer, Mar 24 2013
The rows seem to give (up to sign) the coefficients in the expansion of the integer-valued polynomial binomial(x+n,n)*binomial(x+n,n-1) in the basis made of the binomial(x+i,i). - F. Chapoton, Nov 01 2022
Chapoton's observation above is correct: the precise expansion is binomial(x+n,n)*binomial(x+n,n-1) = Sum_{k = 0..n-1} (-1)^k*T(n-1,n-1-k)*binomial(x+2*n-1-k,2*n-1-k), as can be verified using the WZ algorithm. For example, n = 4 gives binomial(x+4,4)*binomial(x+4,3) = 35*binomial(x+7,7) - 45*binomial(x+6,6) + 15*binomial(x+5,5) - binomial(x+4,4). - Peter Bala, Jun 24 2023

Examples

			n=0: 1;
n=1: 1,  3;
n=2: 1,  8,  10;
n=3: 1, 15,  45,   35;
n=4: 1, 24, 126,  224,   126;
n=5: 1, 35, 280,  840,  1050,   462;
n=6: 1, 48, 540, 2400,  4950,  4752,  1716;
n=7: 1, 63, 945, 5775, 17325, 27027, 21021, 6435;
		

Crossrefs

Programs

  • Maple
    A178301 := proc(n,k)
            binomial(n,k)*binomial(n+k+1,n+1) ;
    end proc: # R. J. Mathar, Mar 24 2013
    R := proc(n) add((-1)^(n+k)*(2*k+1)*orthopoly:-P(k,2*x+1)/(n+1), k=0..n) end:
    for n from 0 to 6 do seq(coeff(R(n), x, k), k=0..n) od; # Peter Luschny, Aug 25 2021
  • Mathematica
    Flatten[Table[Binomial[n,k]Binomial[n+k+1,n+1],{n,0,10},{k,0,n}]] (* Harvey P. Dale, Aug 23 2014 *)
  • Maxima
    create_list(binomial(n,k)*binomial(n+k+1,n+1),n,0,12,k,0,n); /* Emanuele Munarini, Dec 16 2016 */
    
  • PARI
    R(n,x) = sum(k=0,n, (-1)^(n+k) * (2*k+1) * pollegendre(k,2*x+1)) / (n+1); \\ Max Alekseyev, Aug 25 2021

Formula

T(n,k) = A007318(n,k) * A178300(n+1,k+1).
From Peter Bala, Jun 18 2015: (Start)
n-th row polynomial R(n,x) = Sum_{k = 0..n} binomial(n,k)*binomial(n+k+1,n+1)*x^k = Sum_{k = 0..n} (-1)^(n+k)*binomial(n+1,k+1)*binomial(n+k+1,n+1)*(1 + x)^k.
Recurrence: (2*n - 1)*(n + 1)*R(n,x) = 2*(4*n^2*x + 2*n^2 - x - 1)*R(n-1,x) - (2*n + 1)(n - 1)*R(n-2,x) with R(0,x) = 1, R(1,x) = 1 + 3*x.
A182626(n) = -R(n-1,-2) for n >= 1. (End)
From Peter Bala, Jul 20 2015: (Start)
n-th row polynomial R(n,x) = Jacobi_P(n,0,1,2*x + 1).
(1 + x)*R(n,x) gives the row polynomials of A123160. (End)
G.f.: (1+x-sqrt(1-2*x+x^2-4*x*y))/(2*(1+y)*x*sqrt(1-2*x+x^2-4*x*y)). - Emanuele Munarini, Dec 16 2016
R(n,x) = Sum_{k=0..n} (-1)^(n+k)*(2*k+1)*P(k,2*x+1)/(n+1), where P(k,x) is the k-th Legendre polynomial (cf. A100258) and P(k,2*x+1) is the k-th shifted Legendre polynomial (cf. A063007). - Max Alekseyev, Jun 28 2018; corrected by Peter Bala, Aug 08 2021
Polynomial g(n,x) = R(n,-x)/(n+1) delivers the maximum of f(1)^2/(Integral_{x=0..1} f(x)^2 dx) over all polynomials f(x) with real coefficients and deg(f(x)) <= n. This maximum equals (n+1)^2. See dxdy.ru link. - Max Alekseyev, Jun 28 2018

A122899 Triangle with row sums counting directed animals.

Original entry on oeis.org

1, 1, 1, 0, 4, 1, 0, 3, 9, 1, 0, 0, 18, 16, 1, 0, 0, 10, 60, 25, 1, 0, 0, 0, 80, 150, 36, 1, 0, 0, 0, 35, 350, 315, 49, 1, 0, 0, 0, 0, 350, 1120, 588, 64, 1, 0, 0, 0, 0, 126, 1890, 2940, 1008
Offset: 0

Views

Author

Paul Barry, Sep 18 2006

Keywords

Comments

Row sums are A005773(n+1). Product of A007318 and A122899 is A103371.

Examples

			Triangle begins
1,
1, 1,
0, 4, 1,
0, 3, 9, 1,
0, 0, 18, 16, 1,
0, 0, 10, 60, 25, 1,
0, 0, 0, 80, 150, 36, 1,
0, 0, 0, 35, 350, 315, 49, 1,
0, 0, 0, 0, 350, 1120, 588, 64, 1,
0, 0, 0, 0, 126, 1890, 2940, 1008, 81, 1,
0, 0, 0, 0, 0, 1512, 7350, 6720, 1620, 100, 1
		

Crossrefs

Cf. A123160.

Formula

Number triangle T(n,k)=sum{j=0..n, (-1)^(n-j)C(n,j)C(j+1,k+1)C(j,k)}
T(n,k) = C(n,k)*C(k+1,n-k). The columns of this triangle (ignoring leading zeros) give the rows of A123160. - Peter Bala, Jan 24 2008

A306226 Triangle read by rows: T(n,k) = Sum_{i=0..n/2} C(n-i,i)*C(n-i,k-i)*C(n-1,i) (0 <= k <= n).

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 7, 11, 5, 1, 13, 36, 37, 13, 1, 21, 92, 160, 123, 35, 1, 31, 200, 520, 655, 401, 96, 1, 43, 387, 1405, 2575, 2541, 1293, 267, 1, 57, 686, 3325, 8295, 11711, 9492, 4131, 750, 1, 73, 1136, 7112, 23128, 43736, 50148, 34476, 13107, 2123, 1, 91, 1782, 14040, 57708, 140112, 212856, 205332, 122535, 41353, 6046
Offset: 0

Views

Author

Vladimir Kruchinin, Feb 16 2019

Keywords

Examples

			1;
1, 1;
1, 3, 2;
1, 7, 11, 5;
1, 13, 36, 37, 13;
1, 21, 92, 160, 123, 35;
		

Crossrefs

Cf. A123160.
Cf. A005773 (right diagonal).

Programs

  • Maxima
    T(n,k):=sum(binomial(n-i,i)*binomial(n-i,k-i)*binomial(n-1,i),i,0,n/2);

Formula

G.f.: (x*y+x+1)/(2*sqrt((-x*y-x+1)^2-4*x*y*(x*y+x)))+1/2.
Showing 1-5 of 5 results.