A285888 Numbers n such that (1 + n)^n + (-n)^n is prime.

0, 2, 3, 4, 5, 7, 167

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

Table of n, a(n) for n=1..7.

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

Supersequence of A098463.

Cf. A243100, A285886, A285887.

Sequence in context: A259389 A277217 A259384 * A028986 A327324 A063948

Adjacent sequences: A285885 A285886 A285887 * A285889 A285890 A285891

Keyword

nonn,more

Author

Juri-Stepan Gerasimov, Apr 27 2017

Status

approved