cp's OEIS Frontend

A355775 Expansion of 1 + f/(1 + 2*f), where f is the g.f. of the swinging factorial (A056040).

%I A355775 #13 Apr 17 2025 07:04:02
%S A355775 1,1,0,2,-10,30,-84,244,-722,2126,-6252,18364,-53988,158668,-466408,
%T A355775 1370792,-4029154,11842222,-34806972,102303468,-300690348,883782404,
%U A355775 -2597604952,7634834648,-22440207764,65955900268,-193856647736,569780485080,-1674690451976
%N A355775 Expansion of 1 + f/(1 + 2*f), where f is the g.f. of the swinging factorial (A056040).
%F A355775 a(n) = [z^n] (1 + ((1 - 4*z^2)^(3/2) + 4*z^2 - z - 1) / ((1 - 4*z^2)^(3/2) + 2*(4*z^2 - z - 1))).
%F A355775 Let u = 1 - 4*z^2, v = u + z and w = u^(3/2) - v. Then the g.f. can be written as 1 + w/(w - v) and the g.f. of the swinging factorial 1 - w/(w + v).
%p A355775 # Traditional style:
%p A355775 u := 1 - 4*z^2: v := u + z: w := u^(3/2) - v: 1 + w/(w - v):
%p A355775 ser := series(%, z, 30): seq(coeff(%, z, n), n = 0..28);
%p A355775 # Alternative (_Alois P. Heinz_ style, see A355488):
%p A355775 c := n -> n! / iquo(n, 2)!^2:
%p A355775 a := n -> (f -> coeff(series(1 + f/(1 + 2*f), x, n+1), x, n))(add(c(j)*x^j, j=1..n)): seq(a(n), n = 0..28);
%t A355775 c[n_] := n! / Quotient[n, 2]!^2;
%t A355775 a[n_] := With[{f = Sum[c[j]*x^j, {j, 1, n}]}, SeriesCoefficient[1 + f/(1 + 2*f), {x, 0, n}]];
%t A355775 Table[a[n], {n, 0, 28}] (* _Jean-François Alcover_, Apr 17 2025, after _Alois P. Heinz_ *)
%o A355775 (SageMath)
%o A355775 def sw_factorial(n): return factorial(n) // factorial(n // 2)^2
%o A355775 A = QQ[['t']]
%o A355775 f = A([0] + [sw_factorial(n) for n in range(1, 29)]).O(29)
%o A355775 print(list(1 + f / (1 + 2 * f)))  # After _F. Chapoton_ in A355488.
%Y A355775 Cf. A056040, A355488.
%K A355775 sign
%O A355775 0,4
%A A355775 _Peter Luschny_, Jul 22 2022