OFFSET
1,1
PROG
(PARI) listkd(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-p)); ); Vec(list); }
upto(k) = {my(list = List()); for (i=1, k, my(klist = listkd(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 edited by Michel Marcus, Jan 22 2022
STATUS
approved