Solution (source code)

= Solution

Write $Y=\mathbb P^1\setminus\{0,1,\infty\}$. One concrete uniformization uses the image in $\operatorname{PSL}_2(\mathbb R)$ of the level-two <principal congruence subgroup>
$$
\Gamma(2)=\left\{\begin{pmatrix}a&b\\c&d\end{pmatrix}\in\operatorname{SL}_2(\mathbb Z):
a,d\text{ odd},\ b,c\text{ even}\right\}/\{\pm I\}.
$$
It acts on the <complex upper half-plane> by <Möbius transformations>. Integer matrix entries make the group discrete, and the action is <properly discontinuous>. There are no nonidentity elliptic stabilizers: an elliptic element would have integer even trace with absolute value less than two, hence trace zero. But then $d=-a$ gives $bc=-(a^2+1)\equiv2\pmod4$, contradicting the evenness of $b,c$. Thus the quotient is an ordinary <Riemann surface>, with no orbifold points.

A fundamental region is the ideal quadrilateral
$$
\mathcal F=\{\tau\in\mathbb H:-1\leq\operatorname{Re}\tau\leq1,
\ |\tau-\tfrac12|\geq\tfrac12,
\ |\tau+\tfrac12|\geq\tfrac12\}.
$$
Its vertical sides are identified by $T(\tau)=\tau+2$, and its semicircular sides by $U(\tau)=\tau/(1-2\tau)$. The ideal vertices are $-1,0,1,\infty$; $-1$ and $1$ represent one puncture, while zero and infinity give two others. Filling these three punctures gives genus zero: the identified quadrilateral has three vertices, two paired edges and one face, so its compactification has Euler characteristic $3-2+1=2$. Its complex compactification is a sphere, and a <Möbius transformation> puts the punctures at $0,1,\infty$. Hence
$$
\boxed{Y\cong\mathbb H/\Gamma(2).}
$$
This is the <thrice-punctured sphere as a modular quotient>. Equivalently, the <uniformization theorem> applied to its <universal covering map> produces the upper-half-plane model: the sphere cannot cover this noncompact surface, and an entire covering from $\mathbb C$ would contradict <Little Picard theorem> by omitting both zero and one.

The curvature-minus-one metric $ds_{\mathbb H}=|d\tau|/\operatorname{Im}\tau$ is invariant under every real <Möbius transformation> preserving $\mathbb H$. Indeed, for determinant-one coefficients,
$$
\operatorname{Im}\frac{a\tau+b}{c\tau+d}
=\frac{\operatorname{Im}\tau}{|c\tau+d|^2},\qquad
\left|\frac{d}{d\tau}\frac{a\tau+b}{c\tau+d}\right|
=\frac1{|c\tau+d|^2}.
$$
Thus it descends to the <hyperbolic metric> on $Y$. This defines its distance by the infimum of path lengths; a choice of local lift does not affect lengths.

For the <Schwarz-Pick lemma>, first prove the <Schwarz lemma>. If $h:\mathbb D\to\mathbb D$ is a <holomorphic function> and $h(0)=0$, the function $h(z)/z$ extends as a <holomorphic function> across zero. On $|z|=r$ its modulus is at most $1/r$. The <maximum modulus principle> and $r\uparrow1$ yield $|h(z)|\leq|z|$ and $|h'(0)|\leq1$. Equality at a nonzero point or in the derivative makes $h(z)/z$ a constant of modulus one, again by the <maximum modulus principle>, so $h$ is a rotation.

For a general <holomorphic function> $F:\mathbb D\to\mathbb D$, conjugate by the <automorphisms of the unit disk> $A_a(z)=(z-a)/(1-\overline a z)$, putting $a$ and $F(a)$ at zero. Applying the proved <Schwarz lemma> gives
$$
\boxed{\left|\frac{F(z)-F(a)}{1-\overline{F(a)}F(z)}\right|
\leq\left|\frac{z-a}{1-\overline a z}\right|},
\qquad
\boxed{\frac{|F'(a)|}{1-|F(a)|^2}\leq\frac1{1-|a|^2}}.
$$
The disk <hyperbolic metric> is $2|dz|/(1-|z|^2)$, and its distance is twice the inverse hyperbolic tangent of the displayed pseudohyperbolic distance. Both forms of <Schwarz-Pick lemma> therefore express metric contraction. Equality in the derivative at one point makes $F$ a disk automorphism by the equality case already proved.

For $X=\mathbb D\setminus\{0\}$ and a <holomorphic function> $f:X\to Y$, compose with a <universal covering map> of $X$ and lift to a <universal covering map> of $Y$. The simply connected source of the lift is a disk or half-plane. <Schwarz-Pick lemma> there bounds the lifted metric derivative, and the covering maps are local <isometries>, so
$$
f^*ds_Y\leq ds_X,\qquad
\boxed{d_Y(f(z_0),f(z_1))\leq d_X(z_0,z_1)}.
$$
The second inequality follows by applying the length inequality to every joining path and taking the infimum. The <hyperbolic metric on the punctured disk> is its own intrinsic metric,
$$
ds_X=\frac{|dz|}{|z|\log(1/|z|)},
$$
obtained from the covering $z=e^{i\tau}$. It is not merely the ordinary disk metric restricted to the punctured disk.

\b[Yes, equality can occur for distinct points.] The cyclic subgroup $\langle T\rangle$ lies in $\Gamma(2)$. Its quotient is conformally the punctured disk, using $z=e^{\pi i\tau}$, because this identifies precisely $\tau\sim\tau+2$. The subgroup inclusion induces a <covering map> that is a <holomorphic function>
$$
\mathbb D^*\cong\mathbb H/\langle T\rangle
\longrightarrow\mathbb H/\Gamma(2)\cong Y.
$$
It is a local <isometry>. Choose a sufficiently small geodesically convex ball around any target point that is evenly covered. Two distinct sufficiently close points in one lift of that ball are joined by the lift of its minimizing target <geodesic>. The lift has the same length, giving $d_X\leq d_Y$ for this pair, while contraction gives the reverse inequality. Thus equality holds.

For context, if the lifted map is not an automorphism, the derivative inequality is strict at every point; integrating along a minimizing <geodesic> in the complete punctured-disk metric gives a strict distance inequality for every distinct pair. Equality forces the covering-map case, but even a <covering map> need not preserve all distances globally because shorter paths can pass through other sheets. This is the distinction captured by <equality in hyperbolic Schwarz-Pick contraction>.