login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A296443
Numbers k such that 2*10^(2k)-2*10^k+1 is prime.
1
1, 2, 6, 7, 8, 315, 667, 5125, 7301, 10500, 11096
OFFSET
1,2
COMMENTS
Numbers of this form divide 4*10^(4k)+1.
a(8) > 5000. - Jon E. Schoenfield, Dec 16 2017
a(12) > 25000. - Michael S. Branicky, Aug 08 2024
EXAMPLE
181, 19801, 1999998000001, 199999980000001, and 19999999800000001 are prime, while 1998001=277*7213, 199980001=13*41*457*821, and 19999800001=53*5953*63389.
MATHEMATICA
ParallelMap[If[PrimeQ[2*10^(2 #) - 2*10^# + 1], #, Nothing] &, Range@ 4000] (* Robert G. Wilson v, Dec 13 2017 *)
PROG
(PARI) isok(k) = isprime(2*10^(2*k)-2*10^k+1); \\ Michel Marcus, Dec 13 2017
CROSSREFS
See A296444 for 2*10^(2k)+2*10^k+1.
Sequence in context: A196747 A261691 A351088 * A102046 A019913 A379710
KEYWORD
nonn,more
AUTHOR
Patrick A. Thomas, Dec 13 2017
EXTENSIONS
a(6)-a(7) from Michel Marcus, Dec 13 2017
a(8)-a(11) from Michael S. Branicky, Mar 31 2023
STATUS
approved