OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
map(t -> t - 4*10^ilog10(t), select(t -> isprime(t - 4*10^ilog10(t)), [seq(seq(i^2, i=ceil(sqrt(4.1*10^d)) .. floor(sqrt(5*10^d))), d=1..10)])); # Robert Israel, Nov 08 2023
MATHEMATICA
Select[Prime[Range[10^4]], IntegerQ[Sqrt[FromDigits[Join[{4}, IntegerDigits[#]]]]]&] (* G. C. Greubel, Jun 23 2016 *)
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Claudio Meller, Nov 10 2009
EXTENSIONS
More terms from Michel Marcus, Aug 05 2013
STATUS
approved