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-10 of 10 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).
		

Crossrefs

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).
		

Crossrefs

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

A243659 Number of Sylvester classes of 3-packed words of degree n.

Original entry on oeis.org

1, 1, 5, 34, 267, 2279, 20540, 192350, 1853255, 18252079, 182924645, 1859546968, 19127944500, 198725331588, 2082256791048, 21979169545670, 233495834018591, 2494624746580655, 26786319835972799, 288915128642169250, 3128814683222599331, 34007373443388857999
Offset: 0

Views

Author

N. J. A. Sloane, Jun 12 2014

Keywords

Comments

See Novelli-Thibon (2014) for precise definition.

Crossrefs

Column k=3 of A336573.

Programs

  • Maple
    a := proc(n) option remember; if n = 0 then 1 elif n = 1 then 1 else (4*(37604*n^5-158474*n^4+248391*n^3-178459*n^2+58042*n-6720)*a(n-1) - 3*(n-2)*(3*n-4)*(3*n-5)*(119*n^2-85*n+14)*a(n-2) )/ (12*n*(3*n-1)*(3*n+1)*(119*n^2-323*n+218)) fi; end:
    seq(a(n), n = 0..20); # Peter Bala, Sep 08 2024
  • Mathematica
    b[0] = 1; b[n_] := b[n] = 1/n Sum[Sum[2^(j-2i)(-1)^(i-j) Binomial[i, 3i-j] Binomial[i+j-1, i-1], {j, 0, 3i}] b[n-i], {i, 1, n}];
    a[n_] := b[n+1];
    Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Jul 27 2018, after Vladimir Kruchinin *)
  • Maxima
    a(n):=if n=0 then 1 else 1/n*sum(sum(2^(j-2*i)*(-1)^(i-j)*binomial(i,3*i-j)*binomial(i+j-1,i-1),j,0,3*i)*a(n-i),i,1,n); /* Vladimir Kruchinin, Apr 07 2017 */
    
  • PARI
    a(n) = if(n==0, 1, sum(i=1, n, a(n-i)*sum(j=0, 3*i, 2^(j-2*i)*(-1)^(i-j)*binomial(i,3*i-j)*binomial(i+j-1,i-1)))/n); \\ Seiichi Manyama, Jul 26 2020
    
  • PARI
    a(n) = my(A=1+x*O(x^n)); for(i=0, n, A=1-x*A^3*(1-2*A)); polcoeff(A, n); \\ Seiichi Manyama, Jul 26 2020
    
  • PARI
    a(n) = (-1)^n*sum(k=0, n, (-2)^k*binomial(n, k)*binomial(3*n+k+1, n)/(3*n+k+1)); \\ Seiichi Manyama, Jul 26 2020
    
  • PARI
    a(n) = (-1)^n*sum(k=0, n, (-2)^(n-k)*binomial(3*n+1, k)*binomial(4*n-k, n-k))/(3*n+1); \\ Seiichi Manyama, Jul 26 2020

Formula

Novelli-Thibon give an explicit formula in Eq. (182).
a(0) = 1 and a(n) = (1/n) * Sum_{i=1..n} ( Sum_{j=0..3*i} (2^(j-2*i)*(-1)^(i-j) * binomial(i,3*i-j)*binomial(i+j-1,i-1)) *a(n-i) ) for n > 0. - Vladimir Kruchinin, Apr 09 2017
From Seiichi Manyama, Jul 26 2020: (Start)
G.f. A(x) satisfies: A(x) = 1 - x * A(x)^3 * (1 - 2 * A(x)).
a(n) = (-1)^n * Sum_{k=0..n} (-2)^k * binomial(n,k) * binomial(3*n+k+1,n)/(3*n+k+1).
a(n) = ( (-1)^n / (3*n+1) ) * Sum_{k=0..n} (-2)^(n-k) * binomial(3*n+1,k) * binomial(4*n-k,n-k). (End)
a(n) ~ sqrt(24388 + 9221*sqrt(7)) * (316 + 119*sqrt(7))^(n - 1/2) / (sqrt(7*Pi) * n^(3/2) * 2^(n + 3/2) * 3^(3*n + 3/2)). - Vaclav Kotesovec, Jul 31 2021
a(n) = (1/n) * Sum_{k=0..n-1} binomial(n,k) * binomial(4*n-k,n-1-k) for n > 0. - Seiichi Manyama, Aug 08 2023
P-recursive: 12*n*(3*n-1)*(3*n+1)*(119*n^2-323*n+218)*a(n) = 4*(37604*n^5-158474*n^4+248391*n^3-178459*n^2+58042*n-6720)*a(n-1) - (3*n-4)*(3*n-5)*(3*n-6)*(119*n^2-85*n+14)*a(n-2) with a(0) = a(1) = 1. - Peter Bala, Sep 08 2024

Extensions

a(9)-a(21) from Lars Blomberg, Jul 12 2017
a(0)=1 inserted by Seiichi Manyama, Jul 26 2020

A364747 G.f. A(x) satisfies A(x) = 1 + x*A(x)^4 / (1 - x*A(x)).

Original entry on oeis.org

1, 1, 5, 32, 234, 1854, 15490, 134380, 1198944, 10931761, 101412677, 954155059, 9083120975, 87326765375, 846709605539, 8269910074087, 81291388929027, 803592049667495, 7983612883739843, 79671910265120574, 798283229227457304, 8027625597750959053
Offset: 0

Views

Author

Seiichi Manyama, Aug 05 2023

Keywords

Crossrefs

Programs

  • PARI
    a(n) = if(n==0, 1, sum(k=0, n-1, binomial(n, k)*binomial(4*n-3*k, n-1-k))/n);
    
  • PARI
    a(n, r=1, s=1, t=4, u=1) = r*sum(k=0, n, binomial(t*k+u*(n-k)+r, k)*binomial(n+(s-1)*k-1, n-k)/(t*k+u*(n-k)+r)); \\ Seiichi Manyama, Dec 05 2024

Formula

a(n) = (1/n) * Sum_{k=0..n-1} binomial(n,k) * binomial(4*n-3*k,n-1-k) for n > 0.
From Seiichi Manyama, Dec 05 2024: (Start)
G.f. A(x) satisfies A(x) = 1/(1 - x*A(x)^3/(1 - x*A(x))).
If g.f. satisfies A(x) = ( 1 + x*A(x)^(t/r) / (1 - x*A(x)^(u/r))^s )^r, then a(n) = r * Sum_{k=0..n} binomial(t*k+u*(n-k)+r,k) * binomial(n+(s-1)*k-1,n-k)/(t*k+u*(n-k)+r). (End)

A365194 G.f. satisfies A(x) = 1 + x*A(x)^5 / (1 - x*A(x)^6).

Original entry on oeis.org

1, 1, 6, 52, 529, 5889, 69462, 853013, 10791018, 139659604, 1840435530, 24611295075, 333132371248, 4555465710569, 62839303262352, 873363902976309, 12218178082489873, 171918448407833112, 2431415226089290680, 34544425914499450493, 492807213597429920649
Offset: 0

Views

Author

Seiichi Manyama, Aug 25 2023

Keywords

Crossrefs

Programs

  • PARI
    a(n) = sum(k=0, n, binomial(6*n-k+1, k)*binomial(n-1, n-k)/(6*n-k+1));

Formula

a(n) = Sum_{k=0..n} binomial(6*n-k+1,k) * binomial(n-1,n-k)/(6*n-k+1).
a(n) = Sum_{k=0..n} (-1)^(n-k) * binomial(5*n+2*k+1,k) * binomial(n-1,n-k)/(5*n+2*k+1).
a(n) = (1/n) * Sum_{k=0..floor((n-1)/2)} binomial(n,k) * binomial(6*n-k,n-1-2*k) for n > 0. - Seiichi Manyama, Dec 26 2024

A378920 G.f. A(x) satisfies A(x) = 1 + x*A(x)^6/(1 + x*A(x)^2).

Original entry on oeis.org

1, 1, 5, 38, 339, 3308, 34191, 367844, 4076112, 46204209, 533239820, 6244542391, 74016115926, 886276231388, 10704869669941, 130271156244371, 1595708949486866, 19658780721376791, 243429900033986385, 3028086095940468087, 37821457123957529163, 474145963420441744445
Offset: 0

Views

Author

Seiichi Manyama, Dec 11 2024

Keywords

Crossrefs

Programs

  • PARI
    a(n, r=1, s=-1, t=6, u=2) = r*sum(k=0, n, binomial(t*k+u*(n-k)+r, k)*binomial(s*k, n-k)/(t*k+u*(n-k)+r));

Formula

G.f. A(x) satisfies A(x) = 1/(1 - x*A(x)^5/(1 + x*A(x)^2)).
If g.f. satisfies A(x) = ( 1 + x*A(x)^(t/r) * (1 + x*A(x)^(u/r))^s )^r, then a(n) = r * Sum_{k=0..n} binomial(t*k+u*(n-k)+r,k) * binomial(s*k,n-k)/(t*k+u*(n-k)+r).

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

Original entry on oeis.org

1, 1, 5, 33, 250, 2054, 17800, 160183, 1482535, 14022415, 134943095, 1317046306, 13005842030, 129708875695, 1304588594925, 13217663310305, 134775670244250, 1382019265706377, 14242560597119165, 147435736533094415, 1532365596794307010
Offset: 0

Views

Author

Seiichi Manyama, Aug 08 2023

Keywords

Crossrefs

Programs

  • Maple
    A364792 := proc(n)
        if n = 0 then
            1;
        else
            add( binomial(n,k) * binomial(4*n-2*k,n-1-k),k=0..n-1) ;
            %/n ;
        end if ;
    end proc:
    seq(A364792(n),n=0..80); # R. J. Mathar, Aug 10 2023
  • PARI
    a(n) = if(n==0, 1, sum(k=0, n-1, binomial(n, k)*binomial(4*n-2*k, n-1-k))/n);

Formula

a(n) = (1/n) * Sum_{k=0..n-1} binomial(n,k) * binomial(4*n-2*k,n-1-k) for n > 0.
D-finite with recurrence 3*n*(36653*n-48128)*(3*n-1)*(3*n+1)*a(n) +5*(-2160545*n^4 +5139476*n^3 -2463019*n^2 -1385144*n +913296)*a(n-1) +4*(-948403*n^4 +17991137*n^3 -77629283*n^2 +126107767*n -70578450)*a(n-2) +10*(n-3)*(599072*n^3 -5090881*n^2 +13501042*n -11263100)*a(n-3) -50*(6861*n-12886)*(n-3) *(n-4)*(2*n-7)*a(n-4)=0. - R. J. Mathar, Aug 10 2023

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

Original entry on oeis.org

1, 3, 18, 139, 1218, 11511, 114398, 1178421, 12469626, 134734092, 1480317468, 16487870031, 185744716414, 2112756042468, 24230663513604, 279889210974003, 3253295301115290, 38023971948455859, 446603044829013514, 5268557500949993964, 62398899992129490756
Offset: 0

Views

Author

Seiichi Manyama, Dec 11 2024

Keywords

Crossrefs

Programs

  • PARI
    a(n, r=3, s=-1, t=6, u=4) = r*sum(k=0, n, binomial(t*k+u*(n-k)+r, k)*binomial(s*k, n-k)/(t*k+u*(n-k)+r));

Formula

G.f. A(x) satisfies:
(1) A(x) = 1/( 1 - x*A(x)^(5/3)/(1 + x*A(x)^(4/3)) )^3.
(2) A(x) = 1 + x * A(x)^(4/3) * (1 + A(x)^(2/3) + A(x)^(4/3)).
(3) A(x) = B(x)^3 where B(x) is the g.f. of A364765.
If g.f. satisfies A(x) = ( 1 + x*A(x)^(t/r) * (1 + x*A(x)^(u/r))^s )^r, then a(n) = r * Sum_{k=0..n} binomial(t*k+u*(n-k)+r,k) * binomial(s*k,n-k)/(t*k+u*(n-k)+r).

A378919 G.f. A(x) satisfies A(x) = 1 + x*A(x)^6/(1 + x*A(x)).

Original entry on oeis.org

1, 1, 5, 39, 355, 3532, 37206, 407861, 4604493, 53169811, 625067441, 7456004083, 90015754691, 1097834790182, 13505674728174, 167395320811562, 2088350145491232, 26203315734195937, 330460721192844017, 4186559092558049570, 53255890990455126082, 679954025388880445771
Offset: 0

Views

Author

Seiichi Manyama, Dec 11 2024

Keywords

Crossrefs

Programs

  • PARI
    a(n, r=1, s=-1, t=6, u=1) = r*sum(k=0, n, binomial(t*k+u*(n-k)+r, k)*binomial(s*k, n-k)/(t*k+u*(n-k)+r));

Formula

G.f. A(x) satisfies A(x) = 1/(1 - x*A(x)^5/(1 + x*A(x))).
If g.f. satisfies A(x) = ( 1 + x*A(x)^(t/r) * (1 + x*A(x)^(u/r))^s )^r, then a(n) = r * Sum_{k=0..n} binomial(t*k+u*(n-k)+r,k) * binomial(s*k,n-k)/(t*k+u*(n-k)+r).

A378954 G.f. A(x) satisfies A(x) = ( 1 + x*A(x)^3/(1 + x*A(x)^2) )^2.

Original entry on oeis.org

1, 2, 11, 82, 705, 6584, 64902, 664608, 7001006, 75378082, 825810304, 9176278104, 103171720299, 1171558985630, 13416903518301, 154784357304138, 1797153050309355, 20984321920535966, 246252819129444579, 2902768234099178002, 34355158795966317996, 408086199665333171952
Offset: 0

Views

Author

Seiichi Manyama, Dec 12 2024

Keywords

Crossrefs

Programs

  • PARI
    a(n, r=2, s=-1, t=6, u=4) = r*sum(k=0, n, binomial(t*k+u*(n-k)+r, k)*binomial(s*k, n-k)/(t*k+u*(n-k)+r));

Formula

G.f. A(x) satisfies:
(1) A(x) = 1/( 1 - x*A(x)^(5/2)/(1 + x*A(x)^2) )^2.
(2) A(x) = 1 + x * A(x)^2 * (1 + A(x)^(3/2)).
(3) A(x) = B(x)^2 where B(x) is the g.f. of A364765.
If g.f. satisfies A(x) = ( 1 + x*A(x)^(t/r) * (1 + x*A(x)^(u/r))^s )^r, then a(n) = r * Sum_{k=0..n} binomial(t*k+u*(n-k)+r,k) * binomial(s*k,n-k)/(t*k+u*(n-k)+r).
Showing 1-10 of 10 results.