A297422 Numbers k such that 4*10^(2k) + 2*(10^k) + 1 is prime.

0, 1, 3, 5, 10, 196, 251, 310, 415, 1107, 2422

Offset

1,3

Comments

Numbers of this form divide 8*(10^3k)-1.

a(12) > 50000. - Robert Price, Apr 11 2018

Links

Table of n, a(n) for n=1..11.

Example

7, 421, 4002001, and 40000200001 are prime, while 40201=7*5743, 400020001=31*12903871, etc.

Maple

select(k->isprime(4*10^(2*k)+2*(10^k)+1), [$0..1000]); # Muniru A Asiru, Apr 11 2018

Mathematica

Select[Range[1500], PrimeQ[4 (10^(2 #)) + 2 (10^#) + 1] &] (* Michael De Vlieger, Dec 30 2017 *)

Prog

(PARI) isok(k) = isprime(4*10^(2*k)+2*(10^k)+1); \\ Michel Marcus, Dec 30 2017

Crossrefs

Sequence in context: A024329 A354517 A354389 * A064523 A365443 A268034

Adjacent sequences: A297419 A297420 A297421 * A297423 A297424 A297425

Keyword

nonn,more

Author

Patrick A. Thomas, Dec 30 2017

Extensions

a(6)-a(10) from Michael De Vlieger, Dec 30 2017

a(11) from Jon E. Schoenfield, Jan 01 2018

Status

approved