A215624 -id:A215624 - cp's OEIS frontend

A198953 G.f. satisfies A(x) = (1 + xA(x)) (1 + x*A(x)^3).

1, 2, 9, 56, 400, 3095, 25240, 213633, 1859006, 16527544, 149472480, 1370794835, 12718060947, 119158146283, 1125816405458, 10714275588727, 102615375322564, 988302823695146, 9565859385140272, 93000625498797314, 907782305262566776, 8892941663606408172

Offset: 0

Paul D. Hanna, Oct 31 2011

The radius of convergence of g.f. A(x) is r = 0.095007017562450871521918431664620... with A(r) = 1.6228790124092133906198298670423120590101223122... where y=A(r) satisfies 2*y^5 + 6*y^4 - 18*y^3 + 6*y^2 - 3 = 0.

			G.f.: A(x) = 1 + 2*x + 9*x^2 + 56*x^3 + 400*x^4 + 3095*x^5 + 25240*x^6 +...
Related expansions.
A(x)^2 = 1 + 4*x + 22*x^2 + 148*x^3 + 1105*x^4 + 8798*x^5 + 73196*x^6 +...
A(x)^3 = 1 + 6*x + 39*x^2 + 284*x^3 + 2223*x^4 + 18267*x^5 + 155445*x^6 +...
A(x)^4 = 1 + 8*x + 60*x^2 + 472*x^3 + 3878*x^4 + 32948*x^5 + 287300*x^6 +...
where A(x) = 1 + x*(A(x) + A(x)^3) + x^2*A(x)^4.
The logarithm of the g.f. equals the series:
log(A(x)) = (1 + A(x)^2)*x + (1 + 2^2*A(x)^2 + A(x)^4)*x^2/2 +
(1 + 3^2*A(x)^2 + 3^2*A(x)^4 + A(x)^6)*x^3/3 +
(1 + 4^2*A(x)^2 + 6^2*A(x)^4 + 4^2*A(x)^6 + A(x)^8)*x^4/4 +
(1 + 5^2*A(x)^2 + 10^2*A(x)^4 + 10^2*A(x)^6 + 5^2*A(x)^8 + A(x)^10)*x^5/5 +...
more explicitly,
log(A(x)) = 2*x + 14*x^2/2 + 122*x^3/3 + 1118*x^4/4 + 10557*x^5/5 + 101642*x^6/6 + 991916*x^7/7 +...

Vincenzo Librandi, Table of n, a(n) for n = 0..200

Cf. A215623, A215624, A215654, A007863, A036765, A198951, A181734, A073157, A364375.

nmax=20; aa=ConstantArray[0,nmax]; aa[[1]]=2; Do[AGF=1+Sum[aa[[n]]*x^n,{n,1,j-1}]+koef*x^j; sol=Solve[Coefficient[(1 + x*AGF)*(1 + x*AGF^3) - AGF,x,j]==0,koef][[1]];aa[[j]]=koef/.sol[[1]],{j,2,nmax}]; Flatten[{1,aa}] (* Vaclav Kotesovec, Sep 19 2013 *)

a(n):=sum((sum((binomial(2*n+2*k+2,j-k)*binomial(n+2*k,k))/(k+n+1),k,0,j))*(-1)^(n-j)*binomial(2*n-j,n-j),j,0,n); /* Vladimir Kruchinin, Mar 13 2016 */

{a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2*(A+x*O(x^n))^(2*j))*x^m/m))); polcoeff(A, n)}

{a(n)=polcoeff((1/x)*serreverse( 2*x^2*(1+x) / (1 - sqrt(1 - 4*x*(1+x)^2 +x^3*O(x^n)))),n)}

{a(n)=local(A=1+x);for(i=1,n,A=(1 + x*A)*(1 + x*(A+x*O(x^n))^3));polcoeff(A,n)}

G.f. A(x) satisfies:
(1) A(x) = exp( Sum_{n>=1} x^n/n * Sum_{k=0..n} C(n,k)^2 * A(x)^(2*k) ).
(2) A(x) = (1/x)*Series_Reversion( 2*x^2*(1+x) / (1 - sqrt(1 - 4*x*(1+x)^2)) ).
(3) A(x) = G(x*A(x)) where G(x) = A(x/G(x)) is the g.f. of A073157 (Schroeder n-paths containing no FFs).
The formal inverse of g.f. A(x) is (sqrt((1-x^2)^2 + 4*x^3) - (1+x^2)) / (2*x^3).
D-finite with recurrence: 2*n*(n+1)*(2*n+1)*(1275*n^5 - 11696*n^4 + 36827*n^3 - 40618*n^2 - 5828*n + 25368)*a(n) = 6*n*(2*n - 1)*(7650*n^6 - 66351*n^5 + 183953*n^4 - 102147*n^3 - 314787*n^2 + 450754*n - 137760)*a(n-1) - 6*(n-1)*(2*n - 3)*(34425*n^6 - 281367*n^5 + 690471*n^4 - 86579*n^3 - 1831014*n^2 + 2230808*n - 685440)*a(n-2) + 6*(22950*n^8 - 279378*n^7 + 1275447*n^6 - 2461807*n^5 + 518525*n^4 + 5756973*n^3 - 9486182*n^2 + 5962912*n - 1303680)*a(n-3) - 6*(22950*n^8 - 313803*n^7 + 1633059*n^6 - 3736233*n^5 + 1886879*n^4 + 7909228*n^3 - 16107824*n^2 + 11531408*n - 2756544)*a(n-4) + 3*(n-4)*(3*n - 14)*(3*n - 7)*(1275*n^5 - 5321*n^4 + 2793*n^3 + 12437*n^2 - 16992*n + 5328)*a(n-5). - Vaclav Kotesovec, Sep 19 2013
a(n) ~ c*d^n/(sqrt(Pi)*n^(3/2)), where d = 10.5255382776611313... is the root of the equation -27 + 108*d - 108*d^2 + 324*d^3 - 72*d^4 + 4*d^5 = 0 and c = 0.5321376859604656812266678970406658537671... - Vaclav Kotesovec, Sep 19 2013
a(n) = Sum_{j=0..n}((Sum_{k=0..j}((binomial(2*n+2*k+2,j-k)*binomial(n+2*k,k))/(k+n+1)))*(-1)^(n-j)*binomial(2*n-j,n-j)). - Vladimir Kruchinin, Mar 13 2016
a(n) = Sum_{k=0..n} binomial(n+2*k+1,k) * binomial(n+2*k+1,n-k) / (n+2*k+1). - Seiichi Manyama, Jul 19 2023

A215654 G.f. satisfies: A(x) = (1 + xA(x)^2) (1 + x*A(x)^3).

Original entry on oeis.org

1, 2, 11, 81, 684, 6257, 60325, 603641, 6210059, 65272503, 697898849, 7566847547, 82999675563, 919376968734, 10269588489433, 115548651723889, 1308374198000780, 14897993185500455, 170482798370871370, 1959574731164246402, 22614008012647634411, 261915716386286916342

Offset: 0

Paul D. Hanna, Aug 19 2012

nonn

More generally, for fixed parameters p, q, r, and s, if F(x) satisfies:
F(x) = (1 + x^r*F(x)^(p+1)) * (1 + x^(r+s)*F(x)^(p+q+1)), then
F(x) = exp( Sum_{n>=1} x^(n*r)*F(x)^(n*p)/n * [Sum_{k=0..n} C(n,k)^2 * x^(k*s)*F(x)^(k*q)] ).
The radius of convergence of g.f. A(x) is r = 0.08035832347291483065438962031... with A(r) = 1.5393913914574609282262181402132760790902539070... where y=A(r) satisfies 20*y^3 - 38*y^2 + 15*y - 6 = 0.
r = 1/(187/300*17^(2/3) + 119/75*17^(1/3) + 1273/300). - Vaclav Kotesovec, Sep 17 2013
Number of hybrid ternary trees with n internal nodes. [Hong and Park]. - N. J. A. Sloane, Mar 26 2014

			G.f.: A(x) = 1 + 2*x + 11*x^2 + 81*x^3 + 684*x^4 + 6257*x^5 + 60325*x^6 +...
Related expansions.
A(x)^2 = 1 + 4*x + 26*x^2 + 206*x^3 + 1813*x^4 + 17032*x^5 +...
A(x)^3 = 1 + 6*x + 45*x^2 + 383*x^3 + 3519*x^4 + 34023*x^5 +...
A(x)^5 = 1 + 10*x + 95*x^2 + 925*x^3 + 9270*x^4 + 95237*x^5 +...
where A(x) = 1 + x*(A(x)^2 + A(x)^3) + x^2*A(x)^5.
The g.f. also satisfies the series:
A(x) = 1 + 2*x*A(x)^2 + 3*x^2*A(x)^4 + 5*x^3*A(x)^6 + 8*x^4*A(x)^8 + 13*x^5*A(x)^10 + 21*x^6*A(x)^12 + 34*x^7*A(x)^14 +...+ Fibonacci(n+2)*x^n*A(x)^(2*n) +...
and consequently, A( x*(1-x-x^2)^2/(1+x)^2 ) = (1+x)/(1-x-x^2).
The logarithm of the g.f. equals the series:
log(A(x)) = (1 + A(x))*x*A(x) + (1 + 2^2*A(x) + A(x)^2)*x^2*A(x)^2/2 +
(1 + 3^2*A(x) + 3^2*A(x)^2 + A(x)^3)*x^3*A(x)^3/3 +
(1 + 4^2*A(x) + 6^2*A(x)^2 + 4^2*A(x)^3 + A(x)^4)*x^4*A(x)^4/4 +
(1 + 5^2*A(x) + 10^2*A(x)^2 + 10^2*A(x)^3 + 5^2*A(x)^4 + A(x)^5)*x^5*A(x)^5/5 +...
Explicitly,
log(A(x)) = 2*x + 18*x^2/2 + 185*x^3/3 + 2006*x^4/4 + 22412*x^5/5 + 255249*x^6/6 + 2946155*x^7/7 + 34342270*x^8/8 +...+ L(n)*x^n/n +...
where L(n) = [x^n] (1+x)^(2*n)/(1-x-x^2)^(2*n) / 2.

Vincenzo Librandi, Table of n, a(n) for n = 0..200
SeoungJi Hong and SeungKyung Park, Hybrid d-ary trees and their generalization, Bull. Korean Math. Soc. 51 (2014), No. 1, pp. 229-235. See p. 233. - _N. J. A. Sloane_, Mar 26 2014
Sheng-liang Yang and Mei-yang Jiang, Pattern avoiding problems on the hybrid d-trees, J. Lanzhou Univ. Tech., (China, 2023) Vol. 49, No. 2, 144-150. (in Mandarin)

Cf. A007863, A198953, A215624, A198951.

Column k=3 of A245049.

a:= n-> coeff(series(RootOf((1+x*A^2)*(1+x*A^3)-A, A), x, n+1), x, n):
seq(a(n), n=0..33);  # Alois P. Heinz, Apr 04 2019

CoefficientList[Sqrt[1/x*InverseSeries[Series[x*(1-x-x^2)^2/(1+x)^2,{x,0,20}],x]],x] (* Vaclav Kotesovec, Sep 17 2013 *)

a(n):=sum(binomial(2*n+i,i)*binomial(2*n+i+1,n-i),i,0,n)/(2*n+1); /* Vladimir Kruchinin, Apr 04 2019 */

{a(n)=local(A=1+x+x*O(x^n)); for(i=1, n, A=(1 + x*A^2)*(1 + x*A^3)); polcoeff(A, n)}

{a(n)=polcoeff(sqrt((1/x)*serreverse( x*(1-x-x^2)^2/(1+x +x*O(x^n))^2)), n)}
for(n=0,31,print1(a(n),", "))

{a(n)=local(A=1+x+x*O(x^n)); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2*A^j)*x^m*A^m/m))); polcoeff(A, n)}

{a(n)=local(A=1+x+x*O(x^n)); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2/A^j)*x^m*A^(2*m)/m))); polcoeff(A, n)}

{a(n)=polcoeff(((1+x)/(1-x-x^2 +x*O(x^n)))^(2*n+1)/(2*n+1),n)}

G.f. A(x) satisfies:
(1) A(x) = sqrt( (1/x)*Series_Reversion( x*(1-x-x^2)^2/(1+x)^2 ) ).
(2) A(x) = exp( Sum_{n>=1} x^n*A(x)^n/n * Sum_{k=0..n} C(n,k)^2 * A(x)^k ).
(3) A(x) = exp( Sum_{n>=1} x^n*A(x)^(2*n)/n * Sum_{k=0..n} C(n,k)^2 / A(x)^k ).
(4) A(x) = Sum_{n>=0} Fibonacci(n+2) * x^n * A(x)^(2*n).
(5) A(x) = G(x*A(x)) where G(x) = A(x/G(x)) is the g.f. of A007863 (number of hybrid binary trees with n internal nodes).
The formal inverse of g.f. A(x) is (sqrt(1-2*x+5*x^2) - (1+x))/(2*x^3).
a(n) = [x^n] ( (1+x)/(1-x-x^2) )^(2*n+1) / (2*n+1).
Recurrence: 100*(n-1)*n*(2*n-1)*(2*n+1)*(4913*n^3 - 26877*n^2 + 49912*n - 30480)*a(n) = 2*(n-1)*(2*n-1)*(6254249*n^5 - 40468670*n^4 + 99110119*n^3 - 109861414*n^2 + 52822608*n - 8566560)*a(n-1) - 3*(2343501*n^7 - 22194333*n^6 + 87905623*n^5 - 187987155*n^4 + 233161624*n^3 - 166253172*n^2 + 62010112*n - 8952000)*a(n-2) + 6*(n-2)*(2*n-5)*(3*n-8)*(3*n-4)*(4913*n^3 - 12138*n^2 + 10897*n - 2532)*a(n-3). - Vaclav Kotesovec, Sep 17 2013
a(n) ~ 1/1020*sqrt(73695 + 11730*17^(2/3) + 28815*17^(1/3)) * (187/300*17^(2/3) + 119/75*17^(1/3) + 1273/300)^n / (n^(3/2)*sqrt(Pi)). - Vaclav Kotesovec, Sep 17 2013
a(n) = 1/(2*n+1)*Sum_{i=0..n} C(2*n+i,i)*C(2*n+i+1,n-i). - Vladimir Kruchinin, Apr 04 2019

A364338 G.f. satisfies A(x) = (1 + x) * (1 + x*A(x)^5).

Original entry on oeis.org

1, 2, 11, 105, 1140, 13555, 170637, 2235472, 30161255, 416248640, 5848462880, 83378361111, 1203100853951, 17537182300140, 257858115407535, 3819894878557990, 56958234329850060, 854192593184162160, 12875579347191388830, 194963091634569681550, 2964229359714424159370, 45234864131654311730160

Offset: 0

Seiichi Manyama, Jul 19 2023

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..827

Cf. A073157, A364336, A364337, A364339.

Cf. A215624, A234525, A239108, A364331, A364335.

terms = 22; A[] = 0; Do[A[x] = (1+x)(1+x*A[x]^5) + O[x]^terms // Normal, terms]; CoefficientList[A[x], x] (* Stefano Spezia, Mar 24 2025 *)

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

a(n) = Sum_{k=0..n} binomial(5*k+1,k) * binomial(5*k+1,n-k) / (5*k+1).

A215623 G.f. satisfies A(x) = (1 + xA(x)) (1 + x*A(x)^4).

Original entry on oeis.org

1, 2, 11, 89, 836, 8551, 92445, 1039030, 12019135, 142151324, 1711116646, 20894534324, 258195565959, 3222677162409, 40569811695707, 514520507077695, 6567611974106756, 84310605465652750, 1087798325715407703, 14098475168420865396, 183465816241394787196

Offset: 0

Paul D. Hanna, Aug 17 2012

The radius of convergence of g.f. A(x) is r = 0.0712256396327314729661274986100... with A(r) = 1.4248895273944523042559975726479124492235978714420... where y=A(r) satisfies 3*y^7 - 4*y^6 + 16*y^5 - 28*y^4 + 8*y^3 - 4 = 0.

			G.f.: A(x) = 1 + 2*x + 11*x^2 + 89*x^3 + 836*x^4 + 8551*x^5 + 92445*x^6 + ...
Related expansions.
A(x)^4 = 1 + 8*x + 68*x^2 + 652*x^3 + 6750*x^4 + 73544*x^5 + 831078*x^6 + ...
A(x)^5 = 1 + 10*x + 95*x^2 + 965*x^3 + 10350*x^4 + 115507*x^5 + ...
where A(x) = 1 + x*(A(x) + A(x)^4) + x^2*A(x)^5.
The logarithm of the g.f. equals the series:
log(A(x)) = (1 + A(x)^3)*x + (1 + 2^2*A(x)^3 + A(x)^6)*x^2/2 +
  (1 + 3^2*A(x)^3 + 3^2*A(x)^6 + A(x)^9)*x^3/3 +
  (1 + 4^2*A(x)^3 + 6^2*A(x)^6 + 4^2*A(x)^9 + A(x)^12)*x^4/4 +
  (1 + 5^2*A(x)^3 + 10^2*A(x)^6 + 10^2*A(x)^9 + 5^2*A(x)^12 + A(x)^15)*x^5/5 + ...
more explicitly,
log(A(x)) = 2*x + 18*x^2/2 + 209*x^3/3 + 2550*x^4/4 + 32082*x^5/5 + 411705*x^6/6 + 5356416*x^7/7 + ....

Seiichi Manyama, Table of n, a(n) for n = 0..876

Cf. A198953, A215624, A036765, A198951, A181734, A364331, A364336, A364376.

{a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2*(A+x*O(x^n))^(3*j))*x^m/m))); polcoeff(A, n)}

{a(n)=local(A=1+x+x*O(x^n)); for(i=1, n, A=(1+x*A)*(1+x*A^4)+x*O(x^n)); polcoeff(A, n)}
for(n=0,21,print1(a(n),", "))

G.f. satisfies A(x) = exp( Sum_{n>=1} x^n/n * Sum_{k=0..n} C(n,k)^2 * A(x)^(3*k)).
The formal inverse of g.f. A(x) is (sqrt((1-x^3)^2 + 4*x^4) - (1+x^3))/(2*x^4).
a(n) = Sum_{k=0..n} binomial(n+3*k+1,k) * binomial(n+3*k+1,n-k) / (n+3*k+1). - Seiichi Manyama, Jul 19 2023
From  Peter Bala, Sep 10 2024: (Start)
x/series_reversion(x*A(x)) = 1 + 2*x + 7*x^2 + 39*x^3 + 242*x^4 + 1634*x^5 + ..., the g.f. of A364336.
(1/x) * series_reversion(x/A(x)) = 1 + 2*x + 15*x^2 + 163*x^3 + 2070*x^4 + 28698*x^5 + ..., the g.f. of A364331. (End)

A364331 G.f. satisfies A(x) = (1 + xA(x)^2) (1 + x*A(x)^5).

Original entry on oeis.org

1, 2, 15, 163, 2070, 28698, 421015, 6425644, 100977137, 1622885389, 26551709946, 440744175801, 7404449354076, 125657625548824, 2150963575012295, 37094953102567208, 643904274979347286, 11241232087809137759, 197247501440314516840, 3476787208220672891388, 61533794803235280779261

Offset: 0

Seiichi Manyama, Jul 18 2023

Cf. A007863, A069271, A073157, A215654, A215715, A364333.
Cf. A215623, A215624, A239108, A364335, A364338.
Cf. A200719.

A364331 := proc(n)
    add( binomial(2*n+3*k+1,k) * binomial(2*n+3*k+1,n-k)/(2*n+3*k+1),k=0..n) ;
end proc:
seq(A364331(n),n=0..70); # R. J. Mathar, Jul 25 2023

a(n) = sum(k=0, n, binomial(2*n+3*k+1, k)*binomial(2*n+3*k+1, n-k)/(2*n+3*k+1));

a(n) = Sum_{k=0..n} binomial(2*n+3*k+1,k) * binomial(2*n+3*k+1,n-k) / (2*n+3*k+1).

x/series_reversion(x*A(x)) = 1 + 2*x + 11*x^2 + 89*x^3 + 836*x^4 + ..., the g.f. of A215623. - Peter Bala, Sep 08 2024

A216359 G.f. satisfies: A(x) = (1 + xA(x)^2) (1 + x/A(x)).

Original entry on oeis.org

1, 2, 3, 13, 32, 147, 445, 2067, 7019, 32590, 119209, 551551, 2125429, 9795863, 39221165, 180177403, 742575760, 3403131833, 14342166121, 65626369612, 281459129188, 1286834885967, 5596229192396, 25580269950635, 112492633046446, 514323765191879, 2282371511598955

Offset: 0

Paul D. Hanna, Sep 04 2012

The radius of convergence of g.f. A(x) is r = 0.209619875959405379599013693... with A(r) = 2.36951367232829409921688546894691317519410... where y=A(r) satisfies y^7 - 2*y^6 - 4*y^4 + 4*y^3 + 4*y - 2 = 0.

			G.f.: A(x) = 1 + 2*x + 3*x^2 + 13*x^3 + 32*x^4 + 147*x^5 + 445*x^6 +...
Related expansions.
A(x)^2 = 1 + 4*x + 10*x^2 + 38*x^3 + 125*x^4 + 500*x^5 + 1839*x^6 +...
A(x)^3 = 1 + 6*x + 21*x^2 + 83*x^3 + 315*x^4 + 1269*x^5 + 5061*x^6 +...
where A(x) = (1-x^2)*A(x)^2 - x*A(x)^3 - x.
The logarithm of the g.f. equals the series:
log(A(x)) = (1 + 1/A(x)^3)*x*A(x) + (1 + 2^2/A(x)^3 + 1/A(x)^6)*x^2*A(x)^2/2 +
(1 + 3^2/A(x)^3 + 3^2/A(x)^6 + 1/A(x)^9)*x^3*A(x)^3/3 +
(1 + 4^2/A(x)^3 + 6^2/A(x)^6 + 4^2/A(x)^9 + 1/A(x)^12)*x^4*A(x)^4/4 +
(1 + 5^2/A(x)^3 + 10^2/A(x)^6 + 10^2/A(x)^9 + 5^2/A(x)^12 + 1/A(x)^15)*x^5*A(x)^5/5 +...

Robert Israel, Table of n, a(n) for n = 0..1000

Cf. A073157, A198953, A215623, A215624, A364331, A364336.

S:= series(RootOf(x+y+x^2*y^2-y^2+x*y^3, y, 1), x, 41):
seq(coeff(S,x,j),j=0..40); # Robert Israel, Jul 10 2017

nmax=20; aa=ConstantArray[0,nmax]; aa[[1]]=2; Do[AGF=1+Sum[aa[[n]]*x^n,{n,1,j-1}]+koef*x^j; sol=Solve[Coefficient[(1-x^2)*AGF^2 - x*AGF^3 - x - AGF,x,j]==0,koef][[1]];aa[[j]]=koef/.sol[[1]],{j,2,nmax}]; Flatten[{1,aa}] (* Vaclav Kotesovec, Sep 18 2013 *)

{a(n)=local(A=1+x); for(i=1, n, A=(1 + x*A^2)*(1 + x/(A+x*O(x^n)))); polcoeff(A, n)}

{a(n)=local(A=1+x+x*O(x^n)); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2/A^(3*j))*x^m*A^m/m))); polcoeff(A, n)}
for(n=0, 31, print1(a(n), ", "))

G.f. satisfies:
A(x) = exp( Sum_{n>=1} x^n*A(x)^n/n * Sum_{k=0..n} C(n,k)^2 / A(x)^(3*k) ).
The formal inverse of the g.f. A(x) is (sqrt(1 - 2*x^3 + 4*x^4 + x^6) - (1+x^3))/(2*x^2).
Recurrence: n*(n+1)*(1241*n^5 - 21306*n^4 + 135203*n^3 - 381522*n^2 + 435524*n - 104880)*a(n) = 6*n*(1201*n^4 - 19476*n^3 + 114613*n^2 - 287442*n + 255364)*a(n-1) + 2*(12410*n^7 - 237880*n^6 + 1771109*n^5 - 6388366*n^4 + 11032829*n^3 - 6363274*n^2 - 3856020*n + 4157712)*a(n-2) + 6*(2482*n^7 - 51299*n^6 + 419427*n^5 - 1705769*n^4 + 3477465*n^3 - 2797370*n^2 - 637684*n + 1410288)*a(n-3) + 2*(4964*n^7 - 110044*n^6 + 983093*n^5 - 4442260*n^4 + 10160177*n^3 - 8790970*n^2 - 4722180*n + 9233280)*a(n-4) - 6*(2482*n^7 - 58745*n^6 + 553921*n^5 - 2617109*n^4 + 6255337*n^3 - 6022682*n^2 - 1392300*n + 4289616)*a(n-5) + 60*(n-7)*(2*n - 11)*(n^3 - 40*n^2 + 280*n - 552)*a(n-6) + 2*(n-8)*(2*n - 13)*(1241*n^5 - 15101*n^4 + 62389*n^3 - 91339*n^2 - 930*n + 64260)*a(n-7). - Vaclav Kotesovec, Sep 18 2013
a(n) ~ c*d^n/(sqrt(Pi)*n^(3/2)), where d = 4.77053998540509708... is the root of the equation -4 + 12*d^2 - 8*d^3 - 12*d^4 - 20*d^5 + d^7 = 0 and c = 1.27852844884923435863262213680985089152... - Vaclav Kotesovec, Sep 18 2013
In closed form, c = (-4 + (1 + sqrt(1+8/d^2))*d^2) * sqrt((d^3*(1 + sqrt(1+8/d^2) + (4*(4 + d^2*(-3-sqrt(1+8/d^2) + d*(4+d))))/d^6)) / (1 + 1/64*(1 + sqrt(1+8/d^2)-4/d^2)^3*d^3)) / (32*d). - Vaclav Kotesovec, Aug 18 2014
From Peter Bala, Sep 10 2024: (Start)
For n not of the form 3*m + 1, we conjecture that a(n) = Sum_{k = 0..n} binomial(-n+3*k+1, k)*binomial(-n+3*k+1, n-k)/(-n+3*k+1).
Define a sequence operator R: {u(n): n >= 0} -> {v(n): n >= 0} by Sum_{n >= 0} v(n)*x^n = (1/x) * series_reversion(x/Sum_{n >= 0} u(n)*x^n). Then R({a(n)}) = A364336, R^2({a(n)}) = A215623 and R^3({a(n)}) = A364331. Cf. A073157. (End)

cp's OEIS Frontend

A198953 G.f. satisfies A(x) = (1 + x*A(x)) * (1 + x*A(x)^3).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Maxima

PARI

PARI

PARI

Formula

A215654 G.f. satisfies: A(x) = (1 + x*A(x)^2) * (1 + x*A(x)^3).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Maxima

PARI

PARI

PARI

PARI

PARI

Formula

A364338 G.f. satisfies A(x) = (1 + x) * (1 + x*A(x)^5).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A215623 G.f. satisfies A(x) = (1 + x*A(x)) * (1 + x*A(x)^4).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

PARI

Formula

A364331 G.f. satisfies A(x) = (1 + x*A(x)^2) * (1 + x*A(x)^5).

Views

Author

Keywords

Crossrefs

Programs

Maple

PARI

Formula

A216359 G.f. satisfies: A(x) = (1 + x*A(x)^2) * (1 + x/A(x)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

A364335 G.f. satisfies A(x) = (1 + x*A(x)^3) * (1 + x*A(x)^5).

A198953 G.f. satisfies A(x) = (1 + xA(x)) (1 + x*A(x)^3).

A215654 G.f. satisfies: A(x) = (1 + xA(x)^2) (1 + x*A(x)^3).

A215623 G.f. satisfies A(x) = (1 + xA(x)) (1 + x*A(x)^4).

A364331 G.f. satisfies A(x) = (1 + xA(x)^2) (1 + x*A(x)^5).

A216359 G.f. satisfies: A(x) = (1 + xA(x)^2) (1 + x/A(x)).

A364335 G.f. satisfies A(x) = (1 + xA(x)^3) (1 + x*A(x)^5).

A364340 G.f. satisfies A(x) = (1 + xA(x)) (1 + x*A(x)^6).