|
%I #9 Sep 08 2022 08:45:00
%S -1,-8,28,8,5,5,3,3,4,4,3,3,2,3,3,3,2,3,2,2,2,3,2,2,2,3,3,3,2,2,2,2,2,
%T 2,3,3,2,2,3,3,2,2,2,2,2,2,1,2,2,2,2,3,2,2,2,3,3,3,2,2,2,2,2,2,2,3,2,
%U 2,2,3,2,2,1,1,2,2,2,2,1,2,2,2,1,1,2,2,2,2,2,2,2,2,2,3,3,3,2,2,3,3,2,2,1
%N a(n) = round(n/(log(n)-1) - pi(n)).
%H G. C. Greubel, <a href="/A053619/b053619.txt">Table of n, a(n) for n = 1..10000</a>
%H C. K. Caldwell, <a href="http://www.utm.edu/research/primes/howmany.shtml">How Many Primes Are There?</a>
%t Table[Round[n/(Log[n] -1) -PrimePi[n]], {n, 1, 110}] (* _G. C. Greubel_, May 16 2019 *)
%o (PARI) {a(n) = round(n/(log(n)-1) - primepi(n))}; \\ _G. C. Greubel_, May 16 2019
%o (Magma) [Round(n/(Log(n)-1) - #PrimesUpTo(n)): n in [1..110]]; // _G. C. Greubel_, May 16 2019
%o (Sage) [round(n/(log(n)-1) - prime_pi(n)) for n in (1..110)] # _G. C. Greubel_, May 16 2019
%Y Cf. A000720, A052434, A052435.
%K sign
%O 1,2
%A _Henry Bottomley_, Mar 21 2000
%E More terms from _James A. Sellers_, Mar 23 2000
|
