 Solution

Use the curvature $-1$ hyperbolic metric $ds=2|dz|/(1-|z{|}^2)$. For the disk Möbius transformation ${\varphi }_w(z)=(z-w)/(1-\overline{w}z)$, set $r_w(z)=|{\varphi }_w(z)|$. The hyperbolic distance satisfiesThe triangle inequality givesThus convergence of the exponential sum is independent of the base point. Taking the base point zero, it is equivalent to the Blaschke condition ${\sum }_n(1-|z_n|)<\infty$. This condition implies that only finitely many terms lie in any compact subdisk, including finitely many repetitions of any point.

To construct the required holomorphic function, for $a\ne 0$ define the Blaschke factorFor $|z|\le R<1$, writing $a=|a|e^{i\alpha }$ givesConsequently, after the finitely many factors corresponding to small $|a|$ are removed, the sum of $|1-b_a|$ converges uniformly on each compact subdisk. Those factors are uniformly close to one, so their logarithms form an absolutely convergent series there. Their product tail is holomorphic and nonvanishing. Including the finitely many removed factors gives the Blaschke productwith exactly the requested zero multiplicities. Each finite product has modulus at most one in the disk, so the limit does too. Hence $f=B/2$ maps into the open disk and has exactly the prescribed zeros, proving the implication from the summability condition.