|
| |
|
|
A106228
|
|
G.f. satisfies A(x) = 1 + x*A(x)/(1 - x*A(x)^2).
|
|
37
|
|
|
|
1, 1, 2, 6, 21, 80, 322, 1347, 5798, 25512, 114236, 518848, 2384538, 11068567, 51817118, 244370806, 1159883685, 5536508864, 26560581688, 127993221140, 619280193640, 3007251366000, 14651743202152, 71601107803904, 350873710447210, 1723795243004223
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
0,3
|
|
|
COMMENTS
|
Number of paths from (0,0) to (3n-3,0) that stay in the first quadrant (but may touch the horizontal axis), consisting of steps u=(2,1), U=(1,2), or d=(1,-1) and have no tripledescents (ddd). Example: a(3)=6 because we have udud, Uddud, udUdd, UddUdd, uudd and Ududd (the remaining four paths contain the string ddd: uUddd, UdUddd, Uuddd and UUdddd; see A027307). - Emeric Deutsch, Jun 08 2005
a(n) = number of node-labeled ordered trees (A000108) on n vertices, each node labeled with a positive integer <= its outdegree. A node is a non-root non-leaf vertex. Example. a(3)=6 counts the 5 ordered trees on 4 vertices with all labels 1 and the tree
.|.
/ \
with its (one and only) node labeled 2. - David Callan, Jul 14 2006
a(n) = number of Schroeder (n-1)-paths with no triple descents. Example: a(4)=21 counts all 22 Schroeder 3-paths (A006318) except UUUDDD. - David Callan, Jul 14 2006
(1 + 2x + 6x^2 + ...)*(1 + x + 2x^2 + 6x^3) = (1 + 3x + 10x^2 + 37x^3 + ...), where A109081 = (1, 1, 3, 10, 37, ...). - Gary W. Adamson, Nov 15 2011
a(n) = number of Motzkin paths of length 2n-1 with no downsteps in odd position. Example: a(3)=6 counts FFFFF, FFUDF, FUFDF, UDFFF, UDUDF, UFFDF with U an upstep (1,1), F a flatstep (1,0), and D a downstep (1,-1). - David Callan, May 20 2015
Number of permutations of length n that avoid 4123, 4132, and 4213. - Jay Pantone, Oct 01 2015
Conjecturally, the number of sequences (e(1), ..., e(n)), 0 <= e(i) < i, such that there is no triple i < j < k with e(i) > e(j) and e(i) <= e(k). [Martinez and Savage, 2.21] - Eric M. Schmidt, Jul 17 2017
a(n) is the number of permutations of length n avoiding the partially ordered pattern (POP) {1>3, 1>4, 4>2} of length 4. That is, the number of length n permutations having no subsequences of length 4 in which the first element is the largest and the fourth element is larger than the second element. - Sergey Kitaev, Dec 10 2020
a(n) is the number of peakless Motzkin paths of length 2n that do not start with an up edge and where every pair of matching up and down edges occupies positions of the same parity. Equivalently, the number of RNA secondary structures on 2n vertices where the leftmost vertex is not matched and only vertices of the same parity can be matched. - Alexander Burstein, May 17 2021
|
|
|
REFERENCES
|
E. Barcucci, R. Pinzani and R. Sprugnoli, The Motzkin family, P.U.M.A. Ser. A, Vol. 2, 1991, No. 3-4, pp. 249-279.
|
|
|
LINKS
|
N. R. Beaton, M. Bouvel, V. Guerrini, and S. Rinaldi, Slicings of Parallelogram Polyominoes: Catalan, Schröder, Baxter, and Other Sequences, Electr. J. Combin. 26 (2019), P3.13. Also at arXiv:1511.04854 [math.CO], 2015-2019.
Emeric Deutsch, Problem 10658, American Math. Monthly, 107, 2000, 368-370.
|
|
|
FORMULA
|
G.f.: A(x) = (1/x)*series_reversion[x/(1 + x*G001006(x))] and thus G.f. satisfies: A(x) = 1 + x*A(x)*G001006(x*A(x)) where G001006(x) is the g.f. of Motzkin numbers A001006.
G.f.: 1 + x*exp( Sum_{n>=1} A082759(n)*x^n/n ), where A082759(n) = Sum_{k=0..n} binomial(n,k)*trinomial(n,k). - Paul D. Hanna, Nov 02 2012
a(n) = (1/n)sum(binomial(n, j+1)*b(n, j), j=0..n-1), where b(n, j) are the trinomial coefficients [b(n, j)=A027907(n, j)=coefficient of x^j in (1+x+x^2)^n]. - Emeric Deutsch, Jun 08 2005
Given g.f. A(x), then B(x) = x*A(x) satisfies 0 = f(x, B(x)) where f(x, y) = y^3 - (1+y)*x*(y-x). - Michael Somos, Jun 18 2005
a(n+1) = Sum[binomial(2n-2k,n-k)*binomial(n+k,n)/(n+1),{k,0,n}]. - David Callan, Aug 16 2006
For n>0: a(n) = 1/n*sum(binomial(n,j)*sum(binomial(j,i)*binomial(n-j,2*j-n-i-1)*2^(2*n-3*j+2*i+1),i=0..n-1), j=0..n); - Vladimir Kruchinin, Dec 26 2010
a(n) = 1/(n+1)*sum(binomial(n+1,k)*binomial(n+k+1,n-k),k,0,n); - Vladimir Kruchinin, Feb 28 2010
a(n) = upper left term in M^n, M = the production matrix:
1, 1
1, 1, 1
2, 2, 1, 1
3, 3, 2, 1, 1
4, 4, 3, 2, 1, 1
5, 5, 4, 3, 2, 1, 1
...
D-finite with recurrence: 4*n*(2*n+1)*a(n) + 2*(6-5*n-10*n^2)*a(n-1) + 12*(-9*n^2+35*n-33)*a(n-2) - 2*(n-3)*(13*n-28)*a(n-3) - 15*(n-3)*(n-4)*a(n-4) = 0. - R. J. Mathar, Nov 14 2011
a(n) is the sum of top row terms of Q^(n-1), where Q = the following infinite square production matrix:
1, 1, 0, 0, 0, ...
2, 1, 1, 0, 0, ...
3, 2, 1, 1, 0, ...
4, 3, 2, 1, 1, ...
5, 4, 3, 2, 1, ...
... (End)
a(n) = 3_F_2([-n, 1-n, n+1], [1, 3/2], 1/4). - Peter Luschny, Aug 02 2012
A four-term recurrence equation is given in the Maple program. Peter Luschny, Aug 03 2012
a(n) ~ 1/228*sqrt(114)*sqrt((32129+3933*sqrt(57))^(1/3) * ((32129+3933*sqrt(57))^(2/3) + 532 + 38*(32129+3933*sqrt(57))^(1/3))) / ((32129+3933*sqrt(57))^(1/3)) * (((1261+57*sqrt(57))^(2/3) + 112 + 10*(1261+57*sqrt(57))^(1/3)) / (6*(1261+57*sqrt(57))^(1/3)))^n / (sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Sep 16 2013
G.f. satisfies x*F(x)^3 - x*F(x)^2 + (x-1)*F(x) + 1 = 0. - Jay Pantone, Oct 01 2015
a(n) = Sum_{k=0..n} binomial(n,k) * binomial(n/2+3*k/2+1/2,n)/(n+3*k+1). - Seiichi Manyama, Apr 04 2024
|
|
|
EXAMPLE
|
A = 1 + x*A + x^2*A^3 + x^3*A^5 + x^4*A^7 + x^5*A^9 + ...
a(4) = 21 since the top row terms of Q^3 = (10, 7, 3, 1). - Gary W. Adamson, Nov 15 2011
G.f. = 1 + x + 2*x^2 + 6*x^3 + 21*x^4 + 80*x^5 + 322*x^6 + 1347*x^7 + ...
|
|
|
MAPLE
|
a := proc(n) option remember; if n<2 then 1 elif n=2 then 2 else ((380*n^3-840*n^2+496*n-72)*a(n-1)+(76*n^3-282*n^2+302*n-84)*a(n-2)+(57*n^3-297*n^2+402*n-72)*a(n-3))/(76*n^3-54*n^2-46*n) fi end: seq(a(i), i=0..23); # Peter Luschny, Aug 03 2012
|
|
|
MATHEMATICA
|
Flatten[{1, Table[1/n*Sum[Binomial[n, k]*Binomial[n+k, n-k-1], {k, 0, n-1}], {n, 1, 20}]}] (* Vaclav Kotesovec, Sep 16 2013 *)
a[ n_] := If[ n < 0, 0, HypergeometricPFQ[{-n, 1 - n, n + 1}, {1, 3/2}, 1/4]]; (* Michael Somos, May 27 2014 *)
|
|
|
PROG
|
(PARI) {a(n)=local(A=1+x+x*O(x^n)); for(k=1, n, A=1+x*A/(1-x*A^2)); polcoeff(A, n)}
(PARI) {a(n) = local(A); if( n<0, 0, n++; A = (1 + x - sqrt(1 - 2*x - 3*x^2 + x * O(x^n))) / 2; polcoeff( serreverse( x^2 / A), n))}; /* Michael Somos, Jun 18 2005 */
(Sage)
from mpmath import mp
mp.dps = 32; mp.pretty = True
def A106228(n) : return int(mp.hyper([-n, 1-n, n+1], [1, 3/2], 1/4))
(PARI) {a(n) = if( n<1, n==0, polcoeff( serreverse( x / (1 + 2*x + 2*x^2 + x^3) + x * O(x^n)), n))}; /* Michael Somos, Dec 31 2014 */
|
|
|
CROSSREFS
|
|
|
|
KEYWORD
|
nonn
|
|
|
AUTHOR
|
|
|
|
STATUS
|
approved
|
| |
|
|
