login
This site is supported by donations to The OEIS Foundation.
Logo

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A160501 (n+1)^prime(n+1) + n^prime(n). 0
9, 251, 16627, 48844509, 13109522141, 232643574681223, 144347818589843079, 8863082234840576951801, 100000008862938119652501095929, 192043424957750480504146841291811 (list; graph; refs; listen; history; internal format)
OFFSET

1,1

COMMENTS

a(2)=251 is the only prime found for n up to 10000 using The C/Gmp program in the

link which is 17 times faster than the Pari routine.

Here there are divisibility rules: If prime(n) and prime(n+1) do

not differ by 6, then n^2+n+1 is a divisor. So finding primes in

this case will be difficult since 5/6 of the numbers are composite at the on set.

If another prime exists, it is larger than 418977 digits.

LINKS

C. L. Hilliard,Sum of prime indices to respective primes

FORMULA

a(n) = (n+1)^prime(n+1) + n^prime(n) = A062481(n)+A062481(n+1).

EXAMPLE

For n = 3, 4^7 + 3^5 = 16627, the 3rd entry in the sequence.

PROG

(PARI) ppower(n) = { for(x=1, n, y=(x+1)^prime(x+1) + x^prime(x); print1(y", ") ); }

CROSSREFS

Cf. A160491.

Sequence in context: A012072 A007408 A066989 * A075987 A135099 A073427

Adjacent sequences:  A160498 A160499 A160500 * A160502 A160503 A160504

KEYWORD

nonn

AUTHOR

Cino Hilliard (Hillcino368(AT)hotmail.com), May 15 2009

EXTENSIONS

Edited by R. J. Mathar (mathar(AT)strw.leidenuniv.nl), May 30 2009

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified February 15 19:15 EST 2012. Contains 205852 sequences.