A064662 a(n) = n - round( prime(n)/log(prime(n)) ).

-2, -1, 0, 0, 0, 1, 1, 2, 2, 1, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 6, 6, 5, 5, 5, 6, 5, 6, 6, 6, 6, 6, 6, 7, 7, 7, 8, 8, 8, 7, 7, 8, 8, 8, 9, 9, 9, 9, 9, 10, 10, 10, 11, 10, 9, 10, 11, 11, 10, 10, 10, 10, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12

Offset: 1

Author

N. J. A. Sloane, Oct 10 2001

Keywords

• sign

Links

• Harry J. Smith, Table of n, a(n) for n = 1..1000

Crossrefs

Cf. A064659-A064663.

Programs

Magma

[n - Round(NthPrime(n) / Log(NthPrime(n))): n in [1..100]]; // Vincenzo Librandi, May 25 2019

Mathematica

rpl[n_]:=Module[{p=Prime[n]},n-Round[p/Log[p]]]; Array[rpl,90] (* Harvey P. Dale, May 20 2019 *)

PARI

{ for (n=1, 1000, write("b064662.txt", n, " ", n - round( prime(n)/log(prime(n)) )) ) } \\ Harry J. Smith, Sep 21 2009