cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

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

This page as a plain text file.
%I A353122 #34 Aug 08 2024 09:50:05
%S A353122 0,1,2,3,6,9,186,198,8390
%N A353122 Numbers k such that k^k*(k+1) + 1 is prime.
%C A353122 Corresponding primes start 2, 3, 13, 109, 326593, 3874204891, ...
%C A353122 a(9) > 6000. - _Jon E. Schoenfield_, Jun 05 2022
%C A353122 a(10) > 18000. - _Michael S. Branicky_, Aug 08 2024
%e A353122 9 is in the sequence because 9^9*(9+1) + 1 = 3874204891, which is prime.
%t A353122 Join[{0}, Select[Range[200], PrimeQ[#^#*(# + 1) + 1] &]] (* _Amiram Eldar_, Apr 25 2022 *)
%o A353122 (Magma) [n: n in [0..200] | IsPrime(n^n*(n+1) + 1)];
%o A353122 (PARI) isok(k) = ispseudoprime(k^k*(k+1) + 1); \\ _Michel Marcus_, May 16 2022
%Y A353122 Cf. A055897, A081132, A108318, A108879, A191513, A191568, A271718.
%K A353122 nonn,hard,more
%O A353122 1,3
%A A353122 _Juri-Stepan Gerasimov_, Apr 24 2022
%E A353122 a(9) from _Michael S. Branicky_, Dec 22 2023