cp's OEIS Frontend

A374565 Expansion of g.f. A(x) satisfying A(x)^3 = A( x*A(x)^2/(1-x) ).

%I A374565 #21 Jul 23 2024 17:09:44
%S A374565 1,1,1,2,5,11,24,57,141,350,881,2267,5920,15601,41497,111399,301293,
%T A374565 819843,2243058,6167211,17029473,47200752,131270283,366195789,
%U A374565 1024380648,2872770381,8074967031,22745832254,64196912681,181516532273,514107418321,1458407886019,4143318012685
%N A374565 Expansion of g.f. A(x) satisfying A(x)^3 = A( x*A(x)^2/(1-x) ).
%H A374565 Paul D. Hanna, <a href="/A374565/b374565.txt">Table of n, a(n) for n = 1..801</a>
%F A374565 G.f. A(x) = Sum_{n>=1} a(n)*x^n satisfies the following formulas.
%F A374565 (1) A(x) = x + x*Sum_{n>=0} A(x)^(3^n).
%F A374565 (2) A(x) = Series_Reversion( x/(1 + Sum_{n>=0} x^(3^n)) ).
%F A374565 (3) A(x)^3 = A( x*A(x)^2/(1-x) ).
%F A374565 (4) A(x)^9 = A( x*A(x)^8/(1 - x - x*A(x)^2) ).
%F A374565 (5) A(x)^27 = A( x*A(x)^26/(1 - x - x*A(x)^2 - x*A(x)^8) ).
%F A374565 (6) A(x)^(3^n) = A( x*A(x)^(3^n-1) / (1 - x*Sum_{k=0..n-1} A(x)^(3^k-1)) ) for n >= 1.
%F A374565 The radius of convergence r and A(r) satisfy r = 1/(Sum_{n>=0} 3^n*A(r)^(3^n-1)) and A(r) = A( A(r)^2*r/(1-r) )^(1/3), where r = 0.3359879296886914478616860912190963818298151003686099... and A(r) = 0.6985186992950193189255500784091315877737446624401085...
%e A374565 G.f.: A(x) = x + x^2 + x^3 + 2*x^4 + 5*x^5 + 11*x^6 + 24*x^7 + 57*x^8 + 141*x^9 + 350*x^10 + 881*x^11 + 2267*x^12 + 5920*x^13 + 15601*x^14 + 41497*x^15 + ...
%e A374565 where A(x)^3 = A( x*A(x)^2/(1-x) )
%e A374565 and A(x) = x + x*(A(x) + A(x)^3 + A(x)^9 + A(x)^27 + ... A(x)^(3^n) + ...).
%e A374565 RELATED SERIES.
%e A374565 A(x)^3 = x^3 + 3*x^4 + 6*x^5 + 13*x^6 + 33*x^7 + 84*x^8 + 208*x^9 + 522*x^10 + 1341*x^11 + 3479*x^12 + 9078*x^13 + 23907*x^14 + 63560*x^15 + ...
%e A374565 Let B(x) be the series reversion of A(x), B(A(x)) = x, then
%e A374565 B(x) = x - x^2 + x^3 - 2*x^4 + 3*x^5 - 4*x^6 + 6*x^7 - 9*x^8 + 13*x^9 - 20*x^10 + ... + (-1)^(n-1)*A078932(n-1)*x^n + ...
%e A374565 where x/B(x) = 1 + x + x^3 + x^9 + x^27 + x^81 + ... + x^(3^n) + ...
%e A374565 F(x) = A(x/(1+x)) = x + x^4 + 3*x^7 + 13*x^10 + 67*x^13 + 378*x^16 + 2253*x^19 + 13947*x^22 + 88803*x^25 + 577903*x^28 + 3826870*x^31 + 25703868*x^34 + ...
%e A374565 where F(x)^3 = F( x*F(x)^2/(1 - x*F(x)^2) )
%e A374565 and F(x) = x + x*(F(x)^3 + F(x)^9 + F(x)^27 + ... + F(x)^(3^n) + ...).
%e A374565 SPECIFIC VALUES.
%e A374565 A(t) = 2/3 at t = 0.3351780091733165997365854281871805851976265481916...
%e A374565 where 8/27 = A( (4/9)*t/(1-t) )
%e A374565 and t = (2/3)/(1 + Sum_{n>=0} (2/3)^(3^n)).
%e A374565 A(t) = 1/2 at t = 0.3073229277642929985518391822746766756418592443672...
%e A374565 where 1/8 = A( (1/4)*t/(1-t) )
%e A374565 and t = (1/2)/(1 + Sum_{n>=0} (1/2)^(3^n)).
%e A374565 A(1/3) = 0.640317989282342396539425948311398871030928082061168...
%e A374565 where A(1/3)^3 = A( A(1/3)^2/2 ).
%e A374565 A(1/4) = 0.347324237093006237340030053166266719890703533474663...
%e A374565 where A(1/4)^3 = A( A(1/4)^2/3 ).
%e A374565 A(1/5) = 0.254102848699628177600720471035831153854183353627930...
%e A374565 where A(1/5)^3 = A( A(1/5)^2/4 ).
%e A374565 A(1/10) = 0.111264157881789221767410282888976753122883279205707...
%e A374565 where A(1/10)^3 = A( A(1/10)^2/9 ).
%o A374565 (PARI) {a(n) = my(A = serreverse(x/(1 + sum(n=0,ceil(log(n+1)/log(3)), x^(3^n)) + x^3*O(x^n)) )); polcoeff(A,n)}
%o A374565 for(n=1, 40, print1(a(n), ", "))
%o A374565 (PARI) {a(n) = my(A=[1], Ax);
%o A374565 for(i=1, n, A=concat(A, 0); Ax=x*Ser(A);
%o A374565 A[#A] = -polcoeff( Ax^3 - subst(Ax, x, Ax^2*x/(1-x) ), #A+2) ); A[n]}
%o A374565 for(n=1, 40, print1(a(n), ", "))
%Y A374565 Cf. A078932, A075864.
%K A374565 nonn
%O A374565 1,4
%A A374565 _Paul D. Hanna_, Jul 23 2024