OFFSET
0,5
COMMENTS
A Schroeder path of semilength n is a lattice path from (0,0) to (2n,0) consisting of U = (1,1), D = (1,-1) and H = (2,0) steps and never going below the x-axis.
Row n has 1+n^2 terms.
Row sums are the large Schroeder numbers (A006318).
LINKS
Alois P. Heinz, Rows n = 0..32, flattened
FORMULA
G.f.: G(t,z) satisfies G(t,z) = 1 + z*G(t,z) + t*z*G(t,t^2*z)*G(t,z).
Sum_{k>=0} k*T(n,k) = A129180(n).
From Peter Bala, Aug 02 2019: (Start)
O.g.f. as a continued fraction: (t marks the area and z marks the semilength of the path)
G(t,z) = 1/(1 - z - t*z/(1 - t^2*z - t^3*z/(1 - t^4*z - t^5*z/(1 - t^6*z - (...) )))) = 1 + z*(1 + t) + z^2*(1 + 2*t + t^2 + t^3 + t^4) + ....
G(t,z) = 1/(1 - (1 + t)*z/(1 - t^3*z/(1 - (t^2 + t^5)*z/(1 - t^7*z/(1 - (t^4 + t^9)*z/(1 - t^11*z/( (...) ))))))).
O.g.f. as a ratio of q-series: N(t,z)/D(t,z), where N(t,z) = Sum_{n >= 0} (-1)^n*t^(2*n^2+n)*z^n/( (Product_{k = 1..n} 1 - t^(2*k)) * (Product_{k = 1..n+1} 1 - t^(2*k-2)*z) ) and D(t,z) = Sum_{n >= 0} (-1)^n*t^(2*n^2-n)*z^n/( (Product_{k = 1..n} 1 - t^(2*k)) * (Product_{k = 1..n} 1 - t^(2*k-2)*z) ). (End)
Conjecture: T(n, k) = [z^k] R(n, 0) for n >= 0, k >= 0 where R(n, q) = Sum_{j=0..q + (q mod 2) + 1} z^j*R(n-1, j) for n > 0, q >= 0 with R(0, q) = 1 for q >= 0. - Mikhail Kurkov, Aug 03 2023
Conjecture above can be simplified as follows: T(n, k) = [z^k] P(n+1, n+1) for n >= 0, k >= 0 where P(n, k) = P(n-1, k) + z^(n-k)*P(n + ((n-k) mod 2), k-1) for 0 < k <= n with P(n, k) = 0 for k > n, P(n, 0) = 1 for n >= 0. - Mikhail Kurkov, Oct 02 2024
EXAMPLE
T(3,5) = 3 because we have UDUUDD, UUDDUD and UHHD.
Triangle starts:
1;
1,1;
1,2,1,1,1;
1,3,3,3,4,3,2,1,1,1;
1,4,6,7,10,11,10,9,8,7,5,4,3,2,1,1,1;
MAPLE
G:=1/(1-z-t*z*g[1]): for i from 1 to 11 do g[i]:=1/(1-t^(2*i)*z-t^(2*i+1)*z*g[i+1]) od: g[12]:=0: Gser:=simplify(series(G, z=0, 13)): for n from 0 to 11 do P[n]:=sort(coeff(Gser, z, n)) od: for n from 0 to 6 do seq(coeff(P[n], t, j), j=0..n^2) od; # yields sequence in triangular form
# second Maple program:
b:= proc(x, y) option remember; `if`(y>x or y<0, 0,
`if`(x=0, 1, expand(b(x-1, y-1)*z^(y-1/2)
+b(x-2, y)*z^(2*y) +b(x-1, y+1)*z^(y+1/2))))
end:
T:= n-> (p-> seq(coeff(p, z, i), i=0..degree(p)))(b(2*n, 0)):
seq(T(n), n=0..7); # Alois P. Heinz, May 27 2015
MATHEMATICA
b[x_, y_] := b[x, y] = If[y>x || y<0, 0, If[x==0, 1, Expand[b[x-1, y-1]*z^(y-1/2) + b[x-2, y]*z^(2*y) + b[x-1, y+1]*z^(y+1/2)]]]; T[n_] := Function[{p}, Table[ Coefficient[p, z, i], {i, 0, Exponent[p, z]}]][b[2*n, 0]]; Table[T[n], {n, 0, 7}] // Flatten (* Jean-François Alcover, Jun 29 2015, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn,tabf
AUTHOR
Emeric Deutsch, Apr 08 2007
STATUS
approved