The cylinder is given by the equation $x^2 + (y-\frac{a}{2})^2 = (\frac{a}{2})^2$.
The region of the cylinder is given by the limits $0 \le \theta \le \pi$, $0 \le r \le a\sin \theta$ in polar coordinates.
We need to only calculate the surface from a hemisphere and multiply it by two. By implicit functions we have:
$$A=2\iint\frac{\sqrt{\left(\frac{\partial F}{\partial x}\right)^2 + \left(\frac{\partial F}{\partial y}\right)^2 + \left(\frac{\partial F}{\partial z}\right)^2}}{\left|\frac{\partial F}{\partial z} \right|} dA$$
where $F$ is the equation of the sphere.
Plugging in the expressions and simplifying ($z \ge 0)$, we get:
$$A=2a\iint\frac{1}{\sqrt{a^2 - x^2 - y^2}} dxdy$$
Converting to polar coordinates, we have:
$$A = 2a \int_{0}^\pi \int_{0}^{a\sin(\theta)} \frac{r}{\sqrt{a^2 - r^2}} drd\theta$$
Calculating this I get $2\pi a^2$. The answer is $(2\pi - 4)a^2$. Where am I going wrong?