cp's OEIS Frontend

A088491 a(n) = floor(p(n)/p(n-1)), where p(n) = n!/(Product_{j=1..floor(n/2)} A004001(j)).

%I A088491 #13 Mar 28 2022 07:45:34
%S A088491 2,3,4,5,3,7,4,9,3,11,3,13,3,15,4,17,3,19,3,21,3,23,3,25,3,27,3,29,3,
%T A088491 31,4,33,3,35,3,37,3,39,3,41,3,43,3,45,3,47,3,49,3,51,3,53,3,55,3,57,
%U A088491 3,59,3,61,3,63,4,65,3,67,3,69,3,71,3,73,3,75,3,77,3,79,3,81,3,83,3,85,3,87
%N A088491 a(n) = floor(p(n)/p(n-1)), where p(n) = n!/(Product_{j=1..floor(n/2)} A004001(j)).
%H A088491 G. C. Greubel, <a href="/A088491/b088491.txt">Table of n, a(n) for n = 2..5000</a>
%F A088491 a(n) = floor(p(n)/p(n-1)), where p(n) = n!/(Product_{j=1..floor(n/2)} A004001(j)).
%t A088491 Conway[n_]:= Conway[n]= If[n<3, 1, Conway[Conway[n-1]] +Conway[n-Conway[n-1]]];
%t A088491 f[n_]:= f[n]= Product[Conway[i], {i,Floor[n/2]}];
%t A088491 a[n_]:= a[n]= Floor[n*f[n-1]/f[n]];
%t A088491 Table[a[n], {n, 2, 100}] (* modified by _G. C. Greubel_, Mar 27 2022 *)
%o A088491 (Sage)
%o A088491 @CachedFunction
%o A088491 def b(n): # A004001
%o A088491     if (n<3): return 1
%o A088491     else: return b(b(n-1)) + b(n-b(n-1))
%o A088491 def f(n): return product( b(j) for j in (1..(n//2)) )
%o A088491 def A088491(n): return (n*f(n-1)//f(n))
%o A088491 [A088491(n) for n in (2..100)] # _G. C. Greubel_, Mar 27 2022
%Y A088491 Cf. A004001, A088493.
%K A088491 nonn
%O A088491 2,1
%A A088491 _Roger L. Bagula_, Nov 10 2003
%E A088491 Edited by _G. C. Greubel_, Mar 27 2022