The Euclidean algorithm in algebra is applied to a class of gravitational lenses for which the lens equation consists of any set of coupled polynomial equations in the image position. In general, this algorithm allows us to reduce an apparently coupled system to a single polynomial in one variable (say $x$ in Cartesian coordinates) without the other component (say $y$), which is expressed as a function of the first component. This reduction enables us to investigate the lensing properties in an algebraic manner: For instance, we can obtain an analytic expression of the caustics by computing the discriminant of the polynomial equation. To illustrate this Euclidean algorithm, we re-examine a binary gravitational lens and show that the lens equation is reduced to a single real fifth-order equation, in agreement with previous works. We apply this algorithm also to the linearlized Kerr lens and find that the lens equation is reduced to a single real fifth-order one.

It has been recently shown that the lens equation for a binary gravitational lens, which is apparently a coupled system, can be reduced to a real fifth-order (quintic) algebraic equation. Some algebraic properties of the real quintic equation are revealed. We find that the number of images on each side of the separation axis is independent of the mass ratio and separation unless the source crosses the caustics. Furthermore, the discriminant of the quintic equation enables us to study changes in the number of solutions, namely in the number of images. It is shown that this discriminant can be factorized into two parts: One represents the condition that the lens equation can be reduced to a single quintic equation, while the other corresponds to the caustics.