A hybrid SPH/N-body method for star cluster simulations