Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #36 Sep 08 2022 08:46:05
%S 67,113,157,193,257,283,311,337,373,409,419,463,509,599,643,653,661,
%T 743,761,769,797,1013,1031,1039,1103,1129,1193,1237,1301,1381,1399,
%U 1427,1471,1481,1553,1571,1579,1597,1733,1759,1823,1831,1877,2029,2039,2111,2129
%N Primes p where the digital sum of p^2 is equal to 25.
%C From _Bruno Berselli_, Sep 12 2013: (Start)
%C Primes q such that the digital sum of q^2 is 1 < k < 50:
%C k | q
%C ---|------------
%C 4 | 2, 11, 101;
%C 7 | A226803;
%C 9 | 3;
%C 10 | A226802;
%C 13 | A165492;
%C 16 | A165459;
%C 19 | A165493;
%C 22 | 43, 97, 191, 227, 241, 317, 331, 353, ... ;
%C 25 | this sequence;
%C 28 | 163, 197, 233, 307, 359, 397, 431, 467, ... ;
%C 31 | A165502;
%C 34 | 167, 293, 383, 563, 607, 617, 733, 787, ... ;
%C 37 | A165504;
%C 40 | 313, 947, 983, 1303, 1483, 1609, 1663, ... ;
%C 43 | A165504;
%C 46 | 883, 937, 1367, 1637, 2213, 2447, 2683, ... ;
%C 49 | 1667, 2383, 2437, 2617, 2963, 4219, 4457, ... . (End)
%H Vincenzo Librandi, <a href="/A229058/b229058.txt">Table of n, a(n) for n = 1..1000</a>
%t Select[Prime[Range[400]], Total[IntegerDigits[#^2]] == 25 &]
%o (Magma) [p: p in PrimesUpTo(2600) | &+Intseq(p^2) eq 25];
%Y Cf. A165459, A165492, A165493, A165502-A165504, A226802, A226803.
%K nonn
%O 1,1
%A _Vincenzo Librandi_, Sep 12 2013