A074469 Least m such that Sigma-Composite-Harmonic series Sum_{k=1..m} 1/A000203(A002808(k)) >= n.
32, 301, 2123, 13172, 76105, 420007, 2245009, 11719362, 60071831, 303487314, 1515211979
Offset: 1
Crossrefs
Programs
-
Mathematica
c[x_] := FixedPoint[x+PrimePi[ # ]+1&, x] {s=0, s1=0}; Do[s=s+(1/DivisorSigma[1, c[n]]); If[Greater[Floor[s], s1], s1=Floor[s]; Print[{n, Floor[s]}]], {n, 1, 1000000}]
-
PARI
a(n)=my(m,s=0.);for(c=4,(2*n+2)^(n+2),if(isprime(c),next,m++);s+=1/sigma(c);if(s>=n,return(m))) \\ Charles R Greathouse IV, Feb 19 2013
Extensions
a(6)-a(11) from Donovan Johnson, Aug 22 2011