My understanding is that Riemann conceived of his idea in order to devise a coordinate system in which multi-valued functions become single-valued so that the foundation of Complex Analysis of single-valued functions could be applied to multivalued functions. The normal Riemann surface for the square root is a Riemann sphere with a double covering of the complex plane mapped onto it conformally with every point on the sphere equal to ##(z,\sqrt{z})##. And therefore for the function ##f(z)=\sqrt{z}## over this new coordinate system, we have the simple linear function ##f(z,w)=w## which is an entire function zero over any closed contour over its Riemann surface. Contrast this with ##f(z)=\sqrt{z}+\frac{k}{z}##. What do you suppose the integral of this function over a closed contour (winding number of one) about the point ##(0,\infty)## over it's Riemann surface is? ##4k\pi i##? What about ##f(z)=\sqrt[n]{z}+\frac{k}{z}##? ##2nk\pi i##? What might the Riemann surface of ##\sqrt[n]{z}## even look like? I imagine it like one of those cantelopes with deep furrows with each section between the furrows corresponding to one copy of the z-plane. Once this perspective is taken, we can then re-cast the Residue Theorem for multivalued functions over their normal Riemann surfaces, then Laurent's Theorem and then all the other theorems for single-valued functions.
For example, the Laurent series for ##\sqrt{z}## is ##\sqrt{z}##, the series for ##\sqrt{z}+\frac{1}{z}## is ##\frac{1}{\left(\sqrt{z}\right)^2}+\sqrt{z}##. That was easy. What would the Laurent series for ##\frac{1}{\sqrt{z(z-1)(z-2)}}## be about the point ##(0,0)## over its Riemann surface? But if you know that all the theorems of single-valued functions apply to multivalued functions over their Riemann surface, how might you calculate the coefficients for this series expansion? Use Laurent's Theorem.