OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..1560
Eric Weisstein's World of Mathematics, Interprime
EXAMPLE
625 is in this sequence because 625 = 5^4, nontrivial prime power, and 625 = (619+631)/2, with 619 and 631 consecutive primes.
MAPLE
N:= 10^10: # to get all terms <= N
Primes:= select(isprime, [2, seq(i, i=3..isqrt(N), 2)]):
S:= select(t -> t - prevprime(t) = nextprime(t)-t, {seq(seq(p^j, j=2..floor(log[p](N))), p=Primes)}):
sort(convert(S, list)); # Robert Israel, Dec 27 2015
MATHEMATICA
(* version >= 6 *)(#/2 + NextPrime[#]/2) & /@
Select[Prime[Range[5000000]], PrimePowerQ[#/2 + NextPrime[#]/2] &]
(* Wouter Meeussen, Oct 26 2015 *)
PROG
(PARI) {for(i=1, 10^8, if(isprimepower(i)>1&&i==(precprime(i-1)+nextprime(i+1))/2, print1(i, ", ")))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Antonio Roldán, Oct 23 2015
STATUS
approved