OFFSET
1,1
COMMENTS
Corresponding indices are 2, 3, 5, 7, 9, 13, ... with no more up to 12000.
MATHEMATICA
a[0] = 1; a[n_] := a[n] = Sum[Binomial[n, k]*a[n - k], {k, 1, n}]; Select[a@# & /@ Range@ 20, PrimeQ@# &]
PROG
(PARI) lista(nn) = for (n=1, nn, if (isprime(fn=sum(k=1, n, k!*stirling(n, k, 2))), print1(fn, ", ")); ); \\ Michel Marcus, Aug 02 2017
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Jul 28 2017
STATUS
approved