|
%I #22 Sep 12 2018 09:47:51
%S 2,3,5,7,13,19,23,31,43,47,73,113,199,283,467,661,887,1129,1327,1627,
%T 2803,3947,4297,5881,6379,7043,9949,10343,13187,15823,18461,24137,
%U 33647,34763,37663,42863,43067,59753,59797,82619,96017,102679
%N Convex hull primes, that is, prime numbers corresponding to the convex hull of PrimePi, the prime counting function.
%C "Convex hull of PrimePi" is a short wording for "the upper convex hull of the points {p, PrimePi(p)} for p >= 2".
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Convex_hull">Convex hull</a>
%e Prime 83 is not member because there exist two primes from the convex hull, namely 47 and 113, such that (PrimePi(83) - PrimePi(47))/(83 - 47) < (PrimePi(113) - PrimePi(83))/(113 - 83).
%t terms = 42;
%t pMax = 110000;
%t a[1] = 2;
%t a[n_] := a[n] = Module[{}, For[slopeMax = 0; p1 = NextPrime[a[n-1]], p1 <= pMax, p1 = NextPrime[p1], slope = (PrimePi[p1] - PrimePi[a[n-1]])/(p1 - a[n-1]); If[slope > slopeMax, slopeMax = slope; p1Max = p1]]; p1Max];
%t Table[Print["a(", n, ") = ", a[n]]; a[n], {n, 1, 42}]
%Y Cf. A000720, A124661, A167844, A246033 (a subsequence).
%K nonn,more
%O 1,1
%A _Jean-François Alcover_, Sep 11 2018
|
