|
|
A181681
|
|
Primes of the form n^(n+1)+(n-1)^n.
|
|
1
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
The next term of the sequence is a(5)=393^394+392^393 = 15526396580196065…8500849703002799761, which has 1023 decimal digits and is too large to include here.
|
|
LINKS
|
|
|
EXAMPLE
|
For n=5, 16649 is in the sequence because 5^6+4^5=16649 is prime.
|
|
MATHEMATICA
|
lst={}; Do[m=n+1; q=n^m+ (n-1)^n; If[PrimeQ[q], AppendTo[lst, q]], {n, 0, 400}];
lst
|
|
PROG
|
(PARI) A181681(istrt)={for(n=istrt, istrt+400,
if(isprime(n^(n+1)+(n-1)^n), print(n^(n+1)+(n-1)^n))); }
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|