I have tried (b) and had some clue about it. Basically I choose the intersection of the original to map to the intersection in the codomain. That means choosing $\{0, \infty\} \mapsto \{0, 1\}$. After that, I can get only a form like $z \mapsto \frac{az}{az+d}$. And then I give it $x \in \mathbb{R}$ and $iy, y \in \mathbb{R}$ to ensure the two condition is satisfied. With that, I can get $ad \in \mathbb{R},\bar{a}d \in \mathbb{R}$. And I arbitrarily chose a solution and it works.
My question is if there is a general way to do this kind of question? Because in the above, there are too many guessing.