OFFSET
1,1
COMMENTS
The squares are in A225873.
LINKS
Harvey P. Dale and Davin Park, Table of n, a(n) for n = 1..20000 [Terms 1 through 1000 were computed by Harvey P. Dale and terms 1001 through 20000 by Davin Park]
MATHEMATICA
Select[Range[600], PrimeQ[FromDigits[Rest[IntegerDigits[#^2]]]]&] (* Harvey P. Dale, Dec 13 2015 *)
b /. Flatten[Outer[Solve[a + #2*10^#1 == b^2 && 0 <= a < 10^#1 && Sqrt[#2*10^#1] <= b < Sqrt[10^(#1 + 1)] && a \[Element] Primes, {a, b}, Integers] &, Range[0, 5], Range[9]], 2] (* Davin Park, Dec 30 2016 *)
PROG
(PARI) is=(n)->isprime(n^2%10^(#Str(n^2)-1))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
M. F. Hasler, Nov 01 2014
EXTENSIONS
Extended by Davin Park, Dec 30 2016
STATUS
approved