OFFSET
1,2
FORMULA
a(n) = numerator(Sum_{j=1..n} 1/lcm(1..n)).
EXAMPLE
n = 7: LCM values: 1, 2, 6, 12, 60, 60, 420; partial sum = 1 + 1/2 + 1/6 + 1/12 + 1/60 + 1/60 + 1/420 = (420 + 210 + 70 + 35 + 7 + 7 + 1)/420 = 750/420 = 25/14, so a(7) = 25.
MATHEMATICA
q[x_] := Apply[LCM, Table[j, {j, 1, x}]] Table[Numerator[Apply[Plus, Table[1/q[w], {w, 1, m}]]], {m, 1, 30}]
Accumulate[Table[1/LCM@@Range[n], {n, 30}]]//Numerator (* Harvey P. Dale, Aug 08 2021 *)
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Labos Elemer, Oct 08 2001
STATUS
approved