# Strange NSolve failure [duplicate]

In:= NSolve[-LogIntegral + LogIntegral[x] == 119002, x]

During evaluation of In:= NSolve::ifun: Inverse functions are being used by NSolve,
so some solutions may not be found; use Reduce for complete solution information.

Out= {{x -> 1.56751*10^6}}


There are some natural numbers where this fails:

In:= NSolve[-LogIntegral + LogIntegral[x] == 119003, x]

During evaluation of In:= NSolve::nsmet: This system cannot be solved
with the methods available to NSolve.



What is so special about 119003?

In:= NSolve[-LogIntegral + LogIntegral[x] == 119004, x]

During evaluation of In:= NSolve::ifun: Inverse functions are being used by NSolve,
so some solutions may not be found; use Reduce for complete solution information.

Out= {{x -> 1.56754*10^6}}


Some other such integers are 128005, 133003, 137004. What could be the reason?

In case it is version dependent, mine is 11.0.1.0 (on Windows 10)

As Bob says, NSolve has an issue with using machine numbers, so why not use Solve instead? In the following, I also use the InverseFunctions option to avoid the Solve::ifun message, since I do want Solve to use inverse functions, and I include N to convert the Root object output of Solve into a real number:

Map[
{119003, 128005, 133003, 137004}
]


{{{x -> 1.56753*10^6}}, {{x -> 1.6963*10^6}}, {{x -> 1.76809*10^6}}, {{x -> 1.82572*10^6}}}

You need to use arbitrary-precision rather than machine precision.

With machine precision:

Quiet[NSolve[-LogIntegral + LogIntegral[x] == #, x] & /@ {119003, 128005,
133003, 137004}]

x], {{x -> 1.6963*10^6 - 1.62507*10^-15 I}}, {{x ->
1.76809*10^6 + 4.29213*10^-15 I}},

Quiet[NSolve[-LogIntegral + LogIntegral[x] == #, x,