%I #40 Mar 30 2023 01:52:36
%S 1,2,3,4,6,10,12,15,33,55,56,61,154,201,462,570,841,848,1297,1357,
%T 2008,2628,2848,3111,6877,9162,17766
%N Numbers k such that (10^k + 2)/6 is prime.
%C Some of the larger entries may only correspond to probable primes.
%C a(26) > 8000. - _Jon E. Schoenfield_, Dec 24 2017
%t Select[Range[5000], PrimeQ[(10^# + 2) / 6] &] (* _Vincenzo Librandi_, Sep 24 2016 *)
%o (PARI) isok(n) = isprime((10^n + 2)/6); \\ _Michel Marcus_, Aug 28 2013
%o (Magma) [n: n in [0..200] |IsPrime((10^n+2) div 6)]; // _Vincenzo Librandi_, Sep 24 2016
%o (Python)
%o from sympy import isprime, prime, nextprime
%o def afind(limit, startk=1):
%o k, tenk = startk, 10**startk
%o for k in range(startk, limit+1):
%o if isprime((tenk + 2)//6):
%o print(k, end=", ")
%o tenk *= 10
%o afind(1000) # _Michael S. Branicky_, Feb 26 2022
%K nonn,more
%O 1,2
%A _Ed Pegg Jr_, Nov 23 2002
%E More terms from _Jason Earls_, Nov 24 2002
%E a(25) from _Jon E. Schoenfield_, Dec 23 2017
%E a(26) from _Michael S. Branicky_, Feb 27 2022
%E a(27) from _Michael S. Branicky_, Mar 29 2023
|