This paper discusses the effects of cloud turbulence, turbulent entrainment, and entrained cloud condensation nuclei (CCN) activation on the evolution of the cloud droplet size spectrum. We simulate an ensemble of idealized turbulent cloud parcels that are subject to entrainment events modeled as a random process. Entrainment events, subsequent turbulent mixing inside the parcel, supersaturation fluctuations, and the resulting stochastic droplet activation and growth by condensation are simulated using a Monte Carlo scheme. Quantities characterizing the turbulence intensity, entrainment rate, CCN concentration, and the mean fraction of environmental air entrained in an event are all specified as independent external parameters. Cloud microphysics is described by applying Lagrangian particles, the so-called superdroplets. These are either unactivated CCN or cloud droplets that grow from activated CCN. The model accounts for the addition of environmental CCN into the cloud by entraining eddies at the cloud edge. Turbulent mixing of the entrained dry air with cloudy air is described using the classical linear relaxation to the mean model. We show that turbulence plays an important role in aiding entrained CCN to activate, and thus broadening the droplet size distribution. These findings are consistent with previous large-eddy simulations (LESs) that consider the impact of variable droplet growth histories on the droplet size spectra in small cumuli. The scheme developed in this work is ready to be used as a stochastic subgrid-scale scheme in LESs of natural clouds.