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
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
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Mar 09 2007
STATUS
approved