cp's OEIS Frontend

A352701 G.f. (1/x)*Series_Reversion( x*(1-x)*(3 - 2*x + sqrt(1+4*x-4*x^2))/4 ).

%I A352701 #5 Mar 31 2022 09:36:28
%S A352701 1,1,3,7,28,79,350,1075,5020,16180,78023,259417,1278340,4343642,
%T A352701 21740636,75065787,380161308,1328887420,6792111260,23975385148,
%U A352701 123448657904,439228736887,2275311657814,8148868193557,42427160829508
%N A352701 G.f. (1/x)*Series_Reversion( x*(1-x)*(3 - 2*x + sqrt(1+4*x-4*x^2))/4 ).
%C A352701 Essentially an unsigned version of A351771 (after dropping the initial term).
%C A352701 a(2*n) = A352383(n) for n >= 0.
%F A352701 The g.f. A(x) satisfies:
%F A352701 (1) A(x) = (1/x)*Series_Reversion( x*(1-x)*(3 - 2*x + sqrt(1+4*x-4*x^2))/4 );
%F A352701 (2) (A(x) - A(-x))/2 = x*(A(x)^2 + A(-x)^2)/2;
%F A352701 (3) ((A(x) + A(-x))/2)^3 = F(x^2), where F(x) = (1/x)*Series_Reversion( x*(1+x)^3/(1+2*x)^6 );
%F A352701 (4) 1 - x*(A(x) - A(-x))/2  = x/Series_Reversion( x - x*(C(x) + C(-x))/2 ), where C(x) = (1 - sqrt(1-4*x))/2 is the Catalan function (A000108);
%F A352701 (5a) (1/A(x) + 1/A(-x))/2 = ( 1 - x*(A(x) - A(-x))/2 )^2;
%F A352701 (5b) (1/A(x) - 1/A(-x))/2 = (-x)/(1 - 2*x*(A(x) - A(-x))/2);
%F A352701 (6a) (1/A(x)^2 + 1/A(-x)^2)/2 = ( 1 - x*(A(x) - A(-x))/2 )^3.
%F A352701 (6b) (1/A(x)^2 - 1/A(-x)^2)/2 = -2*x*(1 - x*(A(x) - A(-x))/2)^2/( 1 - x*(A(x) - A(-x)) ).
%F A352701 Let B(x) =  Series_Reversion( x*(1-x^2)/(1+x^2)^3 ), then
%F A352701 (7) A(x) = (1 - sqrt(1 - 4*x - 4*x*B(x)^2))/(2*x);
%F A352701 (8) A(x) - x*A(x)^2 = A(-x) + x*A(-x)^2 = 1 + B(x)^2;
%F A352701 (9) 1 - x*(A(x) - A(-x))/2  = 1/(1 + B(x)^2);
%F A352701 (10) 1/A(x) = 1/(1 + B(x)^2)^2 - x*(1 + B(x)^2)/(1 - B(x)^2);
%F A352701 (10a) (1/A(x) + 1/A(-x))/2 = 1/(1 + B(x)^2)^2;
%F A352701 (10b) (1/A(x) - 1/A(-x))/2 = (-x)*(1 + B(x)^2)/(1 - B(x)^2);
%F A352701 (11) 1/A(x)^2 = 1/(1 + B(x)^2)^3 - 2*x/(1 - B(x)^4);
%F A352701 (11a) (1/A(x)^2 + 1/A(-x)^2)/2 = 1/(1 + B(x)^2)^3;
%F A352701 (11b) (1/A(x)^2 - 1/A(-x)^2)/2 = -2*x/(1 - B(x)^4).
%e A352701 G.f.: A(x) = 1 + x + 3*x^2 + 7*x^3 + 28*x^4 + 79*x^5 + 350*x^6 + 1075*x^7 + 5020*x^8 + 16180*x^9 + 78023*x^10 + 259417*x^11 + ...
%e A352701 such that A(x) = (1/x)*Series_Reversion(x*G(x)) and A(x*G(x)) = 1/G(x),
%e A352701 where G(x) = (1-x)*(3 - 2*x + sqrt(1+4*x-4*x^2))/4, which starts
%e A352701 G(x) = 1 - x - x^2 + 3*x^3 - 8*x^4 + 26*x^5 - 92*x^6 + 344*x^7 - 1336*x^8 + 5336*x^9 - 21776*x^10 + ...
%e A352701 Let B(x) =  Series_Reversion( x*(1-x^2)/(1+x^2)^3 ),
%e A352701 B(x) = x + 4*x^3 + 39*x^5 + 496*x^7 + 7180*x^9 + 112236*x^11 + 1846082*x^13 + 31485120*x^15 + ...,
%e A352701 then A(x) = 1 + x*A(x)^2 + B(x)^2, where
%e A352701 B(x)^2 = x^2 + 8*x^4 + 94*x^6 + 1304*x^8 + 19849*x^10 + 320600*x^12 + 5396108*x^14 + 93615864*x^16 + ...
%o A352701 (PARI) {a(n) = my(A = (1/x)*serreverse( x*(1-x)*(3 - 2*x + sqrt(1+4*x-4*x^2 +x*O(x^n) ))/4 ));
%o A352701 polcoeff(A,n)}
%o A352701 for(n=0,30,print1(a(n),", "))
%Y A352701 Cf. A351771, A352383.
%K A352701 nonn
%O A352701 0,3
%A A352701 _Paul D. Hanna_, Mar 29 2022