A307452 Primes p such that the sum of the decimal digits of p^4 is also a prime.

2, 5, 7, 17, 23, 41, 47, 53, 67, 73, 97, 103, 113, 151, 157, 163, 173, 179, 197, 199, 223, 227, 251, 257, 263, 281, 293, 313, 349, 353, 389, 431, 439, 449, 457, 461, 479, 499, 503, 557, 577, 587, 593, 619, 659, 673, 709, 733, 829, 853, 857, 983, 997, 1033

Offset

1,1

Links

Table of n, a(n) for n=1..54.

Example

2^4=16, 1+6=7. 5^4=625, 6+2+5=13. 7^4=2401, 2+4+0+1=7.

Mathematica

Select[Prime@ Range@ 180, PrimeQ@ Total@ IntegerDigits[#^4] &] (* Michael De Vlieger, Jul 08 2019 *)

Prog

(PARI) isok(p) = isprime(p) && isprime(sumdigits(p^4)); \\ Michel Marcus, Apr 09 2019

Crossrefs

Cf. A007953 (sumdigits), A030514 (primes^4).

Sequence in context: A019084 A278487 A176676 * A066028 A066039 A142341

Adjacent sequences: A307449 A307450 A307451 * A307453 A307454 A307455

Keyword

nonn,base,easy

Author

Salvatore Di Guida, Apr 09 2019

Status

approved