A335266 - OEIS

%I #21 Sep 04 2020 15:56:08

%S 1,1,3,3,5,5,5,11,15,7,13,7,13,31,13,9,5,21,3,21,51,17,9,39,15,3,57,

%T 27,3,15,35,11,39,55,13,27,11,35,17,11,15,13,1,13,19,75,5,23,41,9,63,

%U 63,13,25,31,21,11,9,27,11,3,31,21,35,47,33

%N Numbers k such that b(n) = 2*b(n - 1) + k is prime, where b(0) = 1 and k is the smallest positive number.

%C Sequence b, with an offset change, corresponds to A110930. - _Rémy Sigrist_, Jul 18 2020

%t b[0] = 1; b[n_] := b[n] = NextPrime[2*b[n - 1]]; Table[b[n] - 2*b[n - 1], {n, 1, 66}] (* _Amiram Eldar_, Aug 18 2020 *)

%o (PARI) b=1; for(n = 1, 66, b = 2 * b; p = nextprime(b + 1); k = p - b; b = p; print1(k", "))

%Y Cf. A000040, A110930.

%K nonn

%O 0,3

%A _Dimitris Valianatos_, May 29 2020