login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A128466
Primes of the form ((k+1)^k - 1)/k^2 = A060073(k+1).
4
2, 7, 311, 7563707819165039903
OFFSET
1,1
COMMENTS
Corresponding numbers k are listed in A127837.
Terms are the primes in A060073.
Next term has 15850 = 1 + floor((4357*log(4358) - 2*log(4357))/log(10)) digits and is too large to include. - M. F. Hasler, May 22 2007
FORMULA
a(n) = ((A127837(n) + 1)^A127837(n) - 1) / A127837(n)^2.
MATHEMATICA
Select[Table[((n+1)^n-1)/n^2, {n, 500}], PrimeQ] (* Harvey P. Dale, Apr 30 2011 *)
PROG
(PARI) A128466(n)=A060073(A127837(n)+1) /* see there. --- or: */ forprime(p=1, 10^5, if(ispseudoprime(n=((p+1)^p-1)/p^2), print1(n, ", "))); \\ M. F. Hasler, May 22 2007
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Mar 09 2007
STATUS
approved