I think the method is fine. The division algorithm steps can be used to write $1$ as a linear combination of $4a^2+1$ and $2a+1$. Thus $4a^2+1$ and $2a+1$ must be relatively prime.
@Shaun I might have said a few things about transcendence degree that seem to be useful, or were useful in similar problems. Likewise for P^n (i.e. stating whether you know and understand the standard affine cover would be helpful). It's tough (for me) to answer a question in which the asker does not state what they know already.
In general, you have $R/I\otimes_RM\cong M/IM$ for ideals $I$ and $R$-modules $M$, and this allows you to figure out what $R/m\otimes_RmN$ is. (E.g., consider the case $N=R$).
@Wave Briefly, you can localize the exact sequence of $R$-module $0\to p\to R \to R/p\to 0$, using the multiplicative set $S=R\setminus p$. That leads to $0\to pR_p\to R_p\to S^{-1}(R/p)\to 0$. Then, one shows that $S^{-1}(R/p)$ is the field of fractions of $R/p$.
Hold up, it is a direct application of the maximum principle. The real part of $f$ is constant, so $|e^f|$ is constant and a fortiori attains its maximum on the connected open domain. Hence, $e^f$ must be constant. This implies $f$ is discrete-valued and thus locally constant. @Wave