In this paper we derived in QCD the Balitsky-Fadin-Kuraev-Lipatov (BFKL) linear, inhomogeneous equation for the factorial moments of multiplicity distribution () from Le-Mueller-Munier equation. In particular, the equation for the average multiplicity of the color-singlet dipoles () turns out to be the homogeneous BFKL while at small . Second, using the diffusion approximation for the BFKL kernel we show that the factorial moments are equal to which leads to the multiplicity distribution . We also suggest a procedure for finding corrections to this multiplicity distribution which will be useful for descriptions of the experimental data.