%I
%S 11,13,17,19,31,73,97,101,103,107,109,127,251,277,347,367,431,433,457,
%T 499,613,631,677,691,701,709,727,811,877,941,947,1009,1063,1069,1087,
%U 1171,1229,1283,1289,1399,1423,1429,1531,1553,1597,1607,1621,1667,1669,1733
%N Primes p such that p + sum of squares of digits of p is also prime.
%H Vincenzo Librandi, <a href="/A076162/b076162.txt">Table of n, a(n) for n = 1..1000</a>
%e 13 is OK because 13+(1^2+3^2)=23 is a prime.
%t Select[Prime[Range[2,300]],PrimeQ[#+Total[IntegerDigits[#]^2]]&] (* _Vincenzo Librandi_, Apr 13 2013 *)
%K nonn,base
%O 1,1
%A _Zak Seidov_, Nov 01 2002
%E More terms from _Vincenzo Librandi_, Apr 13 2013
