We study the Cauchy problem for the integrable nonlocal focusing nonlinear Schrödingeris the Heaviside step function, and A > 0 and R > 0 are arbitrary constants. Our main aim is to study the large-t behavior of the solution of this problem. We show that for R ∈ (2n−1)π 2A , (2n+1)π
2A, n = 1, 2, . . . , the (x, t) plane splits into 4n+2 sectors exhibiting different asymptotic behavior. Namely, there are 2n+1 sectors where the solution decays to 0, whereas in the other 2n + 1 sectors (alternating with the sectors with decay), the solution approaches (different) constants along each ray x/t = const. Our main technical tool is the representation of the solution of the Cauchy problem in terms of the solution of an associated matrix Riemann-Hilbert problem and its subsequent asymptotic analysis following the ideas of nonlinear steepest descent method.