OFFSET
1,1
COMMENTS
For n >= 1, the corresponding primes are A119029(a(n)) = {2, 1676770323947695709, ...}.
a(15) > 10000. - Lucas A. Brown, Apr 01 2021
LINKS
Lucas A. Brown, A129977.py
MATHEMATICA
Do[ f=Numerator[ Sum[ n^(k-1)/k!, {k, 1, n} ] ]; If[ PrimeQ[f], Print[{n, f}] ], {n, 1, 316} ]
Select[Range[2000], PrimeQ[Numerator[Sum[#^(k-1)/k!, {k, #}]]]&] (* Harvey P. Dale, Jun 15 2019 *)
PROG
(PARI) for( n=1, 1000, if( ispseudoprime( numerator( sum( k=1, n, n^(k-1)/k!))), print1(n", "))) \\ M. F. Hasler, Jun 18 2007
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Alexander Adamchuk, Jun 13 2007
EXTENSIONS
Edited and extended (a(6)..a(8)) by M. F. Hasler, Jun 18 2007
More terms from Ryan Propper, Jan 12 2008
Various sections edited by Petros Hadjicostas, May 12 2020
a(12)-a(14) from Lucas A. Brown, Apr 01 2021
STATUS
approved