OFFSET
1,1
COMMENTS
From Flávio V. Fernandes, May 26 2021: (Start)
Equivalently, primes p such that p is the largest of all prime factors of the numbers in the interval [2*p, nextprime(2*p)-1].
For any prime p, if p is not the largest of all prime factors of the numbers in that interval (i.e., if p is not a term of this sequence), then the largest of all prime factors of the numbers in that interval will be a prime q that occurs in the number 2*q.
For all n, the largest prime < 2*a(n) is a term of A059788. (End)
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
f(precprime(2*p)) = p, where f is the mapping defined by A052248.
EXAMPLE
5 is a term since 7 is the prime preceding 2*5, 11 is the next prime and 5 is the largest of all prime factors of 8, 9 and 10.
MATHEMATICA
Select[Range[300], PrimeQ[#] && NextPrime[2*#] < 2 * NextPrime[#] &] (* Amiram Eldar, Feb 07 2020 *)
PROG
(PARI) {forprime(k=2, 317, p=precprime(2*k); q=nextprime(p+1); m=0; for(j=p+1, q-1, f=factor(j); a=f[matsize(f)[1], 1]; if(m<a, m=a)); if(m==k, print1(k, ", ")))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Klaus Brockhaus, Feb 10 2003
STATUS
approved