A353122 Numbers k such that k^k*(k+1) + 1 is prime.

0, 1, 2, 3, 6, 9, 186, 198, 8390

Offset

1,3

Comments

Corresponding primes start 2, 3, 13, 109, 326593, 3874204891, ...

a(9) > 6000. - Jon E. Schoenfield, Jun 05 2022

a(10) > 18000. - Michael S. Branicky, Aug 08 2024

Links

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

Example

9 is in the sequence because 9^9*(9+1) + 1 = 3874204891, which is prime.

Mathematica

Join[{0}, Select[Range[200], PrimeQ[#^#*(# + 1) + 1] &]] (* Amiram Eldar, Apr 25 2022 *)

Prog

(Magma) [n: n in [0..200] | IsPrime(n^n*(n+1) + 1)];
(PARI) isok(k) = ispseudoprime(k^k*(k+1) + 1); \\ Michel Marcus, May 16 2022

Crossrefs

Cf. A055897, A081132, A108318, A108879, A191513, A191568, A271718.

Sequence in context: A145761 A071714 A077753 * A058911 A249247 A276956

Adjacent sequences: A353119 A353120 A353121 * A353123 A353124 A353125

Keyword

nonn,hard,more

Author

Juri-Stepan Gerasimov, Apr 24 2022

Extensions

a(9) from Michael S. Branicky, Dec 22 2023

Status

approved