|
%I #10 Apr 03 2023 10:36:13
%S 4,5,7,9,13,15,27,29,31,39,41,45,59,61,75,83,87,107,123,153,159,229,
%T 255,333,343,381,399,469,599,671,707,735,737,745,797,981,1115,1121,
%U 1233,1545,1625,1943,2177,2355,2593,2875,2891,6779,6937,7717,7767,9947,12285,13411,15367,15457,15699,17319,18035,20617,33125,37449,56353,61869,80599
%N Numbers k such that k![14]-2 is prime, where k![14] is the fourteen-fold multifactorial.
%C a(66) > 100000.
%C The first 12 primes associated with this sequence: 2, 3, 5, 7, 11, 13, 349, 433, 1579, 10723, 14389, 71143.
%C k = 1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 32, 51, 72, 95, 120.
%H C. Caldwell and H. Dubner (Eds): <a href="https://t5k.org/lists/top_ten/">The top ten prime numbers: from the unpublished collections of R. Ondrejka</a> (May 2001), Table 21 F, p. 75
%H Ken Davis, <a href="http://mfprimes.free-dc.org">Status of Search for Multifactorial Primes</a>.
%t MultiFactorial[n_, k_] := If[n < 1, 1, n*MultiFactorial[n - k, k]];
%t Select[Range[2, 100000], PrimeQ[MultiFactorial[#, 14] - 2] &]
%o (PARI) M(n,k) = if(n<1, 1, n*M(n - k, k));
%o isok(n) = isprime(M(n,14) - 2); \\ _Indranil Ghosh_, Mar 22 2017
%K nonn
%O 1,1
%A _Robert Price_, Mar 22 2017
|
