OFFSET
1,1
COMMENTS
Corresponding primes are 10139, 2026769, 13749310319, 213458046676619, 191898783962510369, ...
a(29) > 50000. - Robert Price, May 08 2015
LINKS
MATHEMATICA
Select[Range[8, 50000], PrimeQ[#!! - 256] &]
PROG
(PARI) is(n)=ispseudoprime(n!! - 2^8) \\ Anders Hellström, Dec 02 2015
(PARI) df(n) = if( n<0, 0, my(E); E = exp(x^2 / 2 + x * O(x^n)); n! * polcoeff( 1 + E * x * (1 + intformal(1 / E)), n));
for(n=1, 1e3, if(ispseudoprime(df(n) - 2^8), print1(n , ", "))) \\ Altug Alkan, Dec 03 2015
CROSSREFS
KEYWORD
hard,more,nonn
AUTHOR
Robert Price, Dec 01 2015
STATUS
approved