Hydrodynamics of particle clusters suspended in viscous fluids is a subject of considerable theoretical and practical importance. Using a multipole expansion of the flow velocity in a series of spherical harmonics, Lamb's fundamental solution of the Stokes flow outside a single sphere is generalized in this work to the case of N nonoverlapping spheres of arbitrary size with slip boundary conditions. The expansion coefficients are found by transforming the boundary conditions to the Lamb form and by transforming the spherical coordinates and solid spherical harmonics centered at different spheres. The problem is reduced to the solution of the linear system of equations for the expansion coefficients, which is carried out numerically. Based on the developed theory, the relation between the hydrodynamic and gyration radius of fractal-like aggregates with different structure is established. In another application, an asymptotic slip-regime dependence of the aggregate hydrodynamic radius on the Knudsen number and the number of particles is found by performing calculations of drag forces acting on the gas-borne fractal-like and straight chain aggregates. A good agreement is shown in comparing predictions of the described theory with available experimental and theoretical results on motion of various small sphere clusters in viscous fluid.