|
|
A128466
|
|
Primes of the form ((k+1)^k - 1)/k^2 = A060073(k+1).
|
|
4
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Corresponding numbers k are listed in A127837.
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
|
|
LINKS
|
|
|
FORMULA
|
|
|
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
|
|
|
STATUS
|
approved
|
|
|
|