OFFSET
1,1
MATHEMATICA
li[n_] := -1/2* 1/ProductLog[-Exp[1 - 1/((2 n - 2))]/(2n-2)]; Table[If[! Round[n*E^(-1 - 1/(2*n))] == Ceiling[li[n]], Print[n]; n], {n, 4, 1000000}]
CROSSREFS
KEYWORD
nonn,hard
AUTHOR
José María Grau Ribas, May 25 2013
STATUS
approved