OFFSET
1,2
COMMENTS
The next term, if it exists, is > 10000. - Hugo Pfoertner, Jan 06 2020
The associated primes are: 13, 37, 881, 4651, 1273609, ...
From Robert Israel, Apr 28 2017: (Start)
All terms other than 0 are primes or powers of 2.
Heuristically, this sequence might be expected to be finite. (End)
LINKS
J. S. Gerasimov, x^(y + 1) - y^x, SeqFan list, Aug 18 2014.
EXAMPLE
4 is in this sequence because (1 + 4)^4 + (-4)^4 = 881 is prime.
MAPLE
N:= 1000: # to get all terms <= N
cands:= select(isprime, {seq(i, i=3..N, 2)}) union {0, seq(2^k, k=1..ilog2(N))}:
select(n -> isprime((1+n)^n + (-n)^n), cands); # Robert Israel, Apr 28 2017
PROG
(Magma) [n: n in [0..170]| IsPrime((n+1)^n + (-n)^n)];
(PARI) is(n)=ispseudoprime((n+1)^n+(-n)^n) \\ Charles R Greathouse IV, Apr 28 2017
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Juri-Stepan Gerasimov, Apr 27 2017
STATUS
approved