For the nonlinear Schrödinger equation (NLSE), in the presence of disorder, exponentially localized stationary states are found. We demonstrate analytically that the localization length is typically independent of the strength of the nonlinearity and is identical to the one found for the corresponding linear equation. The analysis makes use of the correspondence between the stationary NLSE and the Langevin equation as well as of the resulting Fokker-Planck equation. The calculations are performed for the "white noise" random potential, and an exact expression for the exponential growth of the eigenstates is obtained analytically. It is argued that the main conclusions are robust.