Let $C_1$ and $C_2$ be two circles of unequal radius. Circles $C_1$ and $C_2$ intersect at points $A$ and $B$; let $L_1$ be the tangent line to $C_1$ at $A$, and let $L_2$ be the tangent line to $C_2$ at $B$, and let $P$ be the intersection of $L_1$ and $L_2$. Let $M$ and $N$ be points in $C_1$ and $C_2$, respectively, such that $PM$ is tangent to $C_1$, and $PN$ is tangent to $C_2$.
Let $AM$ and $BP$ intersect at $S$, and let $BN$ and $AP$ intersect at $T$. Show that $\square ATBS$ is cyclic.
It was from at list of Power of a Point, so there should be a solution using power of a point.
I tried to show that triangles $\triangle APM$ and $\triangle BPN$ are similar. Also, I tried a lot of things using power of a point, but I did not make any progress.

I tried to solve it with regular geometry and found out that $\overline{AQ}\parallel\overline{BR}$, where $Q=(PB\cap\circ{O1})\setminus{B}$ and $R=(PA\cap\circ{O2})\setminus{A}$.
Usually it is not a bad idea to try inversion, if you try to prove points are cyclic, hence I have performed an inversion with a center at $A$ and the following result is obtained:
I have tried to color-code each line, so it is easier to find the corresponding lines and/or circles, but in generall:
Since points $A,T,B,S$ need to be on a circle it results into $T,B,S$ need to be colinear in the inverted problem. I have labeled most of the angles and the only thing to prove is $\alpha=\beta$ or $\epsilon=\varphi$. Most of the angles can be calculated by inscribed circles or angles between circle segments and tangents (if you have any questions why two angles are equal, feel free to ask, but by doing it on your own you learn a lot).
Now we have from triangles $\triangle SQB$: $$\delta+\gamma+\varphi+\beta=\pi$$ $$\frac{\sin\varphi}{\sin\gamma}=\frac{\overline{SM}}{\overline{QM}}=\frac{\overline{SM}}{\overline{MB}}=\frac{\sin\beta}{\sin\delta}\Rightarrow \frac{\sin\varphi}{\sin\beta}=\frac{\sin\gamma}{\sin\delta}$$
and from $\triangle PNT$ ($\measuredangle PNR = \epsilon$): $$\delta+\gamma+\epsilon+\alpha=\pi$$ $$\frac{\sin\epsilon}{\sin\alpha}=\frac{\overline{PR}}{\overline{NR}}=\frac{\overline{TR}}{\overline{NR}}=\frac{\sin\gamma}{\sin\delta}\Rightarrow \frac{\sin\epsilon}{\sin\alpha}=\frac{\sin\gamma}{\sin\delta}$$
Now we have: $$\delta+\gamma+\varphi+\beta=\pi=\delta+\gamma+\epsilon+\alpha\Rightarrow\varphi+\beta=\epsilon+\alpha$$ $$\frac{\sin\varphi}{\sin\beta}=\frac{\sin\gamma}{\sin\delta}=\frac{\sin\epsilon}{\sin\alpha}\Rightarrow\frac{\sin\varphi}{\sin\beta}=\frac{\sin\epsilon}{\sin\alpha}$$
Modifying the trigonometric equation gives us: $$\sin\alpha\cdot\sin\varphi=\sin\beta\cdot\sin\epsilon$$ $$\frac{1}{2}(\cos(\alpha-\varphi)+\cos(\alpha+\varphi))=\frac{1}{2}(\cos(\beta-\epsilon)+\cos(\beta+\epsilon))$$
Since $\varphi+\beta=\epsilon+\alpha\Rightarrow\varphi-\alpha=\epsilon-\beta\Rightarrow\cos(\varphi-\alpha)=\cos(\epsilon-\beta)$ we get that $$\cos(\alpha+\varphi)=\cos(\beta+\epsilon)$$ $$\cos(\alpha+\varphi)-\cos(\beta+\epsilon)=-2\sin(\frac{\alpha+\varphi+\beta+\epsilon}{2})\sin(\frac{\alpha+\varphi-\beta-\epsilon}{2})=0$$ Following solutions arise:
Hence only the first solution is possible from which follows that $$\alpha=\beta \land \varphi=\epsilon$$
This implies that the points $S,B,T$ must lie on one line and will be on the same circle with the point $A$.