Displaying 1-10 of 14 results found.
G.f. satisfies A(x) = 1 + x*A(x)^2 + x^2*A(x)^4 + x^3*A(x)^6.
+0
11
1, 1, 3, 12, 54, 262, 1337, 7072, 38426, 213197, 1202795, 6879160, 39794416, 232429030, 1368806610, 8118934656, 48458809586, 290832756606, 1754059333738, 10625545472716, 64620970743082, 394409682103262, 2415084675723048, 14832185219521152, 91339478577683664
FORMULA
a(n) = (1/(2*n-1))*Sum_{j=0..2*n-1} binomial(2*n-1,j)*Sum_{i=j..n+j-1} binomial(j,i-j)*binomial(2*n-j-1,3*j-3*n-i+1), n>0.
G.f. A(x) satisfies:
(1) A(x) = sqrt( (1/x)*Series_Reversion( x/(1 + x + x^2 + x^3)^2 ) ).
(2) A( x/(1 + x + x^2 + x^3)^2 ) = 1 + x + x^2 + x^3.
(3) A(x) = G(x*A(x)) where G(x) = A(x/G(x)) = g.f. of A036765 (number of rooted trees with a degree constraint).
(4) a(n) = [x^n] (1 + x + x^2 + x^3)^(2*n+1) / (2*n+1).
(5) A(x) = exp( Sum_{n>=1} x^n*A(x)^n/n * [Sum_{k=0..n} C(n,k)^2 * x^k*A(x)^(2*k)] ).
(6) A(x) = exp( Sum_{n>=1} x^n*A(x)^n/n * [(1-x*A(x)^2)^(2*n+1)*Sum_{k>=0} C(n+k,k)^2*x^k*A(x)^(2*k) )] ).
(End)
a(n) = 1/(2*n + 1)*Sum_{k = 0..floor(n/2)} binomial(2*n + 1,k)*binomial(2*n + 1,n - 2*k).
More generally, the coefficient of x^n in A(x)^r equals r/(2*n + r)*Sum_{k = 0..floor(n/2)} binomial(2*n + r,k)*binomial(2*n + r,n - 2*k) by the Lagrange-Bürmann formula.
O.g.f. A(x) = exp(Sum_{n >= 1} 1/2*b(n)*x^n/n), where b(n) = Sum_{k = 0..floor(n/2)} binomial(2*n,k)*binomial(2*n,n - 2*k). Cf. A036765, A198951, A200731. (End)
Recurrence: 5*n*(5*n - 1)*(5*n + 1)*(5*n + 2)*(5*n + 3)*(13144*n^4 - 57784*n^3 + 90149*n^2 - 59354*n + 13980)*a(n) = 8*(2*n - 1)*(16259128*n^8 - 71478808*n^7 + 108653137*n^6 - 60530902*n^5 - 2811173*n^4 + 12694433*n^3 - 2398482*n^2 - 352503*n + 78570)*a(n-1) + 128*(n-1)*(2*n - 3)*(2*n - 1)*(52576*n^6 - 178560*n^5 + 136156*n^4 + 22938*n^3 - 16067*n^2 - 3138*n - 405)*a(n-2) + 2048*(n-2)*(n-1)*(2*n - 5)*(2*n - 3)*(2*n - 1)*(13144*n^4 - 5208*n^3 - 4339*n^2 + 168*n + 135)*a(n-3). - Vaclav Kotesovec, Nov 17 2017
A(x^2) = (1/x) * series reversion of x/(1 + x^2 + x^4 + x^6). - Peter Bala, Jul 27 2023
EXAMPLE
G.f.: A(x) = 1 + x + 3*x^2 + 12*x^3 + 54*x^4 + 262*x^5 + 1337*x^6 +...
where A(x) = (1 + x*A(x)^2)*(1 + x^2*A(x)^4).
Related expansions:
A(x)^2 = 1 + 2*x + 7*x^2 + 30*x^3 + 141*x^4 + 704*x^5 + 3666*x^6 +...
A(x)^4 = 1 + 4*x + 18*x^2 + 88*x^3 + 451*x^4 + 2392*x^5 + 13022*x^6 +...
A(x)^6 = 1 + 6*x + 33*x^2 + 182*x^3 + 1014*x^4 + 5718*x^5 + 32623*x^6 +...
where A(x) = 1 + x*A(x)^2 + x^2*A(x)^4 + x^3*A(x)^6.
The logarithm of the g.f. A = A(x) equals the series:
log(A(x)) = (1 + x*A^2)*x*A + (1 + 2^2*x*A^2 + x^2*A^4)*x^2*A^2/2 +
(1 + 3^2*x*A^2 + 3^2*x^2*A^4 + x^3*A^6)*x^3*A^3/3 +
(1 + 4^2*x*A^2 + 6^2*x^2*A^4 + 4^2*x^3*A^6 + x^4*A^8)*x^4*A^4/4 +
(1 + 5^2*x*A^2 + 10^2*x^2*A^4 + 10^2*x^3*A^6 + 5^2*x^4*A^8 + x^5*A^10)*x^5*A^5/5 + ...
which involves squares of binomial coefficients. (End)
MAPLE
F:= proc(n) if n::even then
simplify((1/2)*hypergeom([-(1/2)*n, -2*n-1, -(1/2)*n+1/2], [(1/2)*n+1, 3/2+(1/2)*n], -1)*(2*n+2)!/((2*n+1)*((n+1)!)^2))
else
simplify((1/2)*hypergeom([-(1/2)*n, -2*n-1, -(1/2)*n+1/2], [(1/2)*n+1, 3/2+(1/2)*n], -1)*(2*n+2)!/((2*n+1)*((n+1)!)^2))
fi
end proc:
MATHEMATICA
a[n_] := 1/(2n + 1) Sum[Binomial[2n + 1, k] Binomial[2n + 1, n - 2k], {k, 0, n/2}];
(* or: *)
a[n_] := (Binomial[2n + 1, n] HypergeometricPFQ[{-2n - 1, 1/2 - n/2, -n/2}, {n/2 + 1, n/2 + 3/2}, -1])/(2n + 1);
PROG
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=(1+x*A^2)*(1+x^2*A^4)+x*O(x^n)); polcoeff(A, n)} /* Paul D. Hanna */
(PARI) {a(n)=polcoeff(sqrt((1/x)*serreverse(x/(1 + x + x^2 + x^3 +x*O(x^n))^2)), n)} /* Paul D. Hanna */
(PARI) {a(n)=polcoeff( (1 + x + x^2 + x^3+x*O(x^n))^(2*n+1)/(2*n+1), n)} /* Paul D. Hanna */
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, (x*A+x*O(x^n))^m/m*sum(j=0, m, binomial(m, j)^2*x^j*A^(2*j))))); polcoeff(A, n, x)} /* Paul D. Hanna */
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, x^m*A^m/m*(1-x*A^2)^(2*m+1)*sum(j=0, n, binomial(m+j, j)^2*x^j*A^(2*j))))); polcoeff(A, n, x)} /* Paul D. Hanna */
G.f. satisfies A(x) = (1 + x*A(x)^2) * (1 + x^2*A(x)^5).
+0
12
1, 1, 3, 13, 64, 340, 1903, 11053, 65993, 402527, 2497439, 15712220, 100001459, 642719263, 4165537744, 27193644061, 178654643151, 1180282875483, 7836312619243, 52259258911091, 349902441457427, 2351240866736891, 15851508780927739, 107187240225220684, 726784821098903319
COMMENTS
More generally, for fixed parameters p and q, if F(x) satisfies:
F(x) = exp( Sum_{n>=1} x^n * F(x)^(n*p)/n * [Sum_{k=0..n} C(n,k)^2 * x^k * F(x)^(k*q)] ),
then F(x) = (1 + x*F(x)^(p+1))*(1 + x^2*F(x)^(p+q+1)).
FORMULA
G.f. A(x) satisfies:
(1) A(x) = exp( Sum_{n>=1} x^n * A(x)^n/n * [Sum_{k=0..n} C(n,k)^2 * x^k * A(x)^(3*k)] ).
(2) A(x) = exp( Sum_{n>=1} x^n * A(x)^n/n * [(1-x/A(x)^3)^(2*n+1) * Sum_{k>=0} C(n+k,k)^2*x^k * A(x)^(3*k)] ).
Recurrence: 4232*(n-2)*(n-1)*n*(2*n - 3)*(2*n - 1)*(2*n + 1)*(108983978975*n^7 - 1828734495225*n^6 + 13017379495661*n^5 - 50928975062019*n^4 + 118201965098732*n^3 - 162617590602876*n^2 + 122676758610192*n - 39103265134080)*a(n) = 8*(n-2)*(n-1)*(2*n - 3)*(2*n - 1)*(850837923857825*n^9 - 15127768128079400*n^8 + 116088908648008427*n^7 - 502364025369222635*n^6 + 1342887860190877280*n^5 - 2280899268898038065*n^4 + 2433907848768834828*n^3 - 1548429898790214180*n^2 + 521138603722292640*n - 68863424146977600)*a(n-1) - 30*(n-2)*(2*n - 3)*(155302170039375*n^11 - 3227155335853125*n^10 + 29807524885054600*n^9 - 161278340404759950*n^8 + 566950865855228019*n^7 - 1356848300481904461*n^6 + 2250482361655315470*n^5 - 2579665279074165840*n^4 + 1996011605601581864*n^3 - 988803599084885136*n^2 + 280851990522009984*n - 34444332223983360)*a(n-2) + 5*(7288303593953125*n^13 - 187891351713750000*n^12 + 2204843674914291875*n^11 - 15579013461781304250*n^10 + 73867718896175411475*n^9 - 247858726321141236540*n^8 + 604530296941440837821*n^7 - 1082990060568950070282*n^6 + 1421457900098213642392*n^5 - 1345695224728829837040*n^4 + 889319601933492222864*n^3 - 386196670582228097568*n^2 + 97916706472751405568*n - 10797892365692920320)*a(n-3) + 10*(n-3)*(2*n - 7)*(5*n - 18)*(5*n - 14)*(5*n - 12)*(5*n - 11)*(108983978975*n^7 - 1065846642400*n^6 + 4333636082786*n^5 - 9458655747964*n^4 + 11899609166891*n^3 - 8554104592084*n^2 + 3208950812340*n - 473478110640)*a(n-4). - Vaclav Kotesovec, Nov 17 2017
a(n) ~ s * sqrt((r*s*(r*s^3 - 1) - 3) / (7*Pi*(5*r*s*(1 + r*s^3) - 3))) / (2*n^(3/2)*r^n), where r = 0.1385102270697349252376651829944449360743895474888... and s = 1.450646440303399446510765649245639306003224666768... are real roots of the system of equations (1 + r*s^2)*(1 + r^2*s^5) = s, r*s*(2 + 5*r*s^3 + 7*r^2*s^5) = 1. - Vaclav Kotesovec, Nov 22 2017
a(n) = Sum_{k=0..floor(n/2)} binomial(2*n+k+1,k) * binomial(2*n+k+1,n-2*k) / (2*n+k+1). - Seiichi Manyama, Jul 18 2023
EXAMPLE
G.f.: A(x) = 1 + x + 3*x^2 + 13*x^3 + 64*x^4 + 340*x^5 + 1903*x^6 +...
Related expansions:
A(x)^2 = 1 + 2*x + 7*x^2 + 32*x^3 + 163*x^4 + 886*x^5 + 5039*x^6 +...
A(x)^5 = 1 + 5*x + 25*x^2 + 135*x^3 + 765*x^4 + 4481*x^5 + 26920*x^6 +...
A(x)^7 = 1 + 7*x + 42*x^2 + 252*x^3 + 1533*x^4 + 9457*x^5 + 59101*x^6 +...
where A(x) = 1 + x*A(x)^2 + x^2*A(x)^5 + x^3*A(x)^7.
The logarithm of the g.f. A = A(x) equals the series:
log(A(x)) = (1 + x*A^3)*x*A + (1 + 2^2*x*A^3 + x^2*A^6)*x^2*A^2/2 +
(1 + 3^2*x*A^3 + 3^2*x^2*A^6 + x^3*A^9)*x^3*A^3/3 +
(1 + 4^2*x*A^3 + 6^2*x^2*A^6 + 4^2*x^3*A^9 + x^4*A^12)*x^4*A^4/4 +
(1 + 5^2*x*A^3 + 10^2*x^2*A^6 + 10^2*x^3*A^9 + 5^2*x^4*A^12 + x^5*A^15)*x^5*A^5/5 + ...
PROG
(PARI) {a(n)=local(p=1, q=3, A=1+x); for(i=1, n, A=(1+x*A^(p+1))*(1+x^2*A^(p+q+1))+x*O(x^n)); polcoeff(A, n)}
(PARI) {a(n)=local(p=1, q=3, A=1+x); for(i=1, n, A=exp(sum(m=1, n, x^m*(A+x*O(x^n))^(p*m)/m*sum(j=0, m, binomial(m, j)^2*x^j*(A+x*O(x^n))^(q*j))))); polcoeff(A, n, x)}
(PARI) {a(n)=local(p=1, q=3, A=1+x); for(i=1, n, A=exp(sum(m=1, n, x^m*(A+x*O(x^n))^(p*m)/m*(1-x*A^q)^(2*m+1)*sum(j=0, n, binomial(m+j, j)^2*x^j*(A+x*O(x^n))^(q*j))))); polcoeff(A, n, x)}
CROSSREFS
Cf. A200716, A200717, A200718, A200074, A200075, A199876, A199874, A199877, A198951, A198953, A198957, A192415, A198888, A036765, A104545.
G.f. satisfies A(x) = (1 + x*A(x)^2) * (1 + x^2*A(x)^6).
+0
13
1, 1, 3, 14, 75, 433, 2636, 16668, 108399, 720431, 4871555, 33409042, 231817448, 1624503716, 11480658056, 81731416480, 585579734959, 4219179476875, 30552067317233, 222225174139730, 1622894404239115, 11894991079960721, 87472260252499560, 645183802300787356, 4771926560361458884
COMMENTS
More generally, for fixed parameters p and q, if F(x) satisfies:
F(x) = exp( Sum_{n>=1} x^n * F(x)^(n*p)/n * [Sum_{k=0..n} C(n,k)^2 * x^k * F(x)^(k*q)] ),
then F(x) = (1 + x*F(x)^(p+1))*(1 + x^2*F(x)^(p+q+1)).
FORMULA
G.f. A(x) satisfies:
(1) A(x) = sqrt( (1/x)*Series_Reversion( 2*x^5*(1+x)^2/(1 - 2*x^2*(1+x)^2 - sqrt(1 - 4*x^2*(1+x)^2)) ) ).
(2) A(x) = G(x*A(x)^2) where G(x) = A(x/G(x)^2) is the g.f. of A104545 (Motzkin paths of length n having no consecutive (1,0) steps).
(3) A(x) = exp( Sum_{n>=1} x^n * A(x)^n/n * [Sum_{k=0..n} C(n,k)^2 * x^k * A(x)^(4*k)] ).
(4) A(x) = exp( Sum_{n>=1} x^n * A(x)^n/n * [(1-x/A(x)^4)^(2*n+1) * Sum_{k>=0} C(n+k,k)^2*x^k * A(x)^(4*k)] ).
a(n) = Sum_{k=0..floor(n/2)}((binomial(2*n+2*k+1,k)*binomial(2*n+2*k+1,n-2*k))/(2*n+2*k+1)). - Vladimir Kruchinin, Mar 11 2016
EXAMPLE
G.f.: A(x) = 1 + x + 3*x^2 + 14*x^3 + 75*x^4 + 433*x^5 + 2636*x^6 +...
Related expansions:
A(x)^2 = 1 + 2*x + 7*x^2 + 34*x^3 + 187*x^4 + 1100*x^5 + 6784*x^6 +...
A(x)^6 = 1 + 6*x + 33*x^2 + 194*x^3 + 1200*x^4 + 7674*x^5 + 50317*x^6 +...
A(x)^8 = 1 + 8*x + 52*x^2 + 336*x^3 + 2210*x^4 + 14776*x^5 + 100216*x^6 +...
where A(x) = 1 + x*A(x)^2 + x^2*A(x)^6 + x^3*A(x)^8.
The logarithm of the g.f. A = A(x) equals the series:
log(A(x)) = (1 + x*A^4)*x*A + (1 + 2^2*x*A^4 + x^2*A^8)*x^2*A^2/2 +
(1 + 3^2*x*A^4 + 3^2*x^2*A^8 + x^3*A^12)*x^3*A^3/3 +
(1 + 4^2*x*A^4 + 6^2*x^2*A^8 + 4^2*x^3*A^12 + x^4*A^16)*x^4*A^4/4 +
(1 + 5^2*x*A^4 + 10^2*x^2*A^8 + 10^2*x^3*A^12 + 5^2*x^4*A^16 + x^5*A^20)*x^5*A^5/5 + ...
The g.f. of A104545, G(x) = A(x/G(x)^2) where A(x) = G(x*A(x)^2), begins:
G(x) = 1 + x + x^2 + 3*x^3 + 5*x^4 + 11*x^5 + 25*x^6 + 55*x^7 + 129*x^8 +...
MATHEMATICA
a[n_] := Sum[Binomial[2*n + 2*k + 1, k]*Binomial[2*n + 2*k + 1, n - 2*k]/ (2*n + 2*k + 1), {k, 0, n/2}];
PROG
(PARI) {a(n)=polcoeff(sqrt( (1/x)*serreverse( 2*x^5*(1+x)^2/(1 - 2*x^2*(1+x)^2 - sqrt(1 - 4*x^2*(1+x)^2+O(x^(n+6)))) ) ), n)}
(PARI) {a(n)=local(p=1, q=4, A=1+x); for(i=1, n, A=(1+x*A^(p+1))*(1+x^2*A^(p+q+1))+x*O(x^n)); polcoeff(A, n)}
(PARI) {a(n)=local(p=1, q=4, A=1+x); for(i=1, n, A=exp(sum(m=1, n, x^m*(A+x*O(x^n))^(p*m)/m*sum(j=0, m, binomial(m, j)^2*x^j*(A+x*O(x^n))^(q*j))))); polcoeff(A, n, x)}
(PARI) {a(n)=local(p=1, q=4, A=1+x); for(i=1, n, A=exp(sum(m=1, n, x^m*(A+x*O(x^n))^(p*m)/m*(1-x*A^q)^(2*m+1)*sum(j=0, n, binomial(m+j, j)^2*x^j*(A+x*O(x^n))^(q*j))))); polcoeff(A, n, x)}
(Maxima)
a(n):=sum((binomial(2*n+2*k+1, k)*binomial(2*n+2*k+1, n-2*k))/(2*n+2*k+1), k, 0, (n)/2); /* Vladimir Kruchinin, Mar 11 2016 */
CROSSREFS
Cf. A104545, A200716, A200717, A199876, A199877, A198951, A198953, A198957, A192415, A198888, A036765.
G.f. satisfies A(x) = (1 + x^2)*(1 + x*A(x)^2).
+0
10
1, 1, 3, 8, 24, 80, 278, 997, 3670, 13782, 52588, 203314, 794726, 3135540, 12470444, 49942305, 201233170, 815205699, 3318291966, 13565162636, 55669063762, 229257178198, 947142023262, 3924380904498, 16303716754884, 67899954924360, 283425070356740, 1185551594834910
COMMENTS
More generally, for fixed parameters p, q, r, and s, if F(x) satisfies:
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)] ),
then F(x) = (1 + x^r*F(x)^(p+1))*(1 + x^(r+s)*F(x)^(p+q+1)).
FORMULA
G.f. satisfies:
(1) A(x) = (1 - sqrt(1 - 4*x*(1+x^2)^2)) / (2*x*(1+x^2)).
(2) A(x) = exp( Sum_{n>=1} x^n/n * A(x)^n * [Sum_{k=0..n} C(n,k)^2 * x^k / A(x)^(2*k)] ).
(3) A(x) = exp( Sum_{n>=1} (1-x/A(x)^2)^(2*n+1)*[Sum_{k>=0} C(n+k,k)^2*x^k/A(x)^(2*k) )] * x^n*A(x)^n/n ).
(4) A(x) = x / Series_Reversion( x*G(x) ) where G(x) is the g.f. of A200717.
(5) A(x) = G(x/A(x)) where G(x) = A(x*G(x)) is the g.f. of A200717.
Recurrence: (n+1)*a(n) = 2*(2*n-1)*a(n-1) - (n+1)*a(n-2) + 6*(2*n-5)*a(n-3) + 6*(2*n-9)*a(n-5) + 2*(2*n-13)*a(n-7). - Vaclav Kotesovec, Aug 19 2013
a(n) ~ c*d^n/n^(3/2), where d = 4.41997678... is the root of the equation -4-8*d^2-4*d^4+d^5=0 and c = sqrt(d*(8 + 16*d^2 + 8*d^4 + 3*d^5 + d^7) / (Pi*(1 + d^2)^3))/4 = 0.648259186485429075561822659694489853... - Vaclav Kotesovec, Aug 19 2013, updated Oct 11 2018
a(n) = Sum_{i=0..floor(n/2)} C(2*n-4*i+1,i)*C(2*n-4*i+1,n-2*i)/(2*n-4*i+1). - Vladimir Kruchinin, Oct 11 2018
EXAMPLE
G.f.: A(x) = 1 + x + 3*x^2 + 8*x^3 + 24*x^4 + 80*x^5 + 278*x^6 + 997*x^7 +...
Related expansions:
A(x)^2 = 1 + 2*x + 7*x^2 + 22*x^3 + 73*x^4 + 256*x^5 + 924*x^6 + 3414*x^7 +...
where A(x) = 1+x^2 + x*(1+x^2)*A(x)^2.
The logarithm of the g.f. A = A(x) equals the series:
log(A(x)) = (1 + x/A^2)*A*x + (1 + 2^2*x/A^2 + x^2/A^4)*A^2*x^2/2 +
(1 + 3^2*x/A^2 + 3^2*x^2/A^4 + x^3/A^6)*A^3*x^3/3 +
(1 + 4^2*x/A^2 + 6^2*x^2/A^4 + 4^2*x^3/A^6 + x^4/A^8)*A^4*x^4/4 +
(1 + 5^2*x/A^2 + 10^2*x^2/A^4 + 10^2*x^3/A^6 + 5^2*x^4/A^8 + x^5/A^10)*A^5*x^5/5 +
(1 + 6^2*x/A^2 + 15^2*x^2/A^4 + 20^2*x^3/A^6 + 15^2*x^4/A^8 + 6^2*x^5/A^10 + x^6/A^12)*A^6*x^6/6 +...
more explicitly,
log(A(x)) = x + 5*x^2/2 + 16*x^3/3 + 57*x^4/4 + 231*x^5/5 + 938*x^6/6 + 3830*x^7/7 + 15833*x^8/8 +...
MATHEMATICA
nmax=20; aa=ConstantArray[0, nmax]; aa[[1]]=1; Do[AGF=1+Sum[aa[[n]]*x^n, {n, 1, j-1}]+koef*x^j; sol=Solve[Coefficient[(1+x^2)*(1+x*AGF^2)-AGF, x, j]==0, koef][[1]]; aa[[j]]=koef/.sol[[1]], {j, 2, nmax}]; Flatten[{1, aa}] (* Vaclav Kotesovec, Aug 19 2013 *)
CoefficientList[Series[(1 - Sqrt[1 - 4*x*(1 + x^2)^2]) / (2*x*(1 + x^2)), {x, 0, 30}], x] (* Vaclav Kotesovec, Oct 11 2018 *)
Table[Sum[Binomial[2*n - 4*i + 1, i] * Binomial[2*n - 4*i + 1, n - 2*i]/(2*n - 4*i + 1), {i, 0, Floor[n/2]}], {n, 0, 30}] (* Vaclav Kotesovec, Oct 11 2018, after Vladimir Kruchinin *)
PROG
(PARI) {a(n)=polcoeff((1 - sqrt(1 - 4*x*(1+x^2 +x*O(x^n))^2)) / (2*x*(1+x^2 +x*O(x^n))), n)}
for(n=0, 31, print1(a(n), ", "))
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=(1+x*A^2)*(1+x^2)+x*O(x^n)); polcoeff(A, n)}
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2*x^j/A^(2*j))*(x*A+x*O(x^n))^m/m))); polcoeff(A, n, x)}
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, (1-x/A^2)^(2*m+1)*sum(j=0, n, binomial(m+j, j)^2*x^j/A^(2*j))*x^m*A^m/m))); polcoeff(A, n, x)}
G.f. satisfies A(x) = (1 + x*A(x)^2)*(1 + x^2*A(x)^2).
+0
18
1, 1, 3, 10, 37, 147, 611, 2625, 11564, 51953, 237123, 1096420, 5125063, 24178427, 114974387, 550511901, 2651896733, 12843003108, 62494595022, 305400429548, 1498184696271, 7375179807191, 36421312544431, 180383163330765, 895756907248150, 4459095182031675, 22247684478181317
FORMULA
G.f. A(x) satisfies:
(1) A(x) = (1/x)*Series_Reversion( x/(1+x^2) - x^2 ).
(2) A( x*(1-x-x^3)/(1+x^2) ) = (1+x^2)/(1-x-x^3).
(3) a(n) = [x^n] ((1+x^2)/(1-x-x^3))^(n+1) / (n+1).
(4) A(x) = exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^2 * x^k] * x^n*A(x)^n/n ).
(5) A(x) = exp( Sum_{n>=1} [(1-x)^(2*n+1)*Sum_{k>=0} C(n+k,k)^2*x^k )] * x^n*A(x)^n/n ).
Recurrence: 31*(n-1)*n*(n+1)*(85396*n^4 - 902916*n^3 + 3471647*n^2 - 5767203*n + 3503250)*a(n) = 2*(n-1)*n*(6319304*n^5 - 69975436*n^4 + 290875210*n^3 - 559740413*n^2 + 484175751*n - 138985722)*a(n-1) + 2*(n-1)*(2903464*n^6 - 36506072*n^5 + 179801738*n^4 - 439606930*n^3 + 553204983*n^2 - 328951215*n + 67014378)*a(n-2) + 2*(2*n - 5)*(1964108*n^6 - 24695284*n^5 + 123902749*n^4 - 317652203*n^3 + 438313617*n^2 - 307740825*n + 85471038)*a(n-3) - 32*(n-3)*(2*n - 7)*(85396*n^5 - 860218*n^4 + 3249611*n^3 - 5747414*n^2 + 4753791*n - 1471338)*a(n-4) + 8*(n-4)*(n-3)*(2*n - 9)*(85396*n^4 - 561332*n^3 + 1275275*n^2 - 1191073*n + 390174)*a(n-5). - Vaclav Kotesovec, Aug 18 2013
a(n) ~ c*d^n/n^(3/2), where d=5.28245622984... is the root of the equation -16 + 64*d - 92*d^2 - 68*d^3 - 148*d^4 + 31*d^5 = 0 and c = 0.49559010377906722118329... - Vaclav Kotesovec, Aug 18 2013
a(n) = Sum_{k=0..floor(n/2)} binomial(2*n-2*k+1,k) * binomial(2*n-2*k+1,n-2*k) / (2*n-2*k+1). - Seiichi Manyama, Jul 18 2023
EXAMPLE
G.f.: A(x) = 1 + x + 3*x^2 + 10*x^3 + 37*x^4 + 147*x^5 + 611*x^6 +...
where A( x/(1+x^2) - x^2 ) = (1+x^2)/(1-x-x^3).
Related expansions:
A(x)^2 = 1 + 2*x + 7*x^2 + 26*x^3 + 103*x^4 + 428*x^5 + 1838*x^6 +...
A(x)^4 = 1 + 4*x + 18*x^2 + 80*x^3 + 359*x^4 + 1632*x^5 + 7506*x^6 +...
where A(x) = 1 + x*(1+x)*A(x)^2 + x^3*A(x)^4.
The logarithm of the g.f. equals the series:
log(A(x)) = (1 + x)*x*A(x) + (1 + 2^2*x + x^2)*x^2*A(x)^2/2 +
(1 + 3^2*x + 3^2*x^2 + x^3)*x^3*A(x)^3/3 +
(1 + 4^2*x + 6^2*x^2 + 4^2*x^3 + x^4)*x^4*A(x)^4/4 +
(1 + 5^2*x + 10^2*x^2 + 10^2*x^3 + 5^2*x^4 + x^5)*x^5*A(x)^5/5 +
(1 + 6^2*x + 15^2*x^2 + 20^2*x^3 + 15^2*x^4 + 6^2*x^5 + x^6)*x^6*A(x)^6/6 +...
more explicitly,
log(A(x)) = x + 5*x^2/2 + 22*x^3/3 + 101*x^4/4 + 481*x^5/5 + 2330*x^6/6 +...
MATHEMATICA
nmax=20; aa=ConstantArray[0, nmax]; aa[[1]]=1; Do[AGF=1+Sum[aa[[n]]*x^n, {n, 1, j-1}]+koef*x^j; sol=Solve[Coefficient[(1+x*AGF^2)*(1+x^2*AGF^2)-AGF, x, j]==0, koef][[1]]; aa[[j]]=koef/.sol[[1]], {j, 2, nmax}]; Flatten[{1, aa}] (* Vaclav Kotesovec, Aug 18 2013 *)
PROG
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2*x^j)*(x*A+x*O(x^n))^m/m))); polcoeff(A, n, x)}
(PARI) {a(n)=polcoeff((1/x)*serreverse(x/(1+x^2+x*O(x^n))-x^2), n)}
(PARI) {a(n)=polcoeff(((1+x^2)/(1-x-x^3+x*O(x^n)))^(n+1)/(n+1), n)}
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, (1-x)^(2*m+1)*sum(j=0, n, binomial(m+j, j)^2*x^j)*x^m*A^m/m))); polcoeff(A, n, x)}
G.f. satisfies A(x) = (1 + x*A(x)^2)*(1 + x^2*A(x)^3).
+0
12
1, 1, 3, 11, 45, 198, 914, 4367, 21414, 107155, 544987, 2808978, 14640073, 77025373, 408544815, 2182206259, 11727989593, 63373962690, 344109933186, 1876562458845, 10273572074493, 56443282489240, 311097732946200, 1719707775782826, 9531914043637385, 52963938340248863, 294966593345731623
COMMENTS
More generally, for fixed parameters p, q, r, and s, if F(x) satisfies:
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)] ),
then F(x) = (1 + x^r*F(x)^(p+1))*(1 + x^(r+s)*F(x)^(p+q+1)).
FORMULA
G.f.: (1/x)*Series_Reversion( x*(1-x-x^2 + sqrt((1+x+x^2)*(1-3*x+x^2)))/2 ).
a(n) = [x^n] G(x)^(n+1)/(n+1), where 1+x*G(x) is the g.f. of A004148.
G.f. A(x) satisfies:
(1) A(x) = (1/x)*Series_Reversion( x/G(x) ) where 1+x*G(x) is the g.f. of A004148.
(2) A(x) = G(x*A(x)) where G(x) = A(x/G(x)) and 1+x*G(x) is the g.f. of A004148.
(3) A(x) = exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^2 * x^k*A(x)^k] * x^n*A(x)^n/n ).
(4) A(x) = exp( Sum_{n>=1} [(1-x*A(x))^(2*n+1)*Sum_{k>=0} C(n+k,k)^2*x^k*A(x)^k )] * x^n*A(x)^n/n.
Recurrence: 8*n*(2*n+1)*(4*n+1)*(4*n+3)*(1557671*n^7 - 18939961*n^6 + 94817789*n^5 - 252067387*n^4 + 381880748*n^3 - 327052012*n^2 + 145198992*n - 25583040)*a(n) = (2026529971*n^11 - 24640889261*n^10 + 122927623620*n^9 - 322351865586*n^8 + 467303512311*n^7 - 343677276405*n^6 + 61590777290*n^5 + 76066203476*n^4 - 45605627832*n^3 + 4625651136*n^2 + 1916801280*n - 338688000)*a(n-1) + 2*(800642894*n^11 - 10936104295*n^10 + 62803409541*n^9 - 196202081616*n^8 + 357730085364*n^7 - 370711524567*n^6 + 174415015309*n^5 + 25877389846*n^4 - 63266190708*n^3 + 19055552472*n^2 + 1313789760*n - 861840000)*a(n-2) + 6*(308418858*n^11 - 4675368852*n^10 + 30103912361*n^9 - 106665982366*n^8 + 223860428776*n^7 - 274000455628*n^6 + 166116940489*n^5 - 2432493994*n^4 - 54297743044*n^3 + 22033617000*n^2 + 936446400*n - 1315440000)*a(n-3) + 6*(n-2)*(2*n-7)*(3*n-10)*(3*n-8)*(1557671*n^7 - 8036264*n^6 + 13889114*n^5 - 7559372*n^4 - 2491645*n^3 + 2975476*n^2 - 179460*n - 187200)*a(n-4). - Vaclav Kotesovec, Sep 19 2013
a(n) ~ c*d^n/(sqrt(Pi)*n^(3/2)), where d = 1301/1024 + 1/(1024*sqrt(3/(7183147 - (2002819072*2^(2/3))/(3725055779 + 42057117*sqrt(16305))^(1/3) + 1024*(7450111558 + 84114234*sqrt(16305))^(1/3)))) + (1/2)*sqrt(7183147/393216 - (3725055779 + 42057117*sqrt(16305))^(1/3)/(384*2^(2/3)) + 977939/(192*(7450111558 + 84114234*sqrt(16305))^(1/3)) + (1/131072)*(4194454317*sqrt(3/(7183147 - (2002819072*2^(2/3))/(3725055779 + 42057117*sqrt(16305))^(1/3) + 1024*(7450111558 + 84114234*sqrt(16305))^(1/3))))) = 5.89828930084513611... is the root of the equation -108 - 1188*d - 1028*d^2 - 1301*d^3 + 256*d^4 = 0 and c = 0.656947859044624009263362998790812821830934... - Vaclav Kotesovec, Sep 19 2013
a(n) = Sum_{k=0..floor(n/2)} binomial(2*n-k+1,k) * binomial(2*n-k+1,n-2*k) / (2*n-k+1). - Seiichi Manyama, Jul 18 2023
EXAMPLE
G.f.: A(x) = 1 + x + 3*x^2 + 11*x^3 + 45*x^4 + 198*x^5 + 914*x^6 +...
Related expansions:
A(x)^2 = 1 + 2*x + 7*x^2 + 28*x^3 + 121*x^4 + 552*x^5 + 2615*x^6 +...
A(x)^3 = 1 + 3*x + 12*x^2 + 52*x^3 + 237*x^4 + 1122*x^5 + 5463*x^6 +...
A(x)^5 = 1 + 5*x + 25*x^2 + 125*x^3 + 630*x^4 + 3211*x^5 + 16545*x^6 +...
where A(x) = 1 + x*A(x)^2 + x^2*A(x)^3 + x^3*A(x)^5.
The logarithm of the g.f. A = A(x) equals the series:
log(A(x)) = (1 + x*A)*x*A + (1 + 2^2*x*A + x^2*A^2)*x^2*A^2/2 +
(1 + 3^2*x*A + 3^2*x^2*A^2 + x^3*A^3)*x^3*A^3/3 +
(1 + 4^2*x*A + 6^2*x^2*A^2 + 4^2*x^3*A^3 + x^4*A^4)*x^4*A^4/4 +
(1 + 5^2*x*A + 10^2*x^2*A^2 + 10^2*x^3*A^3 + 5^2*x^4*A^4 + x^5*A^5)*x^5*A^5/5 +
(1 + 6^2*x*A + 15^2*x^2*A^2 + 20^2*x^3*A^3 + 15^2*x^4*A^4 + 6^2*x^5*A^5 + x^6*A^6)*x^6*A^6/6 +...
more explicitly,
log(A(x)) = x + 5*x^2/2 + 25*x^3/3 + 129*x^4/4 + 686*x^5/5 + 3713*x^6/6 + 20350*x^7/7 +...
Given G(x) where 1+x*G(x) is the g.f. of A004148, then the coefficients in the powers of G(x) begin:
1: [(1), 1, 2, 4, 8, 17, 37, 82, 185, 423, 978, ...];
2: [1,(2), 5, 12, 28, 66, 156, 370, 882, 2112, ...];
3: [1, 3,(9), 25, 66, 171, 437, 1107, 2790, 7009, ...];
4: [1, 4, 14,(44), 129, 364, 1000, 2696, 7172, 18892, ...];
5: [1, 5, 20, 70,(225), 686, 2015, 5760, 16135, 44500, ...];
6: [1, 6, 27, 104, 363,(1188), 3713, 11214, 32994, 95106, ...];
7: [1, 7, 35, 147, 553, 1932,(6398), 20350, 62734, 188650, ...];
8: [1, 8, 44, 200, 806, 2992, 10460,(34936), 112585, 352560, ...];
9: [1, 9, 54, 264, 1134, 4455, 16389, 57330,(192726), 627406, ...]; ...;
the coefficients in parenthesis form the initial terms of this sequence:
[1/1, 2/2, 9/3, 44/4, 225/5, 1188/6, 6398/7, 34936/8, 192726/9, ...].
The coefficients in the logarithm of the g.f. is also a diagonal in the above table.
MATHEMATICA
CoefficientList[1/x*InverseSeries[Series[x*(1-x-x^2 + Sqrt[(1+x+x^2)*(1-3*x+x^2)])/2, {x, 0, 20}], x], x] (* Vaclav Kotesovec, Sep 19 2013 *)
PROG
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=(1+x*A^2)*(1+x^2*A^3)+x*O(x^n)); polcoeff(A, n)}
(PARI) {a(n)=polcoeff(1/x*serreverse(x*(1-x-x^2 + sqrt((1+x+x^2)*(1-3*x+x^2)+x*O(x^n)))/2), n)}
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2*x^j*A^j)*(x*A+x*O(x^n))^m/m))); polcoeff(A, n, x)}
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, (1-x*A)^(2*m+1)*sum(j=0, n, binomial(m+j, j)^2*x^j*A^j)*x^m*A^m/m))); polcoeff(A, n, x)}
G.f. satisfies A(x) = (1 + x*A(x)^2)*(1 + x^2*A(x)).
+0
12
1, 1, 3, 9, 30, 108, 406, 1577, 6280, 25499, 105169, 439388, 1855636, 7908909, 33975250, 146954693, 639460707, 2797384235, 12295494109, 54272825103, 240480529815, 1069257987503, 4769306203838, 21334400243252, 95687482105807, 430217846136134, 1938651904470374, 8754225470415889
COMMENTS
More generally, for fixed parameters p, q, r, and s, if F(x) satisfies:
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)] ),
then F(x) = (1 + x^r*F(x)^(p+1))*(1 + x^(r+s)*F(x)^(p+q+1)).
FORMULA
G.f. satisfies:
(1) A(x) = exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^2 * x^k*A(x)^(n-k)] * x^n/n ).
(2) A(x) = exp( Sum_{n>=1} [(1-x/A(x))^(2*n+1)*Sum_{k>=0} C(n+k,k)^2*x^k/A(x)^k )] * x^n*A(x)^n/n.
(3) A(x) = x / Series_Reversion( x*G(x) ) where G(x) is the g.f. of A199876.
(4) A(x) = G(x/A(x)) where G(x) = A(x*G(x)) is the g.f. of A199876.
Recurrence: (n+1)*(n+2)*(1241*n^4 - 10636*n^3 + 25417*n^2 - 7382*n - 17136)*a(n) = - 18*(n+1)*(443*n^3 - 3889*n^2 + 9734*n - 5712)*a(n-1) + 4*(6205*n^6 - 53180*n^5 + 115741*n^4 + 64762*n^3 - 370103*n^2 + 246727*n - 25704)*a(n-2) + 6*(2482*n^6 - 24995*n^5 + 76519*n^4 - 36347*n^3 - 185471*n^2 + 293092*n - 140400)*a(n-3) + 2*(4964*n^6 - 57436*n^5 + 228617*n^4 - 276802*n^3 - 361447*n^2 + 956696*n - 320496)*a(n-4) - 6*(2482*n^6 - 32441*n^5 + 140587*n^4 - 173153*n^3 - 266705*n^2 + 677518*n - 291840)*a(n-5) + 12*(n-4)*(2*n - 11)*(11*n^2 + 73*n - 748)*a(n-6) + 2*(n-5)*(2*n - 13)*(1241*n^4 - 5672*n^3 + 955*n^2 + 16508*n - 8496)*a(n-7). - Vaclav Kotesovec, Aug 18 2013
a(n) ~ c*d^n/n^(3/2), where d = 4.770539985405... is the root of the equation -4 + 12*d^2 - 8*d^3 - 12*d^4 - 20*d^5 + d^7 = 0 and c = 0.612892860188927397373456... - Vaclav Kotesovec, Aug 18 2013
a(n) = Sum_{k=0..floor(n/2)} binomial(2*n-3*k+1,k) * binomial(2*n-3*k+1,n-2*k) / (2*n-3*k+1). - Seiichi Manyama, Jul 18 2023
EXAMPLE
G.f.: A(x) = 1 + x + 3*x^2 + 9*x^3 + 30*x^4 + 108*x^5 + 406*x^6 + 1577*x^7 +...
Related expansions:
A(x)^2 = 1 + 2*x + 7*x^2 + 24*x^3 + 87*x^4 + 330*x^5 + 1289*x^6 +...
A(x)^3 = 1 + 3*x + 12*x^2 + 46*x^3 + 180*x^4 + 720*x^5 + 2928*x^6 +...
where A(x) = 1 + x*A(x)^2 + x^2*A(x) + x^3*A(x)^3.
The logarithm of the g.f. A = A(x) equals the series:
log(A(x)) = (A + x)*x + (A^2 + 2^2*x*A + x^2)*x^2/2 +
(A^3 + 3^2*x*A^2 + 3^2*x^2*A + x^3)*x^3/3 +
(A^4 + 4^2*x*A^3 + 6^2*x^2*A^2 + 4^2*x^3*A + x^4)*x^4/4 +
(A^5 + 5^2*x*A^4 + 10^2*x^2*A^3 + 10^2*x^3*A^2 + 5^2*x^4*A + x^5)*x^5/5 +
(A^6 + 6^2*x*A^5 + 15^2*x^2*A^4 + 20^2*x^3*A^3 + 15^2*x^4*A^2 + 6^2*x^5*A + x^6)*x^6/6 +...
more explicitly,
log(A(x)) = x + 5*x^2/2 + 19*x^3/3 + 77*x^4/4 + 331*x^5/5 + 1445*x^6/6 + 6392*x^7/7 + 28565*x^8/8 +...
MAPLE
a:= n-> coeff(series(RootOf(A=(1+x*A^2)*(1+x^2*A), A), x, n+1), x, n):
MATHEMATICA
m = 28; A[_] = 0;
Do[A[x_] = (1 + x A[x]^2)(1 + x^2 A[x]) + O[x]^m, {m}];
PROG
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=(1+x*A^2)*(1+x^2*A^1)+x*O(x^n)); polcoeff(A, n)}
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2*x^j/A^j)*(x*A+x*O(x^n))^m/m))); polcoeff(A, n, x)}
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, (1-x/A)^(2*m+1)*sum(j=0, n, binomial(m+j, j)^2*x^j/A^j)*x^m*A^m/m))); polcoeff(A, n, x)}
G.f. satisfies: A(x) = (1 + x*A(x)^3) * (1 + x^2*A(x)^6).
+0
7
1, 1, 4, 22, 139, 953, 6894, 51796, 400269, 3161262, 25403536, 207043048, 1707345547, 14219399626, 119431172630, 1010495472960, 8604568715969, 73683710894255, 634142349130800, 5482062214763436, 47582484748270453, 414503778412715065, 3622792181209018168, 31758958747482608912
COMMENTS
More generally, for fixed parameters p and q, if F(x) satisfies:
F(x) = exp( Sum_{n>=1} x^n * F(x)^(n*p)/n * [Sum_{k=0..n} C(n,k)^2 * x^k * F(x)^(k*q)] ),
then F(x) = (1 + x*F(x)^(p+1))*(1 + x^2*F(x)^(p+q+1)); here p=2, q=3.
FORMULA
G.f. A(x) satisfies:
(1) a(n) = [x^n] (1 + x + x^2 + x^3)^(3*n+1) / (3*n+1).
(2) A(x) = ( (1/x)*Series_Reversion( x/(1 + x + x^2 + x^3)^3 ) )^(1/3).
(3) A( x/(1 + x + x^2 + x^3)^3 ) = 1 + x + x^2 + x^3.
(4) A(x) = G(x*A(x)^2) where G(x) = A(x/G(x)^2) = g.f. of A036765 (number of rooted trees with a degree constraint).
(5) A(x) = exp( Sum_{n>=1} x^n*A(x)^(2*n)/n * [Sum_{k=0..n} C(n,k)^2 * x^k*A(x)^(3*k)] ).
(6) A(x) = exp( Sum_{n>=1} x^n*A(x)^(2*n)/n * [(1-x*A(x)^2)^(2*n+1)*Sum_{k>=0} C(n+k,k)^2*x^k*A(x)^(3*k) )] ).
a(n) = 1/(3*n + 1)*Sum_{k = 0..floor(n/2)} binomial(3*n + 1,k)*binomial(3*n + 1,n - 2*k).
More generally, the coefficient of x^n in A(x)^r equals r/(3*n + r)*Sum_{k = 0..floor(n/2)} binomial(3*n + r,k)*binomial(3*n + r,n - 2*k) by the Lagrange-Bürmann formula.
O.g.f. A(x) = exp(Sum_{n >= 1} 1/3*b(n)x^n/n), where b(n) = Sum_{k = 0..floor(n/2)} binomial(3*n,k)*binomial(3*n,n - 2*k). Cf. A036765, A186241, A198951. (End)
Recurrence: 128*n*(2*n - 1)*(4*n - 1)*(4*n + 1)*(8*n - 3)*(8*n - 1)*(8*n + 1)*(8*n + 3)*(511073753*n^7 - 4871850365*n^6 + 19478089219*n^5 - 42349790393*n^4 + 54094962928*n^3 - 40605677522*n^2 + 16589611340*n - 2846611200)*a(n) = 3*(3*n - 2)*(3*n - 1)*(3047149994898003*n^13 - 32094344705469618*n^12 + 145743661212727337*n^11 - 373710048777443810*n^10 + 593788894662012231*n^9 - 600683242386376410*n^8 + 377600776651518819*n^7 - 130595257353511374*n^6 + 11334217618972546*n^5 + 8004135084547148*n^4 - 2618300200112616*n^3 + 152383960257264*n^2 + 33025238671680*n - 3264156403200)*a(n-1) - 576*(n-1)*(3*n - 5)*(3*n - 4)*(3*n - 2)*(3*n - 1)*(495741540410*n^10 - 3982082543435*n^9 + 12891395244590*n^8 - 21360691645174*n^7 + 18695904340190*n^6 - 7495052530111*n^5 + 212344193250*n^4 + 656210670544*n^3 - 106487698440*n^2 - 7969373424*n + 1477828800)*a(n-2) + 110592*(n-2)*(n-1)*(3*n - 8)*(3*n - 7)*(3*n - 5)*(3*n - 4)*(3*n - 2)*(3*n - 1)*(511073753*n^7 - 1294334094*n^6 + 979535842*n^5 - 149518418*n^4 - 72732399*n^3 + 16154432*n^2 + 843684*n - 192240)*a(n-3). - Vaclav Kotesovec, Nov 17 2017
a(n) ~ s/(2*sqrt(3*Pi*(4 - 9*r*s^2*(1 + r*s^3)))*n^(3/2)*r^n), where r = 0.1068159753611743655799981945670627355827110854720... and s = 1.345561337338583233012136458010090420775336284226... are real roots of the system of equations (1 + r*s^3)*(1 + r^2*s^6) = s, 3*r*s^2*(1 + 2*r*s^3 + 3*r^2*s^6) = 1. - Vaclav Kotesovec, Nov 22 2017
EXAMPLE
G.f.: A(x) = 1 + x + 4*x^2 + 22*x^3 + 139*x^4 + 953*x^5 + 6894*x^6 +...
where A(x) = (1 + x*A(x)^3)*(1 + x^2*A(x)^6).
Related expansions:
A(x)^3 = 1 + 3*x + 15*x^2 + 91*x^3 + 609*x^4 + 4335*x^5 + 32197*x^6 +...
A(x)^6 = 1 + 6*x + 39*x^2 + 272*x^3 + 1989*x^4 + 15054*x^5 + 116955*x^6 +...
A(x)^9 = 1 + 9*x + 72*x^2 + 570*x^3 + 4545*x^4 + 36639*x^5 + 298662*x^6 +...
where A(x) = 1 + x*A(x)^3 + x^2*A(x)^6 + x^3*A(x)^9.
The logarithm of the g.f. A = A(x) equals the series:
log(A(x)) = (1 + x*A^3)*x*A^2 + (1 + 2^2*x*A^3 + x^2*A^6)*x^2*A^4/2 +
(1 + 3^2*x*A^3 + 3^2*x^2*A^6 + x^3*A^9)*x^3*A^6/3 +
(1 + 4^2*x*A^3 + 6^2*x^2*A^6 + 4^2*x^3*A^9 + x^4*A^12)*x^4*A^8/4 +
(1 + 5^2*x*A^3 + 10^2*x^2*A^6 + 10^2*x^3*A^9 + 5^2*x^4*A^12 + x^5*A^15)*x^5*A^10/5 + ...
which involves squares of binomial coefficients.
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]^3)*(1 + x^2 A[x]^6) + O[x]^(n + 1), x] == 0 /. sol; sol = sol ~Join~ Solve[eq][[1]], {n, 1, nmax}];
sol /. Rule -> Set;
PROG
(PARI) {a(n)=polcoeff( ((1/x)*serreverse(x/(1 + x + x^2 + x^3 +x*O(x^n))^3))^(1/3), n)}
(PARI) {a(n)=polcoeff( (1 + x + x^2 + x^3 +x*O(x^n))^(3*n+1)/(3*n+1), n)}
(PARI) {a(n)=local(p=2, q=3, A=1+x); for(i=1, n, A=(1+x*A^(p+1))*(1+x^2*A^(p+q+1))+x*O(x^n)); polcoeff(A, n)}
(PARI) {a(n)=local(p=2, q=3, A=1+x); for(i=1, n, A=exp(sum(m=1, n, x^m*(A+x*O(x^n))^(p*m)/m*sum(j=0, m, binomial(m, j)^2*x^j*(A+x*O(x^n))^(q*j))))); polcoeff(A, n, x)}
(PARI) {a(n)=local(p=2, q=3, A=1+x); for(i=1, n, A=exp(sum(m=1, n, x^m*(A+x*O(x^n))^(p*m)/m*(1-x*A^q)^(2*m+1)*sum(j=0, n, binomial(m+j, j)^2*x^j*(A+x*O(x^n))^(q*j))))); polcoeff(A, n, x)}
CROSSREFS
Cf. A036765, A200716, A200717, A200718, A200719, A200725, A200074, A200075, A199874, A199876, A199877, A198951, A198953, A198957, A192415, A198888. Cf. A186241.
G.f. satisfies: A(x) = (1 + x*A(x)^4)*(1 + x^2*A(x)^4).
+0
15
1, 1, 5, 31, 222, 1727, 14179, 120930, 1060992, 9514463, 86818391, 803516167, 7524700644, 71169939341, 678877680077, 6523424076116, 63087757216084, 613575943566436, 5997490784042496, 58886692596764215, 580516324380845804, 5743718741275361697, 57017511243375535969
FORMULA
G.f. A(x) satisfies:
(1) A(x) = exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^2 * x^k] * x^n*A(x)^(3*n)/n ).
(2) A(x) = exp( Sum_{n>=1} [(1-x)^(2*n+1)*Sum_{k>=0} C(n+k,k)^2*x^k )] * x^n*A(x)^(3*n)/n.
a(n) ~ s * sqrt((1 + 2*r + 3*r^2*s^4) / (2*Pi*(3 + 3*r + 14*r^2*s^4))) / (2*n^(3/2)*r^n), where r = 0.0940387024218615638441791629908854357421782432118... and s = 1.322930427586092521664829345633697493713415726621... are real roots of the system of equations 1 + r*(1 + r)*s^4 + r^3*s^8 = s, 4*r*s^3*(1 + r + 2*r^2*s^4) = 1. - Vaclav Kotesovec, Nov 22 2017
EXAMPLE
G.f.: A(x) = 1 + x + 5*x^2 + 31*x^3 + 222*x^4 + 1727*x^5 + 14179*x^6 +...
Related expansions:
A(x)^4 = 1 + 4*x + 26*x^2 + 188*x^3 + 1471*x^4 + 12124*x^5 + 103684*x^6 +...
A(x)^8 = 1 + 8*x + 68*x^2 + 584*x^3 + 5122*x^4 + 45792*x^5 + 416196*x^6 +...
where A(x) = 1 + x*(1+x)*A(x)^4 + x^3*A(x)^8.
The logarithm of the g.f. equals the series:
log(A(x)) = (1 + x)*x*A(x)^3 + (1 + 2^2*x + x^2)*x^2*A(x)^6/2 +
(1 + 3^2*x + 3^2*x^2 + x^3)*x^3*A(x)^9/3 +
(1 + 4^2*x + 6^2*x^2 + 4^2*x^3 + x^4)*x^4*A(x)^12/4 +
(1 + 5^2*x + 10^2*x^2 + 10^2*x^3 + 5^2*x^4 + x^5)*x^5*A(x)^15/5 +
(1 + 6^2*x + 15^2*x^2 + 20^2*x^3 + 15^2*x^4 + 6^2*x^5 + x^6)*x^6*A(x)^18/6 +...
more explicitly,
log(A(x)) = x + 9*x^2/2 + 79*x^3/3 + 733*x^4/4 + 7006*x^5/5 + 68229*x^6/6 + 673268*x^7/7 +...
MATHEMATICA
terms = 23; A[_] = 1; Do[A[x_] = (1 + x*A[x]^4)*(1 + x^2*A[x]^4) + O[x]^terms // Normal, terms]; CoefficientList[A[x], x] (* Jean-François Alcover, Jan 09 2018 *)
PROG
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=(1 + x*A^4)*(1 + x^2*A^4)+x*O(x^n)); polcoeff(A, n)}
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2*x^j)*(x*A^3+x*O(x^n))^m/m))); polcoeff(A, n, x)}
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, (1-x)^(2*m+1)*sum(j=0, n, binomial(m+j, j)^2*x^j)*x^m*A^(3*m)/m))); polcoeff(A, n, x)}
G.f. satisfies: A(x) = (1 + x*A(x)^3) * (1 + x^2*A(x)^4).
+0
1
1, 1, 4, 20, 114, 703, 4565, 30752, 212921, 1505916, 10833164, 79018804, 583062388, 4344431508, 32641910199, 247033970128, 1881402836376, 14408753414558, 110897147057354, 857307054338476, 6653979156676983, 51831065993122915, 405060413133136902, 3175019470333290488
COMMENTS
More generally, for fixed parameters p and q, if F(x) satisfies:
F(x) = exp( Sum_{n>=1} x^n * F(x)^(n*p)/n * [Sum_{k=0..n} C(n,k)^2 * x^k * F(x)^(k*q)] ),
then F(x) = (1 + x*F(x)^(p+1))*(1 + x^2*F(x)^(p+q+1)); here p=2 and q=1.
FORMULA
G.f. A(x) satisfies:
(1) A(x) = sqrt( (1/x)*Series_Reversion( x*(1-x-x^3)^2/(1+x^2)^2 ) ).
(2) A( x*(1-x-x^3)^2/(1+x^2)^2 ) = (1+x^2)/(1-x-x^3).
(3) a(n) = [x^n] ((1+x^2)/(1-x-x^3))^(2*n+2) / (n+1).
(4) A(x) = exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^2 * x^k*A(x)^k] * x^n*A(x)^(2*n)/n ).
(5) A(x) = exp( Sum_{n>=1} [(1-x*A(x))^(2*n+1)*Sum_{k>=0} C(n+k,k)^2*x^k*A(x)^k )] * x^n*A(x)^(2*n)/n ).
(6) A(x) = (1/x)*Series_Reversion(x/G(x)) where A(x) = G(x*A(x)) and A(x/G(x)) = G(x) = (1 + x*G(x)^2)*(1 + x^2*G(x)^2) is the g.f. of A199874.
a(n) ~ s * sqrt((1 + 2*r*s + 3*r^2*s^4) / (3*Pi*(1 + 2*r*s + 7*r^2*s^4))) / (2*n^(3/2)*r^n), where r = 0.1194948955213353102456218138370139612914667337222... and s = 1.428770161302757679335810379290625953730830139744... are real roots of the system of equations (1 + r*s^3)*(1 + r^2*s^4) = s, r*s^2*(3 + 4*r*s + 7*r^2*s^4) = 1. - Vaclav Kotesovec, Nov 22 2017
EXAMPLE
G.f.: A(x) = 1 + x + 4*x^2 + 20*x^3 + 114*x^4 + 703*x^5 + 4565*x^6 +...
where A( x*(1-x-x^3)^2/(1+x^2)^2 ) = (1+x^2)/(1-x-x^3).
Related expansions:
A(x)^3 = 1 + 3*x + 15*x^2 + 85*x^3 + 522*x^4 + 3381*x^5 + 22735*x^6 +...
A(x)^4 = 1 + 4*x + 22*x^2 + 132*x^3 + 841*x^4 + 5588*x^5 + 38288*x^6 +...
A(x)^7 = 1 + 7*x + 49*x^2 + 343*x^3 + 2429*x^4 + 17430*x^5 +...
where A(x) = 1 + x*A(x)^3 + x^2*A(x)^4 + x^3*A(x)^7.
The logarithm of the g.f. equals the series:
log(A(x)) = (1 + x*A(x))*x*A(x)^2 + (1 + 2^2*x*A(x) + x^2*A(x)^2)*x^2*A(x)^4/2 +
(1 + 3^2*x*A(x) + 3^2*x^2*A(x)^2 + x^3*A(x)^3)*x^3*A(x)^6/3 +
(1 + 4^2*x*A(x) + 6^2*x^2*A(x)^2 + 4^2*x^3*A(x)^3 + x^4*A(x)^4)*x^4*A(x)^8/4 +
(1 + 5^2*x*A(x) + 10^2*x^2*A(x)^2 + 10^2*x^3*A(x)^3 + 5^2*x^4*A(x)^4 + x^5*A(x)^5)*x^5*A(x)^10/5 +
(1 + 6^2*x*A(x) + 15^2*x^2*A(x)^2 + 20^2*x^3*A(x)^3 + 15^2*x^4*A(x)^4 + 6^2*x^5*A(x)^5 + x^6*A(x)^6)*x^6*A(x)^12/6 +...
more explicitly,
log(A(x)) = x + 7*x^2/2 + 49*x^3/3 + 359*x^4/4 + 2706*x^5/5 + 20767*x^6/6 +...
MATHEMATICA
CoefficientList[Sqrt[1/x * InverseSeries[Series[x*(1 - x - x^3)^2/(1 + x^2)^2, {x, 0, 20}], x]], x] (* Vaclav Kotesovec, Nov 22 2017 *)
PROG
(PARI) {a(n)=polcoeff(sqrt( (1/x)*serreverse( x*(1-x-x^3)^2/(1+x^2+x*O(x^n))^2 ) ), n)}
for(n=0, 30, print1(a(n), ", "))
(PARI) {a(n)=local(p=2, q=1, A=1+x); for(i=1, n, A=(1+x*A^(p+1))*(1+x^2*A^(p+q+1))+x*O(x^n)); polcoeff(A, n)}
(PARI) {a(n)=local(p=2, q=1, A=1+x); for(i=1, n, A=exp(sum(m=1, n, x^m*(A+x*O(x^n))^(p*m)/m*sum(j=0, m, binomial(m, j)^2*x^j*(A+x*O(x^n))^(q*j))))); polcoeff(A, n, x)}
(PARI) {a(n)=local(p=2, q=1, A=1+x); for(i=1, n, A=exp(sum(m=1, n, x^m*(A+x*O(x^n))^(p*m)/m*(1-x*A^q)^(2*m+1)*sum(j=0, n, binomial(m+j, j)^2*x^j*(A+x*O(x^n))^(q*j))))); polcoeff(A, n, x)}
Search completed in 0.013 seconds
|