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

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

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