cp's OEIS Frontend

A267818 Primes whose base-4 representation is a square when read in base 10.

%I A267818 #24 Jun 29 2019 02:03:39
%S A267818 673,1165921,6228961,53378329,128818873,805379161,926711449,
%T A267818 1542173569,2459671993,8365652929,14232980953,26484614209,30755887321,
%U A267818 31937961241,37968797689,78180076129,80974979329,97835803129,168394504129,314911422913,410879716321,529074406681
%N A267818 Primes whose base-4 representation is a square when read in base 10.
%C A267818 Primes in the sequence A267764.
%H A267818 Chai Wah Wu, <a href="/A267818/b267818.txt">Table of n, a(n) for n = 1..100</a>
%e A267818 673 written in base 4 is 22201, and 22201 in base 10 is 149^2.
%t A267818 Select[Prime@ Range[10^6], IntegerQ@ Sqrt@ FromDigits@ IntegerDigits[#, 4] &] (* _Michael De Vlieger_, Jan 24 2016 *)
%o A267818 (PARI) is(n, b=4, c=10)=issquare(subst(Pol(digits(n, b)), x, c))
%o A267818 lista(nn) = forprime(p=2, nn, if (is(p), print1(p, ", "))) \\ _Michel Marcus_, Jan 24 2016
%o A267818 (Python)
%o A267818 from sympy import isprime
%o A267818 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
%Y A267818 In other bases: A267817, A267818, A267819, A267820, A267821, A241246, A267490.
%K A267818 nonn,base
%O A267818 1,1
%A A267818 _Christopher Cormier_, Jan 20 2016
%E A267818 a(8)-a(22) from _Chai Wah Wu_, Feb 22 2016