|
%I #25 Dec 26 2017 18:33:54
%S 4,7,8,11,16,22,28,29,79,91,170,293,392,1898,2597,3641,4099,8236,
%T 11393,17969,18503
%N Numbers k such that 10^k/8 - 1 is prime.
%C Searched up to k=20000. - _Michael B. Porter_, Dec 07 2017
%C Equivalently, numbers k such that 125*10^(k-3) - 1 is prime.
%e 1249, 1249999, and 12499999 are prime, so 4, 7, and 8 are in the sequence.
%t Select[Range[2000], PrimeQ[10^#/8 - 1] &] (* _Michael De Vlieger_, Dec 08 2017 *)
%o (PARI) lista(nn) = for(n=4, nn, if(ispseudoprime(10^n/8 - 1), print1(n, ", "))) \\ _Iain Fox_, Dec 04 2017
%K nonn,more
%O 1,1
%A _Patrick A. Thomas_, Dec 03 2017
%E a(12), a(13), a(15), a(16) from _Jon E. Schoenfield_, Dec 03 2017
%E a(14) and a(17)-a(21) from _Michael B. Porter_, Dec 07 2017
|
