A301632 Numbers k such that 3^k + k + 1 is a prime.

0, 1, 3, 15, 45, 117, 159, 3739, 3777, 9703, 10623, 21459, 86595

Offset

1,3

Links

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

Maple

a:=n->`if`(isprime(3^k+k+1), k, NULL): seq(a(k), k=0..6000); # Muniru A Asiru, Mar 25 2018

Mathematica

Flatten[{0, Select[Range[5000], PrimeQ[3^# + # + 1] &]}] (* Vaclav Kotesovec, Mar 25 2018 *)

Prog

(PARI) for(n=0, 1000, if(isprime(3^n+n+1), print1(n", ")))

Crossrefs

A182327 gives the primes.

Cf. A171058, A173053, A301633.

Sequence in context: A050534 A048099 A030505 * A074355 A201868 A260021

Adjacent sequences: A301629 A301630 A301631 * A301633 A301634 A301635

Keyword

nonn,more,hard

Author

Seiichi Manyama, Mar 25 2018

Extensions

a(8)-a(9) from Muniru A Asiru, Mar 25 2018

a(10)-a(11) from Vaclav Kotesovec, Mar 25 2018

a(12) from Altug Alkan, Apr 01 2018

a(13) from Giovanni Resta, Apr 06 2018

Status

approved