Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A160501
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A160501 (n+1)^prime(n+1) + n^prime(n). +0
1
9, 251, 16627, 48844509, 13109522141, 232643574681223, 144347818589843079, 8863082234840576951801, 100000008862938119652501095929, 192043424957750480504146841291811 (list; graph; listen)
OFFSET

1,1

COMMENT

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.

PROGRAM

(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 this_sequence A075987 A135099 A073427

Adjacent sequences: A160498 A160499 A160500 this_sequence 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

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified December 6 22:55 EST 2009. Contains 170429 sequences.


AT&T Labs Research