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

8, 81, 15625, 5764801, 3138428376721, 3937376385699289, 14063084452067724991009, 37589973457545958193355601, 480250763996501976790165756943041

Offset: 1

Author

Cino Hilliard, Mar 06 2005

Keywords

• easy
• nonn

Comments

Sum of reciprocals rapidly converges to 0.1374098524791901212366977116..

A182938(a(n)) = 0. [Reinhard Zumkeller, Feb 18 2012]

Links

• Reinhard Zumkeller, Table of n, a(n) for n = 1..75

Crossrefs

Cf. A051674, A000040.

Programs

Haskell

a104126 n = p ^ (p + 1) where p = a000040 n
-- Reinhard Zumkeller, Feb 18 2012

Mathematica

#^(#+1)&/@Prime[Range[10]] (* Harvey P. Dale, Dec 12 2021 *)

PARI

ptopp1(n) = { local(x,y,z,sr=0); forprime(x=1,n, y=x^(x+1); z=(x+1)^x; sr+=1./y; print1(y","); ); print(); print(sr) }

Extensions

Offset corrected by Reinhard Zumkeller, Feb 18 2012