In this paper we propose recurrence relations for the dipole densities in QCD, which allows us to find these densities from the solution to the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation. We resolve these relations in the diffusion approximation for the BFKL kernel. Based on this solution, we found the sum of large Pomeron loops. This sum generates the scattering amplitude that decreases at large values of rapidity Y. It turns out that such behavior of the scattering amplitudes is an artifact of diffusion approximation. This approximation leads to the unitarization without saturation both in deep inelastic scattering and in dipole-dipole interaction at high energies.