A076503 Prime numbers whose squares have square digit-sums.

2, 3, 11, 13, 23, 31, 41, 59, 67, 101, 103, 113, 131, 139, 157, 193, 211, 229, 239, 257, 283, 311, 337, 347, 373, 401, 409, 419, 463, 491, 499, 509, 571, 599, 643, 653, 661, 743, 751, 761, 769, 797, 1013, 1021, 1031, 1039, 1103, 1129, 1193, 1201, 1229, 1237

Offset

1,1

Links

Harvey P. Dale, Table of n, a(n) for n = 1..1000

Example

13 is a member because 13 is prime and the digit-sum of its square is 1+6+9=16, which is also square.

Mathematica

Select[Prime[Range[250]], IntegerQ[Sqrt[Total[IntegerDigits[#^2]]]]&] (* Harvey P. Dale, Aug 30 2016 *)

Prog

(PARI) isok(n) = if (! isprime(n), 0, d = digits(n^2); issquare(sum(i=1, #d, d[i]))) \\ Michel Marcus, Jun 20 2013

Crossrefs

Cf. A000040, A061910.

Sequence in context: A089151 A115669 A036956 * A004678 A107715 A090707

Adjacent sequences: A076500 A076501 A076502 * A076504 A076505 A076506

Keyword

nonn,base

Author

Christopher Schloetz (cschloetz(AT)hotmail.com), Nov 09 2002

Extensions

More terms from Michel ten Voorde Jun 13 2003

Status

approved