A Markovian Monte Carlo algorithm for multi-parton production in the high-energy limit is proposed and the matching with unintegrated parton densities is discussed.