OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
a(3)=51 because 51^2 = 2601 is the third perfect square whose binary representation 101000101001 read as the decimal one 101000101001 is prime.
MAPLE
a:=proc(n) if isprime(convert(n^2, binary))=true then n else fi end: seq(a(n), n=1..2400); # Emeric Deutsch, Jul 04 2005
MATHEMATICA
Select[Range[1800], PrimeQ[FromDigits[IntegerDigits[#^2, 2]]]&] (* Harvey P. Dale, May 23 2021 *)
PROG
(PARI) isok(n) = isprime(fromdigits(binary(n^2))); \\ Michel Marcus, Sep 24 2018
CROSSREFS
KEYWORD
easy,nonn,base
AUTHOR
Alexandre Wajnberg, Jun 30 2005
EXTENSIONS
More terms from Emeric Deutsch, Jul 04 2005
STATUS
approved