In this paper, the valuation of European and path-dependent options in foreign exchange (FX) markets is considered when the currency exchange rate evolves according to the Heston model combined with the Cox-Ingersoll-Ross dynamics for the stochastic domestic and foreign short interest rates. The mixed Monte Carlo/PDE method requires that we simulate only the paths of the squared volatility and the two interest rates, while an "inner" Black-Scholes-type expectation is evaluated by means of a PDE. This can lead to a substantial variance reduction and complexity improvements under certain circumstances depending on the contract and the model parameters. In this work, we establish the uniform boundedness of moments of the exchange rate process and its approximation, and prove strong convergence in L p (p ≥ 1) of the latter. Then, we carry out a variance reduction analysis and obtain accurate approximations for quantities of interest. All theoretical contributions can be extended to multi-factor short rates in a straightforward manner. Finally, we illustrate the efficiency of the method for the four-factor Heston-CIR model through a detailed quantitative assessment.