login
A226803
Primes p where the digital sum of p^2 is equal to 7.
4
OFFSET
1,1
COMMENTS
No more terms below 10^9. - Michel Marcus, Nov 02 2013
No more terms below 10^20. - Hiroaki Yamanouchi, Sep 23 2014
Sequence is complete, because A215614 is complete. - Robert Israel, Sep 23 2014
EXAMPLE
5 is in the sequence because 5^2 = 25 and 2 + 5 = 7.
149 is in the sequence because 149^2 = 22201 and 2 + 2 + 2 + 0 + 1 = 7.
MATHEMATICA
Select[Prime[Range[10000]], Total[IntegerDigits[#^2]] == 7 &]
PROG
(Magma) [p: p in PrimesUpTo(6*10^6) | &+Intseq(p^2) eq 7];
(PARI) lista(nn) = {forprime(p=2, nn, if (sumdigits(p^2)==7, print1(p, ", ")); ); } \\ Michel Marcus, Nov 02 2013
CROSSREFS
KEYWORD
nonn,base,bref,fini,full
AUTHOR
Vincenzo Librandi, Jun 24 2013
STATUS
approved