OFFSET
1,1
COMMENTS
Primes in the sequence A267764.
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..100
EXAMPLE
673 written in base 4 is 22201, and 22201 in base 10 is 149^2.
MATHEMATICA
Select[Prime@ Range[10^6], IntegerQ@ Sqrt@ FromDigits@ IntegerDigits[#, 4] &] (* Michael De Vlieger, Jan 24 2016 *)
PROG
(PARI) is(n, b=4, c=10)=issquare(subst(Pol(digits(n, b)), x, c))
lista(nn) = forprime(p=2, nn, if (is(p), print1(p, ", "))) \\ Michel Marcus, Jan 24 2016
(Python)
from sympy import isprime
A267818_list = [int(d, 4) for d in (str(i**2) for i in range(1, 10**6)) if max(d) < '4' and isprime(int(d, 4))] # Chai Wah Wu, Feb 22 2016
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Christopher Cormier, Jan 20 2016
EXTENSIONS
a(8)-a(22) from Chai Wah Wu, Feb 22 2016
STATUS
approved