We discuss the application of the direct Monte Carlo method to the theory of cluster formation. Fractal relationships for the kernels appearing in the Smoluchowski equation are implemented in this method and the scaling behavior of the kernels is investigated using computer simulation. We study the effects of cluster disintegrations and also investigate the effects of "magic" numbers in cluster formation.