cp's OEIS Frontend

A092494 a(n) = Sum_{p prime and p<=n} ceiling(n/p).

%I A092494 #11 Jun 19 2019 14:13:06
%S A092494 0,1,3,4,6,7,10,11,12,13,16,17,20,21,23,25,27,28,31,32,34,36,39,40,42,
%T A092494 43,45,46,49,50,54,55,56,58,60,62,65,66,68,70,73,74,78,79,81,83,86,87,
%U A092494 89,90,92,94,97,98,100,102,104,106,109,110,114,115,117,119,120,122
%N A092494 a(n) = Sum_{p prime and p<=n} ceiling(n/p).
%C A092494 a(n) = A013939(n) + A048865(n).
%H A092494 Robert Israel, <a href="/A092494/b092494.txt">Table of n, a(n) for n = 1..10000</a>
%p A092494 N:= 100: # for a(1)..a(N)
%p A092494 V:= Vector(N):
%p A092494 p:= 0:
%p A092494 do
%p A092494   p:= nextprime(p);
%p A092494   if p > N then break fi;
%p A092494   V[p]:= V[p]+1;
%p A092494   for k from 2 to floor(N/p) do
%p A092494     V[(k-1)*p+1 .. k*p]:= V[(k-1)*p+1 .. k*p] +~ k;
%p A092494   od;
%p A092494   if (k-1)*p+1<=N then V[(k-1)*p+1..N]:= V[(k-1)*p+1..N]+~ k fi
%p A092494 od:
%p A092494 convert(V,list); # _Robert Israel_, Jun 19 2019
%o A092494 (PARI) a(n) = sum(k=1, n, isprime(k)*ceil(n/k)); \\ _Michel Marcus_, Jun 19 2019
%Y A092494 Cf. A006590.
%K A092494 nonn
%O A092494 1,3
%A A092494 _Reinhard Zumkeller_, Apr 05 2004