OFFSET
1,1
FORMULA
prime(a(n)) = A053607(n).
EXAMPLE
Primes 18 and 19 are 61 and 67, and the interval (62, 63, 64, 65, 66) contains the prime-power 64, so 18 is in the sequence.
MATHEMATICA
Select[Range[100], Length[Select[Range[Prime[#]+1, Prime[#+1]-1], PrimePowerQ]]>=1&]
PROG
(Python)
from itertools import count, islice
from sympy import factorint, nextprime
def A377057_gen(): # generator of terms
p, q, k = 2, 3, 1
for k in count(1):
if any(len(factorint(i))<=1 for i in range(p+1, q)):
yield k
p, q = q, nextprime(q)
CROSSREFS
The corresponding primes are A053607.
For no prime-powers we have A377286.
For exactly one prime-power we have A377287.
A000015 gives the least prime-power >= n.
A031218 gives the greatest prime-power <= n.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Oct 25 2024
STATUS
approved