|
|
A025550
|
|
a(n) = ( 1/1 + 1/3 + 1/5 + ... + 1/(2*n-1) )*LCM(1, 3, 5, ..., 2*n-1).
|
|
12
|
|
|
1, 4, 23, 176, 563, 6508, 88069, 91072, 1593269, 31037876, 31730711, 744355888, 3788707301, 11552032628, 340028535787, 10686452707072, 10823198495797, 10952130239452, 409741429887649, 414022624965424, 17141894231615609, 743947082888833412, 750488463554668427
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Or, numerator of 1/1 + 1/3 + ... + 1/(2n-1) up to a(38).
Following similar remark by T. D. Noe in A025547, this coincides with f(n) = numerator of 1 + 1/3 + 1/5 + 1/7 + ... + 1/(2n-1) iff n <= 38. But a(39) = 18048708369314455836683437302413, f(39) = 1640791669937677803334857936583. Note that f(n) = numerator(digamma(n+1/2)/2 + log(2) + euler_gamma/2). - Paul Barry, Aug 19 2005 [See A350669(n-1).]
|
|
LINKS
|
|
|
MAPLE
|
a:= n-> (f-> add(1/p, p=f)*ilcm(f[]))([2*i-1$i=1..n]):
|
|
MATHEMATICA
|
Table[(Total[1/Range[1, 2n-1, 2]])LCM@@Range[1, 2n-1, 2], {n, 30}] (* Harvey P. Dale, Sep 09 2020 *)
|
|
PROG
|
(Haskell)
a025550 n = numerator $ sum $ map (1 %) $ take n [1, 3 ..]
(Magma) [&+[1/d: d in i]*Lcm(i) where i is [1..2*n-1 by 2]: n in [1..21]]; // Bruno Berselli, Apr 16 2015
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy,nice,frac
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|