login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Primes p where the digital sum of p^2 is equal to 25.
2

%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