login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A091367 Primes p such that the sum of the digits raised to the 4th power is prime. 5

%I #14 Feb 24 2024 01:08:28

%S 11,23,29,41,43,47,61,67,83,89,101,113,131,179,191,197,223,269,311,

%T 313,331,353,379,397,401,443,461,601,607,641,661,719,739,809,883,911,

%U 937,971,1013,1019,1031,1033,1091,1097,1103,1109,1181,1301,1303,1367,1433

%N Primes p such that the sum of the digits raised to the 4th power is prime.

%H Michael S. Branicky, <a href="/A091367/b091367.txt">Table of n, a(n) for n = 1..10000</a>

%e a(1) = 11 because 1^4 + 1^4 = 2 which is prime.

%e a(10) = 89 because 8^4 + 9^4 = 10657 which is prime.

%t upto=500;Select[Prime[Range[upto]],PrimeQ[Total[IntegerDigits[#]^4]]&] (* _Paolo Xausa_, Nov 23 2021 *)

%o (Python)

%o from sympy import isprime, primerange

%o def ok(p): return isprime(sum(int(d)**4 for d in str(p)))

%o def aupto(limit): return [p for p in primerange(1, limit+1) if ok(p)]

%o print(aupto(1433)) # _Michael S. Branicky_, Nov 23 2021

%Y Cf. A046704 (primes whose digits sum to a prime), A052034 (primes whose digits squared sum to a prime), A091366 (primes whose digits cubed sum to a prime).

%K base,nonn

%O 1,1

%A Chuck Seggelin (barkeep(AT)plastereddragon.com), Jan 03 2004

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 06:34 EDT 2024. Contains 371265 sequences. (Running on oeis4.)