cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A210767 Numbers whose digit sum as well as sum of the 4th powers of the digits is a prime.

This page as a plain text file.
%I A210767 #19 Apr 01 2019 14:00:01
%S A210767 11,12,14,16,21,23,25,29,32,34,38,41,43,47,52,58,61,67,74,76,83,85,89,
%T A210767 92,98,101,102,104,106,110,111,113,119,120,131,133,140,146,160,164,
%U A210767 166,179,191,197,201,203,205,209,210,223,230,232,250,269,290,296,302
%N A210767 Numbers whose digit sum as well as sum of the 4th powers of the digits is a prime.
%C A210767 This is to the exponent 4 as A182404 is to the exponent 2.
%H A210767 Charles R Greathouse IV, <a href="/A210767/b210767.txt">Table of n, a(n) for n = 1..10000</a>
%F A210767 {n such that A055013(n) and A007953(n) are both primes}.
%e A210767 21 is in the sequence because sum of digits 2+1= 3 is prime, and sum of the 4th powers of the digits 2^4+1^4=17 is a prime.
%t A210767 Select[Range[350],AllTrue[{Total[IntegerDigits[#]],Total[ IntegerDigits[ #]^4]},PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Apr 01 2019 *)
%o A210767 (PARI) dspow(n,b,k)=my(s);while(n,s+=(n%b)^k;n\=b);s
%o A210767 select(n->isprime(sumdigits(n))&&isprime(dspow(n,10,4)), vector(10^3, i, i)) \\ _Charles R Greathouse IV_, May 11 2012
%Y A210767 Cf. A007953, A055013, A182404.
%K A210767 nonn,base,easy
%O A210767 1,1
%A A210767 _Jonathan Vos Post_, May 10 2012