OFFSET
1,2
COMMENTS
There are fewer consecutive terms in the sequence (330 in the first 10000) than expected by chance (p~0.000000005, determined by bootstrap resampling).
LINKS
Christian N. K. Anderson, Table of n, a(n) for n = 1..10000
EXAMPLE
467 is in this sequence because 467^2 = 218089, and 218089 + (2 + 1 + 8 + 0 + 8 + 9)^2 = 218873, which is prime.
PROG
(R) library(gmp); ya=rep(0, 100); cv=as.bigz(0); len=0
while(len<100) if(isprime((cv=cv+1)^2+digsum(cv^2)^2)) ya[(len=len+1)]=as.numeric(cv)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Kevin L. Schwartz and Christian N. K. Anderson, Apr 26 2013
STATUS
approved