|
%I #22 May 08 2020 16:58:54
%S 1,1,1,1,2,2,3,3,12,4,10,10,30,30,105,7,56,56,252,252,1260,60,330,330,
%T 1980,396,2574,286,2002,2002,15015,15015,240240,7280,61880,1768,15912,
%U 15912,151164,3876,38760,38760,406980,406980,4476780,99484,1144066
%N a(n) = lcm{1,2,...,n} / swinging_factorial(n) = A003418(n) / A056040(n).
%C Characterization: Let e_{p}(m) denote the exponent of the prime p in the prime factorization of m and [.] denote the Iverson bracket, then
%C e_{p}(a(n)) = Sum_{k>=1} [floor(n/p^k) is even].
%C This implies, among other things, that no prime > floor(n/2) can divide a(n). The prime exponents e_{2}(a(2n)) give Guy Steele's sequence GS(5,3) A080100.
%C Asymptotics: log a(n) ~ n(1 - log 2). It is conjectured that log a(n) ~ n(1 - log 2) + O(n^{1/2+eps}) for all eps > 0.
%C Bounds: A056040(floor(n/3)) <= a(n) <= A056040(floor(n/2)) if n >= 285.
%H Peter Luschny, <a href="/A180000/b180000.txt">Table of n, a(n) for n = 0..1000</a>
%H Peter Luschny, <a href="/A180000/a180000.pdf">Die schwingende Fakultät und Orbitalsysteme</a>, August 2011.
%F a(n) = 2^(-n)*Product_{1<=k<=n} A014963(k)*(k/2)^((-1)^k).
%p a := proc(n) local A014963, k;
%p A014963 := proc(n) if n < 2 then 1 else numtheory[factorset](n);
%p if 1 < nops(%) then 1 else op(%) fi fi end;
%p mul(A014963(k)*(k/2)^((-1)^k), k=1..n)/2^n end;
%p # Also:
%p A180000 := proc(n) local lcm, sf;
%p lcm := ilcm(seq(i,i=1..n));
%p sf := n!/iquo(n,2)!^2;
%p lcm/sf end;
%t a[0] = 1; a[n_] := LCM @@ Range[n] / (n! / Floor[n/2]!^2); Table[a[n], {n, 0, 46}] (* _Jean-François Alcover_, Jul 23 2013 *)
%o (PARI) L=1; X(n)={ ispower(n, , &n);if(isprime(n),n,1); }
%o Y(n)={ a=X(n); b=if(bitand(1,n),a,a*(n/2)^2); L=(b*L)/n; }
%o A180000_list(n)={ L=1; vector(n,m,Y(m)); } \\ for n>0
%o (Sage)
%o def Exp(m,n) :
%o s = 0; p = m; q = n//p
%o while q > 0 :
%o if is_even(q) :
%o s = s + 1
%o p = p * m
%o q = n//p
%o return s
%o def A180000(n) :
%o A = [1,1,1,1,2,2,3,3,12]
%o if n < 9 : return A[n]
%o R = []; r = isqrt(n)
%o P = Primes(); p = P.first()
%o while p <= n//2 :
%o if p <= r : R.append(p^Exp(p,n))
%o elif p <= n//3 :
%o if is_even(n//p) : R.append(p)
%o else : R.append(p)
%o p = P.next(p)
%o return mul(x for x in R)
%Y Cf. A003418, A014963, A056040.
%K nonn
%O 0,5
%A _Peter Luschny_, Aug 17 2010
|
