login
A226802
Primes p where the digital sum of p^2 is equal to 10.
3
19, 71, 179, 251, 449, 20249, 24499, 100549
OFFSET
1,1
COMMENTS
The next term is > 24154957 (if it exists). - R. J. Mathar, Jul 05 2013
No more terms below 10^12. - Hiroaki Yamanouchi, Sep 23 2014.
No additional terms < 10^15. - Chai Wah Wu, Nov 15 2015
EXAMPLE
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.
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&]
PROG
(Magma) [p: p in PrimesUpTo(5*10^6) | &+Intseq(p^2) eq 10];
CROSSREFS
Sequence in context: A196136 A198002 A093350 * A342737 A142516 A127874
KEYWORD
nonn,base,more
AUTHOR
Vincenzo Librandi, Jun 24 2013
STATUS
approved