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

A001558 Number of hill-free Dyck paths of semilength n+3 and having length of first descent equal to 1 (a hill in a Dyck path is a peak at level 1).

Original entry on oeis.org

1, 3, 10, 33, 111, 379, 1312, 4596, 16266, 58082, 209010, 757259, 2760123, 10114131, 37239072, 137698584, 511140558, 1904038986, 7115422212, 26668376994, 100221202998, 377570383518, 1425706128480, 5394898197448, 20454676622476
Offset: 0

Views

Author

Keywords

Comments

a(n) is also the number of even-length descents to ground level in all Dyck paths of semilength n+2. Example: a(1)=3 because in UDUDUD, UDUU(DD), UU(DD)UD, UUDU(DD) and UUUDDD we have 3 even-length descents to ground level (shown between parentheses). - Emeric Deutsch, Oct 05 2008
Convolution of A000108 with A104629. - Philippe Deléham, Nov 11 2009
The Kn12 triangle sums of A039599 are given by the terms of this sequence. For the definition of this and other triangle sums see A180662. - Johannes W. Meijer, Apr 20 2011

Examples

			a(1)=3 because we have uu(d)ududd, uuu(d)uddd and uu(d)uuddd, where u=(1,1), d=(1,-1) (the first descents are shown between parentheses).
G.f. = 1 + 3*x + 10*x^2 + 33*x^3 + 111*x^4 + 379*x^5 + 1312*x^6 + ...
		

References

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

Crossrefs

Cf. A111301. - Emeric Deutsch, Oct 05 2008

Programs

  • Magma
    m:=30; R:=PowerSeriesRing(Rationals(), m); Coefficients(R!( (1-Sqrt(1-4*x))/(x*(3-Sqrt(1-4*x)))*((1-Sqrt(1-4*x))/(2*x))^3 )); // G. C. Greubel, Feb 12 2019
    
  • Maple
    F:=(1-sqrt(1-4*z))/z/(3-sqrt(1-4*z)): C:=(1-sqrt(1-4*z))/2/z: g:=F*C^3: gser:=series(g,z=0,32): seq(coeff(gser,z,n),n=0..27); # Emeric Deutsch, May 08 2006
  • Mathematica
    CoefficientList[Series[(1-Sqrt[1-4*x])/(x*(3-Sqrt[1-4*x]))*((1-Sqrt[1-4*x])/(2*x))^3, {x, 0, 30}], x] (* Vaclav Kotesovec, Mar 20 2014 *)
  • PARI
    my(x='x+O('x^30)); Vec((1-sqrt(1-4*x))/(x*(3-sqrt(1-4*x)))*((1-sqrt(1-4*x))/(2*x))^3) \\ G. C. Greubel, Feb 12 2019
    
  • Sage
    ((1-sqrt(1-4*x))/(x*(3-sqrt(1-4*x)))*((1-sqrt(1-4*x))/(2*x))^3).series(x, 30).coefficients(x, sparse=False) # G. C. Greubel, Feb 12 2019

Formula

a(n) = A000957(n+4) - A000957(n+3) - A000957(n+2) (A000957 are the Fine numbers). - Emeric Deutsch, May 08 2006
a(n) = A118972(n+3,1). - Emeric Deutsch, May 08 2006
G.f.: F*C^3, where F = (1-sqrt(1-4z))/(z*(3-sqrt(1-4z))) and C = (1-sqrt(1-4z))/(2z) is the Catalan function. - Emeric Deutsch, May 08 2006
a(n) = Sum_{k>=0} k*A111301(n+2,k). - Emeric Deutsch, Oct 05 2008
(n+3)*a(n) = (-(11/2)*n + 21/2)*a(n-3) + ((9/2)*n + 11/2)*a(n-1) + (-(1/2)*n + 9/2)*a(n-2) + (-2n + 5)*a(n-4). - Simon Plouffe, Feb 09 2012
a(n) ~ 11*2^(2*n+4)/(9*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Mar 20 2014

Extensions

Edited by Emeric Deutsch, May 08 2006

A167772 Riordan array (c(x)/(1+x*c(x)), x*c(x)), c(x) the g.f. of A000108.

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 2, 3, 2, 1, 6, 8, 6, 3, 1, 18, 24, 18, 10, 4, 1, 57, 75, 57, 33, 15, 5, 1, 186, 243, 186, 111, 54, 21, 6, 1, 622, 808, 622, 379, 193, 82, 28, 7, 1, 2120, 2742, 2120, 1312, 690, 311, 118, 36, 8, 1, 7338, 9458, 7338, 4596, 2476, 1164, 474, 163, 45, 9, 1
Offset: 0

Views

Author

Philippe Deléham, Nov 11 2009, corrected Nov 12 2009

Keywords

Examples

			Triangle begins:
     1;
     0,    1;
     1,    1,    1;
     2,    3,    2,    1;
     6,    8,    6,    3,   1;
    18,   24,   18,   10,   4,   1;
    57,   75,   57,   33,  15,   5,   1;
   186,  243,  186,  111,  54,  21,   6,  1;
   622,  808,  622,  379, 193,  82,  28,  7,  1;
  2120, 2742, 2120, 1312, 690, 311, 118, 36,  8,  1;
Production matrix begins:
  0, 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;
  ... - _Philippe Deléham_, Mar 03 2013
		

Crossrefs

Programs

  • Haskell
    import Data.List (genericIndex)
    a167772 n k = genericIndex (a167772_row n) k
    a167772_row n = genericIndex a167772_tabl n
    a167772_tabl = [1] : [0, 1] :
                   map (\xs@(:x:) -> x : xs) (tail a065602_tabl)
    -- Reinhard Zumkeller, May 15 2014
    
  • Mathematica
    A065602[n_, k_]:= A065602[n,k]= Sum[(k-1+2*j)*Binomial[2*(n-j)-k-1, n-1]/(2*(n-j) -k-1), {j, 0, (n-k)/2}];
    T[n_, k_]:= If[k==0, A065602[n+1,3] + Boole[n==0], A065602[n+1, k+1]];
    Table[T[n, k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, May 26 2022 *)
  • SageMath
    def A065602(n,k): return sum( (k+2*j-1)*binomial(2*n-2*j-k-1, n-1)/(2*n-2*j-k-1) for j in (0..(n-k)//2) )
    def A167772(n,k):
        if (k==0): return A065602(n+1,3) + bool(n==0)
        else: return A065602(n+1,k+1)
    flatten([[A167772(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, May 26 2022

Formula

Sum_{k=0..n} T(n, k) = A000958(n+1).
From Philippe Deléham, Nov 12 2009: (Start)
Sum_{k=0..n} T(n,k)*2^k = A014300(n).
Sum_{k=0..n} T(n,k)*2^(n-k) = A064306(n). (End)
For n > 0: T(n,0) = A065602(n+1,3), T(n,k) = A065602(n+1,k+1), k = 1..n. - Reinhard Zumkeller, May 15 2014

A167769 Pendular trinomial triangle (p=0), read by rows of 2n+1 terms (n>=0), defined by the recurrence : if 0

Original entry on oeis.org

1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 2, 3, 2, 1, 0, 0, 1, 3, 6, 8, 6, 3, 1, 0, 0, 1, 4, 10, 18, 24, 18, 10, 4, 1, 0, 0, 1, 5, 15, 33, 57, 75, 57, 33, 15, 5, 1, 0, 0, 1, 6, 21, 54, 111, 186, 243, 186, 111, 54, 21, 6, 1, 0, 0, 1, 7, 28, 82, 193, 379, 622, 808, 622, 379, 193, 82, 28, 7, 1, 0, 0
Offset: 0

Views

Author

Philippe Deléham, Nov 11 2009

Keywords

Comments

See A119369 for p=1 and A122445 for p=2. The diagonals may be generated by iterated convolutions of a base sequence B (A000108(n)) with the sequence C (A000957(n+1)) of central terms.

Examples

			Triangle begins :
  1;
  1, 0,  0;
  1, 1,  1,  0,  0;
  1, 2,  3,  2,  1,  0,  0;
  1, 3,  6,  8,  6,  3,  1,  0,  0;
  1, 4, 10, 18, 24, 18, 10,  4,  1, 0, 0,
  1, 5, 15, 33, 57, 75, 57, 33, 15, 5, 1, 0, 0; ...
		

References

  • Kim, Ki Hang; Rogers, Douglas G.; Roush, Fred W. Similarity relations and semiorders. Proceedings of the Tenth Southeastern Conference on Combinatorics, Graph Theory and Computing (Florida Atlantic Univ., Boca Raton, Fla., 1979), pp. 577--594, Congress. Numer., XXIII-XXIV, Utilitas Math., Winnipeg, Man., 1979. MR0561081 (81i:05013) - From N. J. A. Sloane, Jun 05 2012

Crossrefs

Programs

  • Maple
    T:= proc(n, k) option remember;
          if k=0 and n=0 then 1;
        elif k<0 or k>2*(n-1) then 0;
        elif n=2 and k<3 then 1;
        elif kG. C. Greubel, Mar 17 2021
  • Mathematica
    T[n_, k_]:= T[n, k]= If[k==0 && n==0, 1, If[k<0 || k>2*(n-1), 0, If[n==2 && k<3, 1, If[kG. C. Greubel, Mar 17 2021 *)
  • PARI
    T(n, k)=if(k==0 && n==0, 1, if(k>2*n-2 || k<0, 0, if(n==2 && k<=2, 1, if(kPaul D. Hanna, Nov 12 2009
    
  • Sage
    @CachedFunction
    def T(n, k):
        if (k==0 and n==0): return 1
        elif (k<0 or k>2*(n-1)): return 0
        elif (n==2 and k<3): return 1
        elif (kG. C. Greubel, Mar 17 2021

Formula

Sum_{k=0..2*n} T(n,k) = A071724(n) = [n=0] + 3*binomial(2n,n-1)/(n+2) = [n=0] + n*C(n)/(n+2), where C(n) are the Catalan numbers (A000108). - G. C. Greubel, Mar 17 2021

A237619 Riordan array (1/(1+x*c(x)), x*c(x)) where c(x) is the g.f. of Catalan numbers (A000108).

Original entry on oeis.org

1, -1, 1, 0, 0, 1, -1, 1, 1, 1, -2, 2, 3, 2, 1, -6, 6, 8, 6, 3, 1, -18, 18, 24, 18, 10, 4, 1, -57, 57, 75, 57, 33, 15, 5, 1, -186, 186, 243, 186, 111, 54, 21, 6, 1, -622, 622, 808, 622, 379, 193, 82, 28, 7, 1, -2120, 2120, 2742, 2120, 1312, 690, 311, 118, 36, 8, 1
Offset: 0

Views

Author

Philippe Deléham, Feb 10 2014

Keywords

Examples

			Triangle begins:
    1;
   -1,  1;
    0,  0,  1;
   -1,  1,  1,  1;
   -2,  2,  3,  2,  1;
   -6,  6,  8,  6,  3,  1;
  -18, 18, 24, 18, 10,  4, 1;
  -57, 57, 75, 57, 33, 15, 5, 1;
Production matrix begins:
  -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;
		

Crossrefs

Programs

  • Mathematica
    A065602[n_, k_]:= A065602[n, k]= Sum[(k-1+2*j)*Binomial[2*(n-j)-k-1, n-1]/(2*(n - j) -k-1), {j,0,(n-k)/2}];
    T[n_, k_]:= If[k==0, A065602[n, 0], If[n==1 && k==1, 1, A065602[n, k]]];
    Table[T[n, k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, May 27 2022 *)
  • SageMath
    def A065602(n, k): return sum( (k+2*j-1)*binomial(2*n-2*j-k-1, n-1)/(2*n-2*j-k-1) for j in (0..(n-k)//2) )
    def A237619(n, k):
        if (n<2): return (-1)^(n-k)
        elif (k==0): return A065602(n, 0)
        else: return A065602(n, k)
    flatten([[A237619(n, k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, May 27 2022

Formula

Sum_{k=0..n} T(n,k)*x^k = A126983(n), A000957(n+1), A026641(n) for x = 0, 1, 2 respectively.
T(n, k) = A167772(n-1, k-1) for k > 0, with T(n, 0) = A167772(n, 0).
T(n, 0) = A126983(n).
T(n+1, 1) = A000957(n+1).
T(n+2, 2) = A000958(n+1).
T(n+3, 3) = A104629(n) = A000957(n+3).
T(n+4, 4) = A001558(n).
T(n+5, 5) = A001559(n).
T(n, k) = A065602(n, k) for k > 0, with T(n, k) = (-1)^(n-k), for n < 2, and T(n, 0) = A065602(n, 0). - G. C. Greubel, May 27 2022
Showing 1-4 of 4 results.