OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..1615
EXAMPLE
37 and 37+6=43 are both lucky primes, so 37 is in the sequence.
MAPLE
N:= 10^5: # for terms <= N
L:= [seq(i, i=1..N+6, 2)]:
for n from 2 while n < nops(L) do
r:= L[n];
L:= subsop(seq(r*i=NULL, i=1..nops(L)/r), L);
od:
L:= convert(select(isprime, L), set):
A:= L intersect map(`-`, L, 6):
sort(convert(A, list));
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel, Jul 26 2019
STATUS
approved