OFFSET
1,1
COMMENTS
PROG
(PARI) listkq(k) = {my(list = List()); forprime(p=2^(k-1)+1, 2^k, my(q=2^(k+1)-p); if ((q>p) && isprime(q), listput(list, q)); ); Vec(list); }
upto(k) = {my(list = List()); for (i=1, k, my(klist = listkq(i)); if (#klist, for (j=1, #klist, listput(list, klist[j]))); ); Vec(list); }
upto(11) \\ Michel Marcus, Jan 22 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jul 08 2002
EXTENSIONS
Name corrected by Jon E. Schoenfield, Jun 30 2021
STATUS
approved