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

A026378 a(n) = number of integer strings s(0),...,s(n) counted by array T in A026374 that have s(n)=1; also a(n) = T(2n-1,n-1).

Original entry on oeis.org

1, 4, 17, 75, 339, 1558, 7247, 34016, 160795, 764388, 3650571, 17501619, 84179877, 406020930, 1963073865, 9511333155, 46169418195, 224484046660, 1093097083475, 5329784874185, 26018549129545, 127154354598330, 622031993807565
Offset: 1

Views

Author

Keywords

Comments

Number of lattice paths from (0,0) to the line x=n-1 that do not go below the line y=0 and consist of steps U=(1,1), D=(1,-1) and three types of steps H=(1,0) (left factors of 3-Motzkin steps). Example: a(3)=17 because we have UD, UU, 9 HH paths, 3 HU paths and 3 UH paths. - Emeric Deutsch, Jan 22 2004
Also a(n) = number of integer strings s(0), ..., s(n) counted by array U in A026386 that have s(n)=1; a(n) = U(2n-1, n-1).
The Hankel transform of [1,1,4,17,75,339,1558,...] is [1,3,8,21,55,144,377,...] (see A001906). - Philippe Deléham, Apr 13 2007
Number of peaks in all skew Dyck paths of semilength n. A skew Dyck path is a path in the first quadrant which begins at the origin, ends on the x-axis, consists of steps U=(1,1)(up), D=(1,-1)(down) and L=(-1,-1)(left) so that up and left steps do not overlap. The length of the path is defined to be the number of its steps. Example: a(2)=4 because in the 3 (=A002212(2)) skew Dyck paths (UD)(UD), U(UD)D and U(UD)L we have altogether 4 peaks (shown between parentheses). - Emeric Deutsch, Jul 25 2007
Hankel transform of this sequence gives A000012 = [1,1,1,1,1,1,...]. - Philippe Deléham, Oct 24 2007
5th binomial transform of (-1)^n*A000108. - Paul Barry, Jan 13 2009
From Gary W. Adamson, May 17 2009: (Start)
Convolved with A007317, (1, 2, 5, 15, 51, ...) = A026376: (1, 6, 30, 144, ...)
Equals A026375, (1, 3, 11, 45, 195, ...) convolved with A002212 prefaced with
a 1: (1, 1, 3, 10, 36, 137, ...). (End)
From Tom Copeland, Nov 09 2014: (Start)
The array belongs to an interpolated family of arrays associated to the Catalan A000108 (t=1), and Riordan, or Motzkin sums A005043 (t=0), with the interpolating o.g.f. [1-sqrt(1-4x/(1+(1-t)x))]/2 and inverse x(1-x)/[1+(t-1)x(1-x)]. See A091867 for more info on this family. Here the interpolation is t=-4 (mod signs in the results).
Let C(x) = [1 - sqrt(1-4x)]/2, an o.g.f. for the Catalan numbers A000108, with inverse Cinv(x) = x*(1-x) and P(x,t) = x/(1+t*x) with inverse P(x,-t).
O.g.f: G(x) = [-1 + sqrt(1 + 4*x/(1-5x))]/2 = -C[P(-x,5)].
Inverse O.g.f: Ginv(x) = x*(1+x)/[1 + 5x*(1+x)] = -P(Cinv(-x),-5) (signed A039717). (End)

Crossrefs

Half the values of A026387. Bisection of A026380 and A026392.

Programs

  • Maple
    a := n -> (-1)^n*simplify(GegenbauerC(n-2,-n+1,3/2) - GegenbauerC(n-1,-n+1,3/2)): seq(a(n), n=1..23); # Peter Luschny, May 13 2016
  • Mathematica
    CoefficientList[Series[(1/2)/(5*x^2-x)*(1-5*x-(1-6*x+5*x^2)^(1/2)),{x,0,30}],x] (* Vincenzo Librandi, May 13 2012 *)
    Table[Hypergeometric2F1[3/2, 1-n, 2, -4], {n, 1, 20}] (* Vladimir Reshetnikov, Apr 25 2016 *)

Formula

G.f.: (1/2)/(5*x^2-x)*(1-5*x-(1-6*x+5*x^2)^(1/2)). E.g.f.: exp(3*x)*(BesselI(0, 2*x)+BesselI(1, 2*x)). - Vladeta Jovovic, Oct 03 2003
G.f.: [(1-z)/sqrt(1-6z+5z^2)-1]/2 = z + 4z^2 + 17z^3 + ... - Emeric Deutsch, Jan 22 2004
a(n) = coefficient of t^n in (1+t)(1+3t+t^2)^(n-1). - Emeric Deutsch, Jan 30 2004
a(n) = A026380(2n-2). - Emeric Deutsch, Feb 18 2004
a(n) = [2(3n-2)a(n-1) - 5(n-2)a(n-2)]/n for n>=2; a(0)=0, a(1)=1. - Emeric Deutsch, Mar 18 2004
a(n+1) = sum(k=0, n, binomial(n, k)*sum(i=0, k, binomial(k+i, i))). - Benoit Cloitre, Aug 06 2004
a(n+1) = sum(k=0, n, binomial(n, k)*binomial(2*k+1, k+1)). - Benoit Cloitre, Aug 06 2004
a(n) = Sum(k*A126182(n-1,k-1),k=1..n). - Emeric Deutsch, Jul 25 2007
From Paul Barry, Jan 13 2009: (Start)
G.f.: (1/(1-5x))*c(-x/(1-5x)), c(x) the g.f. of A000108;
a(n) = sum{k=0..n, C(n,k)*(-1)^k*A000108(k)*5^(n-k)} (offset 0). (End)
G.f. 1/(1 - 3x - x(1 - x)/(1 - x - x(1 - x)/(1 - x - x(1 - x)/(1 - x - x(1 - x)/(1...(continued fraction). - Aoife Hennessy (aoife.hennessy(AT)gmail.com), Jul 02 2010
a(n) ~ 5^(n-1/2)/sqrt(Pi*n). - Vaclav Kotesovec, Oct 08 2012
a(n) = hypergeom([3/2, 1-n], [2], -4). - Vladimir Reshetnikov, Apr 25 2016
a(n) = (-1)^n*(GegenbauerC(n-2,-n+1,3/2) - GegenbauerC(n-1,-n+1,3/2)). - Peter Luschny, May 13 2016

A280000 Number of free pure symmetric multifunctions in one symbol with n positions.

Original entry on oeis.org

1, 0, 1, 1, 3, 5, 12, 25, 57, 128, 296, 688, 1618, 3839, 9170, 22065, 53370, 129807, 317080, 777887, 1915247, 4731932, 11726476, 29143123, 72614115, 181363151, 453975928, 1138697689, 2861607677, 7204169689
Offset: 1

Views

Author

Gus Wiseman, Dec 24 2016

Keywords

Comments

A free pure symmetric multifunction (PSM) in one symbol x is either (case 1) the symbol x, or (case 2) an expression of the form h[g_1,...,g_k] where h is a PSM in x, each of the g_i for i=1..(k>0) is a PSM in x, and for i < j we have g_i <= g_j under a canonical total ordering such as the Mathematica ordering. The number of positions in a PSM is the number of brackets [...] plus the number of x's.

Examples

			Sequence of free pure symmetric multifunctions (second column) together with their numbers of positions (first column) and j-numbers (third column, see A279944 for details) begins:
1 x            1
3 x[x]         2
4 x[x,x]       8
5 x[x][x]      3
5 x[x[x]]      4
5 x[x,x,x]     128
6 x[x,x][x]    12
6 x[x][x,x]    27
6 x[x,x[x]]    32
6 x[x,x,x,x]   32768
6 x[x[x,x]]    262144
7 x[x][x][x]   5
7 x[x[x]][x]   6
7 x[x][x[x]]   9
7 x[x[x][x]]   16
7 x[x[x[x]]]   64
7 x[x,x,x][x]  145
7 x[x,x][x,x]  1728
7 x[x,x,x[x]]  2048
7 x[x][x,x,x]  2187
7 x[x,x,x,x,x] 2147483648
7 x[x,x[x,x]]  137438953472
7 x[x[x,x,x]]  1378913...3030144
		

Crossrefs

Cf. A005043 (non-symmetric case), A279944.

Programs

  • Mathematica
    multing[t_,n_]:=Array[(t+#-1)/#&,n,1,Times];
    a[n_]:=If[n===1,1,Sum[a[k]*Sum[Product[multing[a[First[s]],Length[s]],{s,Split[p]}],{p,IntegerPartitions[n-k-1]}],{k,1,n-2}]];
    Array[a,15]
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
    seq(n)={my(v=[1]); for(n=2, n, my(t=EulerT(v)); v=concat(v, sum(k=1, n-2, v[k]*t[n-k-1]))); v} \\ Andrew Howroyd, Aug 19 2018

A086810 Triangle obtained by adding a leading diagonal 1,0,0,0,... to A033282.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 1, 5, 5, 0, 1, 9, 21, 14, 0, 1, 14, 56, 84, 42, 0, 1, 20, 120, 300, 330, 132, 0, 1, 27, 225, 825, 1485, 1287, 429, 0, 1, 35, 385, 1925, 5005, 7007, 5005, 1430, 0, 1, 44, 616, 4004, 14014, 28028, 32032, 19448, 4862, 0, 1, 54, 936, 7644, 34398, 91728
Offset: 0

Views

Author

Philippe Deléham, Aug 05 2003

Keywords

Comments

Mirror image of triangle A133336. - Philippe Deléham, Dec 10 2008
From Tom Copeland, Oct 09 2011: (Start)
With polynomials
P(0,t) = 0
P(1,t) = 1
P(2,t) = t
P(3,t) = t + 2 t^2
P(4,t) = t + 5 t^2 + 5 t^3
P(5,t) = t + 9 t^2 + 21 t^3 + 14 t^4
The o.g.f. A(x,t) = {1+x-sqrt[(1-x)^2-4xt]}/[2(1+t)] (see Drake et al.).
B(x,t)= x-t x^2/(1-x)= x-t(x^2+x^3+x^4+...) is the comp. inverse in x.
Let h(x,t) = 1/(dB/dx) = (1-x)^2/(1+(1+t)*x*(x-2)) = 1/(1-t(2x+3x^2+4x^3+...)), an o.g.f. in x for row polynomials in t of A181289. Then P(n,t) is given by (1/n!)(h(x,t)*d/dx)^n x, evaluated at x=0, A = exp(x*h(y,t)*d/dy) y, eval. at y=0, and dA/dx = h(A(x,t),t). These results are a special case of A133437 with u(x,t) = B(x,t), i.e., u_1=1 and (u_n)=-t for n > 1. See A001003 for t = 1. (End)
Let U(x,t) = [A(x,t)-x]/t, then U(x,0) = -dB(x,t)/dt and U satisfies dU/dt = UdU/dx, the inviscid Burgers' equation (Wikipedia), also called the Hopf equation (see Buchstaber et al.). Also U(x,t) = U(A(x,t),0) = U(x+tU,0) since U(x,0) = [x-B(x,t)]/t. - Tom Copeland, Mar 12 2012
Diagonals of A132081 are essentially rows of this sequence. - Tom Copeland, May 08 2012
T(r, s) is the number of [0,r]-covering hierarchies with s segments (see Kreweras). - Michel Marcus, Nov 22 2014
From Yu Hin Au, Dec 07 2019: (Start)
T(n,k) is the number of small Schröder n-paths (lattice paths from (0,0) to (2n,0) using steps U=(1,1), F=(2,0), D=(1,-1) with no F step on the x-axis) that has exactly k U steps.
T(n,k) is the number of Schröder trees (plane rooted tree where each internal node has at least two children) with exactly n+1 leaves and k internal nodes. (End)

Examples

			Triangle starts:
  1;
  0,  1;
  0,  1,  2;
  0,  1,  5,  5;
  0,  1,  9, 21, 14;
  ...
		

Crossrefs

The diagonals (except for A000007) are also the diagonals of A033282.
Row sums: A001003 (Schroeder numbers).

Programs

  • Mathematica
    Table[Boole[n == 2] + If[# == -1, 0, Binomial[n - 3, #] Binomial[n + # - 1, #]/(# + 1)] &[k - 1], {n, 2, 12}, {k, 0, n - 2}] // Flatten (* after Jean-François Alcover at A033282, or *)
    Table[If[n == 0, 1, Binomial[n, k] Binomial[n + k, k - 1]/n], {n, 0, 10}, {k, 0, n}] // Flatten (* Michael De Vlieger, Aug 22 2016 *)
  • PARI
    t(n, k) = if (n==0, 1, binomial(n, k)*binomial(n+k, k-1)/n);
    tabl(nn) = {for (n=0, nn, for (k=0, n, print1(t(n,k), ", ");); print(););} \\ Michel Marcus, Nov 22 2014

Formula

Triangle T(n, k) read by rows; given by [0, 1, 0, 1, 0, 1, ...] DELTA [1, 1, 1, 1, 1, 1, 1, 1, 1, ...] where DELTA is Deléham's operator defined in A084938.
For k>0, T(n, k) = binomial(n-1, k-1)*binomial(n+k, k)/(n+1); T(0, 0) = 1 and T(n, 0) = 0 if n > 0. [corrected by Marko Riedel, May 04 2023]
Sum_{k>=0} T(n, k)*2^k = A107841(n). - Philippe Deléham, May 26 2005
Sum_{k>=0} T(n-k, k) = A005043(n). - Philippe Deléham, May 30 2005
T(n, k) = A108263(n+k, k). - Philippe Deléham, May 30 2005
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A001003(n), A107841(n), A131763(n), A131765(n), A131846(n), A131926(n), A131869(n), A131927(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, respectively. - Philippe Deléham, Nov 05 2007
Sum_{k=0..n} T(n,k)*5^k*(-2)^(n-k) = A152601(n). - Philippe Deléham, Dec 10 2008
Sum_{k=0..n} T(n,k)*(-1)^k*3^(n-k) = A154825(n). - Philippe Deléham, Jan 17 2009
Umbrally, P(n,t) = Lah[n-1,-t*a.]/n! = (1/n)*Sum_{k=1..n-1} binomial(n-2,k-1)a_k t^k/k!, where (a.)^k = a_k = (n-1+k)!/(n-1)!, the rising factorial, and Lah(n,t) = n!*Laguerre(n,-1,t) are the Lah polynomials A008297 related to the Laguerre polynomials of order -1. - Tom Copeland, Oct 04 2014
T(n, k) = binomial(n, k)*binomial(n+k, k-1)/n, for k >= 0; T(0, 0) = 1 (see Kreweras, p. 21). - Michel Marcus, Nov 22 2014
P(n,t) = Lah[n-1,-:Dt:]/n! t^(n-1) with (:Dt:)^k = (d/dt)^k t^k = k! Laguerre(k,0,-:tD:) with (:tD:)^j = t^j D^j. The normalized Laguerre polynomials of 0 order are given in A021009. - Tom Copeland, Aug 22 2016

Extensions

Typo in a(60) corrected by Michael De Vlieger, Nov 21 2019

A091867 Triangle read by rows: T(n,k) = number of Dyck paths of semilength n having k peaks at odd height.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 1, 3, 0, 1, 3, 4, 6, 0, 1, 6, 15, 10, 10, 0, 1, 15, 36, 45, 20, 15, 0, 1, 36, 105, 126, 105, 35, 21, 0, 1, 91, 288, 420, 336, 210, 56, 28, 0, 1, 232, 819, 1296, 1260, 756, 378, 84, 36, 0, 1, 603, 2320, 4095, 4320, 3150, 1512, 630, 120, 45, 0, 1, 1585, 6633, 12760, 15015, 11880, 6930, 2772, 990, 165, 55, 0, 1
Offset: 0

Views

Author

Emeric Deutsch, Mar 10 2004

Keywords

Comments

Number of ordered trees with n edges having k leaves at odd height. Row sums are the Catalan numbers (A000108). T(n,0)=A005043(n). Sum_{k=0..n} k*T(n,k) = binomial(2n-2,n-1).
T(n,k)=number of Dyck paths of semilength n and having k ascents of length 1 (an ascent is a maximal string of consecutive up steps). Example: T(4,2)=6 because we have UdUduud, UduuddUd, uuddUdUd, uudUdUdd, UduudUdd and uudUddUd (the ascents of length 1 are indicated by U instead of u).
T(n,k) is the number of Łukasiewicz paths of length n having k level steps (i.e., (1,0)). A Łukasiewicz path of length n is a path in the first quadrant from (0,0) to (n,0) using rise steps (1,k) for any positive integer k, level steps (1,0) and fall steps (1,-1) (see R. P. Stanley, Enumerative Combinatorics, Vol. 2, Cambridge Univ. Press, Cambridge, 1999, p. 223, Exercise 6.19w; the integers are the slopes of the steps). Example: T(4,1)=4 because we have HU(2)DD, U(2)HDD, U(2)DHD and U(2)DDH, where H=(1,0), U(1,1), U(2)=(1,2) and D=(1,-1). - Emeric Deutsch, Jan 06 2005
T(n,k) = number of noncrossing partitions of [n] containing k singleton blocks. Also, T(n,k) = number of noncrossing partitions of [n] containing k adjacencies. An adjacency is an occurrence of 2 consecutive integers in the same block (here 1 and n are considered consecutive). In fact, the statistics # singletons and # adjacencies have a symmetric joint distribution.
Exponential Riordan array [e^x*(Bessel_I(0,2x)-Bessel_I(1,2x)),x]. - Paul Barry, Mar 03 2011
T(n,k) is the number of ordered trees having n edges and exactly k nodes with one child. - Geoffrey Critzer, Feb 25 2013
From Tom Copeland, Nov 04 2014: (Start)
Summing the coeff. of the partitions in A134264 for a Lagrange inversion formula (see also A249548) containing (h_1)^k = (1')^k gives this triangle, so this array's o.g.f. H(x,t) = x + t * x^2 + (1 + t^2) * x^3 ... is the inverse of the o.g.f. of A104597 with a sign change, i.e., H^(-1)(x,t) = (x-x^2) / [1 + (t-1)(x-x^2)] = Cinv(x)/[1 + (t-1)Cinv(x)] = P[Cinv(x),t-1] where Cinv(x)= x * (1-x) is the inverse of C(x) = [1-sqrt(1-4*x)]/2, an o.g.f. for the Catalan numbers A000108, and P(x,t) = x/(1+t*x) with inverse Pinv(x,t) = -P(-x,t) = x/(1-t*x). Therefore,
O.g.f.: H(x,t) = C[Pinv(x,t-1)] = C[P(x,1-t)] = C[x/(1-(t-1)x)] = {1-sqrt[1-4*x/(1-(t-1)x)]}/2 (for A091867). Reprising,
Inverse O.g.f.: H^(-1)(x,t) = x*(1-x) / [1 + (t-1)x(1-x)] = P[Cinv(x),t-1].
From general arguments in A134264, the row polynomials are an Appell sequence with lowering operator d/dt, having the umbral property (p(.,t)+a)^n=p(n,t+a) with e.g.f. = e^(x*t)/w(x), where 1/w(x)= e.g.f. of first column for the Motzkin numbers in A005043. (Mislabeled argument corrected on Jan 31 2016.)
Cf. A124644 (t-shifted polynomials), A026378 (t=-4), A001700 (t=-3), A005773 (t=-2), A126930 (t=-1) and A210736 (t=-1, a(0)=0, unsigned), A005043 (t=0), A000108 (t=1), A007317 (t=2), A064613 (t=3), A104455 (t=4), A030528 (for inverses).
(End)
The sequence of binomial transforms A126930, A005043, A000108, ... in the above comment appears in A126930 and the link therein to a paper by F. Fite et al. on page 42. - Tom Copeland, Jul 23 2016

Examples

			T(4,2)=6 because we have (ud)uu(ud)dd, uu(ud)dd(ud), uu(ud)(ud)dd, (ud)(ud)uudd, (ud)uudd(ud) and uudd(ud)(ud) (here u=(1,1), d=(1,-1) and the peaks at odd height are shown between parentheses).
Triangle begins:
   1;
   0,   1;
   1,   0,   1;
   1,   3,   0,   1;
   3,   4,   6,   0,  1;
   6,  15,  10,  10,  0,  1;
  15,  36,  45,  20, 15,  0, 1;
  36, 105, 126, 105, 35, 21, 0, 1;
  ...
		

References

  • R. Sedgewick and P. Flajolet, Analysis of Algorithms, Addison and Wesley, 1996, page 254 (first edition)

Crossrefs

Programs

  • Maple
    T := proc(n,k) if k>n then 0 elif k=n then 1 else (binomial(n+1,k)/(n+1))*sum(binomial(n+1-k,j)*binomial(n-k-j-1,j-1),j=1..floor((n-k)/2)) fi end: seq(seq(T(n,k),k=0..n),n=0..12);
    T := (n,k) -> (-1)^(n+k)*binomial(n,k)*hypergeom([-n+k,1/2],[2],4): seq(seq(simplify(T(n, k)), k=0..n), n=0..10); # Peter Luschny, Jul 27 2016
    # alternative Maple program:
    b:= proc(x, y, t) option remember; expand(`if`(x=0, 1,
          `if`(y>0, b(x-1, y-1, 0)*z^irem(t*y, 2), 0)+
          `if`(y (p-> seq(coeff(p, z, i), i=0..n))(b(2*n, 0$2)):
    seq(T(n), n=0..16);  # Alois P. Heinz, May 12 2017
  • Mathematica
    nn=10;cy = ( 1 + x - x y - ( -4x(1+x-x y) + (-1 -x + x y)^2)^(1/2))/(2(1+x-x y)); Drop[CoefficientList[Series[cy,{x,0,nn}],{x,y}],1]//Grid  (* Geoffrey Critzer, Feb 25 2013 *)
    Table[Which[k == n, 1, k > n, 0, True, (Binomial[n + 1, k]/(n + 1)) Sum[Binomial[n + 1 - k, j] Binomial[n - k - j - 1, j - 1], {j, Floor[(n - k)/2]}]], {n, 0, 11}, {k, 0, n}] // Flatten (* Michael De Vlieger, Jul 25 2016 *)

Formula

T(n, k) = [binomial(n+1, k)/(n+1)]*Sum_{j=1..floor((n-k)/2)} binomial(n+1-k, j)*binomial(n-k-j-1, j-1) for kn. G.f.=G=G(t, z) satisfies z(1+z-tz)G^2-(1+z-tz)G+1=0. T(n, k)=r(n-k)*binomial(n, k), where r(n)=A005043(n) are the Riordan numbers.
G.f.: 1/(1-xy-x^2/(1-x-xy-x^2/(1-x-xy-x^2/(1-x-xy-x^2/(1-... (continued fraction). - Paul Barry, Aug 03 2009
Sum_{k=0..n} T(n,k)*x^k = A126930(n), A005043(n), A000108(n), A007317(n), A064613(n), A104455(n) for x = -1,0,1,2,3,4 respectively. - Philippe Deléham, Dec 03 2009
Sum_{k=0..n} (-1)^(n-k)*T(n,k)*x^k = A168491(n), A099323(n+1), A001405(n), A005773(n+1), A001700(n), A026378(n+1), A005573(n), A122898(n) for x = -1, 0, 1, 2, 3, 4, 5, 6 respectively. - Philippe Deléham, Dec 03 2009
E.g.f.: e^(x+xy)*(Bessel_I(0,2x)-Bessel_I(1,2x)). - Paul Barry, Mar 10 2010
From Tom Copeland, Nov 06 2014: (Start)
O.g.f.: H(x,t) = {1-sqrt[1-4x/(1-(t-1)x)]}/2 (shifted index, as given in Copeland's comment, see comp. inverse there).
H(x,t)= x / [1-(C.+(t-1))x] = Sum_{n>=1} (C.+ (t-1))^(n-1)*x^n umbrally, e.g., (a.+b.)^2 = a_0*b_2 + 2 a_1*b1_+ a_0*b_2, where (C.)^n = C_n are the Catalan numbers (1,1,2,5,14,..) of A000108.
This shows directly that the lowering operator for the polynomials is D=d/dt, i.e., D p(n,t)= D(C. + (t-1))^n = n * (C. + (t-1))^(n-1) = n*p(n-1,t), so that the polynomials form an Appell sequence, and that p(n,0) gives a Motzkin sum, or Riordan, number A005043.
(End)
T(n,k) = (-1)^(n+k)*binomial(n,k)*hypergeom([k-n,1/2],[2],4). - Peter Luschny, Jul 27 2016

A358371 Number of leaves in the n-th standard ordered rooted tree.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 2, 3, 2, 2, 2, 3, 2, 3, 3, 4, 1, 3, 2, 3, 2, 3, 3, 4, 3, 3, 3, 4, 3, 4, 4, 5, 2, 2, 3, 4, 2, 3, 3, 4, 3, 3, 3, 4, 3, 4, 4, 5, 2, 4, 3, 4, 3, 4, 4, 5, 4, 4, 4, 5, 4, 5, 5, 6, 2, 3, 2, 3, 3, 4, 4, 5, 3, 3, 3, 4, 3, 4, 4, 5, 2, 4, 3, 4, 3, 4
Offset: 1

Views

Author

Gus Wiseman, Nov 13 2022

Keywords

Comments

We define the n-th standard ordered rooted tree to be obtained by taking the (n-1)-th composition in standard order (graded reverse-lexicographic, A066099) as root and replacing each part with its own standard ordered rooted tree. This ranking is an ordered variation of Matula-Goebel numbers, giving a bijective correspondence between positive integers and unlabeled ordered rooted trees.

Examples

			The standard ordered rooted tree ranking begins:
  1: o        10: (((o))o)   19: (((o))(o))
  2: (o)      11: ((o)(o))   20: (((o))oo)
  3: ((o))    12: ((o)oo)    21: ((o)((o)))
  4: (oo)     13: (o((o)))   22: ((o)(o)o)
  5: (((o)))  14: (o(o)o)    23: ((o)o(o))
  6: ((o)o)   15: (oo(o))    24: ((o)ooo)
  7: (o(o))   16: (oooo)     25: (o(oo))
  8: (ooo)    17: ((((o))))  26: (o((o))o)
  9: ((oo))   18: ((oo)o)    27: (o(o)(o))
For example, the 25th ordered tree is (o,(o,o)) because the 24th composition is (1,4) and the 3rd composition is (1,1). Hence a(25) = 3.
		

Crossrefs

The triangle counting trees by this statistic is A001263, unordered A055277.
The version for unordered trees is A109129, nodes A061775, edges A196050.
The nodes are counted by A358372.
A000081 counts unordered rooted trees, ranked by A358378.
A000108 counts ordered rooted trees.
A358374 ranks ordered identity trees, counted by A032027.
A358375 ranks ordered binary trees, counted by A126120

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join @@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    srt[n_]:=If[n==1,{},srt/@stc[n-1]];
    Table[Count[srt[n],{},{0,Infinity}],{n,100}]

A358372 Number of nodes in the n-th standard ordered rooted tree.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 14 2022

Keywords

Comments

We define the n-th standard ordered rooted tree to be obtained by taking the (n-1)-th composition in standard order (graded reverse-lexicographic, A066099) as root and replacing each part with its own standard ordered rooted tree. This ranking is an ordered variation of Matula-Goebel numbers, giving a bijective correspondence between positive integers and unlabeled ordered rooted trees.

Examples

			The standard ordered rooted tree ranking begins:
  1: o        10: (((o))o)   19: (((o))(o))
  2: (o)      11: ((o)(o))   20: (((o))oo)
  3: ((o))    12: ((o)oo)    21: ((o)((o)))
  4: (oo)     13: (o((o)))   22: ((o)(o)o)
  5: (((o)))  14: (o(o)o)    23: ((o)o(o))
  6: ((o)o)   15: (oo(o))    24: ((o)ooo)
  7: (o(o))   16: (oooo)     25: (o(oo))
  8: (ooo)    17: ((((o))))  26: (o((o))o)
  9: ((oo))   18: ((oo)o)    27: (o(o)(o))
For example, the 25th ordered tree is (o,(o,o)) because the 24th composition is (1,4) and the 3rd composition is (1,1). Hence a(25) = 5.
		

Crossrefs

The triangle counting trees by leaves is A001263, unordered A055277.
The version for unordered trees is A061775, leaves A109129, edges A196050.
The leaves are counted by A358371.
A000081 counts unlabeled rooted trees, ranked by A358378.
A358374 ranks ordered identity trees, counted by A032027.
A358375 ranks ordered binary trees, counted by A126120

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join @@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    srt[n_]:=If[n==1,{},srt/@stc[n-1]];
    Table[Count[srt[n],_,{0,Infinity}],{n,100}]

A057682 a(n) = Sum_{j=0..floor(n/3)} (-1)^j*binomial(n,3*j+1).

Original entry on oeis.org

0, 1, 2, 3, 3, 0, -9, -27, -54, -81, -81, 0, 243, 729, 1458, 2187, 2187, 0, -6561, -19683, -39366, -59049, -59049, 0, 177147, 531441, 1062882, 1594323, 1594323, 0, -4782969, -14348907, -28697814, -43046721, -43046721, 0, 129140163, 387420489, 774840978
Offset: 0

Views

Author

N. J. A. Sloane, Oct 20 2000

Keywords

Comments

Let M be any endomorphism on any vector space, such that M^3 = 1 (identity). Then (1-M)^n = A057681(n)-a(n)*M+z(n)*M^2, where z(0)=z(1)=0 and, apparently, z(n+2)=A057083(n). - Stanislav Sykora, Jun 10 2012
From Tom Copeland, Nov 09 2014: (Start)
This array belongs to an interpolated family of arrays associated to the Catalan A000108 (t=1), and Riordan, or Motzkin sums A005043 (t=0), with the interp. (here t=-2) o.g.f. G(x,t) = x(1-x)/[1+(t-1)x(1-x)] and inverse o.g.f. Ginv(x,t) = [1-sqrt(1-4x/(1+(1-t)x))]/2 (Cf. A005773 and A091867 and A030528 for more info on this family). (End)
{A057681, A057682, A*}, where A* is A057083 prefixed by two 0's, is the difference analog of the trigonometric functions {k_1(x), k_2(x), k_3(x)} of order 3. For the definitions of {k_i(x)} and the difference analog {K_i (n)} see [Erdelyi] and the Shevelev link respectively. - Vladimir Shevelev, Jul 31 2017

Examples

			G.f. = x + 2*x^2 + 3*x^3 + 3*x^4 - 9*x^6 - 27*x^7 - 54*x^8 - 81*x^9 + ...
If M^3=1 then (1-M)^6 = A057681(6) - a(6)*M + A057083(4)*M^2 = -18 + 9*M + 9*M^2. - _Stanislav Sykora_, Jun 10 2012
		

References

  • A. Erdelyi, Higher Transcendental Functions, McGraw-Hill, 1955, Vol. 3, Chapter XVIII.

Crossrefs

Alternating row sums of triangle A030523.

Programs

  • Magma
    I:=[0,1,2]; [n le 3 select I[n] else 3*Self(n-1)-3*Self(n-2): n in [1..45]]; // Vincenzo Librandi, Nov 10 2014
    
  • Maple
    A057682:=n->add((-1)^j*binomial(n,3*j+1), j=0..floor(n/3)):
    seq(A057682(n), n=0..50); # Wesley Ivan Hurt, Nov 11 2014
  • Mathematica
    A[n_] := Array[KroneckerDelta[#1, #2 + 1] - KroneckerDelta[#1, #2] + Sum[KroneckerDelta[#1, #2 -q], {q, n}] &, {n, n}];
    Join[{0,1}, Table[(-1)^(n-1)*Total[CoefficientList[ CharacteristicPolynomial[A[(n-1)], x], x]], {n,2,30}]] (* John M. Campbell, Mar 16 2012 *)
    Join[{0}, LinearRecurrence[{3,-3}, {1,2}, 40]] (* Jean-François Alcover, Jan 08 2019 *)
  • PARI
    {a(n) = sum( j=0, n\3, (-1)^j * binomial(n, 3*j + 1))} /* Michael Somos, May 26 2004 */
    
  • PARI
    {a(n) = if( n<2, n>0, n-=2; polsym(x^2 - 3*x + 3, n)[n + 1])} /* Michael Somos, May 26 2004 */
    
  • SageMath
    b=BinaryRecurrenceSequence(3,-3,1,2)
    def A057682(n): return 0 if n==0 else b(n-1)
    [A057682(n) for n in range(41)] # G. C. Greubel, Jul 14 2023

Formula

G.f.: (x - x^2) / (1 - 3*x + 3*x^2).
a(n) = 3*a(n-1) - 3*a(n-2), if n>1.
Starting at 1, the binomial transform of A000484. - Paul Barry, Jul 21 2003
It appears that abs(a(n)) = floor(abs(A000748(n))/3). - John W. Layman, Sep 05 2003
a(n) = ((3+i*sqrt(3))/2)^(n-2) + ((3-i*sqrt(3))/2)^(n-2). - Benoit Cloitre, Oct 27 2003
a(n) = n*3F2(1/3-n/3,2/3-n/3,1-n/3 ; 2/3,4/3 ; 1) for n>=1. - John M. Campbell, Jun 01 2011
Let A(n) be the n X n matrix with -1's along the main diagonal, 1's everywhere above the main diagonal, and 1's along the subdiagonal. Then a(n) equals (-1)^(n-1) times the sum of the coefficients of the characteristic polynomial of A(n-1), for all n>1 (see Mathematica code below). - John M. Campbell, Mar 16 2012
Start with x(0)=1, y(0)=0, z(0)=0 and set x(n+1) = x(n) - z(n), y(n+1) = y(n) - x(n), z(n+1) = z(n) - y(n). Then a(n) = -y(n). But this recurrence falls into a repetitive cycle of length 6 and multiplicative factor -27, so that a(n) = -27*a(n-6) for any n>6. - Stanislav Sykora, Jun 10 2012
a(n) = A057083(n-1) - A057083(n-2). - R. J. Mathar, Oct 25 2012
G.f.: 3*x - 1/3 + 3*x/(G(0) - 1) where G(k)= 1 + 3*(2*k+3)*x/(2*k+1 - 3*x*(k+2)*(2*k+1)/(3*x*(k+2) + (k+1)/G(k+1)));(continued fraction, 3rd kind, 3-step). - Sergei N. Gladkovskii, Nov 23 2012
G.f.: Q(0,u) -1, where u=x/(1-x), Q(k,u) = 1 - u^2 + (k+2)*u - u*(k+1 - u)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Oct 07 2013
From Vladimir Shevelev, Jul 31 2017: (Start)
For n>=1, a(n) = 2*3^((n-2)/2)*cos(Pi*(n-2)/6);
For n>=2, a(n) = K_1(n) + K_3(n-2);
For m,n>=2, a(n+m) = a(n)*K_1(m) + K_1(n)*a(m) - K_3(n-2)*K_3(m-2), where
K_1 = A057681, K_3 = A057083. (End)

A358373 Triangle read by rows where row n lists the sorted standard ordered rooted tree-numbers of all unlabeled ordered rooted trees with n vertices.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 25, 33, 65, 129, 257, 19, 20, 21, 22, 23, 24, 26, 27, 28, 29, 30, 31, 32, 34, 35, 36, 41, 49, 50, 57, 66, 97, 130, 193, 258, 385, 513, 514, 769, 1025, 2049, 4097, 8193, 16385, 32769, 65537, 131073
Offset: 1

Views

Author

Gus Wiseman, Nov 14 2022

Keywords

Comments

We define the standard ordered rooted tree (SORT)-number of an unlabeled ordered rooted tree to be one plus the standard composition number (A066099) of the SORT-numbers of the branches, or 1 if there are no branches. This ranking is an ordered variation of Matula-Goebel numbers, giving a bijective correspondence between positive integers and unlabeled ordered rooted trees.

Examples

			Triangle begins:
   1
   2
   3   4
   5   6   7   8   9
  10  11  12  13  14  15  16  17  18  25  33  65 129 257
For example, the tree t = ((o,o),o) has branches (o,o) and o with SORT-numbers 4 and 1, and the standard composition number of (4,1) is 17, so t has SORT-number 18 and is found in row 5.
		

Crossrefs

The version for compositions is A000027.
Row-lengths are A000108.
The unordered version (using Matula-Goebel numbers) is A061773.
The version for Heinz numbers of partitions is A215366.
The row containing n is A358372(n).
A000081 counts unlabeled rooted trees, ranked by A358378.
A001263 counts unlabeled ordered rooted trees by nodes and leaves.
A358371 counts leaves in standard ordered rooted trees.

Programs

  • Mathematica
    stcinv[q_]:=Total[2^(Accumulate[Reverse[q]])]/2;
    aotrank[t_]:=If[t=={},1,1+stcinv[aotrank/@t]];
    aot[n_]:=If[n==1,{{}},Join@@Table[Tuples[aot/@c],{c,Join@@Permutations/@IntegerPartitions[n-1]}]];
    Table[Sort[aotrank/@aot[n]],{n,6}]

A358377 Numbers k such that the k-th standard ordered rooted tree is a generalized Bethe tree (counted by A003238).

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 9, 11, 16, 17, 32, 37, 43, 64, 128, 129, 137, 171, 256, 257, 293, 512, 529, 683, 1024, 1025, 2048, 2185, 2341, 2731, 4096, 8192, 10923, 16384, 16913, 18725, 32768, 32769, 32897, 34953, 43691, 65536, 65537, 131072, 131329, 149797, 174763
Offset: 1

Views

Author

Gus Wiseman, Nov 14 2022

Keywords

Comments

We define the n-th standard ordered rooted tree to be obtained by taking the (n-1)-th composition in standard order (graded reverse-lexicographic, A066099) as root and replacing each part with its own standard ordered rooted tree. This ranking is an ordered variation of Matula-Goebel numbers, giving a bijective correspondence between positive integers and unlabeled ordered rooted trees.
A generalized Bethe tree is an unlabeled rooted tree where all branches directly under the same root are equal.

Examples

			The terms together with their corresponding ordered rooted trees begin:
    1: o
    2: (o)
    3: ((o))
    4: (oo)
    5: (((o)))
    8: (ooo)
    9: ((oo))
   11: ((o)(o))
   16: (oooo)
   17: ((((o))))
   32: (ooooo)
   37: (((o))((o)))
   43: ((o)(o)(o))
   64: (oooooo)
  128: (ooooooo)
  129: ((ooo))
  137: ((oo)(oo))
  171: ((o)(o)(o)(o))
		

Crossrefs

These trees are counted by A003238.
The unordered version is A214577, also counted by A003238.
A000081 counts unlabeled rooted trees, ranked by A358378.
A358371 and A358372 count leaves and nodes in standard ordered rooted trees.
A358374 ranks ordered identity trees, counted by A032027.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join @@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    srt[n_]:=If[n==1,{},srt/@stc[n-1]];
    Select[Range[1000],FreeQ[srt[#],[_]?(!SameQ@@#&)]&]

A020474 A Motzkin triangle: a(n,k), n >= 2, 2 <= k <= n, = number of complete, strictly subdiagonal staircase functions.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 0, 2, 4, 0, 0, 1, 5, 9, 0, 0, 0, 3, 12, 21, 0, 0, 0, 1, 9, 30, 51, 0, 0, 0, 0, 4, 25, 76, 127, 0, 0, 0, 0, 1, 14, 69, 196, 323, 0, 0, 0, 0, 0, 5, 44, 189, 512, 835, 0, 0, 0, 0, 0, 1, 20, 133, 518, 1353, 2188, 0, 0, 0, 0, 0, 0, 6, 70, 392, 1422, 3610, 5798, 0, 0, 0, 0
Offset: 2

Views

Author

Keywords

Comments

T(n,k) = number of Dyck n-paths that start UU, contain no DUDU and no subpath of the form UUPDD with P a nonempty Dyck path and whose terminal descent has length n-k+2. For example, T(5,4)=2 counts UUDUUDUDDD, UUUDDUUDDD (each ending with exactly n-k+2=3 Ds). - David Callan, Sep 25 2006

Examples

			Triangle begins:
  1
  0, 1
  0, 1, 2
  0, 0, 2, 4
  0, 0, 1, 5,  9
  0, 0, 0, 3, 12, 21
  0, 0, 0, 1,  9, 30, 51
  0, 0, 0, 0,  4, 25, 76, 127
  0, 0, 0, 0,  1, 14, 69, 196, 323
		

Crossrefs

Main diagonal is A001006.
Other diagonals include A002026, A005322, A005323, A005324, A005325. Row sums are (essentially) A005043.
The triangle version of A062105 has the same recurrence with different initial conditions. - N. J. A. Sloane, Apr 11 2020

Programs

  • Haskell
    a020474 n k = a020474_tabl !! (n-2) !! (k-2)
    a020474_row n = a020474_tabl !! (n-2)
    a020474_tabl = map fst $ iterate f ([1], [0,1]) where
       f (us,vs) = (vs, scanl (+) 0 ws) where
         ws = zipWith (+) (us ++ [0]) vs
    -- Reinhard Zumkeller, Jan 03 2013
    
  • Maple
    M:=16; T:=Array(0..M,0..M,0);
    T[0,0]:=1; T[1,1]:=1;
    for i from 1 to M do T[i,0]:=0; od:
    for n from 2 to M do for k from 1 to n do
    T[n,k]:= T[n,k-1]+T[n-1,k-1]+T[n-2,k-1];
    od: od;
    rho:=n->[seq(T[n,k],k=0..n)];
    for n from 0 to M do lprint(rho(n)); od: # N. J. A. Sloane, Apr 11 2020
  • Mathematica
    a[2,2]=1; a[n_,k_]/;Not[n>2 && 2<=k<=n] := 0; a[n_,k_]/;(n>2 && 2<=k<=n) := a[n,k] = a[n,k-1] + a[n-1,k-1] + a[n-2,k-1]; Table[a[n,k],{n,2,10},{k,2,n}] (* David Callan, Sep 25 2006 *)
  • PARI
    T(n,k)=if(n==0&&k==0,1,if(n<=0||k<=0||nRalf Stephan
    
  • Sage
    @cached_function
    def T(n, k):
        if k<0 or nPeter Luschny, Jun 23 2015

Formula

a(n,k) = a(n,k-1) + a(n-1,k-1) + a(n-2,k-1), n > k >= 2.

Extensions

More terms from James Sellers, Feb 04 2000
Previous Showing 41-50 of 190 results. Next