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!)
A072185 Primes p = x1x2x3...xn in base 10 such that x1^p1 + x2^p2 + x3^p3 + ... + xn^pn is a prime, where pn is the n-th prime. 0
11, 23, 29, 41, 43, 47, 61, 67, 101, 131, 139, 151, 173, 179, 199, 263, 311, 331, 337, 359, 401, 461, 593, 601, 607, 641, 757, 809, 821, 827, 863, 881, 887, 911, 953, 977, 991, 1019, 1091, 1097, 1109, 1123, 1181, 1217, 1301, 1321, 1381, 1451, 1499, 1583 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

LINKS

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

EXAMPLE

23 belongs to the sequence because 2^2 + 3^3 = 31 is a prime.

MAPLE

From R. J. Mathar, Oct 10 2010: (Start)

isA072185 := proc(n) local d ; if isprime(n) then d := convert(n, base, 10) ; sdg := add( op(-i, d)^ithprime(i), i=1..nops(d)) ; isprime(sdg) ; else false; end if; end proc:

for p from 2 to 2000 do if isA072185(p) then printf("%d, ", p) ; fi; end do: (End)

PROG

(UBASIC) 30 for I=10 to 100000 40 L=alen(I) 50 Q=2:A=0 60 A1=val(mid(str(I), 2, 1)) 70 for H=2 to L 80 A=A+val(mid(str(I), H+1, 1))^nxtprm(Q) 90 Q=nxtprm(Q) 100 next 110 P=A+A1^2 120 if prmdiv(P)=P and P<>1 and prmdiv(I)=I and I<>1 then print I 130 next

CROSSREFS

Sequence in context: A164932 A086244 A091939 * A105898 A247228 A243461

Adjacent sequences:  A072182 A072183 A072184 * A072186 A072187 A072188

KEYWORD

easy,base,nonn

AUTHOR

Felice Russo, Jul 01 2002

EXTENSIONS

More terms from R. J. Mathar, Oct 10 2010

STATUS

approved

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 October 1 09:13 EDT 2022. Contains 357140 sequences. (Running on oeis4.)