Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #12 Sep 08 2022 08:45:47
%S 433,613,683,773,827,863,1063,1117,1187,1223,1567,1583,1657,1693,1783,
%T 1907,1997,2017,2087,2141,2143,2161,2267,2357,2393,2467,2557,2593,
%U 2609,2663,2719,2753,2789,2843,2879,2897,2969,2971,3023,3041,3061,3167,3187
%N Primes p with a digits sum of p^2 equal to 37.
%H Vincenzo Librandi, <a href="/A165503/b165503.txt">Table of n, a(n) for n = 1..1000</a>
%F {A000040(i) : A123157(i) = 37} [_R. J. Mathar_, Sep 29 2009]
%e 433 is in the sequence because 433^2=187489 and 1+8+7+4+8+9=37.
%e 1783 is in the sequence because 1783^2=3179089 and 3+1+7+9+0+8+9=37.
%t Select[Prime[Range[500]], Total[IntegerDigits[#^2]]== 37&] (* _Vincenzo Librandi_, Sep 12 2013 *)
%o (Magma) [p: p in PrimesUpTo(6150) | &+Intseq(p^2) eq 37]; // _Vincenzo Librandi_, Sep 12 2013
%Y Cf. primes p where the digital sum of p^2 is equal to: A226803 (7), A165492 (13), A165493 (19), A165502 (31), this sequence (37), A165504 (43).
%K nonn,base
%O 1,1
%A _Vincenzo Librandi_, Sep 21 2009
%E More terms from _R. J. Mathar_, Sep 29 2009