|
%I #33 Sep 08 2022 08:46:05
%S 5,149,1049
%N Primes p where the digital sum of p^2 is equal to 7.
%C No more terms below 10^9. - _Michel Marcus_, Nov 02 2013
%C No more terms below 10^20. - _Hiroaki Yamanouchi_, Sep 23 2014
%C Sequence is complete, because A215614 is complete. - _Robert Israel_, Sep 23 2014
%e 5 is in the sequence because 5^2 = 25 and 2 + 5 = 7.
%e 149 is in the sequence because 149^2 = 22201 and 2 + 2 + 2 + 0 + 1 = 7.
%t Select[Prime[Range[10000]], Total[IntegerDigits[#^2]] == 7 &]
%o (Magma) [p: p in PrimesUpTo(6*10^6) | &+Intseq(p^2) eq 7];
%o (PARI) lista(nn) = {forprime(p=2, nn, if (sumdigits(p^2)==7, print1(p, ", ")););} \\ _Michel Marcus_, Nov 02 2013
%Y Cf. A226802, A165492, A165459, A165493, A215614.
%K nonn,base,bref,fini,full
%O 1,1
%A _Vincenzo Librandi_, Jun 24 2013
| 