|
%I #30 Sep 08 2022 08:46:05
%S 19,71,179,251,449,20249,24499,100549
%N Primes p where the digital sum of p^2 is equal to 10.
%C The next term is > 24154957 (if it exists). - _R. J. Mathar_, Jul 05 2013
%C No more terms below 10^12. - _Hiroaki Yamanouchi_, Sep 23 2014.
%C No additional terms < 10^15. - _Chai Wah Wu_, Nov 15 2015
%e 19 is in the sequence because 19^2=361 and 3+6+1=10.
%e 71 is in the sequence because 71^2=5041 and 5+0+4+1=10.
%p select(p -> isprime(p) and convert(convert(p^2,base,10),`+`)=10, [seq(2*k+1,k=1..100000)]); # _Robert Israel_, Sep 23 2014
%t Select[Prime[Range[70000]], Total[IntegerDigits[#^2]]== 10&]
%o (Magma) [p: p in PrimesUpTo(5*10^6) | &+Intseq(p^2) eq 10];
%Y Cf. A123157, A165492, A165459, A165493. Subset of A135027.
%K nonn,base,more
%O 1,1
%A _Vincenzo Librandi_, Jun 24 2013
| 