OFFSET
1,1
LINKS
Christian N. K. Anderson, Table of n, a(n) for n = 1..10000
EXAMPLE
269 is a prime number, and 2^2+6^2+9^2 = 121 = 11^2.
MATHEMATICA
Select[Prime[Range[2000]], IntegerQ[Sqrt[Total[IntegerDigits[#]^2]]] &] (* T. D. Noe, Apr 05 2013 *)
PROG
(R) ssod<-function(i) sum(as.numeric(strsplit(as.character(i), "")[[1]])^2)
issquare<-function(x) as.integer(sqrt(x))==sqrt(x)
x=as.bigz(c()); i=2
while(length(x)<10000) {if(issquare(ssod(i))) x=c(x, i); i=nextprime(i)}
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Kevin L. Schwartz and Christian N. K. Anderson, Apr 04 2013
STATUS
approved