

A226802


Primes p where the digital sum of p^2 is equal to 10.


3




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


LINKS

Table of n, a(n) for n=1..8.


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

Cf. A123157, A165492, A165459, A165493. Subset of A135027.
Sequence in context: A196136 A198002 A093350 * A342737 A142516 A127874
Adjacent sequences: A226799 A226800 A226801 * A226803 A226804 A226805


KEYWORD

nonn,base,more


AUTHOR

Vincenzo Librandi, Jun 24 2013


STATUS

approved



