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