A052279 Primes such that the sum of the factorials of the digits is a perfect square.
17, 41, 71, 211, 433, 457, 547, 1013, 1031, 1103, 1301, 1489, 2063, 3001, 3011, 4451, 5077, 5441, 5651, 5717, 6203, 6551, 7057, 7507, 7517, 8419, 8941, 10163, 10613, 10631, 16103, 16301, 20369, 20639, 20693, 20873, 20963, 21313, 21661, 23003, 23087, 23131
Offset: 1
Links
- Giovanni Resta, Table of n, a(n) for n = 1..10000
Programs
-
Mathematica
Select[Prime[Range[3000]],IntegerQ[Sqrt[Total[IntegerDigits[#]!]]]&] (* Harvey P. Dale, Apr 11 2016 *)
-
PARI
isok(n) = isprime(n) && (d = digits(n)) && issquare(sum(i=1, #d, d[i]!)); \\ Michel Marcus, Jan 06 2014
Extensions
More terms from Michel Marcus, Jan 06 2014