A052279 - OEIS

%I #12 Apr 11 2016 15:45:11

%S 17,41,71,211,433,457,547,1013,1031,1103,1301,1489,2063,3001,3011,

%T 4451,5077,5441,5651,5717,6203,6551,7057,7507,7517,8419,8941,10163,

%U 10613,10631,16103,16301,20369,20639,20693,20873,20963,21313,21661,23003,23087,23131

%N Primes such that the sum of the factorials of the digits is a perfect square.

%H Giovanni Resta, <a href="/A052279/b052279.txt">Table of n, a(n) for n = 1..10000</a>

%t Select[Prime[Range[3000]],IntegerQ[Sqrt[Total[IntegerDigits[#]!]]]&] (* _Harvey P. Dale_, Apr 11 2016 *)

%o (PARI) isok(n) = isprime(n) && (d = digits(n)) && issquare(sum(i=1, #d, d[i]!)); \\ _Michel Marcus_, Jan 06 2014

%Y Cf. A000040, A000142, A000290, A061602.

%K base,nonn

%O 1,1

%A _G. L. Honaker, Jr._, Feb 05 2000

%E More terms from _Michel Marcus_, Jan 06 2014