OFFSET
1,1
COMMENTS
Corresponding primes are 6563, 6571, 6841, 2504908961, 17961239302561, 81359229958561, 664565853958561, ...
Terms > 68 correspond to probable primes.
a(36) > 50000.
LINKS
Henri & Renaud Lifchitz, PRP Records. Search for n!3+3^8
Joe McLean, Interesting Sources of Probable Primes
OpenPFGW Project, Primality Tester
EXAMPLE
10!3 + 3^4 = 10*7*4*1 + 6561 = 6841 is prime, so 10 is in the sequence.
MATHEMATICA
MultiFactorial[n_, k_] := If[n < 1, 1, If[n < k + 1, n, n*MultiFactorial[n - k, k]]];
Select[Range[0, 50000], PrimeQ[MultiFactorial[#, 3] + 3^8] &]
Select[Range[800], PrimeQ[6561+Times@@Range[#, 1, -3]]&] (* Harvey P. Dale, Mar 08 2023 *)
PROG
(PARI) is(n)=ispseudoprime(n!!! + 3^8) \\ Anders Hellström, Nov 27 2015
(PARI) tf(n) = prod(i=0, (n-1)\3, n-3*i);
for(n=1, 1e4, if(ispseudoprime(tf(n) + 3^8), print1(n , ", "))) \\ Altug Alkan, Dec 03 2015
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Robert Price, Nov 26 2015
STATUS
approved