A123438 - OEIS

%I #22 Dec 16 2024 11:34:38

%S 4,5,6,7,10,11,15,16,34,44,63,66,168,427,575,928,1094,1218,1363,1713,

%T 5278,10814

%N Numbers k such that the denominator of k!/!k (= A000142(k)/A000166(k)) is prime.

%C That is, numbers k such that A053557(k) is prime. - _Michel Marcus_, Aug 28 2013

%t a[n_] := Sum[(-1)^k/k!, {k, 0, n}]; Select[Range[100], PrimeQ[Numerator[a[#]]] &] (* _G. C. Greubel_, Oct 31 2017 *)

%o (PARI) isok(n) = isprime(numerator(sum(k=0, n, (-1)^k/k!))); \\ _Michel Marcus_, Aug 28 2013

%Y Cf. A000142, A000166, A053557.

%K nonn,more

%O 1,1

%A _Ed Pegg Jr_, Jul 11 2008

%E a(16)-a(20) from _G. C. Greubel_, Nov 01 2017

%E a(21)-a(22) from _Michael S. Branicky_, Dec 16 2024