|
%I #9 Jan 08 2019 09:12:41
%S 2,3,11,13,23,31,41,59,67,101,103,113,131,139,157,193,211,229,239,257,
%T 283,311,337,347,373,401,409,419,463,491,499,509,571,599,643,653,661,
%U 743,751,761,769,797,1013,1021,1031,1039,1103,1129,1193,1201,1229,1237
%N Prime numbers whose squares have square digit-sums.
%H Harvey P. Dale, <a href="/A076503/b076503.txt">Table of n, a(n) for n = 1..1000</a>
%e 13 is a member because 13 is prime and the digit-sum of its square is 1+6+9=16, which is also square.
%t Select[Prime[Range[250]],IntegerQ[Sqrt[Total[IntegerDigits[#^2]]]]&] (* _Harvey P. Dale_, Aug 30 2016 *)
%o (PARI) isok(n) = if (! isprime(n), 0, d = digits(n^2); issquare(sum(i=1, #d, d[i]))) \\ _Michel Marcus_, Jun 20 2013
%Y Cf. A000040, A061910.
%K nonn,base
%O 1,1
%A Christopher Schloetz (cschloetz(AT)hotmail.com), Nov 09 2002
%E More terms from _Michel ten Voorde_ Jun 13 2003
|
