[go: up one dir, main page]

login
A074633
a(n) is the smallest index m such that Sum_{k=1..m} 1/PrimePi(k) >= n, where PrimePi()=A000720().
4
2, 4, 8, 12, 18, 27, 37, 51, 68, 89, 116, 147, 186, 232, 287, 352, 428, 518, 623, 745, 887, 1050, 1240, 1456, 1704, 1987, 2309, 2674, 3090, 3557, 4087, 4684, 5353, 6105, 6949, 7892, 8944, 10121, 11431, 12885, 14502, 16298, 18286, 20485, 22917, 25607
OFFSET
1,1
EXAMPLE
a(85) = 927685 because 927686 is the smallest m such that Sum_{k=1..m} 1/PrimePi(k) >= 85.
MATHEMATICA
{s=0, s1=0}; Do[s=s+(1/PrimePi[n]); If[Greater[Floor[s], s1], s1=Floor[s]; Print[{n, Floor[s]}]], {n, 2, 1000000}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Aug 28 2002
EXTENSIONS
Edited by Jon E. Schoenfield, Apr 04 2023
STATUS
approved