login
a(n) = Sum_{k=1..n} round(prime(k) / log(prime(k))).
1

%I #11 Jan 21 2019 19:01:55

%S 3,6,9,13,18,23,29,35,42,51,60,70,81,92,104,117,131,146,162,179,196,

%T 214,233,253,274,296,318,341,364,388,414,441,469,497,527,557,588,620,

%U 653,687,722,757,793,830,867,905,944,985,1027,1069,1112,1156,1200,1245,1291

%N a(n) = Sum_{k=1..n} round(prime(k) / log(prime(k))).

%C Partial sums of A002036. - _Robert Israel_, Jul 05 2015

%H Robert Israel, <a href="/A092847/b092847.txt">Table of n, a(n) for n = 1..10000</a>

%p N:= 1000: # to get a(1) to a(N)

%p M:= ithprime(N):

%p Primes:= select(isprime,[2,seq(2*i+1,i=1..(M-1)/2)]):

%p A002036:= map(t -> round(t/log(t)),Primes):

%p ListTools:-PartialSums(%); # _Robert Israel_, Jul 05 2015

%Y Cf. A002036.

%K easy,nonn

%O 1,1

%A _Jorge Coveiro_, Apr 15 2004