OFFSET
1,2
COMMENTS
Some of the larger entries may only correspond to probable primes.
MATHEMATICA
e (* from A000043 *) = {2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607, 1279, 2203, 2281, 3217, 4253, 4423, ...}; PrevPrim[n_] := Block[{k = n - 1}, While[ !PrimeQ[k], k-- ]; k]; Do[p = 2^(e[[n]] - 1)*(2^e[[n]] - 1); Print[p - PrevPrim[p]], {n, 1, 20}]
#-NextPrime[#, -1]&/@PerfectNumber[Range[15]] (* The program generates the first 15 terms of the sequence. To generate more, increase the Range constant but the program may take a long time to run. *) (* Harvey P. Dale, Nov 22 2023 *)
CROSSREFS
KEYWORD
hard,more,nonn
AUTHOR
Jason Earls, Dec 01 2002
EXTENSIONS
Edited and extended by Robert G. Wilson v, Dec 03 2002
a(19)-a(26) from Donovan Johnson, Feb 22 2008
STATUS
approved