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.
%I A087147 #26 Feb 16 2025 08:32:51
%S A087147 0,3,7,9,67,291,1343,6984,12861
%N A087147 Numbers k such that k! + (k+1)! + 1 is prime.
%C A087147 291 is in the sequence and also is in the sequence A087146, thus (291!+292!-1,291!+292!+1) is a twin pair of primes. Any additional terms are greater than 1800 with the next prime having more than 5086 digits.
%C A087147 Next term is greater than 4200. - Gabriel Cunningham (gcasey(AT)mit.edu), Sep 09 2003
%C A087147 a(10) > 25000. - _Robert Price_, Aug 26 2015
%C A087147 k+1 is not prime because (p-1)! + p! + 1 == 0 mod p and (p-1)! + p! + 1 > p where p is prime. - _Seiichi Manyama_, Mar 22 2018
%D A087147 H. Dubner, Factorial and primorial primes, J. Rec. Math., 19 (No.3, 1987)
%H A087147 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/WilsonsTheorem.html">Wilson's Theorem</a>
%e A087147 3 is in the sequence because 3!+4!+1=31 is prime.
%t A087147 v={}; Do[If[PrimeQ[n!+(n+1)!+1], v=Append[v, n]; Print[v]], {n, 1800}]; v
%t A087147 Select[Range[0,25000],PrimeQ[#!+(#+1)!+1]&] (* _Robert Price_, Aug 26 2015 *)
%o A087147 (PARI) isok(k) = ispseudoprime(k!+(k+1)!+1); \\ _Altug Alkan_, Mar 22 2018
%Y A087147 Cf. A087146, A051856.
%Y A087147 Primes in A118913. [From _Dmitry Kamenetsky_, Oct 21 2008]
%K A087147 more,nonn
%O A087147 1,2
%A A087147 _Farideh Firoozbakht_, Aug 19 2003
%E A087147 a(8)-a(9) from _Robert Price_, Aug 26 2015