OFFSET
1,1
COMMENTS
Denominator of a second order prime harmonic number.
FORMULA
a(n) is the denominator of Sum_{j=1..n} Sum_{i=1..j} 1/prime(i).
EXAMPLE
1/2, 4/3, 71/30, 124/35, 11111/2310, 92402/15015, 257189/34034, ...
MATHEMATICA
Table[Sum[(n - k + 1)/Prime[k], {k, 1, n}], {n, 1, 20}] // Denominator
PROG
(Python)
from fractions import Fraction
from sympy import prime, primerange
def a(n): return sum(Fraction(n-i, p) for i, p in enumerate(primerange(1, prime(n)+1))).denominator
print([a(n) for n in range(1, 21)]) # Michael S. Branicky, Jun 09 2022
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Ilya Gutkovskiy, Jun 09 2022
STATUS
approved