OFFSET
1,1
COMMENTS
Corresponding primes are 1047511, 4129751, 24285271, 2504843351, 126757680265156951, ... .
a(41) > 50000.
LINKS
Henri & Renaud Lifchitz, PRP Records. Search for n!3-59049.
Joe McLean, Interesting Sources of Probable Primes
EXAMPLE
19!3 - 3^10 = 19*16*13*10*7*4*1 - 59049 = 1047511 is prime, so 19 is in the sequence.
MATHEMATICA
MultiFactorial[n_, k_] := If[n < 1, 1, If[n < k + 1, n, n*MultiFactorial[n - k, k]]];
Select[Range[17, 50000], PrimeQ[MultiFactorial[#, 3] - 3^10] &]
PROG
(PARI) tf(n) = prod(i=0, (n-1)\3, n-3*i);
for(n=1, 1e4, if(ispseudoprime(tf(n) - 3^10), print1(n , ", "))) \\ Altug Alkan, Dec 04 2015
CROSSREFS
KEYWORD
hard,more,nonn
AUTHOR
Robert Price, Dec 04 2015
STATUS
approved