|
%I #12 Sep 27 2013 10:57:25
%S 1,4,10,16,31,32,40,41,43,62,71,76,94,95,97,98,121,142,158,163,164,
%T 166,179,188,208,211,214,227,229,259,260,265,284,301,313,317,320,328,
%U 331,340,352,355,356,365,380,382,386,392,397,401,418,424,425,431,436,439
%N Numbers n such that n^2+sum-of-digits(n^2) is prime.
%C This is the sequence of indices of prime numbers in A171613.
%C The Ulam spiral for this sequence is a near-perfect line y=-x (see links).
%H Christian N. K. Anderson, <a href="/A224966/b224966.txt">Table of n, a(n) for n = 1..10000</a>
%H Christian N. K. Anderson, The <a href="/A224966/a224966.jpg">Ulam Spiral</a> for the prime numbers derived from this sequence, i.e. a(n)^2+sum of digits(a(n)^2)
%e a(12)=76 because 76^2=5776, and 5776+(5+7+7+6)=5801, which is prime.
%o (R) library(gmp); digsum<-function(x) sum(as.numeric(unlist(strsplit(as.character(x),split=""))))
%o ans=as.bigz(rep(0,100)); n=1; i=as.bigz(1)
%o while(n<=100) {
%o if(isprime((w=i^2+digsum(i^2)))) ans[(n=n+1)-1]=i
%o i=i+1
%o }; ans
%Y Cf. A048521.
%Y Cf. numbers of the form n^2+sum-of-digits(n^2) A171613, and subsets A171614, A171615.
%Y Cf. A062028.
%K nonn,base
%O 1,2
%A _Kevin L. Schwartz_ and _Christian N. K. Anderson_, Apr 21 2013
|
