A120055 a(n) is the least semiprime s such that the Sum_{semiprime i <= s} 1/i exceeds n.
4, 35, 871, 43217, 5296623, 2011783478
Offset: 0
Examples
a(0)=4 because 1/4 > 0. a(1)=35 because 1/4 + 1/6 + 1/9 + 1/10 + 1/14 + 1/15 + 1/21 + 1/22 + 1/25 + 1/26 + 1/33 + 1/34 + 1/35 = 15708817/15315300 > 1.
Programs
-
Mathematica
s = 0; k = 1; Do[ While[s <= n, If[ Plus @@ Last /@ FactorInteger@k == 2, s = N[ s + 1/k, 20]]; k++ ]; Print[{k - 1, s}]; k ++, {n, 0, 5}]
Comments