login
Least common multiple of {1,3,5,...,2n-1}.
35

%I #51 May 15 2024 11:04:18

%S 1,3,15,105,315,3465,45045,45045,765765,14549535,14549535,334639305,

%T 1673196525,5019589575,145568097675,4512611027925,4512611027925,

%U 4512611027925,166966608033225,166966608033225,6845630929362225,294362129962575675,294362129962575675

%N Least common multiple of {1,3,5,...,2n-1}.

%C This sequence coincides with the sequence f(n) = denominator of 1 + 1/3 + 1/5 + 1/7 + ... + 1/(2n-1) iff n <= 38. But a(39) = 6414924694381721303722858446525, f(39) = 583174972216520118520259858775. - _T. D. Noe_, Aug 04 2004 [See A350670(n-1).]

%C Coincides for n=1..42 with the denominators of a series for Pi*sqrt(2)/4 and then starts to differ. See A127676.

%C a(floor((n+1)/2)) = gcd(a(n), A051426(n)). - _Reinhard Zumkeller_, Apr 25 2011

%C A051417(n) = a(n+1)/a(n).

%H T. D. Noe, <a href="/A025547/b025547.txt">Table of n, a(n) for n = 1..200</a>

%H Yue-Wu Li and Feng Qi, <a href="https://doi.org/10.3390/axioms13050317">A New Closed-Form Formula of the Gauss Hypergeometric Function at Specific Arguments</a>, Axioms (2024) Vol. 13, Art. No. 317. See p. 11 of 24.

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/JeepProblem.html">Jeep Problem</a>, <a href="http://mathworld.wolfram.com/Pi.html">Pi</a>, <a href="http://mathworld.wolfram.com/PiContinuedFraction.html">Pi Continued Fraction</a>, <a href="http://mathworld.wolfram.com/LeastCommonMultiple.html">Least Common Multiple</a>

%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Least_common_multiple">Least common multiple</a>

%H <a href="/index/Lc#lcm">Index entries for sequences related to lcm's</a>

%p A025547:=proc(n) local i,t1; t1:=1; for i from 1 to n do t1:=lcm(t1,2*i-1); od: t1; end;

%p f := n->denom(add(1/(2*k-1),k=0..n)); # a different sequence!

%t a = 1; Join[{1}, Table[a = LCM[a, n], {n, 3, 125, 2}]] (* _Zak Seidov_, Jan 18 2011 *)

%t nn=30;With[{c=Range[1,2*nn,2]},Table[LCM@@Take[c,n],{n,nn}]] (* _Harvey P. Dale_, Jan 27 2013 *)

%o (Haskell)

%o a025547 n = a025547_list !! (n-1)

%o a025547_list = scanl1 lcm a005408_list

%o -- _Reinhard Zumkeller_, Oct 25 2013, Apr 25 2011

%o (PARI) a(n)=lcm(vector(n,k,2*k-1)) \\ _Charles R Greathouse IV_, Nov 20 2012

%o (Python) # generates initial segment of sequence

%o from math import gcd

%o from itertools import accumulate

%o def lcm(a, b): return a * b // gcd(a, b)

%o def aupton(nn): return list(accumulate((2*i+1 for i in range(nn)), lcm))

%o print(aupton(23)) # _Michael S. Branicky_, Mar 28 2022

%Y Cf. A007509, A025550, A075135. The numerators are in A074599.

%Y Cf. A003418 (LCM of {1..n}).

%Y Cf. A005408, A350670.

%K easy,nice,nonn

%O 1,2

%A _Clark Kimberling_