OFFSET
1,1
COMMENTS
The sequence consists of the prime numbers p that are champions (left to right maxima) of the function |psi(p^3)-p^3|, where psi(p) is the Chebyshev psi function.
LINKS
M. Planat and P. Solé, Efficient prime counting and the Chebyshev primes, arXiv:1109.6489 [math.NT], 2011.
MATHEMATICA
ChebyshevPsi[n_] := Range[n] // MangoldtLambda // Total;
Reap[For[max=0; p=2, p < 1000, p = NextPrime[p], f = Abs[ChebyshevPsi[p^3] - p^3]; If[f > max, max = f; Print[p]; Sow[p]]]][[2, 1]] (* Jean-François Alcover, Dec 03 2018 *)
PROG
(Perl) use ntheory ":all"; my($max, $f)=(0); forprimes { $f=abs(chebyshev_psi($_**3)-$_**3); if ($f > $max) { say; $max=$f; } } 1000; # Dana Jacobsen, Dec 28 2015
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Michel Planat, Oct 11 2011
EXTENSIONS
More terms from Dana Jacobsen, Dec 28 2015
STATUS
approved