OFFSET
1,1
COMMENTS
That is, numbers k such that A053557(k) is prime. - Michel Marcus, Aug 28 2013
MATHEMATICA
a[n_] := Sum[(-1)^k/k!, {k, 0, n}]; Select[Range[100], PrimeQ[Numerator[a[#]]] &] (* G. C. Greubel, Oct 31 2017 *)
PROG
(PARI) isok(n) = isprime(numerator(sum(k=0, n, (-1)^k/k!))); \\ Michel Marcus, Aug 28 2013
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Ed Pegg Jr, Jul 11 2008
EXTENSIONS
a(16)-a(20) from G. C. Greubel, Nov 01 2017
a(21)-a(22) from Michael S. Branicky, Dec 16 2024
STATUS
approved