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”).

A191513
Numbers n such that n*(n-1)^n-1 is prime.
1
3, 5, 10, 11, 18, 127, 286, 560, 1025
OFFSET
1,1
EXAMPLE
a(1)=3 because 3*2^3-1=23 is prime, a(2)=5 because 5*4^5-1=5119 is prime, a(3)=10 because 10*9^10-1=34867844009 is prime, a(4)=11 because 11*10^11-1=1099999999999 is prime.
MAPLE
isA191513 := proc(n) n*(n-1)^n-1 ; isprime(%) ; end proc:
for n from 1 do if isA191513(n) then print(n); end if; end do: # R. J. Mathar, Jun 30 2011
PROG
(Magma) [n: n in [1..600]| IsPrime(n*(n-1)^n-1)]; // Vincenzo Librandi, Jun 13 2011
(PARI) for(n=1, 1e4, if(ispseudoprime(n*(n-1)^n-1), print1(n", "))) \\ Charles R Greathouse IV, Jun 14 2011
CROSSREFS
Sequence in context: A332049 A113858 A101130 * A254540 A326419 A102309
KEYWORD
nonn
AUTHOR
EXTENSIONS
Corrected and extended by Vincenzo Librandi and Charles R Greathouse IV, Jun 14 2011
STATUS
approved