%I #48 Jun 28 2024 14:00:53
%S 1,3,59,361139,43922730588128390
%N a(n) = smallest k such that Sum_{ i = 1..k } 1/prime(i) exceeds n.
%C The corresponding primes prime(a(n)) are in A016088.
%C Index m for which the prime harmonic number p[ m ] := Sum[ 1/Prime[ k ],{k,1,m} ] >= n.
%H E. Bach, D. Klyve, and J. P. Sorenson, <a href="http://dx.doi.org/10.1090/S0025-5718-09-02249-2">Computing prime harmonic sums</a>, Math. Comp. 78 (2009) 2283-2305.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PrimeNumber.html">Prime Number</a>.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/HarmonicSeriesofPrimes.html">Harmonic Series of Primes</a>.
%F From _Jonathan Sondow_, Apr 17 2013: (Start)
%F a(n) = A000720(A016088(n)) = A000720(A096232(n))+1.
%F a(n) = e^(e^(n + O(1))), see comment in A223037. [corrected by _Charles R Greathouse IV_, Aug 22 2013] (End)
%F a(n) = A103591(2*n). - _Michel Marcus_, Aug 22 2013
%t Table[m = 1; s = 0; While[(s = s + 1/Prime[m]) <= n, m++];
%t m, {n, 0, 4}] (* _Robert Price_, Mar 27 2019 *)
%o (PARI) a(n)=my(t); forprime(p=2,, t+=1./p; if(t>n, return(primepi(p)))) \\ _Charles R Greathouse IV_, Apr 29 2015
%Y Cf. A004080, A016088, A096232, A223037.
%Y Cf. A024451/A002110 for Sum_{i = 1..n} 1/prime(i).
%K nonn,more,nice
%O 0,2
%A _Eric W. Weisstein_
%E a(4) found by Tomás Oliveira e Silva (tos(AT)det.ua.pt), using the fourth term of A016088. - Dec 14 2005
%E a(0) from _Jonathan Sondow_, Apr 16 2013