OFFSET
1,1
COMMENTS
Observation: All the terms in this sequence, after a(1), are the larger of the two consecutive primes which have positive difference either 2^3 or 4^3.
Superset of A031927 as the sequence contains for example numbers like 89753, 107441, 288647,.. (with gaps of 4^3...) that are not in A031927. - R. J. Mathar, Jun 06 2014
LINKS
K. D. Bajpai, Table of n, a(n) for n = 1..10000
EXAMPLE
97 is prime and appears in the sequence because 97 - 89 = 8 = 2^3.
397 is prime and appears in the sequence because 397 - 389 = 8 = 2^3.
MAPLE
MATHEMATICA
n = 0; Do[t = Prime[k] - Prime[k - 1]; If[IntegerQ[t^(1/3)], n++; Print[n, " ", Prime[k]]], {k, 2, 15*10^4}]
PROG
(PARI) s=[]; forprime(p=3, 4000, if(ispower(p-precprime(p-1), 3), s=concat(s, p))); s \\ Colin Barker, Jun 03 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
K. D. Bajpai, May 31 2014
STATUS
approved