Hybrid integral equation/Monte Carlo approach to complexation thermodynamics