I use the property that bilinear transformation maps inverse points w. r. t a circle to inverse points w. r. t. the image circle
Suppose that $w = \frac{az+b}{cz+d }$ is a transformation mapping $|z|=R$ onto $|w|=R$
Now $w=0$ , $w=\infty$ are inverse points for $|w|=R$ and they are transforms of
$z=-\frac{b}{a}, z=-\frac{d}{c}$
Respectively
$\implies -\frac{b}{a}, -\frac{d}{c}$
are inverse point for $|z|=R$
if we write $\alpha = -\frac{b}{a}$
inverse point of any point $\alpha$ W. R. T circle is $|z|=R$ is $\frac{R^2}{\bar{\alpha }}$
$ \implies \frac{-d}{c}= \frac{R^2}{\bar{\alpha }}$
So we rewrite $ w$
$w= \frac{a(z+b/a) }{c(z+d/c) }$
Using above relations
$w = \frac{a}{c} \frac{z-\alpha }{z-
R^2/ \bar{\alpha }}$
$w= \frac{a\bar{\alpha }}{c}. \frac{z-\alpha }{\bar{\alpha }z-R^2}$
Let $K =\frac{a\bar{\alpha }}{c}.$
$$w = K \frac{z-\alpha }{\bar{\alpha }z-R^2}$$
Using fact $|w|=R, |z|=R$ you can easily verify $|K|=R^2$