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

A219537 G.f. satisfies A(x) = 1 + x*(A(x)^2 - A(x)^3 + A(x)^4).

Original entry on oeis.org

1, 1, 3, 13, 66, 366, 2148, 13115, 82449, 530095, 3469401, 23037642, 154820262, 1050999343, 7196493255, 49644745965, 344704716018, 2407157839593, 16895247295947, 119121868831235, 843306880720218, 5992060655349521, 42718501097385207, 305476181765843358
Offset: 0

Views

Author

Paul D. Hanna, Nov 21 2012

Keywords

Comments

a(n) is the number of noncrossing partial matchings on points 1, 2, ... , 3*n where point 1 is unmatched if n>0 and only points congruent modulo 3 can be matched. See Example 57 on p. 47 of the Burstein-Shapiro reference. - Alexander Burstein, Jun 03 2022

Examples

			G.f.: A(x) = 1 + x + 3*x^2 + 13*x^3 + 66*x^4 + 366*x^5 + 2148*x^6 +...
Related expansions:
A(x)^2 = 1 + 2*x + 7*x^2 + 32*x^3 + 167*x^4 + 942*x^5 + 5593*x^6 +...
A(x)^3 = 1 + 3*x + 12*x^2 + 58*x^3 + 312*x^4 + 1794*x^5 + 10794*x^6 +...
A(x)^4 = 1 + 4*x + 18*x^2 + 92*x^3 + 511*x^4 + 3000*x^5 + 18316*x^6 +...
A(x)^5 = 1 + 5*x + 25*x^2 + 135*x^3 + 775*x^4 + 4651*x^5 + 28845*x^6 +...
A(x)^6 = 1 + 6*x + 33*x^2 + 188*x^3 + 1116*x^4 + 6852*x^5 + 43204*x^6 +...
where A(x) = 1 + x*(A(x)^2 - A(x)^3 + A(x)^4),
and A(x)^2 = 1 + x*(A(x)^2 + A(x)^5),
and A(x)^3 = 1 + x*(A(x)^2 + A(x)^4 + A(x)^6),
and A(x)^4 = 1 + x*(A(x)^2 + A(x)^4 + A(x)^5 + A(x)^7),
and A(x)^5 = 1 + x*(A(x)^2 + A(x)^4 + A(x)^5 + A(x)^6 + A(x)^8), etc.
The g.f. satisfies A(x) = F(x*A(x)^2) and F(x) = A(x/F(x)^2) where
F(x) = 1 + x + x^2 + 2*x^3 + 4*x^4 + 9*x^5 + 21*x^6 + 51*x^7 +...+ A001006(n-1)*x^n +...
is a g.f. of the Motzkin numbers (A001006, shifted right 1 place).
The g.f. satisfies A(x) = G(x*A(x)) and G(x) = A(x/G(x)) where
G(x) = 1 + x + 2*x^2 + 6*x^3 + 21*x^4 + 80*x^5 + 322*x^6 +...+ A106228(n)*x^n +...
satisfies G(x) = 1 + x*G(x)/(1 - x*G(x)^2).

Links

Crossrefs

Cf. A219534, A219535, A219536, A219538, A363573, A366400.
Cf. A001006, A106228, A271469, A364765.

Programs

  • Maple
    rec := {(36*n^4+126*n^3+126*n^2+36*n)*a(n)+(-276*n^4-1548*n^3-3198*n^2-2898*n-972)*a(n+1)+(940*n^4+7090*n^3+19916*n^2+24650*n+11316)*a(n+2)+(-845*n^4-9000*n^3-34159*n^2-53004*n-26136)*a(n+3)+(-260*n^4-5200*n^3-37454*n^2-116538*n-133128)*a(n+4)+(459*n^4+9774*n^3+77955*n^2+276012*n+366060)*a(n+5)+(-54*n^4-1242*n^3-10686*n^2-40758*n-58140)*a(n+6), a(0) = 1, a(1) = 1, a(2) = 3, a(3) = 13, a(4) = 66, a(5) = 366}:
f:= gfun:-rectoproc(rec,a(n),remember):
map(f, [$0..50]); # Robert Israel, Feb 25 2018
  • Mathematica
    nmax = 23; sol = {a[0] -> 1};
Do[A[x_] = Sum[a[k] x^k, {k, 0, n}] /. sol; eq = CoefficientList[A[x] - (1 + x (A[x]^2 - A[x]^3 + A[x]^4)) + O[x]^(n + 1), x] == 0 /. sol; sol = sol ~Join~ Solve[eq][[1]], {n, 1, nmax}];
sol /. Rule -> Set;
a /@ Range[0, nmax] (* Jean-François Alcover, Nov 01 2019 *)
  • PARI
    /* Formula A(x) = 1 + x*(A(x)^2 - A(x)^3 + A(x)^4): */
{a(n)=local(A=1);for(i=1,n,A=1+x*(A^2-A^3+A^4) +x*O(x^n));polcoeff(A,n)}
for(n=0,25,print1(a(n),", "))
  • PARI
    /* Formula using Series Reversion involving Motzkin numbers: */
{a(n)=local(A=1);A=(1+x-sqrt(1-2*x-3*x^2+x^3*O(x^n)))/(2*x); polcoeff(sqrt(1/x*serreverse(x/A^2)), n)}
for(n=0,25,print1(a(n),", "))

Formula

G.f. A(x) satisfies [from Paul D. Hanna, Mar 21 2016]: (Start)
(1) A(x)^2 = 1 + x*(A(x)^2 + A(x)^5).
(2) A(x)^3 = 1 + x*(A(x)^2 + A(x)^4 + A(x)^6).
Let F(x) = (1+x - sqrt(1 - 2*x - 3*x^2)) / (2*x), then g.f. A(x) satisfies:
(3) A(x) = sqrt( (1/x)*Series_Reversion(x/F(x)^2) ),
(4) A(x) = F(x*A(x)^2) and F(x) = A(x/F(x)^2),
where F(x) = 1 + x*M(x) such that M(x) = 1 + x*M(x) + x^2*M(x)^2 is the g.f. of the Motzkin numbers (A001006).
Let G(x) = 1 + x*G(x)/(1 - x*G(x)^2), then g.f. A(x) satisfies:
(5) A(x) = (1/x)*Series_Reversion(x/G(x)),
(6) A(x) = G(x*A(x)) and G(x) = A(x/G(x)).
where G(x) is the g.f. of A106228. (End)
Recurrence: 3*n*(3*n-1)*(3*n+1)*(5*n-11)*(5*n-8)*(5*n-6)*a(n) = 6*(5*n-11)*(900*n^5 - 3870*n^4 + 6033*n^3 - 4165*n^2 + 1238*n - 120)*a(n-1) - 2*(n-2)*(5*n-1)*(950*n^4 - 5510*n^3 + 11199*n^2 - 9207*n + 2430)*a(n-2) + 6*(n-3)*(n-2)*(2*n-5)*(5*n-6)*(5*n-3)*(5*n-1)*a(n-3). - Vaclav Kotesovec, Aug 19 2013
a(n) ~ sqrt(300+75*10^(2/3)+30*10^(1/3))/90 * (5/9*10^(2/3)+10/9*10^(1/3)+8/3)^n / (sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Aug 19 2013
Recurrence: 18*n*(2*n+1)*(n+2)*(n+1)*a(n)-(6*(n+1))*(n+2)*(46*n^2+120*n+81)*a(n+1)+(2*(n+2))*(470*n^3+2605*n^2+4748*n+2829)*a(n+2)-(n+3)*(845*n^3+6465*n^2+14764*n+8712)*a(n+3)-(2*(n+4))*(130*n^3+2080*n^2+10407*n+16641)*a(n+4)+(3*(n+5))*(153*n^3+2493*n^2+13520*n+24404)*a(n+5)-(6*(n+5))*(3*n+17)*(3*n+19)*(n+6)*a(n+6) = 0. - Robert Israel, Feb 25 2018
G.f. A(x) satisfies: A(-x*A(x)^5) = 1/A(x). - Alexander Burstein, Jun 03 2022
a(n) = (1/n) * Sum_{k=0..n-1} binomial(n,k) * binomial(2*n+k,n-1-k) for n > 0. - Seiichi Manyama, Aug 05 2023
a(n) = (1/n) * Sum_{k=0..n-1} (-1)^k * binomial(n,k) * binomial(4*n-2*k,n-1-k) for n > 0. - Seiichi Manyama, Aug 06 2023
G.f.: A(x) = sqrt(B(x)) where B(x) is the g.f. of A366400. - Seiichi Manyama, Mar 31 2024
a(n) = (1/n) * Sum_{k=0..floor((n-1)/2)} binomial(n,k) * binomial(3*n-k,n-1-2*k) for n > 0. - Seiichi Manyama, Apr 01 2024
a(n) = Sum_{k=0..n} binomial(n,k) * binomial(n+3*k/2+1/2,n)/(2*n+3*k+1). - Seiichi Manyama, Apr 04 2024
G.f.: Sum_{k>=0} binomial(5*k/2 + 1/2, k)*x^k/((5*k + 1)*(1 - x)^((5*k + 1)/2)). - Miles Wilson, Feb 02 2025

A271469 G.f. satisfies A(x) = 1 + x*(A(x)^3 - A(x)^4 + A(x)^5).

Original entry on oeis.org

1, 1, 4, 23, 155, 1142, 8910, 72350, 605056, 5175866, 45077560, 398348733, 3562916317, 32192775763, 293410452560, 2694283228653, 24902681767987, 231496130358758, 2162985033344112, 20301976721356134, 191336242071696514, 1809916398759630481, 17178063381786563194, 163536967014934201972, 1561247114394683682834, 14943175106109268856975
Offset: 0

Views

Author

Paul D. Hanna, Apr 08 2016

Keywords

Examples

			G.f.: A(x) = 1 + x + 4*x^2 + 23*x^3 + 155*x^4 + 1142*x^5 + 8910*x^6 +...
Related expansions:
A(x)^2 = 1 + 2*x + 9*x^2 + 54*x^3 + 372*x^4 + 2778*x^5 + 21873*x^6 +...
A(x)^3 = 1 + 3*x + 15*x^2 + 94*x^3 + 663*x^4 + 5025*x^5 + 39970*x^6 +...
A(x)^4 = 1 + 4*x + 22*x^2 + 144*x^3 + 1041*x^4 + 8016*x^5 + 64470*x^6 +...
A(x)^5 = 1 + 5*x + 30*x^2 + 205*x^3 + 1520*x^4 + 11901*x^5 + 96850*x^6 +...
A(x)^6 = 1 + 6*x + 39*x^2 + 278*x^3 + 2115*x^4 + 16848*x^5 + 138816*x^6 +...
A(x)^7 = 1 + 7*x + 49*x^2 + 364*x^3 + 2842*x^4 + 23044*x^5 + 192325*x^6 +...
where A(x) = 1 + x*(A(x)^3 - A(x)^4 + A(x)^5),
and A(x)^2 = 1 + x*(A(x)^3 + A(x)^6),
and A(x)^3 = 1 + x*(A(x)^3 + A(x)^5 + A(x)^7),
and A(x)^4 = 1 + x*(A(x)^3 + A(x)^5 + A(x)^6 + A(x)^8),
and A(x)^5 = 1 + x*(A(x)^3 + A(x)^5 + A(x)^6 + A(x)^7 + A(x)^9), etc.
The g.f. satisfies A(x) = F(x*A(x)^3) and F(x) = A(x/F(x)^3) where
F(x) = 1 + x + x^2 + 2*x^3 + 4*x^4 + 9*x^5 + 21*x^6 + 51*x^7 +...+ A001006(n-1)*x^n +...
is a g.f. of the Motzkin numbers (A001006, shifted right 1 place).
The g.f. satisfies A(x) = G(x*A(x)) and G(x) = A(x/G(x)) where
G(x) = 1 + x + 3*x^2 + 13*x^3 + 66*x^4 + 366*x^5 + 2148*x^6 +...+ A219537(n)*x^n +...
satisfies G(x) = 1 + x*(G(x)^2 - G(x)^3 + G(x)^4).

Links

Crossrefs

Cf. A363573, A001006, A370474.
Cf. A001006, A106228, A219537, A364765, A365194.
Cf. A001764, A370476.

Programs

  • Mathematica
    CoefficientList[(1/x*InverseSeries[Series[8*x^4/(1 + x - Sqrt[1 - 2*x - 3*x^2])^3, {x, 0, 20}], x])^(1/3), x] (* Vaclav Kotesovec, Apr 16 2016 *)
  • PARI
    /* Formula A(x) = 1 + x*(A(x)^3 - A(x)^4 + A(x)^5): */
{a(n)=local(A=1); for(i=1, n, A=1+x*(A^3-A^4+A^5) +x*O(x^n)); polcoeff(A, n)}
for(n=0, 25, print1(a(n), ", "))
  • PARI
    /* Formula using Series Reversion involving Motzkin numbers: */
{a(n)=local(A=1); A=(1+x-sqrt(1-2*x-3*x^2+x^3*O(x^n)))/(2*x); polcoeff( (1/x*serreverse(x/A^3))^(1/3), n)}
for(n=0, 25, print1(a(n), ", "))

Formula

G.f. A(x) satisfies:
(1) A(x)^2 = 1 + x*(A(x)^3 + A(x)^6).
(2) A(x)^3 = 1 + x*(A(x)^3 + A(x)^5 + A(x)^7).
Let F(x) = (1+x - sqrt(1 - 2*x - 3*x^2)) / (2*x), then g.f. A(x) satisfies:
(3) A(x) = ( (1/x)*Series_Reversion(x/F(x)^3) )^(1/3),
(4) A(x) = F(x*A(x)^3) and F(x) = A(x/F(x)^3),
where F(x) = 1 + x*M(x) such that M(x) = 1 + x*M(x) + x^2*M(x)^2 is the g.f. of the Motzkin numbers (A001006).
Let G(x) = 1 + x*(G(x)^2 - G(x)^3 + G(x)^4), then g.f. A(x) satisfies:
(5) A(x) = (1/x)*Series_Reversion(x/G(x)),
(6) A(x) = G(x*A(x)) and G(x) = A(x/G(x)),
where G(x) is the g.f. of A219537.
a(n) ~ sqrt((34 + (34102 - 8262*sqrt(17))^(1/3) + (34102 + 8262*sqrt(17))^(1/3)) / 1632) * ((28 + (513243 - 4131*sqrt(17))^(1/3)/3 + (19009 + 153*sqrt(17))^(1/3)) / 8)^n / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Apr 16 2016
D-finite recurrence: 8*n*(2*n-1)*(4*n-1)*(4*n+1)*(204*n^4 - 1341*n^3 + 3191*n^2 - 3286*n + 1242)*a(n) = 12*(45696*n^8 - 391776*n^7 + 1376164*n^6 - 2580579*n^5 + 2808064*n^4 - 1797694*n^3 + 651566*n^2 - 119476*n + 8160)*a(n-1) - 6*(n-2)*(29376*n^7 - 237168*n^6 + 760044*n^5 - 1236774*n^4 + 1082233*n^3 - 496791*n^2 + 108530*n - 8400)*a(n-2) + 9*(n-3)*(n-2)*(3*n-8)*(3*n-4)*(204*n^4 - 525*n^3 + 392*n^2 - 111*n + 10)*a(n-3). - Vaclav Kotesovec, Apr 16 2016
From Seiichi Manyama, Aug 06 2023: (Start)
a(n) = (1/n) * Sum_{k=0..n-1} binomial(n,k) * binomial(3*n+k,n-1-k) for n > 0.
a(n) = (1/n) * Sum_{k=0..n-1} (-1)^k * binomial(n,k) * binomial(5*n-2*k,n-1-k) for n > 0. (End)
G.f.: A(x) = sqrt(B(x)) where B(x) is the g.f. of A370474. - Seiichi Manyama, Mar 31 2024
a(n) = (1/n) * Sum_{k=0..floor((n-1)/2)} binomial(n,k) * binomial(4*n-k,n-1-2*k) for n > 0. - Seiichi Manyama, Apr 01 2024
a(n) = Sum_{k=0..n} binomial(n,k) * binomial(3*n/2+3*k/2+1/2,n)/(3*n+3*k+1). - Seiichi Manyama, Apr 04 2024
Showing 1-2 of 2 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.