OFFSET
1,1
MATHEMATICA
f[n_] := IntegerExponent[n + 2^IntegerExponent[n, 2], 2]; f/@ (Select[Table[2^t-2^k-1, {t, 1, 20}, {k, 1, t-1}] // Flatten // Union, PrimeQ] + 1) (* Amiram Eldar, Dec 17 2018 after Jean-François Alcover at A181741 *)
PROG
(PARI) listt(nn) = {for (n=3, nn, forstep(k=n-1, 1, -1, if (isprime(2^n-2^k-1), print1(n, ", ")); ); ); } \\ Michel Marcus, Dec 17 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, Nov 08 2010
EXTENSIONS
Corrected and extended by Michel Marcus, Dec 17 2018
STATUS
approved