%I #30 Jul 05 2024 11:49:21
%S 0,1,4,7,8,9,13,16,28,54,86,129,190,351,466,697,938,1510,2748,2878,
%T 3396,4057,4384,5534,7069,10364
%N Numbers k such that 4*k! + 1 is prime.
%C a(25) > 6311. - _Jinyuan Wang_, Feb 06 2020
%H Factor Database, <a href="http://factordb.com/index.php?query=4*z%21%2B1&perpage=200">Factors of the numbers 4z!+1</a>
%e k = 7 is a term because 4*7! + 1 = 20161 is prime.
%t Select[Range[5000],PrimeQ[4#!+1]&] (* _Harvey P. Dale_, Mar 23 2011 *)
%o (PARI) is(k) = ispseudoprime(4*k!+1); \\ _Jinyuan Wang_, Feb 06 2020
%Y Cf. m*n!-1 is a prime: A076133, A076134, A099350, A099351, A180627-A180631.
%Y Cf. m*n!+1 is a prime: A051915, A076679-A076683, A178488, A180626, A126896.
%Y Cf. n!/m-1 is a prime: A002982, A082671, A139056, A139199-A139205.
%Y Cf. n!/m+1 is a prime: A002981, A082672, A089085, A139061, A139058, A139063, A139065, A151913, A137390, A139071.
%K nonn,more
%O 1,3
%A Phillip L. Poplin (plpoplin(AT)bellsouth.net), Oct 25 2002
%E Corrected (added missed terms 2748, 2878) by _Serge Batalov_, Feb 24 2015
%E a(24) from _Jinyuan Wang_, Feb 06 2020
%E a(25)-a(26) from _Michael S. Branicky_, Jul 04 2024