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”).
%I #19 Jun 11 2021 09:45:18
%S 2,7,311,7563707819165039903
%N Primes of the form ((k+1)^k - 1)/k^2 = A060073(k+1).
%C Corresponding numbers k are listed in A127837.
%C Terms are the primes in A060073.
%C 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
%F a(n) = ((A127837(n) + 1)^A127837(n) - 1) / A127837(n)^2.
%t Select[Table[((n+1)^n-1)/n^2,{n,500}],PrimeQ] (* _Harvey P. Dale_, Apr 30 2011 *)
%o (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
%Y Cf. A127837, A037205, A060072, A060073, A058128, A128456
%K nonn
%O 1,1
%A _Alexander Adamchuk_, Mar 09 2007