login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

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