OFFSET
1,1
COMMENTS
I have searched through prime powers up to 2^8532. It is very unlikely that there are any other values between the ones listed here, but no prime has been proved to be absent from this sequence. - David Wasserman, Mar 31 2005
EXAMPLE
3 = 128 - 125 = 2^7 - 5^3; 7 = 16 - 9 = 32768 - 32761; 17 = 49 - 32 = 81 - 64 = 529 - 512; 4241 = 528529 - 524288 = 727^2 - 2^19.
MATHEMATICA
pp = Sort[ Flatten[ Table[ Prime[n]^i, {n, 1, PrimePi[ Sqrt[10^16]]}, {i, 1, Log[ Prime[n], 10^16]}]]]; l = Length[pp]; d = Sort[Take[pp, -l + 1] - Take[pp, l - 1]]; a = {}; Do[ If[ PrimeQ[ d[[n]]], a = Append[a, d[[n]]]], {n, 1, l - 1}]; Union[a] a = {}; Do[ If[ PrimeQ[ pp[[n + 1]] - pp[[n]]], a = Append[a, pp[[n + 1]] - pp[[n]]]], {n, 1, Length[pp] - 1}]; Union[a]
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov and Robert G. Wilson v, Oct 29 2002
EXTENSIONS
More terms from David Wasserman, Mar 31 2005
STATUS
approved