%I #22 Sep 08 2022 08:45:47
%S 7,29,47,61,79,151,349,389,461,601,1051,1249,1429,1451,1789,2239,3499,
%T 4799,5011,5101,5501,5749,10501,11149,14249,50101,51001,55001,100501,
%U 100649,101149,150001,318751,448999,501001,1009951,5000011,5000101,5001001,50000101
%N Primes p where the digital sum of p^2 is equal to 13.
%H Vincenzo Librandi and Hiroaki Yamanouchi, <a href="/A165492/b165492.txt">Table of n, a(n) for n = 1..61</a> (terms a(1)-a(44) from Vincenzo Librandi)
%F {A000040(i) : A123157(i) = 13}. [_R. J. Mathar_, Sep 29 2009]
%e 29 is in the sequence because 29^2 = 841 and 8+4+1 = 13.
%e 1789 is in the sequence because 1789^2 = 3200521 and 3+2+0+0+5+2+1 = 13.
%t Select[Prime[Range[80000]],Total[IntegerDigits[#^2]]==13&] (* _Harvey P. Dale_, Apr 08 2011 *)
%o (Magma) [p: p in PrimesUpTo(5*10^7) | &+Intseq(p^2) eq 13]; // _Vincenzo Librandi_, Sep 26 2015
%Y Cf. A165502 - A165504.
%K nonn,base
%O 1,1
%A _Vincenzo Librandi_, Sep 21 2009
%E More terms from _R. J. Mathar_, Sep 29 2009
%E a(37)-a(40) from _Hiroaki Yamanouchi_, Sep 23 2014
|