|
%I #17 Sep 08 2022 08:44:56
%S 10,11,12,14,20,21,23,27,32,34,40,41,43,47,72,74,100,101,103,110,111,
%T 113,130,131,157,175,236,263,289,298,300,301,310,311,326,348,355,362,
%U 384,438,445,454,458,483,485,507,517,535,544,548,553,570,571,584,623
%N Sum of digits of a(n) raised to its digits powers is prime.
%C If p is prime, then R_p (see A002275) is a term, so the sequence is infinite. - _Marius A. Burtea_, Sep 18 2019
%H Marius A. Burtea, <a href="/A046431/b046431.txt">Table of n, a(n) for n = 0..5000</a>
%e 10 is a term: 10 -> 1^1 + 0^0 = 2 (prime).
%e 289 is a term: 289 -> 2^2 + 8^8 + 9^9 = 404197709 (prime).
%p q:= n-> (l-> isprime(add(i^i, i=l)))(convert(n, base, 10)):
%p select(q, [$0..1000])[]; # _Alois P. Heinz_, Sep 18 2019
%t sopQ[n_] := PrimeQ[Total[#^# & /@ (IntegerDigits[n] /. {0 -> 1})]]; Select[Range[625], sopQ[#] &] (* _Jayanta Basu_, Jun 25 2013 *)
%o (Magma) [k:k in [1..650]| IsPrime( &+[a[s]^(a[s]): s in [1..# a]]) where a is Intseq(k) ]; // _Marius A. Burtea_, Sep 18 2019
%Y Cf. A046432-A046441.
%K nonn,base
%O 0,1
%A _Patrick De Geest_, Jul 15 1998
|
