[proofplan] We realize each point $z \in \mathbb{H}$ as the elliptic curve $\mathbb{C}/(\mathbb{Z}z+\mathbb{Z})$ with the marked $N$-torsion point represented by $1/N$. The key computation is that a change of lattice basis by a matrix in $SL_2(\mathbb{Z})$ preserves this marked point exactly when the matrix lies in $\Gamma_1(N)$. Surjectivity follows by choosing a lattice basis adapted to a given torsion point of exact order $N$, and injectivity follows by lifting an isomorphism of complex tori to multiplication by a nonzero complex number and reading off the induced basis change. [/proofplan]

[step:Verify that the standard marked point has exact order $N$] Fix $z \in \mathbb{H}$. The lattice $\Lambda_z := \mathbb{Z}z+\mathbb{Z}$ is a rank-two discrete subgroup of $\mathbb{C}$, and $E_z := \mathbb{C}/\Lambda_z$ is the associated complex elliptic curve. Define the marked point $P_z \in E_z$ by \begin{align*} P_z := \frac{1}{N}+\Lambda_z. \end{align*} Since \begin{align*} N P_z = 1+\Lambda_z = 0 \in E_z, \end{align*} we have $P_z \in E_z[N]$. If $1 \leq k \leq N$ and $kP_z = 0$, then $k/N \in \Lambda_z$. Thus there exist $m,n \in \mathbb{Z}$ such that \begin{align*} \frac{k}{N} = mz+n. \end{align*} Because $\operatorname{Im}(z)>0$, the complex number $z$ is not real, so equality with the real number $k/N$ forces $m=0$. Hence $k/N=n \in \mathbb{Z}$, and therefore $N \mid k$. For $1 \leq k \leq N$, this gives $k=N$. Thus $P_z$ has exact order $N$. [/step]

[step:Compute precisely when a modular change of basis preserves the marked point] Let \begin{align*} \gamma = \begin{pmatrix} a & b \\ c & d \end{pmatrix} \in SL_2(\mathbb{Z}) \end{align*} and define $\gamma z \in \mathbb{H}$ by \begin{align*} \gamma z := \frac{az+b}{cz+d}. \end{align*} Since $\operatorname{Im}(z)>0$ and $\gamma \in SL_2(\mathbb{Z})$, one has $\gamma z \in \mathbb{H}$. Define the complex-[linear map](/page/Linear%20Map) \begin{align*} \varphi_\gamma: \mathbb{C} &\to \mathbb{C} \\ u &\mapsto (cz+d)u. \end{align*} The identities \begin{align*} (cz+d)(\gamma z) &= az+b, \\ (cz+d)\cdot 1 &= cz+d \end{align*} show that \begin{align*} (cz+d)\Lambda_{\gamma z} &= \mathbb{Z}(az+b)+\mathbb{Z}(cz+d) \\ &= \mathbb{Z}z+\mathbb{Z} \\ &= \Lambda_z, \end{align*} because the integer matrix with columns $(a,b)$ and $(c,d)$ has determinant $ad-bc=1$. Hence $\varphi_\gamma$ descends to an isomorphism of complex elliptic curves \begin{align*} \overline{\varphi}_\gamma: E_{\gamma z} &\to E_z \\ u+\Lambda_{\gamma z} &\mapsto (cz+d)u+\Lambda_z. \end{align*} The image of the standard marked point is \begin{align*} \overline{\varphi}_\gamma(P_{\gamma z}) &= \frac{cz+d}{N}+\Lambda_z \\ &= \frac{c}{N}z+\frac{d}{N}+\Lambda_z. \end{align*} Therefore $\overline{\varphi}_\gamma(P_{\gamma z})=P_z$ if and only if \begin{align*} \frac{c}{N}z+\frac{d-1}{N}\in \Lambda_z. \end{align*} Since $z$ and $1$ are $\mathbb{R}$-linearly independent, this condition is equivalent to \begin{align*} c \equiv 0 \pmod N, \qquad d \equiv 1 \pmod N. \end{align*} Together with $ad-bc=1$, these congruences imply $a \equiv 1 \pmod N$. Thus $\overline{\varphi}_\gamma(P_{\gamma z})=P_z$ exactly for $\gamma \in \Gamma_1(N)$. It follows that if $z$ and $z'$ lie in the same $\Gamma_1(N)$-orbit, then $(E_z,P_z)$ and $(E_{z'},P_{z'})$ are isomorphic. Hence the displayed map from $\Gamma_1(N)\backslash \mathbb{H}$ to isomorphism classes is well-defined. [/step]

[step:Choose a lattice basis adapted to any point of exact order $N$] Let $(E,P)$ be a pair consisting of a complex elliptic curve $E$ and a point $P \in E[N]$ of exact order $N$. By the definition of a complex elliptic curve as a one-dimensional complex torus, there is a lattice $\Lambda \subset \mathbb{C}$ such that $E \cong \mathbb{C}/\Lambda$. Choose a $\mathbb{Z}$-basis $(\omega_1,\omega_2)$ of $\Lambda$. Since $P \in E[N]$, there exists $\lambda \in \Lambda$ such that \begin{align*} P = \frac{\lambda}{N}+\Lambda. \end{align*} Write $\lambda = r\omega_1+s\omega_2$ with $r,s \in \mathbb{Z}$. The exact order of $P$ is $N$, so the class of $(r,s)$ in $(\mathbb{Z}/N\mathbb{Z})^2$ has exact order $N$. Equivalently, \begin{align*} \gcd(r,s,N)=1. \end{align*} Choose integers $u,v$ such that \begin{align*} u \equiv r \pmod N, \qquad v \equiv s \pmod N, \qquad \gcd(u,v)=1. \end{align*} We justify the choice explicitly. If $r=0$, then $\gcd(s,N)=1$, so we may take $u=N$ and $v=s$. If $r \neq 0$, take $u=r$ and choose an integer $t$ so that $v:=s+tN$ is not divisible by any prime divisor of $r$. For a prime $p \mid r$ with $p \mid N$, the condition $\gcd(r,s,N)=1$ gives $p \nmid s$, hence $p \nmid s+tN$ for every $t$. For a prime $p \mid r$ with $p \nmid N$, exactly one residue class of $t$ modulo $p$ makes $s+tN \equiv 0 \pmod p$; avoiding these finitely many forbidden residue classes gives $p \nmid v$ for every prime $p \mid r$. Therefore $\gcd(u,v)=\gcd(r,v)=1$, and the congruences modulo $N$ hold by construction. Define \begin{align*} \lambda_2 := u\omega_1+v\omega_2 \in \Lambda. \end{align*} Then $\lambda_2 \equiv \lambda \pmod{N\Lambda}$, so \begin{align*} P = \frac{\lambda_2}{N}+\Lambda. \end{align*} Because $\gcd(u,v)=1$, Bezout's identity gives integers $x,y \in \mathbb{Z}$ such that \begin{align*} xv-yu=1. \end{align*} Define \begin{align*} \lambda_1 := x\omega_1+y\omega_2 \in \Lambda. \end{align*} The change-of-basis matrix from $(\omega_1,\omega_2)$ to $(\lambda_1,\lambda_2)$ has determinant $xv-yu=1$, so $(\lambda_1,\lambda_2)$ is a $\mathbb{Z}$-basis of $\Lambda$. If necessary, replace $\lambda_1$ by $-\lambda_1$; this keeps $(\lambda_1,\lambda_2)$ a lattice basis and makes \begin{align*} z := \frac{\lambda_1}{\lambda_2} \end{align*} belong to $\mathbb{H}$. Multiplication by $\lambda_2^{-1}$ gives an isomorphism \begin{align*} \psi: \mathbb{C}/\Lambda &\to \mathbb{C}/\Lambda_z \\ u+\Lambda &\mapsto \frac{u}{\lambda_2}+\Lambda_z, \end{align*} and under this isomorphism \begin{align*} \psi(P) = \frac{1}{N}+\Lambda_z = P_z. \end{align*} Thus every isomorphism class of pairs $(E,P)$ occurs in the image. [/step]

[step:Show that isomorphic marked elliptic curves give the same $\Gamma_1(N)$-orbit] Suppose $z,w \in \mathbb{H}$ and suppose there is an isomorphism of marked elliptic curves \begin{align*} \Psi: E_z &\to E_w \end{align*} such that $\Psi(P_z)=P_w$. We use the standard lifting property for homomorphisms of complex tori: a holomorphic group homomorphism $\mathbb{C}/\Lambda_z \to \mathbb{C}/\Lambda_w$ lifts through the universal covering map $\mathbb{C} \to \mathbb{C}/\Lambda_z$ to a complex-linear map $\mathbb{C} \to \mathbb{C}$ carrying $\Lambda_z$ into $\Lambda_w$. Since $\Psi$ is an isomorphism, the lifted linear map is multiplication by some $\alpha \in \mathbb{C}^{\times}$ and carries $\Lambda_z$ bijectively onto $\Lambda_w$. Thus $\Psi$ lifts to \begin{align*} \widetilde{\Psi}: \mathbb{C} &\to \mathbb{C} \\ u &\mapsto \alpha u \end{align*} with $\alpha\Lambda_z=\Lambda_w$. Because $\alpha z$ and $\alpha$ form a $\mathbb{Z}$-basis of $\Lambda_w=\mathbb{Z}w+\mathbb{Z}$, there are integers $a,b,c,d \in \mathbb{Z}$ such that \begin{align*} \alpha z &= aw+b, \\ \alpha &= cw+d. \end{align*} Multiplication by the nonzero complex number $\alpha$ preserves orientation on the real plane $\mathbb{C}$, and the ordered bases $(z,1)$ and $(w,1)$ have the same orientation because $\operatorname{Im}(z)>0$ and $\operatorname{Im}(w)>0$. Therefore \begin{align*} ad-bc=1, \end{align*} so \begin{align*} \gamma := \begin{pmatrix} a & b \\ c & d \end{pmatrix} \in SL_2(\mathbb{Z}). \end{align*} Dividing the identity $\alpha z=aw+b$ by $\alpha=cw+d$ gives \begin{align*} z = \frac{aw+b}{cw+d} = \gamma w. \end{align*} The condition $\Psi(P_z)=P_w$ means \begin{align*} \frac{\alpha}{N}+\Lambda_w = \frac{1}{N}+\Lambda_w. \end{align*} Equivalently, \begin{align*} \frac{\alpha-1}{N}\in \Lambda_w. \end{align*} Since $\alpha=cw+d$, this is \begin{align*} \frac{c}{N}w+\frac{d-1}{N}\in \mathbb{Z}w+\mathbb{Z}. \end{align*} The numbers $w$ and $1$ are $\mathbb{R}$-linearly independent, so \begin{align*} c \equiv 0 \pmod N, \qquad d \equiv 1 \pmod N. \end{align*} Using $ad-bc=1$, we also obtain $a \equiv 1 \pmod N$. Hence $\gamma \in \Gamma_1(N)$, and since $z=\gamma w$, the points $z$ and $w$ determine the same element of $\Gamma_1(N)\backslash\mathbb{H}$. This proves injectivity. [/step]

[step:Conclude the moduli interpretation] The construction $z \mapsto (E_z,P_z)$ is well-defined on $\Gamma_1(N)$-orbits, every pair $(E,P)$ with $P$ of exact order $N$ is isomorphic to some $(E_z,P_z)$, and two points $z,w \in \mathbb{H}$ give isomorphic marked elliptic curves exactly when they lie in the same $\Gamma_1(N)$-orbit. Therefore the map \begin{align*} \Gamma_1(N)\backslash \mathbb{H} &\longrightarrow \left\{ \text{isomorphism classes of pairs } (E,P) \right\} \\ [z] &\longmapsto [E_z,P_z] \end{align*} is a bijection. This is precisely the claimed moduli interpretation of $Y_1(N)(\mathbb{C})$. [/step]