A226802 Primes p where the digital sum of p^2 is equal to 10.
19, 71, 179, 251, 449, 20249, 24499, 100549
Offset: 1
Examples
19 is in the sequence because 19^2=361 and 3+6+1=10. 71 is in the sequence because 71^2=5041 and 5+0+4+1=10.
Programs
-
Magma
[p: p in PrimesUpTo(5*10^6) | &+Intseq(p^2) eq 10];
-
Maple
select(p -> isprime(p) and convert(convert(p^2,base,10),`+`)=10, [seq(2*k+1,k=1..100000)]); # Robert Israel, Sep 23 2014
-
Mathematica
Select[Prime[Range[70000]], Total[IntegerDigits[#^2]]== 10&]
Comments